FMT2000 ctsm53 control (f.e30_alpha04a.FMT2000BgcCrop.CTRL.lmwg_dev_76)

[slevis-lmwg]

Description:
Investigating the dry bias in the Amazon seen in cesm_dev 116 (NCAR/cesm_dev#19) and discussed in #75.

Similar to NCAR/cesm_dev#19 and NCAR/amwg_dev#614 but as F2000BgcCrop.

• modern code base beta04 tag (whatever was used for cesm_dev 116)
• Leung2023 dust; tuning turned on and other CAM mods as in cesm_dev 116 (or similar); Adam H can advise.
• CTSM5.3 using CLM6-BGC
• initial conditions from CTSM 5.3 Spinup and HIST on ne30 grid (2nd Attempt - with li2024crujra fire ON this time and using filled Antarctica datm data) #68
• Run 1 (this issue) with CTSM5.3(PPE) parameter file vs. Run 2 (FMT2000 ctsm53 with 52 params (f.e30_alpha04a.FMT2000BgcCrop.52params.lmwg_dev_77) #77) with CTSM5.x parameter file that either goes back to a CTSM5.1 parameter file OR has more target changes to investigate tropical precipitation and clouds? Change kmax back to ctsm5.2 values. Linnia will advise.
• Submit 40-yr simulation

---

Case directory:
Locally:
./create_newcase --compset 2000_CAM70%MT_CLM60%BGC-CROP_CICE%PRES_DOCN%DOM_MOSART_SGLC_SWAV_SESP --res ne30pg3_ne30pg3_mg17 --case ~slevis/cases_LMWG_dev/f.e30_alpha04a.FMT2000BgcCrop.CTRL.lmwg_dev_76 --run-unsupported

---

Sandbox:
Locally:
/glade/work/hannay/cesm_tags/cesm3_0_alpha04a

---

user_nl_ changes:
CLM

hist_avgflag_pertape = 'A','A'
hist_nhtfrq = 0,0
hist_mfilt  = 1,1
hist_ndens  = 2,2
hist_dov2xy = .true.,.false.
hist_fincl2 = 'TLAI', 'TSA','TREFMNAV','TREFMXAV', 'BTRANMN', 'Vcmx25Z', 'FSH', 'VEGWP', 'FCTR', 'FCEV', 'FGEV', 'FIRE', 'FSR', 'FIRA', 'FSA', 'GSSUNLN', 'GSSHALN', 'TSKIN'

reseed_dead_plants = .true.

! turn off shifting cultivation
do_grossunrep = .false.

CAM

! file h0 h1 h2 h3 h4
mfilt = 0
nhtfrq = 0
ndens = 2
interpolate_output = .true.
interpolate_nlat = 192
interpolate_nlon = 288

empty_htapes = .true.

! interpolated output
! output for ADF monthly-interpolated (h0)
fincl1 = 'AODDUST', 'AODVIS', 'AODDUSTdn','AODVISdn','AODBCdn','AODSOAdn','AODSSdn','AODPOMdn','AODSO4dn',
'BURDENBC', 'BURDENDUST', 'BURDENPOM', 'BURDENSEASALT',
'BURDENSO4', 'BURDENSOA', 'CDNUMC', 'CLDHGH', 'CLDICE', 'CLDLIQ', 'CLDLOW', 'CLDMED', 'CLDTOT', 'CLOUD', 'FLDS', 'FLNS', 'FLNSC', 'FLNT', 'FLNTC', 'FLUT',
'FSNS', 'FSNSC', 'FSNT', 'FSNTC', 'FSNTOA', 'ICEFRAC', 'LHFLX', 'LWCF','OMEGA500', 'PBLH','PRECC', 'PRECL', 'PRECSC', 'PRECSL', 'PRECT', 'PS', 'PSL','Q',
'QFLX', 'RELHUM', 'SHFLX', 'SOLIN', 'SST', 'SWCF', 'T', 'TAUX', 'TAUY', 'TGCLDIWP', 'TGCLDLWP', 'TMQ', 'TREFHT', 'TS', 'U', 'U10', 'V', 'LANDFRAC'

rad_climate = 'A:Q:H2O', 'N:O2:O2',
'A:CO2:CO2', 'N:ozone:O3',
'A:N2O:N2O', 'A:CH4:CH4', 'N:CFC11STAR:CFC11', 'A:CFC12:CFC12',
'M:mam4_mode1:/glade/campaign/cesm/cesmdata/inputdata/atm/cam/physprops/mam4_mode1_rrtmg_aeronetdust_sig1.6_dgnh.48_c140304.nc',
'M:mam4_mode2:/glade/campaign/cesm/cesmdata/inputdata/atm/cam/physprops/mam4_mode2_rrtmg_aitkendust_c141106.nc',
'M:mam4_mode3:/glade/campaign/cesm/cesmdata/inputdata/atm/cam/physprops/mam4_mode3_rrtmg_aeronetdust_c141106.nc',
'M:mam4_mode4:/glade/campaign/cesm/cesmdata/inputdata/atm/cam/physprops/mam4_mode4_rrtmg_c130628.nc',
'N:VOLC_MMR1:/glade/campaign/cesm/cesmdata/inputdata/atm/cam/physprops/volc_camRRTMG_byradius_sigma1.6_mode1_c210211.nc',
'N:VOLC_MMR2:/glade/campaign/cesm/cesmdata/inputdata/atm/cam/physprops/volc_camRRTMG_byradius_sigma1.6_mode2_c210211.nc',
'N:VOLC_MMR3:/glade/campaign/cesm/cesmdata/inputdata/atm/cam/physprops/volc_camRRTMG_byradius_sigma1.2_mode3_c210211.nc'

ubc_file_path= "/glade/campaign/cesm/cesmdata/inputdata/atm/cam/chem/ubc/b.e21.BWHIST.f09_g17.CMIP6-historical-WACCM.ensAvg123.cam.h0zm.H2O.1849-2014_c240604.nc"

se_statefreq=144

ncdata = '/glade/u/home/slevis/cases_LMWG_dev/f.e30_alpha04a.FMT2000BgcCrop.CTRL.lmwg_dev_76/ERA5_x_ne30np4_L93_rgC2_WO_aero.2000-01-01-00000.nc'

use_gw_front = .true.
taubgnd  = 1.25D-3

seasalt_emis_scale = 1.0D0
clubb_c8 = 4.4
dust_emis_method = 'Leung_2023'
dust_emis_fact = 2.3

radiation_scheme = 'rrtmgp'

---

SourceMods:
None

---

Diagnostics:
Diags (if available)
https://webext.cgd.ucar.edu/FMT2000/f.e30_alpha04a.FMT2000BgcCrop.CTRL.lmwg_dev_76/lnd/f.e30_alpha04a.FMT2000BgcCrop.CTRL.lmwg_dev_76.11-20-ctsm53n04ctsm52028_f09_hist_2004_2023/setsIndex.html

---

Output:
Output (if still available):
Initially:
/glade/derecho/scratch/slevis/f.e30_alpha04a.FMT2000BgcCrop.CTRL.lmwg_dev_76/run
Later:
/glade/derecho/scratch/slevis/archive/f.e30_alpha04a.FMT2000BgcCrop.CTRL.lmwg_dev_76
Finally at:
/glade/campaign/cgd/tss/people/slevis/f_cases/f.e30_alpha04a.FMT2000BgcCrop.CTRL.lmwg_dev_76

---

Contacts:
@slevis-lmwg

---

Extra details:
Regarding user_nl_cam:

• Due to lack of permissions, I had to change ncdata from
  /glade/campaign/cgd/amp/juliob/ERA5/ne30np4/L93_cdf5/ERA5_x_ne30np4_L93_rgC2_WO_aero.2000-01-01-00000.nc
  to
  /glade/u/home/slevis/cases_LMWG_dev/f.e30_alpha04a.FMT2000BgcCrop.CTRL.lmwg_dev_76/ERA5_x_ne30np4_L93_rgC2_WO_aero.2000-01-01-00000.nc

[slevis-lmwg]

@wwieder the case is ready to start, so please have a look when you get a chance.

[wwieder]

Hey Sam. Thanks for getting this set up.
There are a number of differences in the atm_in files between this case and the F_008 case on AMWG_dev.

Presumably these are mostly related to the start date (2000 here vs 1850 transient on AMWG_dev). The most notable diffs are below, with a quick (but very naive comment from me).

diff atm_in /glade/campaign/cesm/cesmdata/cseg/runs/cesm2_0/f.e30_alpha04a.FMTHIST.ne30_L93_tuning.008/CaseDocs/atm_in | more

<  clubb_c8		= 4.4.       *THIS SEEMS RIGHT TO ME, based on the B_116 case*
---
>  clubb_c8		= 4.95

<  radiation_scheme		= 'rrtmg'
---
>  radiation_scheme		= 'rrtmgp'.     *I THINK WE WANT THIS, again in both F008 and B116 cases*

[wwieder]

I think after you modify the radiation scheme this seems good to try running.

[slevis-lmwg]

@wwieder I should have asked while we were talking:
What's the minimum list of fields that we need in the h1 file?
You said TLAI, but anything else?

[wwieder]

I've been using this list that Keith typically has in his cases:

! h1 stream (monthly average, finest sub-grid)
! Emon, Lmon
hist_mfilt(2) = 1
hist_dov2xy(2) = .false.
hist_nhtfrq(2) = 0
hist_type1d_pertape(2) = ' '
hist_fincl2 += 'TLAI', 'TSA','TREFMNAV','TREFMXAV', 'BTRANMN', 'Vcmx25Z', 'FSH', 'VEGWP', 'FCTR', 'FCEV', 'FGEV', 'FIRE', 'FSR', 'FIRA', 'FSA', 'GSSUNLN', 'GSSHALN', 'TSKIN'

[wwieder]

@slevis-lmwg comparing this run with #77 will be helpful, but I'd also like to see results with run 116 / cesm_dev#19 to see if we can replicate the precip and cloud biases we're seeing in the B-case.

FWIW, @olyson had posted diagnostics comparing two recent B-cases here

[slevis-lmwg]

@olyson we could go over this on Weds morning or sooner if necessary:
I tried using your existing ANALYSIS/climo/...b116_14_23 data to compare with #76, but so far I have failed. I don't want to rerun diagnostics for that case from scratch because the regridding phase will take forever :-) so I'm hoping to troubleshoot with you.

[wwieder]

Hey, not sure where we are on diagnostics, but would it be hard to compare this run vs #69 (likely years 2000-2019)?

[slevis-lmwg]

Let's confirm in our meeting coming up.

UPDATE
Started diagnostics #76 vs. #69
Wait to start NCAR/cesm_dev#25 vs. #76 on Friday p.m.

[slevis-lmwg]

Land diags for #76 vs. #69 are now posted above.

[wwieder]

Sorry @slevis-lmwg where is the diagnostics linked above?

[slevis-lmwg]

In the "diags" section of the issue's introductory post