Case Name
f.e23.FWmaHIST.ne16pg3_mg17_L135_tag160_beres0.40_YL_CF12
Keywords
2deg, ne16, WACCM, 135L, tag160, derecho, cam_dev, gw_convect.F90 (SourceMods), CF (20 --> 12), fixed bug in vertical_diffusion.F90 (Adding a bug fix on the vertical sponge layer diffusion Adam found (https://github.com/ESCOMP/CAM/issues/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.40
Decrease conversion factor (CF) from 20 (default) to 12
ssh derecho
cd /glade/work/hannay/cesm_tags/cam6_3_160/cime/scripts
Case Name f.e23.FWmaHIST.ne16pg3_mg17_L135_tag160_beres0.40_YL_CF12 Keywords 2deg, ne16, WACCM, 135L, tag160, derecho, cam_dev, gw_convect.F90 (SourceMods), CF (20 --> 12), fixed bug in vertical_diffusion.F90 (Adding a bug fix on the vertical sponge layer diffusion Adam found (https://github.com/ESCOMP/CAM/issues/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.40 Decrease conversion factor (CF) from 20 (default) to 12
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.40_YL_CF12 --run-unsupported --project P93300607
cd /glade/u/home/mijeong/cesm2/cases/f.e23.FWmaHIST.ne16pg3_mg17_L135_tag160_beres0.40_YL_CF12
./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 = 12._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.40D0 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