f.e23.FWmaHIST.ne16pg3_mg17_L135_tag160_beres0.85_PHscale0.25

[mijeong135]

Case Name f.e23.FWmaHIST.ne16pg3_mg17_L135_tag160_beres0.85_PHscale0.25 Keywords 2deg, ne16, WACCM, 135L, tag160, derecho, cam_dev, beres=0.85, parcel_hscale factor=0.25

---

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

sensitivity test of varying parcel hdepth (zmconvr.F90) with help from Adam H.

modify the "parcel_hscale" in zmconvr.F90 (see, here)

real(kind_phys), parameter :: parcel_hscale = 0.25_kind_phys (0.5 --> 0.25)

!+++arh !!!! parcel_dz(i) = max(zi(i,pver),60._kind_phys) ! do not use this--it hardwires parcel_dz parcel_dz(i) = max(zi(i,pver),parcel_hscale pbl_dz(i)) ! PBL mixing depth [parcel_hscaleBoundary, but no thinner than zi(i,pver)]

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 clone f.e23.FWmaHIST.ne16pg3_mg17_L135_tag160_beres0.85_ZM2andHB

---

ssh derecho

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

./create_clone --case /glade/u/home/mijeong/cesm2/cases/f.e23.FWmaHIST.ne16pg3_mg17_L135_tag160_beres0.85_PHscale0.25 --clone /glade/u/home/sglanvil/cases/f.e23.FWmaHIST.ne16pg3_mg17_L135_tag160_beres0.85_ZM2andHB --cime-output-root /glade/derecho/scratch/mijeong/

cd /glade/u/home/mijeong/cesm2/cases/f.e23.FWmaHIST.ne16pg3_mg17_L135_tag160_beres0.85_PHscale0.25

./case.setup

Compare SourceMods, env_run.xml, env_mach_pes.xml, and env_build.xml to original case.

;modify the "parcel_hscale" in zmconvr.F90 (see, here)

real(kind_phys), parameter :: parcel_hscale = 0.25_kind_phys

!+++arh !!!! parcel_dz(i) = max(zi(i,pver),60._kind_phys) ! do not use this--it hardwires parcel_dz parcel_dz(i) = max(zi(i,pver),parcel_hscalepbl_dz(i)) ! PBL mixing depth [parcel_hscaleBoundary, but no thinner than zi(i,pver)]

---

;edit user_nl_cam 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 gw_qbo_hdepth_scaling = 0.25D0 effgw_cm = 1.0D0 effgw_rdg_beta = 1.2D0 effgw_rdg_beta_max = 1.2D0 gw_top_taper = .false. frontgfc = 1.75D-15 taubgnd = 1.25D-3 nlte_limit_co2=.true.

se_hypervis_subcycle_sponge = 2 se_nu_top = 1E6 se_hypervis_subcycle = 9 se_rsplit = 2 se_nsplit = 2 se_large_courant_incr = .true. se_nu_div = 6.0E15 se_nu = 6.0E15 se_sponge_del4_lev = 20 se_sponge_del4_nu_div_fac = 8

---

./xmlchange DOUT_S_ROOT='/glade/campaign/acom/acom-climate/mijeong/archive/\$CASE' ./xmlchange DOUT_S_SAVE_INTERIM_RESTART_FILES=FALSE ./xmlchange DOUT_S=TRUE ./xmlchange STOP_OPTION=nmonths ./xmlchange STOP_N=12 ./xmlchange RESUBMIT=9 ./xmlchange CONTINUE_RUN=FALSE

./case.build --clean-all ./case.build

;check clocks ./xmlchange JOB_WALLCLOCK_TIME=12:00:00 --subgroup case.run ./xmlchange JOB_WALLCLOCK_TIME=04:00:00 --subgroup case.st_archive

./case.submit