Scheduled Downtime
On Tuesday 24 October 2023 @ 5pm MT the forums will be in read only mode in preparation for the downtime. On Wednesday 25 October 2023 @ 5am MT, this website will be down for maintenance and expected to return online later in the morning.
Normal Operations
The forums are back online with normal operations. If you notice any issues or errors related to the forums, please reach out to help@ucar.edu

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

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.

 
Top