how to set scale factors for SpinupStability.ncl in regional simulation

Status
Not open for further replies.

xiaoxiaokuishu

Ru Xu
Member
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...
1748253947755.png
 

oleson

Keith Oleson
CSEG and Liaisons
Staff member
We haven't formally offered a method for setting these thresholds, it's up to the user. See here for some possibly useful discussion on this topic.

 

xiaoxiaokuishu

Ru Xu
Member
We haven't formally offered a method for setting these thresholds, it's up to the user. See here for some possibly useful discussion on this topic.

Hi, Oleson
The provided link is very helpful, but I have 2 further questions,
(1) the scale factor is set as A/B, A=land pixel number of my region, B=land pixel number of globel,
not sure how to get the land numbers.. my res is 0.5deg.
(2)From your shared link, I find her output from the SpinupStability.ncl (under clm5.0/tools/contrib )is a regional map, my is still global..
do you know what happened...


1748424996639.png
 
Status
Not open for further replies.
Back
Top