f.cam6_3_109.FMTHIST_v0b.ne30.tuningF.001

[cecilehannay]

Description: FMTHIST with f.cam6_3_109 + tuning "F"

compset: FMTHIST_v0b add cosp output use clm5.2 dataset

Tuning:

clubb_l_predict_upwp_vpwp=.true. 
clubb_l_mono_flux_lim_um   = .true.  
clubb_l_mono_flux_lim_vm   = .true.  
clubb_c_uu_shr = 0.1
clubb_c7=0.1

dust_emis_fact         = 0.80D0

micro_mg_dcs=0.00025
micro_mg_vtrmi_factor=0.5
microp_aero_wsubi_scale=2.5
seasalt_emis_scale=1.5

user_nl_cam

mfilt           =       0,       5,     20,      40,      12,       120,      1,   1
nhtfrq              =       0,     -24,    -24,      -3,       0,       -2,      0,  -8760
ndens               =       2,       2,      2,       2,       2,       1,      2,   1
interpolate_output  =  .true.,  .true., .true., .false., .false., .true.,  .true.
interpolate_nlat    =     192,     192,    192,     192,     192,     192,   192
interpolate_nlon    =     288,     288,    288,     288,     288,     288,   288 

empty_htapes = .true.

fincl1 = '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', 
'CLD_MISR','FISCCP1_COSP','CLD_CAL','CLD_MISR','CLDTOT_CAL','CLDHGH_CAL', 'CLDMED_CAL','CLDLOW_CAL','CLMODIS', 'AODVISdn',
'CCN3', 'CDNUMC'

fincl3 = 'PRECT', 'PRECC', 'FLUT', 'U850', 'U200', 'V850', 'V200', 'OMEGA', 'PSL'

fincl4 =  'PRECC','PRECL'

fincl5 = 'Uzm','Vzm','Wzm','THzm', 'VTHzm','WTHzm','UVzm','UWzm'
phys_grid_ctem_nfreq=-6
phys_grid_ctem_zm_nbas=120
phys_grid_ctem_za_nlat=90

fincl6 =  'LGHTNG_CLD2GRND'

fincl7= 'AQSO4_H2O2','AQSO4_O3', 'bc_a1', 'bc_a4', 'dst_a1', 'dst_a2', 'dst_a3', 'ncl_a1',
'ncl_a1', 'ncl_a2', 'ncl_a3', 'pom_a1', 'pom_a4', 'so4_a1', 'so4_a2', 'so4_a3',
'soa_a1', 'num_a1', 'num_a2', 'num_a3', 'num_a4',
'bc_a1SFWET', 'bc_a4SFWET', 'dst_a1SFWET', 'dst_a2SFWET', 'dst_a3SFWET', 'ncl_a1SFWET',
'ncl_a2SFWET', 'ncl_a3SFWET', 'pom_a1SFWET', 'pom_a4SFWET', 'so4_a1SFWET', 'so4_a2SFWET', 'so4_a3SFWET', 'soa_a1SFWET',
'soa_a2SFWET', 'bc_c1SFWET', 'bc_c4SFWET', 'dst_c1SFWET', 'dst_c2SFWET', 'dst_c3SFWET', 'ncl_c1SFWET', 'ncl_c2SFWET',
'ncl_c3SFWET', 'pom_c1SFWET', 'pom_c4SFWET', 'so4_c1SFWET', 'so4_c2SFWET', 'so4_c3SFWET', 'soa_c1SFWET', 'soa_c2SFWET',
'bc_a1DDF', 'bc_a4DDF', 'dst_a1DDF', 'dst_a2DDF', 'dst_a3DDF', 'ncl_a1DDF', 'ncl_a2DDF', 'ncl_a3DDF',
'pom_a1DDF', 'pom_a4DDF', 'so4_a1DDF', 'so4_a2DDF', 'so4_a3DDF', 'soa_a1DDF', 'soa_a2DDF',
'so4_a1_CLXF', 'so4_a2_CLXF', 'SFbc_a4', 'SFpom_a4', 'SFso4_a1', 'SFso4_a2',
'so4_a1_sfgaex1', 'so4_a2_sfgaex1', 'so4_a3_sfgaex1', 'soa_a1_sfgaex1', 'soa_a2_sfgaex1',
'SFdst_a1','SFdst_a2', 'SFdst_a3', 'SFncl_a1', 'SFncl_a2', 'SFncl_a3',
'num_a2_sfnnuc1', 'so4_a5', 'num_a5', 'SFSO2', 'OCN_FLUX_DMS', 'SAD_SULFC', 'SAD_TROP', 'SAD_AERO'

