N1850OCBDRDDMS_f19_tn14_250119#

Data storage#

The data is stored on NIRD @ sigma2

/projects/NS2345K/DO_NOT_TOUCH/NS2345K/noresm/cases/N1850OCBDRDDMS_f19_tn14_250119

Path to case directory#

on Vilje @ sigma2

/home/metno/oyvinds/NorESM2old/noresmdevcases/N1850OCBDRDDMS_f19_tn14_250119

Path to diagnostics#

http://ns2345k.web.sigma2.no/diagnostics/noresm/common/N1850OCBDRDDMS_f19_tn14_250119/

Summary of simulation#

  • To reduce the net radiation imbalance @TOM (top of model) the sea-salt emissions were increased by 10%

  • Increased (x2) error tolerance in energy conservation test in CICE

  • Modifications to the convection code included as SourceMod

  • Continued to use the previous namelist changes compared to repository for CAM6-Nor, MICOM and CLM5

  • New namelist additions compared to repository for CAM6-Nor:

    • aerotab_table_dir = ‘/cluster/shared/noresm/inputdata/noresm-only/atm/cam/camoslo/AeroTab_20oct16’ -> ‘/cluster/shared/noresm/inputdata/noresm-only/atm/cam/camoslo/AeroTab_8jun17’

  • Namelist changes compared to repository for MICOM:

    • bkopal = 1.e-6 -> 5.e-6

    • rcalc = 48. -> 35.

    • ropal = 35. -> 45.

For all user name list specifics, see bottom of this page

Simulation specifics#

CESM parent

CESM2.0.0

Parent

N1850OCBDRDDMS_f19_tn14_201218

Run type

branch

Branch time from parent

01-01-0091

Simulated years

01-01-0091 - 31-12-0645

Compset

1850_CAM60%PTAERO_CLM50%BGC-CROP_CICE_MICOM%ECO_MOSART_SGLC_SWAV_BGC%BDRDDMS

Git branch

featureCAM5-OsloDevelopment_trunk2.0-6

Git commit

9aa22b2

Resolution

f19_tn14

Machine

Vilje

Node allocation#

 <entry id="NTASKS">
      <type>integer</type>
      <values>
        <value compclass="ATM">768</value>
        <value compclass="CPL">768</value>
        <value compclass="OCN">186</value>
        <value compclass="WAV">300</value>
        <value compclass="GLC">768</value>
        <value compclass="ICE">504</value>
        <value compclass="ROF">8</value>
        <value compclass="LND">256</value>
        <value compclass="ESP">1</value>
      </values>
      <desc>number of tasks for each component</desc>
    </entry>

Code modifications (SourceMods)#

iHAMOCC modifications#

In components/micom/hamocc/beleg_bgc.F90

Line 209

bkopal = 1.e-6    !i.e. 1.0 mmol Si/m3

to

bkopal = 5.e-6    !i.e. 1.0 mmol Si/m3

Line 269 from

#elif defined(WLIN) 
      rcalc = 48.  ! calcium carbonate to organic phosphorous production ratio
      ropal = 35.  ! opal to organic phosphorous production ratio

to

#elif defined(WLIN) 
      rcalc = 35.  ! calcium carbonate to organic phosphorous production ratio
      ropal = 45.  ! opal to organic phosphorous production ratio

1.1 x seasalt emissions#

Line 176 in components/cam/src/chemistry/oslo_aero/seasalt_model.F90

!Eqn. 9 in Salter et al. (2015)
       numberFlux(:ncol,n) = 1.1_r8*whitecapAreaFraction(:ncol)*                                                     &
                            ( coeffA(n)*(sst(:ncol)-273.15_r8)*(sst(:ncol)-273.15_r8)*(sst(:ncol)-273.15_r8)  &
                            + coeffB(n)*(sst(:ncol)-273.15_r8)*(sst(:ncol)-273.15_r8)                         &
                            + coeffC(n)*(sst(:ncol)-273.15_r8)                                                &
                            + coeffD(n) )

Moist convection in CAM#

Moist convection modifications (“zmst” modifications) in

components/cam/src/NorESM/zm_convF90:

Increased error tolerance in energy conservation test in CICE#

ferr = energy conservation error (W m-2)

Line 2390 in /components/cice/src/source/ice_therm_vertical.F90

changed from

if (ferr > ferrmax) then

to

if (ferr > 2*ferrmax) then

Includes the long wave aod error#

Information about the bug: The aerosol long wave calculations used information from the aerosol shortwave interpolation on aerosol size. The result was that aerosol longwave forcing was not included during night. A first estimate based on estimates from AMIP simulation is + 0.03 W/m2. The forcing is not evenly distributed, but mostly focused on Sahara including downstream and the Arabian peninsula. The numbers here are around 1-2 W/m2.

Note this bug was fixed in N1850OCBDRDDMS_f19_tn14_07052019

User name lists#

user_nl_cam#

&dyn_fv_inparm
 fv_am_correction= .true.
 fv_am_diag      = .true.
 fv_am_fix_lbl   = .true.
 fv_am_fixer     = .true.

&phys_ctl_nl
 dme_energy_adjust = .true.


&zmconv_nl
 zmconv_c0_lnd          =  0.0200D0
 zmconv_c0_ocn          =  0.0200D0
 zmconv_ke              =  8.0E-6

&micro_mg_nl
 micro_mg_dcs             = 5.0e-4

&clubb_params_nl
 clubb_gamma_coef = 0.258

&gw_drag_nl
 tau_0_ubc                = .true.

&cldfrc_nl
 cldfrc_iceopt          =  4
 
&phys_ctl_nl
aerotab_table_dir =
'/work/shared/noresm/inputdata/noresm-only/atm/cam/camoslo/AeroTab_8jun17'

user_nl_clm#

Reset snow: Remove infiltration excess water as runoff if the temperature of the surface water pool is below freezing.

finidat = '/work/shared/noresm/inputdata/cesm2_init/b.e20.B1850.f09_g17.pi_control.all.297/0308-01-01/b.e20.B1850.f09_g17.pi_control.all.297.clm2.r.0308-01-01-00000.nc'
use_init_interp = .true.
reset_snow = .true.

Time series of spinup#

LMspinup2 NorESM2-LM spinup simulations
Left column (from top to bottom): Globally and annually averaged Net radiation @ top of model, Surface (2m) air temperature, Sea surface temperature (SST), global and volume averaged ocean temperature, Atlantic meridional oveturning circulation (AMOC) @ 26.5N.
Right column (from top to bottom): Globally and annually sum of Sea salt surface emissions, DMS (dimethylsulfide) surface emissions, globally and annually averaged vertically-integrated total cloud cover, shortwave cloud forcing and longwave cloud forcing.