New dust emissions Test

[tilmes]

Purpose of the run

Test new dust emissions in the new tag cam running a FMTHIST.

Description of the run

FMTHIST starting in 2000 using the latest CAM tag: cam6_4_043. We can then compare to issue 600 run.

Suffix in the casename

Suffix

Namelist modifications

Standard output and add AOD output diagnostics:
AODDUSTnd, AOD

! 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','AODBCdn','AODSOAdn','AODSSdn','AODPOMdn','AODSO4dn',
'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',
AODDUSTdn','AODVISdn','AODBCdn','AODSOAdn','AODSSdn','AODPOMdn','AODSO4dn',
'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',
'SFncl_a1', 'SFncl_a2', 'SFncl_a3'
'ncl_a1', 'ncl_a1DDF', 'ncl_a1SFWET',
'ncl_a2', 'ncl_a2DDF', 'ncl_a2SFWET', 
'ncl_a3', 'ncl_a3DDF', 'ncl_a3SFWET', 
'ncl_c1', 'ncl_c1DDF', 'ncl_c1SFWET', 'ncl_c2', 'ncl_c2DDF', 'ncl_c2SFWET', 'ncl_c3', 'ncl_c3DDF', 'ncl_c3SFWET',
'SFbc_a4', 'bc_a1', 'bc_a1DDF', 'bc_a1SFWET', 'bc_a4', 'bc_a4_CLXF', 'bc_a4DDF', 'bc_a4SFWET',
'bc_c1', 'bc_c1DDF', 'bc_c1SFWET', 'bc_c4', 'bc_c4DDF', 'bc_c4SFWET', 
'SFdst_a1','SFdst_a2', 'SFdst_a3'
'dst_a1', 'dst_a1DDF', 'dst_a1SFWET', 'dst_a2', 'dst_a2DDF',
'dst_a2SFWET', 'dst_a3', 'dst_a3DDF', 'dst_a3SFWET', 'dst_c1', 'dst_c1DDF', 'dst_c1SFWET', 'dst_c2',
'dst_c2DDF', 'dst_c2SFWET', 'dst_c3', 'dst_c3DDF', 'dst_c3SFWET'
'SFpom_a4', 'pom_a1', 'pom_a1DDF', 'pom_a1SFWET', 'pom_a4', 'pom_a4DDF', 'pom_a4SFWET', 
'pom_c1','pom_c1DDF', 'pom_c1SFWET', 'pom_c4', 'pom_c4DDF', 'pom_c4SFWET'
'so4_a1', 'so4_a1_CHMP', 'so4_a1_CLXF', 'so4_a1DDF', 'so4_a1_sfgaex1',
         'so4_a1SFWET', 'so4_c1', 'so4_c1AQH2SO4', 'so4_c1AQSO4', 'so4_c1DDF', 'so4_c1SFWET'
'so4_a2', 'so4_a2_CHMP', 'so4_a2_CLXF', 'so4_a2DDF', 'so4_a2_sfgaex1', 'so4_a2_sfnnuc1', 'so4_a2SFWET',
          'so4_c2', 'so4_c2AQH2SO4', 'so4_c2AQSO4', 'so4_c2DDF', 'so4_c2SFWET'      
'so4_a3', 'so4_a3DDF', 'so4_a3_sfgaex1', 'so4_a3SFWET', 'so4_a5SFWET',
          'so4_c3', 'so4_c3AQH2SO4', 'so4_c3AQSO4', 'so4_c3DDF', 'so4_c3SFWET', 
'soa_a1', 'soa_a1_CHML', 'soa_a1DDF', 'soa_a1_sfgaex1', 'soa_a1SFWET', 
'soa_a2', 'soa_a2_CHML', 'soa_a2DDF', 'soa_a2_sfgaex1', 'soa_a2SFWET', 
soa_c1', 'soa_c1DDF',  'soa_c1SFWET', 'soa_c2', 'soa_c2DDF', 'soa_c2SFWET',
'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
!ne30 value
phys_grid_ctem_za_nlat=90 

