f.e30_cam6_4_036.FMTHIST.ne120_L93.001

[cecilehannay]

Purpose:

cam7 high resolution run: ne120 L93

---

Description:

• pe_layout

./xmlchange NTASKS=6144

• rrtmgp

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

• start in 2000

./xmlchange RUN_STARTDATE=2000-01-01

• user_nl_cam

! file                      h0      h1     h2       h3          h4        h5       h6       h7        h8      h9 
mfilt		        =       0,      30,    0,        100,       30,       4,       80,      40,        4,     48
nhtfrq              =       0,     -24,    0,        -24,       -24,      -6,      -3,      -3,       -6,     -1
ndens               =       2,       2,     2,       2,          2,        2,       2,       2,        2,      2
interpolate_output  =  .true.,  .true., .false., .false.,   .false.,     .false.,  .false., .false., .false.,  .false.
interpolate_nlat    =     768,      768,    192,     192,      192,       192,     192,      192,     192,   192
interpolate_nlon    =     1152,    1152,    288,     288,      288,       288,      288,     288,     288,   288

empty_htapes = .true.

! interpolated output
! output for ADF  monthly-interpolated (h0)    
fincl1 = 'AODDUST', 'AODVIS', 'AODDUSTdn','AODVISdn','BURDENBC', 'BURDENDUST', 'BURDENPOM', 'BURDENSEASALT', 
'BURDENSO4', 'BURDENSOA', 'CDNUMC', 'CLDHGH', 'CLDICE', 'CLDLIQ', 'CLDLOW', 'CLDMED', 'CLDTOT', 'CLOUD', 'FLDS', 'FLNS', 'FLNSC', 'FLNT', 'FLNTC', 'FLUT', 
'FSNS', 'FSNSC', 'FSNT', 'FSNTC', 'FSNTOA', 'ICEFRAC', 'LHFLX', 'LWCF','OMEGA500', 'PBLH','PRECC', 'PRECL', 'PRECSC', 'PRECSL', 'PRECT', 'PS', 'PSL','Q', 
'QFLX', 'RELHUM', 'SHFLX', 'SOLIN', 'SST', 'SWCF', 'T', 'TAUX', 'TAUY', 'TGCLDIWP', 'TGCLDLWP', 'TMQ', 'TREFHT', 'TS', 'U', 'U10', 'V'

!output for MDTF (h1)
fincl2='PRECT', 'PRECC', 'FLUT', 'U850', 'U200', 'V850', 'V200', 'OMEGA500', 'TS', 'SST', 'PSL'

! native grid
!monthly (h2)
fincl3='ACTNI', 'ACTNL', 'ACTREI', 'ACTREL', 'AODDUST', 'AODVIS', 'AODVISdn','BURDENBC', 'BURDENDUST', 'BURDENPOM', 'BURDENSEASALT', 
'BURDENSO4', 'BURDENSOA', 'CAPE', 'CCN3', 'CDNUMC', 'CH4', 'CLDHGH', 'CLDICE', 'CLDLIQ', 'CLDLOW', 'CLDMED', 'CLDTOT', 'CLOUD', 'CMFMC_DP', 
'CT_H2O', 'DCQ', 'DQCORE', 'DTCOND', 'DTCORE', 'DTV', 'EVAPPREC', 'EVAPSNOW', 'FCTI', 'FCTL', 'FICE', 'FLDS', 'FLNS', 'FLNSC', 'FLNT', 'FLNTC', 'FLUT', 
'FREQZM', 'FSDS', 'FSDSC', 'FSNS', 'FSNSC', 'FSNT', 'FSNTC', 'FSNTOA', 'ICEFRAC', 'LANDFRAC', 'LHFLX', 'LWCF', 'MPDICE', 'MPDLIQ', 'MPDQ', 'MPDT', 
'OCNFRAC', 'OMEGA', 'OMEGA500', 'PBLH', 'PHIS', 'PINT', 'PMID', 'PRECC', 'PRECL', 'PRECSC', 'PRECSL', 'PRECT', 'PS', 'PSL', 'PTEQ', 'PTTEND', 'Q', 
'QFLX', 'QRL', 'QRS', 'QTGW', 'RCMTEND_CLUBB', 'RELHUM', 'RVMTEND_CLUBB', 'SHFLX', 'SOLIN', 'SST', 'STEND_CLUBB', 'SWCF', 
'T', 'TAUX', 'TAUY', 'TFIX', 'TGCLDIWP', 'TGCLDLWP', 'TMQ', 'TREFHT', 'TS', 'TTGW', 'U', 'U10', 'UBOT', 'UTGWORO', 'UTGW_TOTAL', 
'V', 'VBOT', 'VTGWORO', 'VTGW_TOTAL', 'WPRTP_CLUBB', 'WPTHLP_CLUBB', 'Z3', 'ZMDQ', 'ZMDT', 'N2O', 'CO2','CFC11','CFC12',
'AODVISdn','CCN3', 'CDNUMC', 'H2O', 'NUMICE', 'NUMLIQ','OMEGA500',
'TROP_P','TROP_T','TROP_Z',
'num_a1','num_a2','num_a3','num_a4',
'dgnd_a01','dgnd_a02','dgnd_a03','dgnd_a04',
'ncl_a1','ncl_a2','ncl_a3',
'bc_a1','bc_a4',
'dst_a1','dst_a3',
'pom_a1','pom_a4',
'so4_a1','so4_a2','so4_a3',
'soa_a1','soa_a2',
'CCN3','CCN4',
'BUTGWSPEC',  'BVTGWSPEC','BVTGWSPEC','TTGWORO','BTTGWSPEC','TAUGWX', 'TAUGWY','TAUBLJX','TAUBLJY', 'UTEND_PHYSTOT', 'UTEND_GWDTOT', 'UTEND_DCONV', 'UTEND_CLUBB', 'UTEND_CORE',
'VTEND_PHYSTOT', 'VTEND_GWDTOT', 'VTEND_DCONV', 'VTEND_CLUBB', 'VTEND_CORE'

