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

Description: Investigating dry bias in Amazon seen in cesm_dev 116 (https://github.com/NCAR/cesm_dev/issues/19) as discussed in #75.

Try replicating Adam's case https://github.com/NCAR/amwg_dev/issues/614 but as F2000BgcCrop.

[x] modern code base beta04 tag
[x] Leung2023 dust; tuning turned on and other CAM mods as in cesm_dev 116 (or similar); Adam H can advise.
[x] CTSM5.3 using CLM6-BGC
[x] initial conditions from #68
[x] Run 1 (this issue): CTSM5.3(PPE) parameter file vs.
Run 2 (open another issue for this): 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 may advise.
[ ] Submit 40-yr simulation: resolving issue with ncdata setting

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

usernl changes: CLM

hist_fincl1 = 'QICE_FORC','TSRF_FORC','TOPO_FORC'
hist_fincl2 = 'PCT_LANDUNIT','ATM_TOPO','TSA','LWup','LWdown',
             'ALBD','ALBI','ALBGRD','ALBGRI','SWdown','SWup',
             'RAIN_ICE','SNOW_ICE','QSNOMELT_ICE','QICE_MELT',
             'QSNOFRZ_ICE','QRUNOFF_ICE','QICE','QSOIL_ICE',
             'SNOW_DEPTH','SNOWDP','H2OSNO','FSNO','QICE_FRZ',
             'QSNOMELT','SOILWATER_10CM','TWS','QFLX_EVAP_TOT',
             'H2OSOI','TSOI','TSOI_10CM','TSOI_ICE','H2OSNO_ICE'
             'QRUNOFF','SNOW','RAIN','FSH','TREFMNAV','TREFMXAV'                 
hist_avgflag_pertape(1) = 'A'
hist_avgflag_pertape(2) = 'A'
hist_nhtfrq = 0,-24
hist_mfilt  = 1,365
hist_ndens  = 2,2

reseed_dead_plants = .true.

! turn off shifting cultivation
do_grossunrep = .false.

! do this for isotopes, not critical at this point, but in general I'd like these to be on for coupled model runs
use_c13 = .true.
use_c14 = .true.
use_c13_timeseries = .true.
use_c14_bombspike = .true.

CAM

! file h0 h1 h2 h3 h4 h5 h6 h7 h8 h9
mfilt = 0, 30, 0, 100, 30, 4, 80, 40, 4, 48
nhtfrq = 0, -24, 0, -24, -24, -6, -3, -3, -6, -1
ndens = 2, 2, 2, 2, 2, 2, 2, 2, 2, 2
interpolate_output = .true., .true., .false., .false., .false., .false., .false., .false., .false., .false.
interpolate_nlat = 192, 192, 192, 192, 192, 192, 192, 192, 192, 192
interpolate_nlon = 288, 288, 288, 288, 288, 288, 288, 288, 288, 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'

!output for MDTF (h1)
fincl2='PRECT', 'PRECC', 'FLUT', 'U850', 'U200', 'V850', 'V200', 'OMEGA500', 'TS', 'SST', 'PSL', 'Z500'

! native grid
!monthly (h2)
fincl3='ACTNI', 'ACTNL', 'ACTREI', 'ACTREL', 'AODDUST', 'AODVIS', 'AODVISdn',
'AODDUSTdn','AODVISdn','AODBCdn','AODSOAdn','AODSSdn','AODPOMdn','AODSO4dn',
'BURDENBC', 'BURDENDUST', 'BURDENPOM', 'BURDENSEASALT',
'BURDENSO4', 'BURDENSOA', 'CAPE', 'CCN3', 'CDNUMC', 'CH4', 'CLDHGH', 'CLDICE', 'CLDLIQ', 'CLDLOW', 'CLDMED', 'CLDTOT', 'CLOUD', 'CMFMC_DP',
'CT_H2O', 'DCQ', 'DQCORE', 'DTCOND', 'DTCORE', 'DTV', 'EVAPPREC', 'EVAPSNOW', 'FCTI', 'FCTL', 'FICE', 'FLDS', 'FLNS', 'FLNSC', 'FLNT', 'FLNTC', 'FLUT',
'FREQZM', 'FSDS', 'FSDSC', 'FSNS', 'FSNSC', 'FSNT', 'FSNTC', 'FSNTOA', 'ICEFRAC', 'LANDFRAC', 'LHFLX', 'LWCF', 'MPDICE', 'MPDLIQ', 'MPDQ', 'MPDT',
'OCNFRAC', 'OMEGA', 'OMEGA500', 'PBLH', 'PHIS', 'PINT', 'PMID', 'PRECC', 'PRECL', 'PRECSC', 'PRECSL', 'PRECT', 'PS', 'PSL', 'PTEQ', 'PTTEND', 'Q',
'QFLX', 'QRL', 'QRS', 'QTGW', 'RCMTEND_CLUBB', 'RELHUM', 'RVMTEND_CLUBB', 'SHFLX', 'SOLIN', 'SST', 'STEND_CLUBB', 'SWCF',
'T', 'TAUX', 'TAUY', 'TFIX', 'TGCLDIWP', 'TGCLDLWP', 'TMQ', 'TREFHT', 'TS', 'TTGW', 'U', 'U10', 'UBOT', 'UTGWORO', 'UTGW_TOTAL',
'V', 'VBOT', 'VTGWORO', 'VTGW_TOTAL', 'WPRTP_CLUBB', 'WPTHLP_CLUBB', 'Z3', 'ZMDQ', 'ZMDT', 'N2O', 'CO2','CFC11','CFC12',
'AODVISdn','CCN3', 'CDNUMC', 'H2O', 'NUMICE', 'NUMLIQ','OMEGA500',
'TROP_P','TROP_T','TROP_Z',
'num_a1','num_a2','num_a3','num_a4',
'dgnd_a01','dgnd_a02','dgnd_a03','dgnd_a04',
'ncl_a1','ncl_a2','ncl_a3',
'SFncl_a1', 'SFncl_a2', 'SFncl_a3'
'ncl_a1', 'ncl_a1DDF', 'ncl_a1SFWET',
'ncl_a2', 'ncl_a2DDF', 'ncl_a2SFWET',
'ncl_a3', 'ncl_a3DDF', 'ncl_a3SFWET',
'ncl_c1', 'ncl_c1DDF', 'ncl_c1SFWET', 'ncl_c2', 'ncl_c2DDF', 'ncl_c2SFWET', 'ncl_c3', 'ncl_c3DDF', 'ncl_c3SFWET',
'SFbc_a4', 'bc_a1', 'bc_a1DDF', 'bc_a1SFWET', 'bc_a4', 'bc_a4_CLXF', 'bc_a4DDF', 'bc_a4SFWET',
'bc_c1', 'bc_c1DDF', 'bc_c1SFWET', 'bc_c4', 'bc_c4DDF', 'bc_c4SFWET',
'SFdst_a1','SFdst_a2', 'SFdst_a3'
'dst_a1', 'dst_a1DDF', 'dst_a1SFWET', 'dst_a2', 'dst_a2DDF',
'dst_a2SFWET', 'dst_a3', 'dst_a3DDF', 'dst_a3SFWET', 'dst_c1', 'dst_c1DDF', 'dst_c1SFWET', 'dst_c2',
'dst_c2DDF', 'dst_c2SFWET', 'dst_c3', 'dst_c3DDF', 'dst_c3SFWET'
'SFpom_a4', 'pom_a1', 'pom_a1DDF', 'pom_a1SFWET', 'pom_a4', 'pom_a4DDF', 'pom_a4SFWET',
'pom_c1','pom_c1DDF', 'pom_c1SFWET', 'pom_c4', 'pom_c4DDF', 'pom_c4SFWET'
'so4_a1', 'so4_a1_CHMP', 'so4_a1_CLXF', 'so4_a1DDF', 'so4_a1_sfgaex1',
'so4_a1SFWET', 'so4_c1', 'so4_c1AQH2SO4', 'so4_c1AQSO4', 'so4_c1DDF', 'so4_c1SFWET'
'so4_a2', 'so4_a2_CHMP', 'so4_a2_CLXF', 'so4_a2DDF', 'so4_a2_sfgaex1', 'so4_a2_sfnnuc1', 'so4_a2SFWET',
'so4_c2', 'so4_c2AQH2SO4', 'so4_c2AQSO4', 'so4_c2DDF', 'so4_c2SFWET'
'so4_a3', 'so4_a3DDF', 'so4_a3_sfgaex1', 'so4_a3SFWET', 'so4_c3', 'so4_c3AQH2SO4', 'so4_c3AQSO4', 'so4_c3DDF', 'so4_c3SFWET',
'soa_a1', 'soa_a1_CHML', 'soa_a1DDF', 'soa_a1_sfgaex1', 'soa_a1SFWET',
'soa_a2', 'soa_a2_CHML', 'soa_a2DDF', 'soa_a2_sfgaex1', 'soa_a2SFWET',
'soa_c1', 'soa_c1DDF', 'soa_c1SFWET', 'soa_c2', 'soa_c2DDF', 'soa_c2SFWET',
'CCN3','CCN4',
'BUTGWSPEC', 'BVTGWSPEC','BVTGWSPEC','TTGWORO','BTTGWSPEC','TAUGWX', 'TAUGWY','TAUBLJX','TAUBLJY', 'UTEND_PHYSTOT', 'UTEND_GWDTOT', 'UTEND_DCONV', 'UTEND_CLUBB', 'UTEND_CORE',
'VTEND_PHYSTOT', 'VTEND_GWDTOT', 'VTEND_DCONV', 'VTEND_CLUBB', 'VTEND_CORE'

!monthly (h3)
fincl4 = 'Uzm','Vzm','Wzm','THzm', 'VTHzm','WTHzm','UVzm','UWzm'
phys_grid_ctem_nfreq=-6
phys_grid_ctem_zm_nbas=120
!ne30 value
phys_grid_ctem_za_nlat=90

!daily - average (h4)
fincl5 = 'SHFLX', 'LHFLX', 'PBLH',
'PRECT', 'PRECC', 'FLUT', 'U850', 'U200', 'V850', 'V200', 'OMEGA500', 'TS', 'SST', 'PSL'

!6 hourly - instantaneous with input from Julio (h5)
fincl6 = 'PS:I','U:I','V:I','T:I','OMEGA:I', 'PMID:I', 'PSDRY:I'

!3 hourly - average (h6)
fincl7= 'PRECC:A', 'PRECL:A'

!3 hourly - instantaneous (h7)
fincl8= 'PRECC:I','PRECL:I',
'U850:I','V850:I','T500:I','T200:I','Z500:I','Z300:I','PSL:I','PS:I','U10:I','UBOT:I','VBOT:I',
'FLUT:I','TMQ:I','OMEGA500:I','uIVT:I', 'vIVT:I'

!6 hourly - instantaneous with input from Ben (h8)
fincl9 = 'um:I','upwp:I','up2_tp:I'
clubb_vars_zt='um'
clubb_vars_zm='upwp','up2_tp'
clubb_history=.true.

!1 hourly - instantaneous for ming (h9)
fincl10 = 'Z500:I','PSL:I','U10:I','PRECT:A','PRECC:A','uIVT:I','vIVT:I','U850:I','U200:I','V850:I','V200:I','Q850:I','T850:I','FLUT:I'

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/campaign/cgd/amp/juliob/ERA5/ne30np4/L93_cdf5/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

SourceMods:

Diagnostics: Diags (if available) Look at the last 20 yrs. https://webext.cgd.ucar.edu/I1850/$CASE/lnd/

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

Contacts: @slevis-lmwg

Extra details: Regarding user_nl_cam:

Due to lack of permissions, I had to change to this path for ncdata: /glade/derecho/scratch/slevis/f.e30_alpha04a.FMT2000BgcCrop.CTRL.lmwg_dev_76/run/ERA5_x_ne30np4_L93_rgC2_WO_aero.2000-01-01-00000.nc

Due to this error, I had to remove these variables from the fincl statements:

FLDLST: Z850 in fincl(7, 8) not found
FLDLST: UPWP_UG_CLUBB in fincl(4, 9) not found
FLDLST: UPWP_UG_CT_CLUBB in fincl(5, 9) not found
FLDLST: DUDZ_CLUBB in fincl(6, 9) not found
FLDLST: UPWP_UG_MEAN_LEV in fincl(7, 9) not found
FLDLST: UPWP_UG_MIN_LEV in fincl(8, 9) not found
FLDLST: UPWP_UG_MAX_LEV in fincl(9, 9) not found
FLDLST: UPWP_UG_MEAN_HGT in fincl(10, 9) not found
FLDLST: Z250 in fincl(4, 10) not found

NCAR / LMWG_dev

F2000BGCCROP ctsm53 control (f.e30_alpha04a.FMT2000BgcCrop.CTRL.lmwg_dev_76) #76