ctsm5.4.029_f19_wetlandCH4_157_ [AD, SASU, pSASU, HIST]

[wwieder]

Description:

• Latest updates to wetland branch + parameter changes to increase both FCH4 and unsaturated oxidation.
• Turning off isotopes for speed.
• This is a crop case (mainly for CLM6 evaluation), but likely should use noCrop for WIEMIP.
• Starting from a clean cold start, AD, SASU, etc (tried shortcutting this in 156).

Parameter changes from the generic crop PR

f_ch4=0.16
porosmin=1.e-6
unsat_aere_ratio=0.01
vmax_oxid_unsat=1.25e-05

---

Case directory:
Locally:
/glade/work/wwieder/ctsm5.4.0/ctsm5.4.029_f19_wetlandCH4_157*

---

Sandbox:
Locally:
/glade/u/home/wwieder/wetland_update
(commit 3af9c708ca86dfdd7bee3c65371b5060d4fa3baf (HEAD -> 5.4.024, swensosc/wetland_update))

./create_newcase --case /glade/work/wwieder/ctsm5.4.0/ctsm5.4.029_f19_wetlandCH4_157_AD --compset I1850Clm60BgcCropCrujra --res f19_g17 --run-unsupported

./create_newcase --case /glade/work/wwieder/ctsm5.4.0/ctsm5.4.029_f19_wetlandCH4_157_HIST --compset HIST_DATM%CRUJRA2024_CLM60%BGC_SICE_SOCN_MOSART_SGLC_SWAV --res f19_g17 --run-unsupported

---

user_nl_ changes:

paramfile='/glade/work/wwieder/modified_inputs/ctsm60_params.genericCropFix.ch4cal_c.c260410.nc'
use_c13 = .false.
use_c14 = .false.

---

1. AD

./xmlchange MOSART_MODE=NULL
./xmlchange CLM_ACCELERATED_SPINUP=on,CLM_FORCE_COLDSTART=on
./xmlchange STOP_N=120,STOP_OPTION=nyears,REST_N=60
./xmlchange CONTINUE_RUN=FALSE,RESUBMIT=2

2. SASU

export REFCASE=/ctsm5.4.029_f19_wetlandCH4_157_AD 

./xmlchange MOSART_MODE=NULL
./xmlchange RUN_TYPE=hybrid
./xmlchange STOP_N=120,REST_N=60,STOP_OPTION=nyears
./xmlchange CLM_ACCELERATED_SPINUP=sasu,CLM_FORCE_COLDSTART=off
./xmlchange CONTINUE_RUN=FALSE,RESUBMIT=2
./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
./xmlchange --subgroup case.run JOB_WALLCLOCK_TIME=12:00:00

3. pSASU

export REFCASE=/ctsm5.4.029_f19_wetlandCH4_157_SASU

/xmlchange MOSART_MODE=NULL
./xmlchange STOP_N=120,STOP_OPTION=nyears,REST_N=60
./xmlchange CONTINUE_RUN=FALSE,RESUBMIT=1
./xmlchange CLM_ACCELERATED_SPINUP=off,CLM_FORCE_COLDSTART=off
./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/$USER/archive/$REFCASE/rest/0361-01-01-00000
./xmlchange --subgroup case.run JOB_WALLCLOCK_TIME=12:00:00

4. HIST
   user_nl_clm

paramfile='/glade/work/wwieder/modified_inputs/ctsm60_params.genericCropFix.ch4cal_c.c260410.nc'
use_c13 = .false.
use_c14 = .false.


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) = ' '

finidat='/glade/derecho/scratch/wwieder/archive/ctsm5.4.029_f19_wetlandCH4_157_pSASU/rest/0241-01-01-00000/ctsm5.4.029_f19_wetlandCH4_157_pSASU.clm2.r.0241-01-01-00000.nc'

1850-1901

./xmlchange RUN_TYPE=startup
./xmlchange RUN_STARTDATE=1850-01-01
./xmlchange STOP_OPTION=nyears,CONTINUE_RUN=FALSE
./xmlchange RESUBMIT=0,STOP_N=51
./xmlchange DATM_YR_ALIGN=1901,DATM_YR_START=1901
./xmlchange DATM_YR_END=1920
./xmlchange --subgroup case.run JOB_WALLCLOCK_TIME=06:00:00

echo "finidat='/glade/derecho/scratch/wwieder/archive/ctsm5.4.029_f19_wetlandCH4_157_pSASU/rest/0241-01-01-00000/ctsm5.4.029_f19_wetlandCH4_157_pSASU.clm2.r.0241-01-01-00000.nc'" >> user_nl_clm

1901-2024

./xmlchange DATM_YR_END=2023
./xmlchange DATM_SKIP_RESTART_READ=TRUE
./xmlchange CONTINUE_RUN=True
./xmlchange STOP_N=123,RESUBMIT=1
./xmlchange --subgroup case.run JOB_WALLCLOCK_TIME=12:00:00

---

SourceMods:

---

Diagnostics:
Diags (if available)

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

---

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

---

Contacts:
@wwieder

---

Extra details:

[wwieder]

@linniahawkins I'm doing a:

• BGC crop run that's based on this case:
  /glade/work/wwieder/ctsm5.4.0/ctsm5.4.029_f19_wetlandCH4_157_cropHIST

• a parallel noCrop case using the new generic crop parameter file is here:
  /glade/work/wwieder/ctsm5.4.0/ctsm5.4.029_f19_wetlandCH4_157_HIST

If I forget to run LDF on this remind me next week.

[wwieder]

@linniahawkins diagnostics for crop vs. noCrop runs are here.

These are f19 cases that were each initialized with an AD, SASU, pSASU case from a crop case.

[wwieder]

Seems like global GPP looks as expected (higher in the crop case).

Drilling in a little closer to look at seasonal patterns produces a bit more of a nuanced look at C & water fluxes (crop case with higher GPP, but not much difference in LAI or ET for the Central US, below). On ballance I'm not seeing anything overly concerning. What do you think @linniahawkins?

[linniahawkins]

Thanks - I agree, looking through these nothing is jumping out so I think it seems reasonable. I'll plot the PFT level output tomorrow just out of curiosity.

One note (and it doesn't matter for these runs) but it looks like your CO2 stream might be recycling 2014 for the 2014-2023 period.

[wwieder]

Also very astute, @linniahawkins! I forgot to add this to the user_nl_datm_streams files in one of these recent runs so opted to consistently make the error in all of them. Hopefully we can get this ironed out soon, @ekluzek is working on our compset definitions now.