!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

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=144

ncdata = '/glade/campaign/cgd/amp/juliob/ERA5/ne30np4/L93_cdf5/ERA5_x_ne30np4_L93_rgC2_WO_aero.2000-01-01-00000.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  


fsurdat = '/glade/derecho/scratch/aherring/restart/ctsm5.2.028/tools/mksurfdata_esmf/surfdata_ne30pg3_SSP3-7.0_1979_78pfts_c241014.nc'
flanduse_timeseries = '/glade/derecho/scratch/aherring/restart/ctsm5.2.028/tools/mksurfdata_esmf/landuse.timeseries_ne30pg3_SSP3-7.0_1979-2026_78pfts_c241014.nc' 
finidat = '/glade/campaign/cesm/cesmdata/inputdata/lnd/clm2/initdata_esmf/ctsm5.3/ctsm53n04ctsm52028_ne30pg3t232_BgcCrop_exice_hist.clm60.r.2000-01-01.nc'
use_init_interp = .true.
use_excess_ice = .false.

Source modifications

SourceMods: no source mods..

Sandbox

cam6_4_043

Contact info

@cecilehannay @adamrher @PeterHjortLauritzen

Any other relevant information

Enter relevant info

[cecilehannay]

@tilmes: The run crashed but i think it is just an incompatibility between the new land and the user_nl_cam. I will restart removing datasets from user_nl_cam

fsurdat = '/glade/derecho/scratch/aherring/restart/ctsm5.2.028/tools/mksurfdata_esmf/surfdata_ne30pg3_SSP3-7.0_1979_78pfts_c241014.nc'
flanduse_timeseries = '/glade/derecho/scratch/aherring/restart/ctsm5.2.028/tools/mksurfdata_esmf/landuse.timeseries_ne30pg3_SSP3-7.0_1979-2026_78pfts_c241014.nc' 
finidat = '/glade/campaign/cesm/cesmdata/inputdata/lnd/clm2/initdata_esmf/ctsm5.3/ctsm53n04ctsm52028_ne30pg3t232_BgcCrop_exice_hist.clm60.r.2000-01-01.nc'

[adamrher]

We need to be running the alpha04a tag. None of the cam tags have Danny's tuning that came into ctsm5.3.005 (cam6_4_043 uses ctsm5.3.002, cesm3_0_alpha04a uses ctsm5.3.007).

And yes, please remove those user_nl_clm entries pointing to those clm surface datasets (not sure about the finidat, I'd see if it works first).

[cecilehannay]

Thanks for your feedback. I can switch to the cesm3_0_alpha04a. There was some discussions at the end of the cam7 meeting and I think the reason cam6_4_043 was picked is that we thought the lastest tag included Danny's tuning. But it is doesn't we certainly want to go to alpha04a

[jedwards4b]

Yes -- please use cesm3_0_alpha04a. Let us know if there are any issues.

[cecilehannay]

@tilmes : I removed: so4_a5SFWET. It is not available

[cecilehannay]

@adamrher and @jedwards4b: this tag is happily running. Thanks for the help.

[olyson]

@adamrher , Per the request at the CESM project meeting, I'm happy to run land diagnostics on this. I don't see case name/location in this issue. A link to the ADF diagnostics would be useful as well so I know what to compare to (issue #600 ?), thanks.

[cecilehannay]

These are: #605 and #604
f.e30_alpha04a.FMTHIST.ne30_L93.002
f.e30_alpha04a.FMTHIST.ne30_L93.001

[adamrher]

@olyson I want to see the impact of ctsm5.2->ctsm5.3, so can we compare #604 to this baseline #600? Both of these use the Zender dust scheme, but RESTOM goes up by 0.2 W/m2 in alpha04a.

The comparison of #605 to #604 would not be as interesting to me, but that would show you the impact of going from Zender to Leung, using the same code base.

[cecilehannay]

Started as:
#605 and #604