!monthly  (h3)
fincl4 = 'Uzm','Vzm','Wzm','THzm', 'VTHzm','WTHzm','UVzm','UWzm'    
phys_grid_ctem_nfreq=-6
phys_grid_ctem_zm_nbas=120
phys_grid_ctem_za_nlat=180

!daily - average (h4)
fincl5 = 'SHFLX', 'LHFLX', 'PBLH',      
'PRECT', 'PRECC', 'FLUT', 'U850', 'U200', 'V850', 'V200', 'OMEGA500', 'TS', 'SST', 'PSL'

!6 hourly - instantaneous with input from Julio (h5)
fincl6 = 'PS:I','U:I','V:I','T:I','OMEGA:I', 'PMID:I', 'PSDRY:I'

!3 hourly - average (h6)
fincl7= 'PRECC:A', 'PRECL:A'	

!3 hourly - instantaneous (h7)
fincl8= 'PRECC:I','PRECL:I',
'U850:I','V850:I','T500:I','T200:I','Z850:I','Z500:I','Z300:I','PSL:I','PS:I','U10:I','UBOT:I','VBOT:I',
'FLUT:I','TMQ:I','OMEGA500:I','uIVT:I', 'vIVT:I' 

!6 hourly - instantaneous with input from Ben (h8)
fincl9 = 'um:I','upwp:I','up2_tp:I','UPWP_UG_CLUBB:I','UPWP_UG_CT_CLUBB:I','DUDZ_CLUBB:I','UPWP_UG_MEAN_LEV:I','UPWP_UG_MIN_LEV:I','UPWP_UG_MAX_LEV:I','UPWP_UG_MEAN_HGT:I'
clubb_vars_zt='um'
clubb_vars_zm='upwp','up2_tp'
clubb_history=.true.

!1 hourly - instantaneous for ming (h9)
fincl10 = 'Z500:I','PSL:I','U10:I','Z250:I','PRECT:A','PRECC:A','uIVT:I','vIVT:I','U850:I','U200:I','V850:I','V200:I','Q850:I','T850:I','FLUT:I'

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'
seasalt_emis_scale             = 0.75D0

seasalt_emis_scale             = 0.75D0

clubb_c8		=  4.95 

ubc_file_path=  "/glade/campaign/cesm/cesmdata/inputdata/atm/cam/chem/ubc/b.e21.BWHIST.f09_g17.CMIP6-historical-WACCM.ensAvg123.cam.h0zm.H2O.1849-2014_c240604.nc"

se_statefreq=576  

ncdata = '/glade/campaign/cgd/amp/aherring/restart/inic/ne120/ERA5_x_ne120np4_L93_rgC2_WO.2000-01-01-00000.nc'
bnd_topo =
'/glade/campaign/cgd/amp/aherring/restart/TopoCESM2/cube_to_target/output/ne120pg3_gmted2010_modis_bedmachine_nc3000_Laplace0025_noleak_20240913.nc'

use_gw_front        = .true.
taubgnd        = 1.5D-6

• user_nl_clm

hist_fincl1 = 'QICE_FORC','TSRF_FORC','TOPO_FORC'                      
hist_fincl2 = 'PCT_LANDUNIT','ATM_TOPO','TSA','LWup','LWdown',         
             'ALBD','ALBI','ALBGRD','ALBGRI','SWdown','SWup',          
             'RAIN_ICE','SNOW_ICE','QSNOMELT_ICE','QICE_MELT',         
             'QSNOFRZ_ICE','QRUNOFF_ICE','QICE','QSOIL_ICE',           
             'SNOW_DEPTH','SNOWDP','H2OSNO','FSNO','QICE_FRZ', 
             'QSNOMELT','SOILWATER_10CM','TWS','QFLX_EVAP_TOT',
             'H2OSOI','TSOI','TSOI_10CM','TSOI_ICE','H2OSNO_ICE' 
             'QRUNOFF','SNOW','RAIN','FSH','TREFMNAV','TREFMXAV'                             
