search

NCAR / amwg_dev

Repo to store model sandboxes and cases used for CAM development
9 stars 2 forks source link

b.e23_alpha14a.BLT1850.ne30_t232.032 #283

Open cecilehannay opened 1 year ago

cecilehannay commented 1 year ago

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

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:

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/](https://webext.cgd.ucar.edu/BLT1850/b.e23_alpha14a.BLT1850.ne30_t232.032/atm/)

**Contacts:**
@cecilehannay, @gustavo-marques, @JulioTBacmeister
klindsay28 commented 1 year ago

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 commented 1 year ago

@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 commented 1 year ago

@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 commented 1 year ago

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

klindsay28 commented 1 year ago

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 commented 1 year ago

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

adamrher commented 1 year ago

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 commented 1 year ago

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 https://github.com/NCAR/amwg_dev/issues/283#issuecomment-1629828334, 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.