4x5 BGC_crop cases for CUPID tutorial

[wwieder]

Description:
We decided that the cupid tutorial will run more smoothly if students just run a new case from the current CTSM tag to compare with a the 5.4 baseline. I also decided that diagnostics and CUPiD plots would be better if this was compared to a CLM5 physics case with GSWP3 forcing.

To generate these case I decided to do a full spinup and historical run (maybe overkill), but mainly because we don't have any 4x5 cases for comparision and no recent runs for CLM5 initial conditions.

---

Case directory:
Locally:
/glade/work/wwieder/ctsm5.4.0/ctsm5.4.002_clm6_BGCcrop_crujra_4x5_AD
/glade/work/wwieder/ctsm5.4.0/ctsm5.4.002_clm5_BGCcrop_gswp3_4x5_AD

---

Sandbox:
Locally:
/glade/u/home/wwieder/CTSM (ctsm5.4.002)

CLM6 physics

./create_newcase --case /glade/work/wwieder/ctsm5.4.0/ctsm5.4.002_clm6_BGCcrop_crujra_4x5_AD --compset I1850Clm60BgcCropCrujraG --res f45_g37 --run-unsupported

./create_clone --case /glade/work/wwieder/ctsm5.4.0/ctsm5.4.002_clm6_BGCcrop_crujra_4x5_SASU --clone /glade/work/wwieder/ctsm5.4.0/ctsm5.4.002_clm6_BGCcrop_crujra_4x5_AD --keepexe

./create_clone --case /glade/work/wwieder/ctsm5.4.0/ctsm5.4.002_clm6_BGCcrop_crujra_4x5_pSASU --clone /glade/work/wwieder/ctsm5.4.0/ctsm5.4.002_clm6_BGCcrop_crujra_4x5_AD --keepexe

./create_newcase --case /glade/work/wwieder/ctsm5.4.0/ctsm5.4.002_clm6_BGCcrop_crujra_4x5_HIST --compset IHistClm60BgcCropCrujra --res f45_g37 --run-unsupported

and for the CLM5 physics

./create_newcase --case /glade/work/wwieder/ctsm5.4.0/ctsm5.4.002_clm5_BGCcrop_gswp3_4x5_AD --compset I1850Clm50BgcCropG --res f45_g37 --run-unsupported

./create_clone --clone /glade/work/wwieder/ctsm5.4.0/ctsm5.4.002_clm5_BGCcrop_gswp3_4x5_AD --case /glade/work/wwieder/ctsm5.4.0/ctsm5.4.002_clm5_BGCcrop_gswp3_4x5_SASU --keepexe

./create_clone --clone /glade/work/wwieder/ctsm5.4.0/ctsm5.4.002_clm5_BGCcrop_gswp3_4x5_AD --case /glade/work/wwieder/ctsm5.4.0/ctsm5.4.002_clm5_BGCcrop_gswp3_4x5_pSASU --keepexe

./create_newcase --case /glade/work/wwieder/ctsm5.4.0/ctsm5.4.002_clm5_BGCcrop_gswp3_4x5_HIST --compset IHistClm50BgcCrop --res f45_g37 --run-unsupported

requests 5 nodes & getting 180 year / day in AD mode

---

user_nl_ changes:

1. AD

./xmlchange MOSART_MODE=NULL
./xmlchange STOP_N=80,STOP_OPTION=nyears
./xmlchange CLM_ACCELERATED_SPINUP=on,CLM_FORCE_COLDSTART=on
./xmlchange CONTINUE_RUN=FALSE,RESUBMIT=3

2. SASU

export REFCASE=ctsm5.4.002_clm6_BGCcrop_crujra_4x5_AD
or
export REFCASE=ctsm5.4.002_clm5_BGCcrop_gswp3_4x5_AD

