b.e23_alpha14a.BLT1850.ne30_t232.032

[cecilehannay]

Description:
Coupled simulation with updated components as of May 26, 2023

• Tag: cesm2_3_alpha14a
• Compset: BLT1850_v0c
• Resolution: ne30pg3_t232

The run initial result were advertised at co-chair meeting on May 30, 2023

Case directory:
Locally (if still available):
/glade/p/cesmdata/cseg/runs/cesm2_0/b.e23_alpha14a.BLT1850.ne30_t232.032

On github:
https://github.com/NCAR/amwg_dev/tree/b.e23_alpha14a.BLT1850.ne30_t232.032

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

On github:

git clone https://github.com/ESCOMP/cesm cesm2_3_alpha14a
cd cesm2_3_alpha14a
git checkout  cesm2_3_alpha14a
./manage_externals/checkout_externals -o

Namelists:

• user_nl_cam
  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
seasalt_emis_scale=1.5

micro_mg_dcs=0.00025
micro_mg_vtrmi_factor=0.5
microp_aero_wsubi_scale=2.5

use_gw_convect_dp = .true.
use_gw_front = .true.
tau_0_ubc = .true.

• output specified below

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','OMEGA500'

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

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', 'SFSO2', 'OCN_FLUX_DMS', 'SAD_SULFC', 'SAD_TROP', 'SAD_AERO'

• user_nl_clm

 use_init_interp = .true. 
 fsurdat = '/glade/work/slevis/git/mksurfdata_toolchain/tools/mksurfdata_esmf/surfdata_ne30np4.pg3_SSP5-8.5_78pfts_CMIP6_1850-2100_c230227.nc'
 finidat = '/glade/p/cesmdata/cseg/inputdata/cesm2_init/b.cesm3_cam058_mom_e.B1850MOM.f09_L32_t061.cam6_cice5.016/0058-01-01/b.cesm3_cam058_mom_e.B1850MOM.f09_L32_t061.cam6_cice5.016.clm2.r.0058-01-01-00000.nc'

• user_nl_cice

f_fresh = 'mxxxx'
f_fresh_ai = 'mxxxx'
ice_ic = '/glade/p/cesmdata/cseg/inputdata/cesm2_init/gmom.e23.GJRAv4.TL319_t061_zstar_N65.tx2_3v2.001/0061-01-01/gmom.e23.GJRAv4.TL319_t061_zstar_N65.tx2_3v2.001.cice.r.0061-01-01-00000.nc'

• user_nl_mom

MEKE_GEOMETRIC_ALPHA = 0.10

Source Mods and config options:

• CAM_CONFIG_OPTS
  ./xmlchange --append CAM_CONFIG_OPTS=-cosp

In /glade/p/cesmdata/cseg/runs/cesm2_0/b.e23_alpha14a.BLT1850.ne30_t232.032/SourceMods

• src.clm:
  clm_varpar.F90 surfrdMod.F90

• src.mom:
  diag_table

Diagnostics:
AMWG diags (if available)
https://webext.cgd.ucar.edu/BLT1850/b.e23_alpha14a.BLT1850.ne30_t232.032/atm/

Contacts:
@cecilehannay, @gustavo-marques, @JulioTBacmeister

[klindsay28]

The casename in SourceMods/src.mom/diag_table has a typo in it.
It has the substring alpha13b, instead of alpha14a.
So the short-term archiver is not working, and neither are my plotting scripts.
@cecilehannay , could you correct the casename before submitting additional years?

A next step would be to rename the mom6 output files in RUNDIR and manually run the short-term archiver,
so that the existing years of output can be handled with existing scripts.

[adamrher]

@cecilehannay @JulioTBacmeister can we run diagnostics on this B-case run compared to the equivalent F-case? I'd like to see exactly how coupling to the ocean lowers RESTOM by ~2 W/m2.

I'm actually not seeing a clean F case equivalent with a post 110 tag (110 and earlier has double counting via the external H2O forcing file; this B-case here uses tag 111). Maybe one of the gw tuning runs? I'd suspect the gw tunings don't materially effect the cloud forcing.

[cecilehannay]

@klindsay28 and @gmarques: I rename the mom6 output files and I moved them to the archive. I hope your scripts are working now. Please let me know if they are more issues.
I didn't extent the run as it was crashing with the error:

2121:WARNING from PE   321: Extreme surface sfc_state detected: i= 477 j= 478 lon=  71.532 lat=  66.371 x=  30.667 y=  87.645 D= 5.1124E+01 SSH= 9.3342E+00 SST= 8.2610E-02 SSS=-6.3630E-30 U-= 2.3987E-03 U+= 0.0000E+00 V-= 0.0000\
E+00 V+=-1.2492E-03
2121:
2121:
2121:FATAL from PE   321: There were a total of         1 locations detected with extreme surface values!
2121:

@adamrher: I asked @justin-richling to run this case vs f.cam6_3_112.FLTHIST_v0c.ne30.non-ogw-ubcT.001. I would add that a 2 W/m2 RESTOM difference between the F case and the B case is very standard. It has been that way for as long as I remember.

[justin-richling]

ADF vs f.cam6_3_112.FLTHIST_v0c.ne30.non-ogw-ubcT.001 here

[klindsay28]

MOM is aborting because SSS is going negative, -6e-30. In the preceding month, Dec-0016, the daily SSS values are decreasing from 8.6e-5 to 7.1e-12. So it's not surprising that it's dipping below zero.

I'm not sure what process is driving SSS negative at this point, which is a coastal point off northern Russia.

That said, @dabail10, I seem to recall that you had ideas to cap/reduce salt uptake by CICE from MOM when SSS is low. What is the status of those ideas being implemented?

[gustavo-marques]

@cecilehannay, can you please add MIN_SALINITY = 1.0E-06 in user_nl_mom and submit the run again?

[adamrher]

Does anyone happen to have know whether this 25 yr run produced reasonable lab sea ice and AMOC mixed layer depths? I'm just trying to get a feel for the coupled climate when those aspects look OK.

[JulioTBacmeister]

I recall that this run produced a frozen Lab Sea.

…

On Mon, Jul 10, 2023 at 4:43 PM Adam Herrington ***@***.***> wrote: Does anyone happen to have know whether this 25 yr run produced reasonable lab sea ice and AMOC mixed layer depths? I'm just trying to get a feel for the coupled climate when those aspects look OK. — Reply to this email directly, view it on GitHub <#283 (comment)>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/ACGLMTWTJU3DPC3KDMMKKVLXPSARJANCNFSM6AAAAAAYUQ3PKI> . You are receiving this because you were mentioned.Message ID: ***@***.***>

-- *My working day may not be your working day. Please do not feel obliged to reply to this message outside of your normal working hours.*