f.e30.FWscHIST.ne30pg3_mg17_L135_cam6_4_038_beres0.70

[mijeong135]

Case Name
f.e30.FWscHIST.ne30pg3_mg17_L135_cam6_4_038_beres0.70
Keywords
1deg, ne30, WACCM, 135L, cam6_4_038, derecho, cam7

---

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

---

Tag cam6_4_038

cd /glade/u/home/mijeong/cesm_tags
git clone https://github.com/ESCOMP/CAM.git cam6_4_038
cd cam6_4_038
git checkout cam6_4_038
./bin/git-fleximod update

Set up a WACCM run (135L), Tag: cam6_4_038
use specified chemistry (FWscHIST)

---

ssh derecho

cd /glade/u/home/mijeong/cesm_tags/cam6_4_038/cime/scripts

./create_newcase --compset FWscHIST --res ne30pg3_ne30pg3_mg17 --case /glade/u/home/mijeong/cesm3/cases/f.e30.FWscHIST.ne30pg3_mg17_L135_cam6_4_038_beres0.70 --run-unsupported --project P93300607

cd /glade/u/home/mijeong/cesm3/cases/f.e30.FWscHIST.ne30pg3_mg17_L135_cam6_4_038_beres0.70

./case.setup

Edit env_build.xml:
./xmlchange CAM_CONFIG_OPTS="-phys cam7 -nlev 135 -age_of_air_trcs -chem waccm_sc_mam4 -cppdefs -Dwaccm_debug -pcols 9"

Turn on RRTMGP
./xmlchange --append CAM_CONFIG_OPTS="-rad rrtmgp"

---

Edit user_nl_cam:

ncdata='/glade/campaign/cgd/amp/pel/inic/L135_ne30pg3_ne30pg3_mg17.cam.i.1979-01-11-00000.nc'
bnd_topo= '/glade/work/juliob/bndtopo/ne30pg3_gmted2010_modis_bedmachine_nc3000_Laplace0100_noleak_20240720.nc'

interpolate_output = .true.,.true.,.true.,.true.
interpolate_nlat = 192,192,192,192
interpolate_nlon = 288,288,288,288

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','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 = -6
phys_grid_ctem_za_nlat = 90
phys_grid_ctem_zm_nbas = 120

inithist = 'YEARLY'
se_statefreq = 144
effgw_beres_dp = 0.70D0
se_nsplit=4

seasalt_emis_scale = 0.75D0
clubb_c8 = 4.95

---

XMLchanges:
./xmlchange RUN_REFDATE=1980-01-01
./xmlchange RUN_STARTDATE=1980-01-01
./xmlchange STOP_OPTION=nmonths
./xmlchange STOP_N=6
./xmlchange CONTINUE_RUN=FALSE
./xmlchange RESUBMIT=2
./xmlchange DOUT_S_SAVE_INTERIM_RESTART_FILES=TRUE
./xmlchange DOUT_S_ROOT=/glade/campaign/acom/acom-climate/mijeong/archive/$CASE
./xmlchange DOUT_S=TRUE

;env_mach_pes.xml (ntasks=-20 ; 20*128=2160)
./xmlchange NTASKS=-20

./case.build

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

./preview_namelists
./preview_run
./case.submit

[PeterHjortLauritzen]

If you want to match Cecile's setup (see NCAR/amwg_dev#583) a couple of things need to be added:

./xmlchange --append CAM_CONFIG_OPTS="-rad rrtmgp"

and

 seasalt_emis_scale             = 0.75D0
clubb_c8		=  4.95 

to user_nl_cam.

I am not sure if this is needed in WACCM (in user_nl_cam):

 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'

Also, Cecile does not have:

dust_emis_fact = 1.3

Maybe @tilmes can comment on that?

[tilmes]

Regarding the dust factor, I would do what is used in Cecile's run for now.

[PeterHjortLauritzen]

Here are code mods for an increased del4 sponge:

/glade/u/home/pel/src/for-mijeong2/global_norms_mod.F90

Added code to increase del4 in top 5 layers:

If you are using se_nsplit=4 then you need to set

se_hypervis_subcycle=17

If you are trying se_nsplit=3 then

se_hypervis_subcycle=23

(=17*4/3)

[mijeong135]

Adding this to my run:

cp /glade/u/home/pel/src/for-mijeong2/global_norms_mod.F90 ./SourceMods/src.cam

!In user_nl_cam:

se_nsplit=4
se_hypervis_subcycle=17

[mijeong135]

Hi Peter, Francis said remove the 'rad_climate =...' part in user_nl_cam. In WACCM, CO2 is not advected and without including this, the model will use the default setting. We will see if this helps.

…

-Mijeong

On Tue, Oct 15, 2024 at 1:44 PM Peter Hjort Lauritzen < ***@***.***> wrote: Here are code mods for an increased del4 sponge: /glade/u/home/pel/src/for-mijeong2/global_norms_mod.F90 Added code to increase del4 in top 5 layers: Screenshot.2024-10-15.at.1.41.57.PM.png (view on web) <https://github.com/user-attachments/assets/b062926b-e0b7-47be-9800-efe3b26c9064> If you are using se_nsplit=4 then you need to set se_hypervis_subcycle=17 If you are trying se_nsplit=3 then se_hypervis_subcycle=23 (=17*4/3) — Reply to this email directly, view it on GitHub <#63 (comment)>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/BCCRGOLMA6LEVOSDVDPVHCDZ3VWCTAVCNFSM6AAAAABPTIVQI6VHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDIMJUHA3DMMZSGM> . You are receiving this because you authored the thread.Message ID: ***@***.***>

[mijeong135]

Adding a bug fix on the vertical sponge layer diffusion Adam found (ESCOMP/CAM#1182)

download vertical_diffusion.F90 (here)

cp vertical_diffusion.F90 ./SourceMods/src.cam
change the line explained below

Below is an email from Adam
;----------
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)

This pertains to line 1106 in the most recent cam tag, here. You can either change your source tree vertical_diffusion.F90, or you can copy the original into SourceMods/src.cam/ of an individual case, and implement the bug fix for an individual case.