Case Name
f.e23.FWmaHIST.ne16pg3_mg17_L135_tag160_beres0.85_gw_HYBRID
Keywords
2deg, ne16, WACCM, 135L, tag160, derecho, cam_dev, beres=0.85, CLUBB + ZM heating
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
;key points
using CLUBB & ZM heating to drive the convective gravity wave scheme with help from Adam H.
! set hdepth to max of clubb and zm
hdepth = MAX(hdepth,hdepth_clubb)
! ...and adjust indices
do i = 1, ncol
topi(i) = MIN(topi(i),topic(i))
boti(i) = pver
end do
! set to clubb heating rate
!q0 = q0_clubb
! set to clubb + zm heating rate
q0 = q0 + q0_clubb
!limit heating to 25 K/day
*q0 = MIN(25._r8/(24._r83600._r8),q0)**
add PRECT & FLUT to the 6-hourly output (h2)
also add FREQZM, CAPE, CMFMC_DP, ZMDT, PRECC, PRECL to the monthly output (h0)
use tag 160 (2 deg, 135L)
beres = 0.85
make sure HB modification is in effect
HB Diffusion on but with Modified Mixing (https://github.com/NCAR/amwg_dev/issues/537)
cp /glade/u/home/pel/src/for-unstable-Ri-hb/pbl_utils.F90 to SourceMods/src.cam
ssh derecho
cd /glade/work/hannay/cesm_tags/cam6_3_160/cime/scripts
Case Name f.e23.FWmaHIST.ne16pg3_mg17_L135_tag160_beres0.85_gw_HYBRID Keywords 2deg, ne16, WACCM, 135L, tag160, derecho, cam_dev, beres=0.85, CLUBB + ZM heating
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
;key points using CLUBB & ZM heating to drive the convective gravity wave scheme with help from Adam H.
;copy cd SourceMods/src.cam cp /glade/derecho/scratch/aherring/for_mijeoug/gw_drag.F90 . cp /glade/derecho/scratch/aherring/for_mijeoug/clubb_intr.F90 . cp /glade/derecho/scratch/aherring/for_mijeoug/gw_convect.F90 .
;Settings in gw_convect.F90
! set hdepth to clubb !hdepth = hdepth_clubb
! set hdepth to max of clubb and zm hdepth = MAX(hdepth,hdepth_clubb)
! ...and adjust indices do i = 1, ncol topi(i) = MIN(topi(i),topic(i)) boti(i) = pver end do
! set to clubb heating rate !q0 = q0_clubb
! set to clubb + zm heating rate q0 = q0 + q0_clubb
!limit heating to 25 K/day *q0 = MIN(25._r8/(24._r83600._r8),q0)**
add PRECT & FLUT to the 6-hourly output (h2) also add FREQZM, CAPE, CMFMC_DP, ZMDT, PRECC, PRECL to the monthly output (h0)
use tag 160 (2 deg, 135L) beres = 0.85 make sure HB modification is in effect HB Diffusion on but with Modified Mixing (https://github.com/NCAR/amwg_dev/issues/537) cp /glade/u/home/pel/src/for-unstable-Ri-hb/pbl_utils.F90 to SourceMods/src.cam
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.85_gw_HYBRID --run-unsupported --project P93300607
cd /glade/u/home/mijeong/cesm2/cases/f.e23.FWmaHIST.ne16pg3_mg17_L135_tag160_beres0.85_gw_HYBRID
./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 from Simone) from here: cp /glade/u/home/sglanvil/cases/f.e23.FWmaHIST.ne16pg3_mg17_L135_tag160_beres0.65/user_nl_clm .
Compare SourceMods, env_run.xml, env_mach_pes.xml, and env_build.xml to original case.
;edit 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 = '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 = -24 phys_grid_ctem_za_nlat = 90 phys_grid_ctem_zm_nbas = 90
se_statefreq = 144 effgw_beres_dp = 0.85D0
Edit/verify namelist settings (compare to Mijeong): ./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 clocks ./xmlchange JOB_WALLCLOCK_TIME=12:00:00 --subgroup case.run ./xmlchange JOB_WALLCLOCK_TIME=04:00:00 --subgroup case.st_archive
./xmlchange NTASKS=-8 ./case.setup --reset ./case.build --clean-all ./case.build ./case.submit