./xmlchange MOSART_MODE=NULL
./xmlchange RUN_TYPE=hybrid
./xmlchange STOP_N=60,STOP_OPTION=nyears
./xmlchange CLM_ACCELERATED_SPINUP=sasu,CLM_FORCE_COLDSTART=off
./xmlchange CONTINUE_RUN=FALSE,RESUBMIT=5
./xmlchange GET_REFCASE=TRUE
./xmlchange RUN_REFDATE=0361-01-01
./xmlchange DRV_RESTART_POINTER=rpointer.cpl.0361-01-01-00000
./xmlchange RUN_REFCASE=$REFCASE
./xmlchange RUN_REFDIR=/glade/derecho/scratch/wwieder/archive/$REFCASE/rest/0361-01-01-00000

3. pSASU

export REFCASE=ctsm5.4.002_clm6_BGCcrop_crujra_4x5_SASU
or
export REFCASE=ctsm5.4.002_clm5_BGCcrop_gswp3_4x5_SASU


./xmlchange MOSART_MODE=NULL
./xmlchange CLM_ACCELERATED_SPINUP=off,CLM_FORCE_COLDSTART=off
./xmlchange CONTINUE_RUN=FALSE,RESUBMIT=0
./xmlchange STOP_N=180,STOP_OPTION=nyears,REST_N=60
./xmlchange RUN_TYPE=hybrid
./xmlchange GET_REFCASE=TRUE
./xmlchange RUN_REFDATE=0361-01-01
./xmlchange DRV_RESTART_POINTER=rpointer.cpl.0361-01-01-00000
./xmlchange RUN_REFCASE=$REFCASE
./xmlchange RUN_REFDIR=/glade/derecho/scratch/wwieder/archive/$REFCASE/rest/0361-01-01-00000

4. HIST
   user_nl_datm_streams

!! manually update datm_streams to get through 2023
co2tseries.20tr:year_last = 2024
co2tseries.20tr:datafiles = /glade/campaign/cesm/cesmdata/inputdata/atm/datm7/CO2/fco2_datm_globalSSP3-7.0_simyr_1750-2501_CMIP6_c201101.nc

user_nl_clm

finidat = '/glade/derecho/scratch/wwieder/archive/ctsm5.4.002_clm5_BGCcrop_gswp3_4x5_pSASU/rest/0181-01-01-00000/ctsm5.4.002_clm5_BGCcrop_gswp3_4x5_pSASU.clm2.r.0181-01-01-00000.nc'

use_init_interp = .true.
init_interp_fill_missing_urban_with_HD = .true.
ndep_taxmode = 'extend'

! "Survivability" output
hist_dov2xy(1) = .true.
hist_dov2xy(2) = .false.
hist_fields_list_file = .false.
hist_fincl2 = 'TLAI', 'TSA', 'TREFMNAV', 'TREFMXAV', 'BTRANMN', 'Vcmx25Z', 'FSH', 'VEGWP', 'FCTR', 'FCEV', 'FGEV',
        'FIRE', 'FSR', 'FIRA', 'FSA', 'GSSUNLN', 'GSSHALN', 'TSKIN', 'GPP', 'NPP', 'AGNPP', 'TOTVEGC',
        'NPP_NUPTAKE', 'AR', 'HR', 'HTOP'
hist_mfilt(1) = 1
hist_mfilt(2) = 1
hist_nhtfrq(1) = 0
hist_nhtfrq(2) = 0
hist_type1d_pertape(1) = ' '
hist_type1d_pertape(2) = ' '

1850-1901

./xmlchange RUN_TYPE=startup
./xmlchange RUN_STARTDATE=1850-01-01
./xmlchange CONTINUE_RUN=FALSE,RESUBMIT=0
./xmlchange DATM_YR_ALIGN=1901,DATM_YR_START=1901
./xmlchange DATM_YR_END=1920
./xmlchange STOP_N=51,STOP_OPTION=nyears
./case.setup

1901-2023

./xmlchange DATM_YR_END=2023
./xmlchange DATM_SKIP_RESTART_READ=TRUE
./xmlchange STOP_N=122
./xmlchange RESUBMIT=0

---

SourceMods:

---

Diagnostics:
Diags (if available)

https://webext.cgd.ucar.edu/I1850/$CASE/lnd/

---

Output:
Output (if still available):
/glade/scratch/$USER/archive/$CASE/

---

Contacts:
@olyson

---

Extra details: