Skip to content

b4b-dev: Update SpinupStability_BGC_v12_SE.ncl to work for FATES and ne16pg3 #3776

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
slevis-lmwg merged 3 commits into from
Mar 12, 2026
+38 −16
Merged

b4b-dev: Update SpinupStability_BGC_v12_SE.ncl to work for FATES and ne16pg3 #3776

Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
9e87625
Update SpinupStability_BGC_v12_SE.ncl to work for FATES and ne16pg3
slevis-lmwg Feb 19, 2026
f243530
Merge branch 'b4b-dev' into spinupstab_script_forfates
ekluzek Mar 12, 2026
83ba68b
Merge branch 'b4b-dev' into spinupstab_script_forfates
slevis-lmwg Mar 12, 2026
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
54 changes: 38 additions & 16 deletions tools/contrib/SpinupStability_BGC_v12_SE.ncl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -2,8 +2,9 @@
; SpinupStability_BGC_v12_SE.ncl
; Script to examine stability of spinup simulation.
; This version operates on either monthly mean or multi-annual mean multi-variable history files
; and supports ne30 grids only
; it supports "ne" grids only, and it supports FATES simulations
; Keith Oleson, July 2025
; Sam Levis, Feb 2026
;
; ARH mods, April 2020 - modified to work for unstructured grids via remapping to FV 1 deg. grid.
; KO mods, March 2025 - regrid history file using ncremap to get area and landfrac instead of
Expand Down Expand Up @@ -35,24 +36,27 @@ begin
;=======================================================================;

; GLOBAL EXAMPLE
caseid = "ctsm53n02ctsm52027_ne30pg3t232_AD"
username = "oleson"
annual_hist = True
caseid = "ctsm60_fates_nocomp_adspinup"
username = "afoster"
annual_hist = False
region = "Global" ; Global
subper = 20 ; Subsampling period in years
paleo = False ; Use paleo map ; UNTESTED IN THIS VERSION
fates = True ; FATES simulation?
hist_ext = "h0a" ; Valid options: h0a, h0

;--- ARH mods ---
map_method = "conserve" ;bilinear,conserve or patch
dout_s = False ;short term storage (archive) directory structure
dout_s = True ;short term storage (archive) directory structure

;srcGrid = "ne120pg3"
;srcGridName="/glade/campaign/cesm/cesmdata/cseg/inputdata/atm/cam/coords/ne120pg3_scrip_c170628.nc"
srcGrid = "ne30pg3"
srcGridName="/glade/campaign/cesm/cesmdata/cseg/inputdata/atm/cam/coords/ne30pg3_scrip_170604.nc"
;srcGrid = "ne30pg3"
;srcGridName="/glade/campaign/cesm/cesmdata/cseg/inputdata/atm/cam/coords/ne30pg3_scrip_170604.nc"
;srcGrid = "ne30pg2"
;srcGridName="/glade/campaign/cesm/cesmdata/cseg/inputdata/atm/cam/coords/ne30pg2_scrip_c170608.nc"
srcGrid = "ne16pg3"
srcGridName="/glade/campaign/cesm/cesmdata/cseg/inputdata/atm/cam/coords/ne16pg3_scrip_170725.nc"

dstGrid = "f09"
dstGridName="/glade/campaign/cesm/cesmdata/cseg/inputdata/share/scripgrids/0.9x1.25_SCRIP_desc.181018.nc"
Expand Down Expand Up @@ -144,9 +148,15 @@ begin
if (isfilevardim(data[0],"landfrac","lndgrid")) then
totecosysc = data[:]->TOTECOSYSC(year,:)
totsomc = data[:]->TOTSOMC(year,:)
totvegc = data[:]->TOTVEGC(year,:)
tlai = data[:]->TLAI(year,:)
gpp = data[:]->GPP(year,:)
if (fates) then
totvegc = data[:]->FATES_VEGC(year,:) * 1000 ; kg -> g
tlai = data[:]->FATES_LAI(year,:)
gpp = data[:]->FATES_GPP(year,:) * 1000 ; kg -> g
else
totvegc = data[:]->TOTVEGC(year,:)
tlai = data[:]->TLAI(year,:)
gpp = data[:]->GPP(year,:)
end if
tws = data[:]->TWS(year,:)
if (isfilevar(data[0],"H2OSNO")) then
h2osno = data[:]->H2OSNO(year,:)
Expand All @@ -157,9 +167,15 @@ begin
else
totecosysc = data[:]->TOTECOSYSC(year,:,:)
totsomc = data[:]->TOTSOMC(year,:,:)
totvegc = data[:]->TOTVEGC(year,:,:)
tlai = data[:]->TLAI(year,:,:)
gpp = data[:]->GPP(year,:,:)
if (fates) then
totvegc = data[:]->FATES_VEGC(year,:,:) * 1000 ; kg -> g
tlai = data[:]->FATES_LAI(year,:,:)
gpp = data[:]->FATES_GPP(year,:,:) * 1000 ; kg -> g
else
totvegc = data[:]->TOTVEGC(year,:,:)
tlai = data[:]->TLAI(year,:,:)
gpp = data[:]->GPP(year,:,:)
end if
tws = data[:]->TWS(year,:,:)
if (isfilevar(data[0],"H2OSNO")) then
h2osno = data[:]->H2OSNO(year,:,:)
Expand All @@ -171,9 +187,15 @@ begin
else
totecosysc = month_to_annual(data[:]->TOTECOSYSC,1)
totsomc = month_to_annual(data[:]->TOTSOMC,1)
totvegc = month_to_annual(data[:]->TOTVEGC,1)
tlai = month_to_annual(data[:]->TLAI,1)
gpp = month_to_annual(data[:]->GPP,1)
if (fates) then
totvegc = month_to_annual(data[:]->FATES_VEGC,1) * 1000 ; kg -> g
tlai = month_to_annual(data[:]->FATES_LAI,1)
gpp = month_to_annual(data[:]->FATES_GPP,1) * 1000 ; kg -> g
else
totvegc = month_to_annual(data[:]->TOTVEGC,1)
tlai = month_to_annual(data[:]->TLAI,1)
gpp = month_to_annual(data[:]->GPP,1)
end if
if (isfilevar(data[0],"TWS")) then
tws = month_to_annual(data[:]->TWS,1)
else
Expand Down
Loading