clubb_l_predict_upwp_vpwp=.true. 
clubb_l_mono_flux_lim_um   = .true.  
clubb_l_mono_flux_lim_vm   = .true.  
clubb_c_uu_shr = 0.1
clubb_c7=0.1

dust_emis_fact         = 0.80D0

micro_mg_dcs=0.00025
micro_mg_vtrmi_factor=0.5
microp_aero_wsubi_scale=2.5
seasalt_emis_scale=1.5

user_nl_clm

use_init_interp = .true.
check_finidat_year_consistency = .false.
fsurdat = '/glade/work/slevis/git/mksurfdata_toolchain/tools/mksurfdata_esmf/surfdata_ne30np4.pg3_SSP5-8.5_78pfts_CMIP6_1850-2100_c230227.nc'
flanduse_timeseries = '/glade/work/slevis/git/mksurfdata_toolchain/tools/mksurfdata_esmf/landuse.timeseries_ne30np4.pg3_SSP5-8.5_78_CMIP6_1850-2100_c230227.nc'

Case directory: Locally (if still available): /glade/p/cesmdata/cseg/runs/cesm2_0/f.cam6_3_109.FMTHIST_v0b.ne30.tuningF.001

On github: https://github.com/NCAR/amwg_dev/tree/f.cam6_3_109.FMTHIST_v0b.ne30.tuningF.001

Sandbox: Locally (if still available): /glade/work/hannay/cesm_tags/cam6_3_109

