n1850.ne30_tn14.hybrid_fatessp.202401007

[mvertens]

Purpose:
Exactly like n1850.ne30_tn14.hybrid.20240829

• use of FATES-SP and spun up initial conditions for FATES-SP using offline spinup
• NOTE: a memory leak is present in the run and hits agains the node limit when running past 11 years.
  See memory leak in n1850 case NorESM#550
• NOTE: the goal is to run 20 years. When that is done - the code will be moved to /nird/datalake.

---

Description:

The following should be done from login3. Before you submit you should issue the command
ml swap hpcx/2.20

• create_newcase commnad

 ./create_newcase --case <user_case_path>/n1850.ne30_tn14.hybrid_fatessp.202401007 --compset 1850_CAM%DEV%LT%NORESM%CAMoslo_CLM60%FATES-SP_CICE_BLOM%ECO_MOSART_DGLC%NOEVOLVE_SWAV_SESP --res ne30pg3_tn14 --project nn9039k --run-unsupported --mach betzy --compiler intel-oneapi

• user_nl_cam

 use_aerocom = .true.
 interpolate_nlat   = 192
 interpolate_nlon   = 288
 interpolate_output = .true.
 history_aerosol = .true. 
 zmconv_c0_lnd  =  0.0075D0
 zmconv_c0_ocn  =  0.0300D0
 zmconv_ke =  5.0E-6
 zmconv_ke_lnd =  1.0E-5
 clim_modal_aero_top_press =  1.D-4
 bndtvg = '/cluster/shared/noresm/inputdata/atm/cam/ggas/noaamisc.r8.nc'

• user_nl_blom
  ICFILE = "/cluster/shared/noresm/inputdata/restart/NOIIAJRAOC20TR_TL319_tn14_ppm_20240816/rest/1775-01-01-00000/NOIIAJRAOC20TR_TL319_tn14_ppm_20240816.blom.r.1775-01-01-00000.nc"

• user_nl_cice
  ice_ic = "/cluster/shared/noresm/inputdata/restart/NOIIAJRAOC20TR_TL319_tn14_ppm_20240816/rest/1775-01-01-00000/NOIIAJRAOC20TR_TL319_tn14_ppm_20240816.cice.r.1775-01-01-00000.nc"

• user_nl_clm
  Add the following to user_nl_clm

finidat='/cluster/shared/noresm/inputdata/lnd/clm2/initdata_esmf/ctsm5.3/testing/1850_fates_spinup_new_dust.clm2.r.0081-01-01-00000.nc'
z0param_method='ZengWang2007'

• xmlchange commands

./xmlchange NTASKS=1920
./xmlchange NTASKS_OCN=256
./xmlchange ROOTPE=0
./xmlchange ROOTPE_OCN=1920
./xmlchange BLOM_VCOORD=cntiso_hybrid,BLOM_TURBULENT_CLOSURE=
./xmlchange STOP_OPTION=nyears
./xmlchange STOP_N=7
./xmlchange REST_N=1
./xmlchange REST_OPTION=nyears
./xmlchange RESUBMIT=3
./xmlchange --subgroup case.run JOB_WALLCLOCK_TIME=48:00:00
./xmlchange --subgroup case.st_archive JOB_WALLCLOCK_TIME=03:00:00
./xmlchange --append  CLM_BLDNML_OPTS="-ignore_warnings"

---

Compset long name
1850_CAM%DEV%LT%NORESM%CAMoslo_CLM60%BGC-CROP_CICE_BLOM%ECO_MOSART_DGLC%NOEVOLVE_SWAV_SESP

---

Model grid long name
a%ne30np4.pg3_l%ne30pg4.pg3_oi%tnx1v4_r%r05_w%null_z%null_g%gris4_m%tnx1v4

---

Case directory:

• Locally on betzy:
  /cluster/projects/nn9560k/mvertens/cases/n1850.ne30_tn14.hybrid_fatessp.202401007
• On github:

---

Code version on github:

• https://github.com/NorESMhub/NorESM/releases/tag/noresm2_5_alpha06

---

Diagnostics:

Atmosphere:

• https://ns2345k.web.sigma2.no/datalake/diagnostics/ADF/plots/n1850.ne30_tn14.hybrid_fatessp.202401007_5_14_vs_Obs/website/index.html
• https://ns2345k.web.sigma2.no/datalake/diagnostics/ADF/plots/n1850.ne30_tn14.hybrid_fatessp.202401007_15_25_vs_Obs/website/index.html

Ocean / sea ice:

• https://ns2345k.web.sigma2.no/datalake/diagnostics/noresm/adagj/n1850.ne30_tn14.hybrid_fatessp.202401007/

AeroVal:

• https://aeroval.met.no/pages/evaluation/?project=noresm&experiment=Autumn24cases&model=HYBRID-FATESSP.202401007#

---

Output:

• Raw output: /nird/datalake/NS9560K/noresm3/cases/ n1850.ne30_tn14.hybrid_fatessp.20241007
• CAM Climos: TODO

---

Contacts:
@mvertens, @adagj

---

Extra details:

[mvertens]

@adagj - please give this a try.

[adagj]

@adagj - please give this a try.

Will do!

[adagj]

To make this simulation run we needed to include (in addition to what Mariana already listed):

in the machine settings: $SRCROOT/ccs_config/machines/betzy/config_machines.xml
<env name="OMPI_MCA_coll_hcoll_enable">0</env>

and in the case direcory, before building
./xmlchange CLM_FORCE_COLDSTART=off
./xmlchange --append CLM_BLDNML_OPTS='-ignore_warnings'

I didn't include
ml swap hpcx/2.20

[oyvindseland]

I looked into some of the aerosol output / aerosol diagnostics for the simulation.
For many of the fields the results looks reasonable or at least as expected from earlier simulations. However this simulation does not have DMS. There seem to be a calculation of the surface flux, possibly coming from prescribed concentration fields but the output of odms (dms concentration in the upper ocean) is 0.

One large difference is that isoprene emissions from Fates-Megan are much lower than from CLM -Megan almost to the point of being negligible, at least in the tropics. Since the Fates output is not on a lat-lon grid I was not able to do a simple comparison, but the highest value in the land model output was comparable to the value in CAM. Was this large change in BVOC an expected outcome from Fates?

One thing I found a bit worrisome was that the grid peak negative emissions number for monoterpene was quite high (1 % of the highest "positive" emissions). I suppose this is just a matter of grid conversion.

When looking at DMS I also noticed the that there is a strong positive trend in Antarctic sea-ice over the simulation period but I presume that I am not the only one looking into that.

[rosiealice]

Hi @oyvindseland. I was -just- also looking into this isoprene issue. But your observation is interesting as I have been running with an updated branch where I modify the mapping to the FATES PFTs, and this code runs with the default. So it is interesting to see that the tiny isoprene fluxes also show up there...

[rosiealice]

I figured out the issue with the isoprene fluxes- it's that we aren't correctly passing the leaf internal CO2 from FATES into MEGAN. There are some residual differences in carene, thujene_a and methanol that I have not got to the bottom of yet, but they are nothing like as bad as the isoprene problem.

[rosiealice]

@oyvindseland the differences in the carene etc were from a mapping error. Just need to make the model pass internal leaf co2 into Megan and it should be ok. (This is slightly non trivial so I need to think about the right way to do it)

[maritsandstad]

@mvertens did this case output move/get deleted?

[mvertens]

@maritsandstad - the output is at /nird/datalake/NS9560K//noresm3/cases/n1850.ne30_tn14.hybrid_fatessp.202401007/