f.e23.FWmaHIST.ne16pg3_mg17_L135_tag160_beres0.65_YL_CF16

[mijeong135]

Case Name
f.e23.FWmaHIST.ne16pg3_mg17_L135_tag160_beres0.65_YL_CF16
Keywords
2deg, ne16, WACCM, 135L, tag160, derecho, cam_dev, gw_convect.F90 (SourceMods), CF (20 --> 16), fixed bug in vertical_diffusion.F90 (Adding a bug fix on the vertical sponge layer diffusion Adam found (ESCOMP/CAM#1182))

---

Case Dir
/glade/u/home/mijeong/cesm2/cases/$CASE
Run Dir
/glade/derecho/scratch/mijeong/$CASE/run
Archive Dir
/glade/campaign/acom/acom-climate/mijeong/archive/$CASE

---

Tag
tag160 (Sasha's run)
Apply HB mod (HB Diffusion on but with Modified Mixing (NCAR/amwg_dev#537))
add Yuanpu's SourceMods (gw_convect.F90)

Fix a bug in vertical_diffusion.F90
In src/physics/cam/vertical_diffusion.F90 you just need to change the hard coded 1 in the kvm arrays to k:
!!!kvm(:ncol,1) = kvm(:ncol,1)+kvm_sponge(k) !+++bug!
kvm(:ncol,k) = kvm(:ncol,k)+kvm_sponge(k)

Changing effgw_beres_dp = 0.65
Decrease conversion factor (CF) from 20 (default) to 16

---

ssh derecho

cd /glade/work/hannay/cesm_tags/cam6_3_160/cime/scripts

./create_newcase --compset FWmaHIST --res ne16pg3_ne16pg3_mg17 --case /glade/u/home/mijeong/cesm2/cases/f.e23.FWmaHIST.ne16pg3_mg17_L135_tag160_beres0.65_YL_CF16 --run-unsupported --project P93300607

cd /glade/u/home/mijeong/cesm2/cases/f.e23.FWmaHIST.ne16pg3_mg17_L135_tag160_beres0.65_YL_CF16

./case.setup

Edit env_build.xml
./xmlchange CAM_CONFIG_OPTS="-phys cam_dev -nlev 135 -age_of_air_trcs -chem waccm_ma_mam5 -cppdefs -Dwaccm_debug -pcols 9"

Copy user_nl_clm:
cp /glade/u/home/mijeong/cesm2/cases/f.e23.FWmaHIST.ne16pg3_mg17_L135_tag132_cam6_beres0.65/user_nl_clm .

cp /glade/u/home/pel/src/for-unstable-Ri-hb/pbl_utils.F90 ./SourceMods/src.cam
cp ~/sourcemods/gw_convect.F90 ./SourceMods/src.cam

cp /glade/work/hannay/cesm_tags/cam6_3_160/src/physics/cam/vertical_diffusion.F90 ./SourceMods/src.cam
!!!kvm(:ncol,1) = kvm(:ncol,1)+kvm_sponge(k) !+++bug!
kvm(:ncol,k) = kvm(:ncol,k)+kvm_sponge(k)

edit gw_convect.F90 (line 115)
real(r8), parameter :: CF = 20._r8 --> CF = 16._r8

Edit/verify namelist settings (user_nl_cam):

ncdata = '/glade/derecho/scratch/pel/acom_ne16pg3_ne16pg3_mg17_1536_long2/run/acom_ne16pg3_ne16pg3_mg17_1536_long2.cam.i.1980-01-01-00000.nc'
bnd_topo = '/glade/derecho/scratch/pel/ne16pg3_gmted2010_modis_bedmachine_nc3000_Laplace0200_20230202.nc'
drydep_srf_file='/glade/u/home/nadavis/input_data/atmsrf_ne16pg3_110920.nc'
flbc_file = '/glade/campaign/cesm/cesmdata/inputdata/atm/waccm/lb/LBC_17500116-20150116_CMIP6_0p5degLat_c180227.nc'
flbc_list = 'CCL4', 'CF2CLBR', 'CF3BR', 'CFC11', 'CFC113', 'CFC12', 'CH3BR', 'CH3CCL3', 'CH3CL', 'CH4',
'CO2', 'H2', 'HCFC22', 'N2O', 'CFC114', 'CFC115', 'HCFC141B', 'HCFC142B', 'H2402', 'OCS', 'SF6', 'CFC11eq'

interpolate_output = .true.,.true.,.true.,.true.
empty_htapes=.true.

fincl1 = 'U','V','Q','T','PSL','OMEGA','PS','PRECT','OMEGA500','OMEGA850',
'FREQZM','CAPE','CMFMC_DP','ZMDT','PRECC','PRECL'
fincl2 = 'U','V','Q','T','PSL','OMEGA','PS','PRECT','OMEGA500','OMEGA850',
'UTGWORO','VTGWORO','UTGWSPEC','BUTGWSPEC','O3','CO2','N2O','H2O'
fincl3 = 'ZMDT','DTCOND','NETDT', 'HDEPTH', 'MAXQ0', 'PRECT', 'FLUT'
fincl8 = 'Uzm','Vzm','Wzm','THzm', 'VTHzm','WTHzm','UVzm','UWzm'

mfilt = 0, 5, 20, 30
avgflag_pertape = 'A','A','I','A'
nhtfrq = 0,-24,-6,-24

phys_grid_ctem_nfreq = -6
phys_grid_ctem_zm_nbas = 120
phys_grid_ctem_za_nlat = 90

se_statefreq = 144
effgw_beres_dp = 0.65D0
se_nsplit=3

---

Edit/verify namelist settings:
./xmlchange RUN_REFDATE=1980-01-01
./xmlchange RUN_STARTDATE=1980-01-01
./xmlchange STOP_OPTION=nmonths
./xmlchange STOP_N=12
./xmlchange CONTINUE_RUN=FALSE
./xmlchange DOUT_S_ROOT='/glade/campaign/acom/acom-climate/mijeong/archive/$CASE'
./xmlchange RESUBMIT=9

Check env_mach_pes.xml for settings
cp /glade/u/home/mijeong/cesm2/cases/f.e23.FWmaHIST.ne16pg3_mg17_L135_tag132_cam6_beres0.65/env_mach_pes.xml .

Check env_workflow.xml for clock settings
./xmlchange JOB_WALLCLOCK_TIME=12:00:00 --subgroup case.run

./case.build --clean-all
./case.build
./preview_namelists
./preview_run
./case.submit