On github: hash: 962921a (https://github.com/ESCOMP/CAM/tree/cam6_3_109)

Diagnostics: ADF diags (if available) https://webext.cgd.ucar.edu/FMTHIST/f.cam6_3_109.FMTHIST_v0b.ne30.tuningF.001/atm/

Contacts:

[cecilehannay]

The MT/L93 gives an error:

  trcdata_init: data type: INTERP_MISSING_MONTHS file: tracer_cnst_halons_3D_L70_
 1849-2015_CMIP6ensAvg_c180927.nc
 (GETFIL): attempting to find local file 
 tracer_cnst_halons_3D_L70_1849-2015_CMIP6ensAvg_c180927.nc
 (GETFIL): using 
 /glade/p/cesmdata/cseg/inputdata/atm/cam/tracer_cnst/tracer_cnst_halons_3D_L70_
 1849-2015_CMIP6ensAvg_c180927.nc
 open_trc_datafile: 
 /glade/p/cesmdata/cseg/inputdata/atm/cam/tracer_cnst/tracer_cnst_halons_3D_L70_
 1849-2015_CMIP6ensAvg_c180927.nc
 trcdata_init: file%has_ps =  T
 HALONS                           is not an invariant
 ERROR: tracer_cnst_init

[cecilehannay]

@fvitt or @tilmes: do you have guidance about that error?

[fvitt]

@cecilehannay The case is using trop_mam4 chemistry rather than ghg_mam4.

[cecilehannay]

Should we change the FMTHIST to use ghg_mam4 ?

[tilmes]

@fvitt I think the question would be if we should use trop_mam4 (as it is now) for both FMTHIST and FLHIST. Is there a reason why the ghg_mam4 does not use the same setting as trop_mam4?

Differences are: trop_mam4 includes HALONS as GHGs, and loss of methane and loss of N2O, etc. The specific additional reactions maybe not be relevant for the low-top model, but they may also not hurt to have, in order to have consistency between the two cases.

trop_mam4 has 4 extra tracers.

[fvitt]

I could be mistaken, but I thought the plan was to use ghg_mam4. FLHIST uses ghg_mam4.

[tilmes]

@cecilehannay can you add 'OCN_FLUX_DMS' for FMTHIST?

[tilmes]

@fvitt not sure who decided that, but sure. I am not sure what the differences are, maybe we should ask Dan Marsh? Not sure what his Github name is.

[fvitt]

@dan800

[adamrher]

I could be mistaken, but I thought the plan was to use ghg_mam4. FLHIST uses ghg_mam4.

I also thought that. Since FLTHIST has CAM%DEV%LT%GHGMAM4 I would've expected FMTHIST to be identical except with an "MT", but instead it's this CAM%DEV%MT. Maybe we decided to punt due to uncertainties in UBC and gravity wave settings at MT? But those are just namelists; I think we should've set the -chem flag to that in FLTHIST, as that's been settled, right?

[tilmes]

@adamrher I am not sure if I was part of the discussion what chemistry to use. Can someone point me to that discussion? Sorry, maybe there was a long github issue and I should refresh my mind. @cecilehannay this chemistry does not include _a5 aerosol, since it is not interactive.

[adamrher]

%GHG_MAM4 turns on -chem ghg_mam4. I think most of the discussions is in here https://github.com/ESCOMP/CAM/issues/762 and here https://github.com/ESCOMP/CAM/issues/765.

[tilmes]

@dan800 as stated before: Differences are: trop_mam4 includes HALONS as GHGs, and loss of methane and loss of N2O, etc. The specific additional reactions maybe not be relevant for the low-top model, but they may also not hurt to have, in order to have consistency between the two cases.

The advantage of using at least the production of H2O would be that we don't have to produce the external forcing files for water vapor from the full chemistry version. I don't see a discussion on the other tracers like N2O and CFCs. Would you want to comment on that?

@fvitt @cecilehannay @adamrher I am fine to use the -chem ghg_mam4 version for now (but would be great to get Dan's answer).

[dan800]

All physics configurations not controlled by a namelist variable (inc. chemistry) are identical to LT for MT. To begin with just use the LT namelist with a different IC (ncdata) file with MT levels. The default namelist should be pretty much identical to LW except:

1) Non-orogragraphic GWs: Specific tuning of the GW settings for MT vs LT should be discussed here: https://github.com/NCAR/amwg_dev/discussions/248#discussioncomment-5722514

2) The flux boundary condition on H2O should be should be changed to zero flux.

3) Other stuff I may have forgotten about that will show up when we try to run :-)

[dan800]

I could be mistaken, but I thought the plan was to use ghg_mam4. FLHIST uses ghg_mam4.

I also thought that. Since FLTHIST has CAM%DEV%LT%GHGMAM4 I would've expected FMTHIST to be identical except with an "MT", but instead it's this CAM%DEV%MT. Maybe we decided to punt due to uncertainties in UBC and gravity wave settings at MT? But those are just namelists; I think we should've set the -chem flag to that in FLTHIST, as that's been settled, right?

No punting! - we start from LT configuration for MT for all physics modules. Not changes for surface boundary conditions, emissions, anything except changing UBC conditions for H2O (zero flux) and tau_0_ubc = True

[adamrher]

No punting! - we start from LT configuration for MT for all physics modules. Not changes for surface boundary conditions, emissions, anything except changing UBC conditions for H2O (zero flux) and tau_0_ubc = True

OK! Then let's plan to discuss how to get this on the trunk quickly at the SE priorities mtg tmrw a.m. @cacraigucar.

Ill note separately that the BMT/BLT compsets we are currently making, are setting GHG_MAM4 for both https://github.com/ESCOMP/CESM/pull/223#pullrequestreview-1417118514. Looking like this B case PR will probably go into a cesm tag today.

[cacraigucar]

As was just pointed out in https://github.com/ESCOMP/CESM/pull/223#pullrequestreview-1417118514, @JulioTBacmeister, @cecilehannay and I have put a hold on making a CESM tag until there are use_cases for BLT1850_v0b and BMT1850_v0b. The discussion for the use_case settings is here: https://github.com/NCAR/amwg_dev/discussions/261

[cecilehannay]

Issue with the FMTHIST compset in cam6_3_109. This run was replaced by #263.