search

NCAR / amwg_dev

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

b.cesm3_cam058_mom_c.B1850WscMOM.f09_L58_t061.012 #140

Open cecilehannay opened 2 years ago

cecilehannay commented 2 years ago

Description: New coupled simulation similar to 009 (issue https://github.com/NCAR/amwg_dev/issues/136) but with the FV dycore.

Uses the FV tuning for clubb an dust

CAM_CONFIG_OPTS = -phys cam_dev -microphys mg2 -chem waccm_sc_mam4 -nlev 58

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

empty_htapes = .true.

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'

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

fincl4 = 'PRECC','PRECL'

tracer_cnst_specifier = 'O3','OH','NO3','HO2','HALONS'

co2_flag=.false.

bnd_topo = '/glade/work/juliob/bndtopo/fv_0.9x1.25_gmted2010_modis_bedmachine_nc3000_Laplace0100_20220708.nc'

clubb_c14 = 1.6D0 clubb_gamma_coef = 0.270 dust_emis_fact = 0.8D0


- user_nl_cice

update_ocn_f = .true. f_fresh = 'mxxxx' f_fresh_ai = 'mxxxx' saltflux_option = 'prognostic'

- user_nl_mom

ENTHALPY_FROM_COUPLER = True MEKE_GEOMETRIC_ALPHA = 0.10

- user_nl_clm

use_init_interp = .true. 



**Case directory:**
Locally (if still available):
/glade/p/cesmdata/cseg/runs/cesm2_0/b.cesm3_cam058_mom_c.B1850WscMOM.f09_L58_t061.012

On github:
https://github.com/NCAR/amwg_dev/tree/b.cesm3_cam058_mom_c.B1850WscMOM.f09_L58_t061.012

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

On github:
https://github.com/ESCOMP/CESM/tree/cesm3_cam6_3_058_MOM
hash: 228563a

**Diagnostics:**
AMWG diags (if available)
[https://webext.cgd.ucar.edu/BWsc1850MOM/b.cesm3_cam058_mom_c.B1850WscMOM.f09_L58_t061.012/atm/](https://webext.cgd.ucar.edu/BWsc1850MOM/b.cesm3_cam058_mom_c.B1850WscMOM.f09_L58_t061.012/atm/)

**Contacts:**
@cecilehannay, @gustavo-marques, @JulioTBacmeister
cecilehannay commented 2 years ago

@dabail10 I noticed that the variable CICE_CONFIG_OPTS in env_run.xml comes out as: CICE_CONFIG_OPTS= "-phys cice5" whether the model is using cice5 or cice6

For instance, in /glade/p/cesmdata/cseg/runs/cesm2_0/b.cesm3_cam058_mom_c.B1850WscMOM.f09_L58_t061.012

The README.case says: 2022-07-12 13:16:48: Component ICE is Sea ICE (cice) model version 6

However, in the env_run.xml CICE_CONFIG_OPTS= "-phys cice5"

Is it a bug? Is it something we should have modified manually? We didn't do that.

Please could you advice. Thanks

cecilehannay commented 2 years ago

This run is now 12 years and it looks like the Labrador is starting to freeze. @gustavo-marques : could you look at the convection. @justin-richling : could you plot the ADF compared to 009 including the timeseries. Thanks.

gustavo-marques commented 2 years ago

Convection in the Lab Sea stopped at year 10 and has not recovered so far.

cecilehannay commented 2 years ago

It is also what we see in the icefrac from the cam file.

@dabail10 have better diags but it seems clear that Labrador Deniz'i Dondu.

justin-richling commented 2 years ago

The ADF diags w/ time series are now available here

The categories may be a bit misleading and not easy to search, just for reference here is the list if you are looking for a specific variable but can't find it in the drop downs:

           'Clouds': {'ACTNI', 'ACTNL', 'ACTREI', 'ACTREL', 'ADRAIN', 'ADSNOW',
                       'AREI', 'AREL', 'CCN3', 'CDNUMC', 'CLDHGH', 'CLDICE',
                       'CLDLIQ', 'CLDLOW', 'CLDMED', 'CLDTOT', 'CLOUD', 'CONCLD',
                       'EVAPPREC', 'EVAPSNOW', 'FCTI', 'FCTL', 'FICE', 'FREQI',
                       'FREQL', 'FREQR', 'FREQS', 'MPDQ', 'PRECC', 'PRECL',
                       'PRECSC', 'PRECSL', 'PRECT', 'TGCLDIWP', 'TGCLDLWP'},
            'Deep Convection': {'CAPE', 'CMFMC_DP', 'FREQZM', 'ZMDQ', 'ZMDT'},
            'COSP': {'CLDTOT_ISCCP', 'CLIMODIS', 'CLTMODIS', 'CLWMODIS',
                     'FISCCP1_COSP', 'ICE_ICLD_VISTAU', 'IWPMODIS',
                     'LIQ_ICLD_VISTAU', 'LWPMODIS', 'MEANCLDALB_ISCCP',
                     'MEANPTOP_ISCCP', 'MEANTAU_ISCCP', 'MEANTB_ISCCP',
                     'MEANTBCLR_ISCCP', 'PCTMODIS', 'REFFCLIMODIS', 'REFFCLWMODIS',
                     'SNOW_ICLD_VISTAU', 'TAUTMODIS', 'TAUWMODIS',
                     'TOT_CLD_VISTAU', 'TOT_ICLD_VISTAU'},
            'Budget': {'DCQ', 'DQCORE', 'DTCORE', 'MPDICE', 'MPDLIQ', 'PTEQ'},
            'Radiation': {'FLNS', 'FLNSC', 'FLNT', 'FLNTC', 'FLUT', 'FSDS',
                          'FSDSC', 'FSNS', 'FSNSC', 'FSNT', 'FSNTC', 'FSNTOA',
                          'LHFLX', 'LWCF', 'QRL', 'QRS', "RESTOM", 'SHFLX', 'SWCF'},
            'State': {'OMEGA', 'OMEGA500', 'PINT', 'PMID', 'PS', 'PSL', 'Q',
                      'RELHUM', 'T', 'U', 'V', 'Z3', 'Wind'},
            'Surface': {'PBLH', 'QFLX', 'TAUX', 'TAUY', 'TREFHT', 'U10',
                        'Surface_Wind_Stress'},
            'GW': {'QTGW', 'UGTW_TOTAL', 'UTGWORO', 'VGTW_TOTAL', 'VTGWORO'},
            'CLUBB': {'RVMTEND_CLUBB', 'STEND_CLUBB', 'WPRTP_CLUBB', 'WPTHLP_CLUBB'}
justin-richling commented 2 years ago

I should also add the naming of the climo years is 0001-0010, but this is an error, the actual data climo years do go from 0001-0012