Hi, all
I run a regional simulation in spinup mode, I first do a 400-years AD spinup, then a post-AD ( 400 years),
then I evalute is the model balanced at post-AD stage, as I do a regional simulation, I need to give the
scale fators for the threshold..
I set it as scale_factor=2640/(360*720*0.29) here, 2640 is all the pixel numbers of my simulation (incluing land and ocean),
360*720 is my res is 0.5deg. Not sure if it is the right way to set the threshold.
glob_thresh_totecosysc = 0.02*scale_factor ; global threshold for TOTECOSYSC equilibrium (delta PgC / yr)
glob_thresh_totsomc = 0.02*scale_factor ; global threshold for TOTSOMC equilibrium (delta PgC / yr)
glob_thresh_totvegc = 0.02*scale_factor ; global threshold for TOTVEGC equilibrium (delta PgC / yr)
glob_thresh_tlai = 0.02 ; global threshold for TLAI equilibrium (delta m2/m2 / yr)
glob_thresh_gpp = 0.02*scale_factor ; global threshold for GPP equilibrium (delta PgC / yr)
glob_thresh_tws = 0.001 ; global threshold for TWS equilibrium (delta m / yr)
glob_thresh_h2osno = 1.0 ; global threshold for H2OSNO equilibrium (delta mm / yr)
glob_thresh_area = 3.0*scale_factor ; global threshold percent area with TOTECOSYSC disequilibrium gt 1 g C/m2/yr
totecosysc_thresh = 1. ; disequilibrium threshold for individual gridcells (gC/m2/yr)
end if
after I run the sripts, it present like below, it seems my carbon is balanced? but land surface is not balanced?
not sure how can I do next...

I run a regional simulation in spinup mode, I first do a 400-years AD spinup, then a post-AD ( 400 years),
then I evalute is the model balanced at post-AD stage, as I do a regional simulation, I need to give the
scale fators for the threshold..
I set it as scale_factor=2640/(360*720*0.29) here, 2640 is all the pixel numbers of my simulation (incluing land and ocean),
360*720 is my res is 0.5deg. Not sure if it is the right way to set the threshold.
glob_thresh_totecosysc = 0.02*scale_factor ; global threshold for TOTECOSYSC equilibrium (delta PgC / yr)
glob_thresh_totsomc = 0.02*scale_factor ; global threshold for TOTSOMC equilibrium (delta PgC / yr)
glob_thresh_totvegc = 0.02*scale_factor ; global threshold for TOTVEGC equilibrium (delta PgC / yr)
glob_thresh_tlai = 0.02 ; global threshold for TLAI equilibrium (delta m2/m2 / yr)
glob_thresh_gpp = 0.02*scale_factor ; global threshold for GPP equilibrium (delta PgC / yr)
glob_thresh_tws = 0.001 ; global threshold for TWS equilibrium (delta m / yr)
glob_thresh_h2osno = 1.0 ; global threshold for H2OSNO equilibrium (delta mm / yr)
glob_thresh_area = 3.0*scale_factor ; global threshold percent area with TOTECOSYSC disequilibrium gt 1 g C/m2/yr
totecosysc_thresh = 1. ; disequilibrium threshold for individual gridcells (gC/m2/yr)
end if
after I run the sripts, it present like below, it seems my carbon is balanced? but land surface is not balanced?
not sure how can I do next...
