noresm3_0_041_cam6_4_121: Misc. bug fixes and additions

[gold2718]

Summary: Changes to radiative CO2 and diagnostic additions and bug fixes.

Contributors: @gold2718, @DirkOlivie

Reviewers: @DirkOlivie

Purpose of changes:

Github PR URL: #291

Changes made to build system: Adjustments to some use cases (see details)

Changes made to the namelist: None

Changes to the defaults for the boundary datasets: None

Substantial timing or memory changes: None

Details:

• Remove rad_climate from use cases.
  • Allow build_namelist to compute this
• Update sim_year to 2022 for cases with extended forcing files
• Update Oslo Aero tag with some additions and bug fixes

Testing:

• aux_cam_noresm on Olivia:
• aux_cam_noresm on Betzy:

Issues addressed by this PR: [ For each issue include a GitHub issue entry, one per line. ]

[DirkOlivie]

Hi @gold2718 , this looks fine to me. I however think that maybe we also have to change sim_year in namelist_defaults_cam.xml on various places - otherwise the 1850-2022 will nowhere be recognized.

[gold2718]

Hi @gold2718 , this looks fine to me. I however think that maybe we also have to change sim_year in namelist_defaults_cam.xml on various places - otherwise the 1850-2022 will nowhere be recognized.

@DirkOlivie, do you have a suggestion for what to do with files whose data ends in 2015? I have found these, there may be others:

atm/cam/ggas/emissions-cmip6_CO2_anthro_surface_175001-201512_fv_0.9x1.25_c20181011.nc
atm/cam/ggas/emissions-cmip6_CO2_anthro_ac_175001-201512_fv_0.9x1.25_c20181011.nc

[gold2718]

@DirkOlivie,

Removing rad_climate from the use cases resulted in a rad_climate which asks for N:VOLC_MMR, however, the pbuf for VOLC_MMR is not getting created. Do we want this species in our runs?

[DirkOlivie]

Hi @gold2718, related to the updated ozone file (this is the same info as in the google-document (overview document):

• new ozone fields are available : they are on derecho (I assume) in the directory : /glade/campaign/acom/acom-climate/shawnh/Forcing_Files/
• the file name is : ozone_strataero_b.e30_alpha08o.BHISTC_MTt4s.ne30_t232_wgx3.330_branch_1850_1850_2021_transient_zm_5day_c20260625.nc
• the file should be put in the directory (Olivia and Betzy) : inputdata/atm/cam/ozone_strataero
• in the use-case hist_camnor_lt_osloaero.xml the variable prescribed_ozone_file should be changed

[DirkOlivie]

Hi @gold2718, related to the CO2 surface and aircraft emission files which end in 2015 (this is the same info as in the google-document (overview document). These files can now be replaced by more recent files:

• these files are currently on Betzy in the directory /cluster/shared/noresm/inputdata/atm/cam/ggas
• NorESM3 aircraft file : emissions-cmip7_CO2_anthro_ac_CMIP_CEDS-CMIP-2025-04-18_185001-202312_0.5x0.5_c20250119.nc
• NorESM3 surface emissions : emissions-cmip7_CO2_anthro_surface_CMIP_CEDS-CMIP-2025-04-18_185001-202312_0.5x0.5_c20250119.nc
• These files are listed in cam/bld/namelist_files/namelist_defaults_cam.xml
• Simultaneously we will have to update co2flux_fuel_year_last and aircraft_co2_year_last: they can be put on 2023
• Simultaneously we will have to update the meshfiles : co2flux_fuel_meshfile and aircraft_co2_meshfile. There is currently no 0.5x0.5 meshfile on Olivia, however on Betzy there is currently (possibly one of these can be used)
  - 360x720_120830_ESMFmesh_c20210507_cdf5.nc
  - 360x720_120830_ESMFmesh_c20210507.nc
  - 360x720_120830_ESMFmesh_tweaked_latlons_c20260108.nc

[gold2718]

@DirkOlivie,

Please have a look at my recent changes. For the rad_climate issue, I added them back to the use cases but also added a case for when co2_cycle is on. This leads to the following rad_climate namelist values for our compsets:

NF1850:  'A:Q:H2O', 'N:O2:O2', 'N:CO2:CO2', 'N:ozone:O3', 'N:N2O:N2O', 'N:CH4:CH4', 'N:CFC11:CFC11', 'N:CFC12:CFC12'
NF1850esm:  'A:Q:H2O', 'N:O2:O2', 'A:CO2:CO2', 'N:ozone:O3', 'N:N2O:N2O', 'N:CH4:CH4', 'N:CFC11:CFC11', 'N:CFC12:CFC12'
NF1850esmGhg:  'A:Q:H2O', 'N:O2:O2', 'A:CO2:CO2', 'N:ozone:O3', 'A:N2O:N2O', 'A:CH4:CH4', 'N:CFC11STAR:CFC11',
NF1850ghg:  'A:Q:H2O', 'N:O2:O2', 'A:CO2:CO2', 'N:ozone:O3', 'A:N2O:N2O', 'A:CH4:CH4', 'N:CFC11STAR:CFC11',
NF2000:  'A:Q:H2O', 'N:O2:O2', 'N:CO2:CO2', 'N:ozone:O3', 'N:N2O:N2O', 'N:CH4:CH4', 'N:CFC11:CFC11', 'N:CFC12:CFC12'
NFHIST:  'A:Q:H2O', 'N:O2:O2', 'N:CO2:CO2', 'N:ozone:O3', 'N:N2O:N2O', 'N:CH4:CH4', 'N:CFC11:CFC11', 'N:CFC12:CFC12'

Do these look correct?

[gold2718]

Other changes I would like reviewed:

• Changes to forcing files in bld/namelist_files/use_cases/hist_camnor_lt_osloaero.xml
• Changes to co2flux_xxx and aircraft_co2_xxx in bld/namelist_files/namelist_defaults_cam.xml

[DirkOlivie]

@gold2718 I think it is a good idea to bring rad_climate back into the use-cases, and use the co2_cycle to distinguish between N:CO2:CO2 an A:CO2:CO2. When I looked into use-case files from the PR, it looks all fine (just in the comment section above were you list the 6 rad_climate changes, there is no CFC12 for N1850esmGhg and N1850ghg - but the code is correct).

[DirkOlivie]

@gold2718 the changes to the aircraft and CO2 emissions files, and to the forcings (historical ozone) are also fine.

[DirkOlivie]

Hi Steve, it all looks fine:

• the rad_climate changes
• the co2 emissions forcing changes
• the ozone field changes

[gold2718]

there is no CFC12 for N1850esmGhg and N1850ghg - but the code is correct).

Sorry, that was just a cut-and-past error from the atm_in files.

Just to double-check, is it correct that CO2 is radiatively active for GHG chemistry, even for concentration-driven runs?

[DirkOlivie]

Hi Steve, yes that is correct (CO2 is radiatively active for GHG chemistry, even for concentration-driven runs).