noresm2_3_develop : Introduction of possibility to run with full atmospheric chemistry

[DirkOlivie]

Summary: Introduction of the ability to run with more complex atmospheric gas-phase chemistry.

Contributors: Dirk Olivié

Reviewers: Steve Goldhaber

Purpose of changes: #100 Add trop+strat chemistry scheme. The new scheme is a combination of the aerosols from OSLO_AERO and the gas-phase species from the TS1 chemistry (pp_trop_strat_mam4_vbs).

Github PR URL:

Changes made to build system: None

Changes made to the namelist:

1. the namelist used by the new compsets contains elements from the standard aerosol-only NorESM namelist, but some external forcings have been removed as not used anymore (ozone climatology, ozone climatology, volcanic stratospheric AOD).
2. New things are : more emissions files are used, more GHGs/ODSs are prescribed, and the lightning NOx scaling factor is now used.
3. Concerning MEGAN fluxes from the land-model : more species are mapped into monoterpene (to make it more in line with CESM).
4. 3D-chemical species are taken into account in the radiation

Changes to the defaults for the boundary datasets:

1. there are additional emission files
2. some climatologies are not used anymore (ozone, oxidants)

Substantial timing or memory changes: an N1850_f19_tn14 simulation on 1024 CPUs is much slower. The throughput is 8.67 simulated years/day (compared to 19.84 simulated years/day in aerosol-only version)

Detailed description of changes :
General files :
bld/build-namelist
bld/config_files/definition.xml [new chemical package "tropstrat_mam_oslo]
bld/configure
bld/namelist_files/use_cases/1850_cam6_noresm_tropstratchem.xml [4 new use cases]
bld/namelist_files/use_cases/hist_cam6_noresm_tropstratchem.xml
bld/namelist_files/use_cases/ssp370_cam6_noresm_tropstratchem.xml
bld/namelist_files/use_cases/ssp585_cam6_noresm_tropstratchem.xml
cime_config/config_component.xml
cime_config/config_compsets.xml
cime_config/testdefs/testlist_cam.xml

Files in specific chemistry directory pp_tropstrat_mam_oslo (generated by chemical preprocessor):
src/chemistry/pp_tropstrat_mam_oslo/chem_mech.doc
src/chemistry/pp_tropstrat_mam_oslo/chem_mech.in
src/chemistry/pp_tropstrat_mam_oslo/chem_mods.F90
src/chemistry/pp_tropstrat_mam_oslo/m_rxt_id.F90
src/chemistry/pp_tropstrat_mam_oslo/m_spc_id.F90
src/chemistry/pp_tropstrat_mam_oslo/mo_adjrxt.F90
src/chemistry/pp_tropstrat_mam_oslo/mo_exp_sol.F90
src/chemistry/pp_tropstrat_mam_oslo/mo_imp_sol.F90
src/chemistry/pp_tropstrat_mam_oslo/mo_indprd.F90
src/chemistry/pp_tropstrat_mam_oslo/mo_lin_matrix.F90
src/chemistry/pp_tropstrat_mam_oslo/mo_lu_factor.F90
src/chemistry/pp_tropstrat_mam_oslo/mo_lu_solve.F90
src/chemistry/pp_tropstrat_mam_oslo/mo_nln_matrix.F90
src/chemistry/pp_tropstrat_mam_oslo/mo_phtadj.F90
src/chemistry/pp_tropstrat_mam_oslo/mo_prod_loss.F90
src/chemistry/pp_tropstrat_mam_oslo/mo_rxt_rates_conv.F90
src/chemistry/pp_tropstrat_mam_oslo/mo_setrxt.F90
src/chemistry/pp_tropstrat_mam_oslo/mo_sim_dat.F90

Test run: ERP_D_Ln9.f19_f19_mtn14.NF1850norbc_tropstratchem.betzy_intel.cam-outfrq9s

closes #100

[gold2718]

This mostly looks good but I have some questions and concerns.

[gold2718]

Please update the tag in Externals_CAM.cfg to be oslo_aero_2_3a2.

[DirkOlivie]

Hi Steve (@gold2718), I have now updated the OSLO_AERO tag in Externals_CAM.cfg to oslo_aero2_3a2.

[gold2718]

Looks good now!

[gold2718]

Are there any expected answer or namelist changes when not running the new chemistry scheme?

[DirkOlivie]

Hi Steve (@gold2718), no there should not be any answer changes. I will try to look into it now and do some tests.

[DirkOlivie]

Hi Steve (@gold2718),
I have now done two simulations for the N1850 compset and they gave identical results. The simulations are using :

• my latest version of CAM and OSLO_AERO I had on Betzy, where the full-chemistry code was in
• CAM (commit 3db5e69) and OSLO_AERO (commit 8b3e7dde) versions from before the full-chemistry code was merged in.
  Do these comparisons make sense, or did I do something wrong?

[gold2718]

I believe the OSLO_AERO commit should be tag = oslo_aero_2_3a1 which is NorESMhub/OSLO_AERO@1a03f8c
The CAM hash should be 3db5e69 .
It looks like your CAM commit is correct but your OSLO_AERO commit is different.

[DirkOlivie]

Hi Steve (@gold2718), I have now checked out for OSLO_AERO the oslo_aero_2_3a1 tag. It differs in what I have in oslo_aero_depos.F90. The oslo_aero_2_3a1 tag does not have the lines :

! Increase scavenging efficiency for large soluble particles.
              if ((lphase==1).and.((mm==l_ss_a2).or.(mm==l_ss_a3).or.(mm==l_so4_a2))) then
                 sol_factic=1.0_r8
                 f_act_conv=1.0_r8
              end if

Might this be related to the difference?

[gold2718]

I think they could be. These lines come from a commit by @oyvindseland on 4. jul. as PR NorESMhub/OSLO_AERO#30

Commit message:
Increased in-cloud convective scavenging coefficient for sea-salt and aqueous phase produced sulfate

I need to make a new baseline with that change, thank you for finding this.

[gold2718]

Good news, using the correct hash for OSLO_AERO, this PR is BFB for the previously existing tests.