hist_avgflag_pertape(1) = 'A'                                           
hist_avgflag_pertape(2) = 'A'                                           
hist_nhtfrq = 0,-24                                                     
hist_mfilt  = 1,365                                                     
hist_ndens  = 2,2  

flanduse_timeseries = '/glade/derecho/scratch/aherring/restart/ctsm5.2.028/tools/mksurfdata_esmf/landuse.timeseries_ne120pg3_SSP3-7.0_1979-2026_78pfts_c240925.nc' 
fsurdat = '/glade/derecho/scratch/aherring/restart/ctsm5.2.028/tools/mksurfdata_esmf/surfdata_ne120pg3_SSP3-7.0_1979_78pfts_c240925.nc'
finidat = '/glade/campaign/cesm/cesmdata/inputdata/lnd/clm2/initdata_map/clmi.F2000.2000-01-01.ne120pg3_mt13_simyr2000_c200728.nc'
use_init_interp = .true.
use_excess_ice = .false.

• SourceMods
  add Z850, Ben's output

/glade/campaign/cesm/cesmdata/cseg/runs/cesm2_0/f.e30_cam6_4_036.FMTHIST.ne120_L93.001/SourceMods/src.cam/

---

Case directory:

• Locally (if still available):
  /glade/campaign/cesm/cesmdata/cseg/runs/cesm2_0/f.e30_cam6_4_036.FMTHIST.ne120_L93.001

• On github:
  https://github.com/NCAR/amwg_dev/tree/f.e30_cam6_4_036.FMTHIST.ne120_L93.001

---

Sandbox:

• Locally (if still available):
  /glade/work/hannay/cesm_tags/cam6_4_036

• On github:
  https://github.com/ESCOMP/CAM/tree/cam6_4_036

---

Diagnostics:

• ADF diags (if available)
  https://webext.cgd.ucar.edu/FMTHIST/f.e30_cam6_4_036.FMTHIST.ne120_L93.001/atm/

---

Output:

• Raw output on derecho (if still available)
  /glade/derecho/scratch/hannay/archive/f.e30_cam6_4_036.FMTHIST.ne120_L93.001/
• CAM Climos: /glade/campaign/cgd/amp/amwg/climo/f.e30_cam6_4_036.FMTHIST.ne120_L93.001

---

Contacts:
@JulioTBacmeister, @adamrher, @islasimpson, @swrneale, @PeterHjortLauritzen, @cmcclus, @bstephens82, @megandevlan, @jedwards4b, @billsack, @wwieder, @ekluzek, @shieldsca, @slevis-lmwg

---

Extra details:

[cecilehannay]

From @JulioTBacmeister

First JFM and JJA zonal mean U and T

Memory of the initial conditions here seems remarkable, i.e., our zonal means for the 3-month JFM 2000 look a lot like the ERA JFM 2000 mean.

JJA U shows clearly that summertime stratospheric jet is not closing off. This is also evident in JFM. This is almost certainly caused by the fact that we have shut off the frontal gravity waves. Not too worried about it. More worried about the warming that seems to be taking place at the top of the model. Maybe it's just seasonal, but are we sure have turned on radiation above 1Pa?

---

@brianpm confirmed the radiaton change went in with cam6_4_016 (ESCOMP/CAM#1100).

@JulioTBacmeister thinks we might be just seeing seasonal variation. He will keep an eye on the model top T.

[cecilehannay]

Note we restarted the simulation clean => this is not relevant anymore.

change in gw parameters during the run

• from: 2000-01-01 to 2003-05

use_gw_front        = .true.

• from: 2003-05-01

use_gw_front        = .true.
taubgnd        = 1.5D-6

@JulioTBacmeister

[cecilehannay]

@PeterHjortLauritzen: The run is crashing with the error:

dec2012.hsn.de.hpc.ucar.edu 5720:  ERROR: negative moist layer thickness.  timestep or remap time too large

Should I adjust the sponge?

[cecilehannay]

From: @PeterHjortLauritzen:
We are getting this warning from the estimates

can you try se_hypervis_subcycle_sponge=9

[cecilehannay]

I successfully run with se_hypervis_subcycle_sponge=9 for a couple of years but the run crashed again with:

 ERROR: negative moist layer thickness.  timestep or remap time too large

@PeterHjortLauritzen @adamrher : please advice if I can change the sponge further.

[adamrher]

It's crashing at layer 8 with 250+ m/s winds. @PeterHjortLauritzen should we increase se_nsplit or se_rsplit?

k,dp_star_moist,u,v,T    8  -0.015685742200    253.63     91.88    185.27