f.cam6_3_100.FWscHIST.ne30_L58.ebudget.001

[cecilehannay]

Description: Science evaluate "new energy" code (incl. condensates in pressure and having energy fixer use an energy consistent with the dynamical core) and compare to cam6_3_100 baseline (#227)

mfilt           =  0,   5,   20,  40, 1,  240,  365, 73
nhtfrq         =  0, -24,  -6,  -3,  0,   1,    -24, -24
fexcl1 = ' '  

fincl1 = 'ACTNI', 'ACTNL', 'ACTREI', 'ACTREL', 'AODDUST', 'AODVIS', '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',
'WV_phBF:A','WL_phBF:A','WI_phBF:A','SE_phBF:A','KE_phBF:A',
'WV_phBP:A','WL_phBP:A','WI_phBP:A','SE_phBP:A','KE_phBP:A',
'WV_phAP:A','WL_phAP:A','WI_phAP:A','SE_phAP:A','KE_phAP:A',
'WV_phAM:A','WL_phAM:A','WI_phAM:A','SE_phAM:A','KE_phAM:A'

fincl2 = 'OMEGA', 'PMID', 'PS', 'Q', 'QRL', 'QRS', 'T', 'TROP_P', 'TROP_T', 'U', 'V', 'Z3'

fincl4 =  'PRECC','PRECL'

history_aerosol = .true.
thermo_budget_history = .true.
thermo_budget_histfile_num = 5
avgflag_pertape                = 'A', 'A', 'A', 'A', 'N', 'A', 'A', 'A', 'I'

ncdata = '/glade/p/cesm/amwg_dev/juliob/FWsc_ne30pg3_58L_GRID_48_taperstart10km_lowtop_BL10_v3_beta1p75_Top_43km.nc'

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

se_rsplit            = 3
se_nsplit           = 2
se_hypervis_subcycle =  4

bnd_topo  = '/glade/p/cgd/amp/juliob/bndtopo/latest/ne30pg3_gmted2010_modis_bedmachine_nc3000_Laplace0100_20230105.nc'

use_gw_convect_dp = .false.
use_gw_front = .false.

ubc_specifier = 'Q:H2O->UBC_FILE'
ubc_file_path = '/glade/p/cesmdata/inputdata/atm/cam/chem/ubc/f.e21.FWHISTBgcCrop.f09_f09_mg17.CMIP6-AMIP-WACCM.ensAvg123.cam.h0zm.UBC.195001-201412_c220322.nc'
ubc_file_input_type = 'SERIAL'

co2_flag = .false.

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

Case directory: Locally (if still available): /glade/p/cesmdata/cseg/runs/cesm2_0/f.cam6_3_100.FWscHIST.ne30_L58.ebudget.001

On github: https://github.com/NCAR/amwg_dev/tree/f.cam6_3_100.FWscHIST.ne30_L58.ebudget.001

SourceMods: Contains SourceMods for:

• Thomas bugfix (@tto061)
• Keith bugfix (https://github.com/NCAR/amwg_dev/issues/227#issuecomment-1490669430)

Sandbox: Locally (if still available): /glade/work/hannay/cesm_tags/f.cam6_3_100.FWscHIST.ne30_L58.ebudget.001

On github: hash: 80b076d

git clone https://github.com/jtruesdal/CAM-1.git cam6_3_100_ebudget_dev_update
cd cam6_3_100_ebudget_dev_update
git checkout ebudget_dev_update
./manage_externals/checkout_externals

Diagnostics: AMWG diags (if available) https://webext.cgd.ucar.edu/FWscHIST/f.cam6_3_100.FWscHIST.ne30_L58.ebudget.001/atm/

Contacts: @PeterHjortLauritzen, @jtruesdal, @cecilehannay

[cecilehannay]

@justin-richling: could you also run the ADF for f.cam6_3_100.FWscHIST.ne30_L58.ebudget.001 versus f.cam6_3_100.FWscHIST.ne30_L58.001 asap. Please use 5-year climos.

[justin-richling]

Just saw this, I'm on it

[justin-richling]

ADF diags are here

[adamrher]

@justin-richling the FLNT values in the table are 0.2 W/m2 different, but the global mean in the lat-lon maps are 237.82 vs. 237.83, or a 0.01 W/m2 difference.

[adamrher]

From Julio

FYI for Justin none of the FLNT numbers on the maps match with the table.

[adamrher]

Just going to post some zonal mean plots here.

That T anomaly in the Northern high-latitude stratosphere is from the dynamics:

I have a hard time connecting that change to these new energy mods; it could alternatively be due to inter-annual variability, since there's a corresponding reduction in U winds in the same location (see ADF diags).

Moving on to the total physics tendency:

My understanding is that we should expect changes in this field as the physics tendencies will be a bit larger in regions with water species. But then I think there will also be a change in the global energy fixer, since phys and dyn energies are more consistent in the new run, and so potentially a smaller contribution from the energy fixer (?). Let me know if you can think of other variables for me to plot to get at these differences more directly.

[PeterHjortLauritzen]

Agreed that it is hard to interpret the results! Some thoughts below:

Temperature increment 1

For a fixed heating increment dQ the temperature increment with the new energy code should be smaller:

dT = dQ/(rho x cp_d). (control)

dT = dQ/(rho_d x cp_d+rho_wv x cp_wv+condensate terms). (new energy run)

where rho is density (incl. all forms of water), rho_wv is density of water vapor, cp_d specific heat at constant pressure for dry air and cp_wv specific heat at constant pressure for water vapor. Since cp_wv>cp_d the temeprature increment (all else being equal) should be smaller.

Temperature increment 2

I expect the energy fixer T tendency (uniform everywhere) to have changed mainly because the dry-mass adjustment has changed from ~0.3Wm/^2 to ~0.6W/m^2. The energy dissipitation of the dynamical core is ~ -0.08 W/m^2. (these numbers are 1 months average values from the last month of the run). Hence the energy fixer should remove around ~0.5W/m^2 compared to ~0.2W/m^2 in terms of a negative temperature tendency.

Sponge changes

The new energy run does not have frictional heating in the sponge.

[justin-richling]

@justin-richling the FLNT values in the table are 0.2 W/m2 different, but the global mean in the lat-lon maps are 237.82 vs. 237.83, or a 0.01 W/m2 difference.

@JulioTBacmeister @adamrher

I apologize, this has been an open Issue for the ADF but hasn't be resolved in a PR yet. I've put a current set of ADF diagnostics that should reflect the plot values in the tables:

https://webext.cgd.ucar.edu/FWscHIST/f.cam6_3_100.FWscHIST.ne30_L58.ebudget.001/atm/f.cam6_3_100.FWscHIST.ne30_L58.ebudget.001_1979_1985_vs_f.cam6_3_100.FWscHIST.ne30_L58.001_1979_1985/

[justin-richling]

Updated ADF diags for 1979-1991 are here