Add support for v3 1950 compsets #6053

wlin7 · 2023-11-06T00:08:21Z

This PR adds supporting files for v3 1950 compsets.

Volcanic SO2 emission is 1950 eruption plus the historical average
Other surface and upper level forcings set to cycle over year 1950 of the historical forcings

[BFB] Non-1950 tests are not affected.

Notes:

The 1950 ne30pg2 surfdata with TOP file currently contains place holder for variables related to TOP scheme.
1950 for ne120 still needs IC file for atm.
With this PR, F1950 can run with L80 v3 config. B-case 1950 is expected to also be able to run.

github-actions · 2023-11-06T00:10:03Z

PR Preview Action v1.4.7
🚀 Deployed preview to https://E3SM-Project.github.io/E3SM/pr-preview/pr-6053/
on branch `gh-pages` at 2024-02-18 05:07 UTC

wlin7 · 2023-11-06T00:20:38Z

@hwangacme , you created forcings for use by v2 1950 compsets. The forcings seem to be a subset of historical forcings (1950-1960). Should we still use the same for v3? The v3 1850 compsets set to cycle over a single year of 1850 for forcings. I am using the same format for 1950 for now.

@keziming , regarding the volcanic SO2 emission, I am using the average of 1850-2014 as you created for piControl. The time variables are modified to align with 1950 compsets (to allow a single cycling year for all elevated forcings). Do you feel there is a need to use an average around 1950 instead? In V2, it was set to cycle over a single year of 1950 from a specified 1950-1960 file created by Hailong.

hwangacme · 2023-11-06T04:14:40Z

@hwangacme , you created forcings for use by v2 1950 compsets. The forcings seem to be a subset of historical forcings (1950-1960). Should we still use the same for v3? The v3 1850 compsets set to cycle over a single year of 1850 for forcings. I am using the same format for 1950 for now.

@keziming , regarding the volcanic SO2 emission, I am using the average of 1850-2014 as you created for piControl. The time variables are modified to align with 1950 compsets (to allow a single cycling year for all elevated forcings). Do you feel there is a need to use an average around 1950 instead? In V2, it was set to cycle over a single year of 1950 from a specified 1950-1960 file created by Hailong.

@wlin7 I didn't make the decision of which forcing files to be used for the 1950 compset, but I think it makes sense to take the single year 1950 forcing from the historical 1850-2014 forcing files.

keziming · 2023-11-10T22:47:19Z

@hwangacme , you created forcings for use by v2 1950 compsets. The forcings seem to be a subset of historical forcings (1950-1960). Should we still use the same for v3? The v3 1850 compsets set to cycle over a single year of 1850 for forcings. I am using the same format for 1950 for now.
@keziming , regarding the volcanic SO2 emission, I am using the average of 1850-2014 as you created for piControl. The time variables are modified to align with 1950 compsets (to allow a single cycling year for all elevated forcings). Do you feel there is a need to use an average around 1950 instead? In V2, it was set to cycle over a single year of 1950 from a specified 1950-1960 file created by Hailong.

@wlin7 I didn't make the decision of which forcing files to be used for the 1950 compset, but I think it makes sense to take the single year 1950 forcing from the historical 1850-2014 forcing files.

@wlin7 there was a volcanic eruption during 1950s (1956 I guess). In general, the volcanic activities in 1950s are similar to 1850-2014 average, not too strong or weak. Do you mean you plan to create a new SO2 external emission file? If so, it would be OK to use piControl VOLC sector plus all other SO2 external emission sectors in 1950s. If I remember correctly, there is no eruption in year 1950 alone.

keziming

For file cmip6_1850_2014_avg_so2_elev_strat_1950_c20230310.nc, I need a clarification of its sectors. All sectors but volcanic in this file should be 1950s average, right? Thanks.

keziming · 2023-11-10T22:55:32Z

components/eam/bld/namelist_files/use_cases/1950_eam_CMIP6_chemUCI-Linoz-mam5-vbs.xml

+<ext_frc_cycle_yr     >1950</ext_frc_cycle_yr>
+<no2_ext_file         >atm/cam/chem/trop_mozart_aero/emis/chem_gases/2degrees/emissions-cmip6_NO2_aircraft_vertical_1750-2015_1.9x2.5_c20170608.nc </no2_ext_file>
+<!--  Use 1850-2014 averaged volcanic SO2 emission for 1950 -->
+<so2_ext_file         >atm/cam/chem/trop_mozart_aero/emis/DECK_ne30/cmip6_1850_2014_avg_so2_elev_strat_1950_c20230310.nc </so2_ext_file>


@wlin7 It may need some clarifications. SO2 external forcing contains several sectors of emission, volcanic emission is only a part of it. In this file, ideally, it should be 1950s emissions from other sectors in addition to 1850-2014 averaged volcanic emission.

Thanks for carefully reviewing the settings, @keziming . In v2, the SO2 external forcing for 1950 compset appears to be cycling over a single year 1950 instead of an average. (file atm/cam/chem/trop_mozart_aero/emis/DECK_ne120/cmip6_mam4_so2_elev_1950-1960_c180203.nc). We can also use single year 1950 or an average around 1950 for non-volcanic SO2 emissions. Can you please help to update those other sectors? We should have some discussion on 1950 forcings in general in next atm group meeting.

tangq

The changes look good to me.

darincomeau · 2023-11-13T21:50:39Z

Thanks @wlin7 for putting this together - greatly appreciated!

Passes SMS.ne30pg2_EC30to60E2r2.WCYCL1950.anvil_intel and SMS.ne30pg2_ECwISC30to60E2r1.CRYO1950.anvil_intel

I tried a tri-grid configuration, but SMS.ne30pg2_r05_ECwISC30to60E3r2.WCYCL1950.anvil_intel failed on setup with

ERROR : CLM build-namelist::ELMBuildNamelist::add_default() : No default value found for fsurdat.

Is tri-grid meant to be supported? The same grid but with WCYCL1850 works.

So far I've only run the short SMS tests, but if you'd like a longer test where we can look at output, happy to do that.

darincomeau · 2023-11-13T21:51:09Z

Also, here's the issue this is fixing: #6016

wlin7 · 2023-11-13T22:01:05Z

@darincomeau , thanks for testing and nice to see the bigrid test passes. The data for tri-grid r05, simyr1950 has never been established. I just added it to the list of files needing to be updated/created.

darincomeau · 2023-11-13T22:03:32Z

Ok great, good to know. Thanks again @wlin7 !

rljacob · 2023-12-05T00:01:11Z

@wlin7 how is the creation of those files going? This PR needs to be done be the v3.0.0 tag is made.

rljacob · 2023-12-07T18:22:06Z

Waiting on @keziming to provide forcing files.

rljacob · 2024-01-04T18:28:53Z

@keziming are the forcing files ready?

keziming · 2024-01-04T19:17:24Z

@keziming are the forcing files ready?

I plan to finished it before this weekend. Sorry for the delay.
update:
I created a new file [so2_elev_strat_1950.nc] under /lcrc/group/e3sm/data/inputdata/atm/cam/chem/trop_mozart_aero/emis/DECK_ne30
What I did:
I grabed 1949-1951 period SO2 external emission from [cmip6_mam4_so2_elev_1850-2014_c180205.nc] (V2 SO2 historical emission file)
Then I added historical averaged eruptive volcano emission (used for piControl simulation) and modified the vertical layer coordinates.
if our simulation use "cycle 1950", it should serve the purpose.
please give a test run with the 1950 compset.

rljacob · 2024-01-19T18:58:40Z

@wlin7 can you check that this works with the new files @keziming made?

hwangacme

Note that the tracer_cnst_cycle_yr is cycling the year 1955 input. I cannot recall why it was set this way.

rljacob · 2024-01-29T17:52:48Z

@wlin7 please confirm this is ready.

rljacob · 2024-02-01T18:35:53Z

Notes: @wlin7 will test new files.

Add support for v3 1950 compsets This PR adds supporting files for v3 1950 compsets. Volcanic SO2 emission is 1950 eruption plus the historical average Other surface and upper level forcings set to cycle over year 1950 of the historical forcings Fixes #6016. [BFB] Non-1950 tests are not affected. * origin/wlin/atm/v3_1950_compset: Update elevated so2 forcing for lnd config for 1950 compsets Correct 1950 DMS file Add support for v3 1950 compsets

wlin7 · 2024-02-22T01:51:16Z

@keziming , this compset is now using the updated SO2 emission for 1950 you produced (so2_elev_strat_1950.nc). Thanks.

Merged to next.

wlin7 · 2024-02-22T20:04:27Z

Merged to master.

keziming · 2024-02-22T20:17:36Z

@wlin7 great

zzzeamon · 2024-05-14T03:04:15Z

@keziming , this compset is now using the updated SO2 emission for 1950 you produced (so2_elev_strat_1950.nc). Thanks.

Merged to next.

could you please upload so2_elev_strat_1950.nc to "https://web.lcrc.anl.gov/public/e3sm/inputdata/atm/cam/chem/trop_mozart_aero/emis/DECK_ne30/" ?

keziming · 2024-05-15T22:18:41Z

@keziming , this compset is now using the updated SO2 emission for 1950 you produced (so2_elev_strat_1950.nc). Thanks.
Merged to next.

could you please upload so2_elev_strat_1950.nc to "https://web.lcrc.anl.gov/public/e3sm/inputdata/atm/cam/chem/trop_mozart_aero/emis/DECK_ne30/" ?

just copied file to https://web.lcrc.anl.gov/public/e3sm/inputdata/atm/cam/chem/trop_mozart_aero/emis/DECK_ne30/

zzzeamon · 2024-05-16T01:26:43Z

@keziming , this compset is now using the updated SO2 emission for 1950 you produced (so2_elev_strat_1950.nc). Thanks.
Merged to next.

could you please upload so2_elev_strat_1950.nc to "https://web.lcrc.anl.gov/public/e3sm/inputdata/atm/cam/chem/trop_mozart_aero/emis/DECK_ne30/" ?

just copied file to https://web.lcrc.anl.gov/public/e3sm/inputdata/atm/cam/chem/trop_mozart_aero/emis/DECK_ne30/

thank you.

Add support for v3 1950 compsets

6fce80d

wlin7 added Atmosphere compset BFB PR leaves answers BFB Atmospheric chemistry labels Nov 6, 2023

wlin7 requested review from darincomeau, hwangacme and keziming November 6, 2023 00:09

rljacob assigned wlin7 Nov 6, 2023

Correct 1950 DMS file

a2443d6

wlin7 mentioned this pull request Nov 7, 2023

Fix tests incompatitle with L80 and restore tentative config options for SSP #6054

Merged

wlin7 requested a review from tangq November 10, 2023 17:56

keziming reviewed Nov 10, 2023

View reviewed changes

tangq approved these changes Nov 11, 2023

View reviewed changes

darincomeau approved these changes Nov 13, 2023

View reviewed changes

rljacob added this to the v3.0.0-rc milestone Jan 8, 2024

hwangacme approved these changes Jan 22, 2024

View reviewed changes

keziming approved these changes Feb 1, 2024

View reviewed changes

This was referenced Feb 11, 2024

Add new SOwISC12to30E3r2 ocean and sea-ice mesh E3SM-Ocean-Discussion/E3SM#72

Closed

Add prescribed tidal effect on ice shelf basal melting E3SM-Ocean-Discussion/E3SM#76

Closed

Update elevated so2 forcing for lnd config for 1950 compsets

067211b

wlin7 merged commit 9dcf200 into master Feb 22, 2024
10 checks passed

wlin7 deleted the wlin/atm/v3_1950_compset branch February 22, 2024 20:03

rljacob modified the milestones: v3.0.0-rc1, v3.0.0-rc2, v3.0.0 Feb 29, 2024

Provide feedback

Saved searches

Use saved searches to filter your results more quickly

Add support for v3 1950 compsets #6053

Add support for v3 1950 compsets #6053

wlin7 commented Nov 6, 2023 •

edited

Loading

github-actions bot commented Nov 6, 2023 •

edited

Loading

wlin7 commented Nov 6, 2023

hwangacme commented Nov 6, 2023

keziming commented Nov 10, 2023

keziming left a comment

keziming Nov 10, 2023

wlin7 Nov 13, 2023 •

edited

Loading

tangq left a comment

darincomeau commented Nov 13, 2023

darincomeau commented Nov 13, 2023

wlin7 commented Nov 13, 2023 •

edited

Loading

darincomeau commented Nov 13, 2023

rljacob commented Dec 5, 2023

rljacob commented Dec 7, 2023 •

edited

Loading

rljacob commented Jan 4, 2024

keziming commented Jan 4, 2024 •

edited

Loading

rljacob commented Jan 19, 2024

hwangacme left a comment

rljacob commented Jan 29, 2024

rljacob commented Feb 1, 2024

wlin7 commented Feb 22, 2024

wlin7 commented Feb 22, 2024

keziming commented Feb 22, 2024

zzzeamon commented May 14, 2024

keziming commented May 15, 2024

zzzeamon commented May 16, 2024

Add support for v3 1950 compsets #6053

Add support for v3 1950 compsets #6053

Conversation

wlin7 commented Nov 6, 2023 • edited Loading

github-actions bot commented Nov 6, 2023 • edited Loading

wlin7 commented Nov 6, 2023

hwangacme commented Nov 6, 2023

keziming commented Nov 10, 2023

keziming left a comment

Choose a reason for hiding this comment

keziming Nov 10, 2023

Choose a reason for hiding this comment

wlin7 Nov 13, 2023 • edited Loading

Choose a reason for hiding this comment

tangq left a comment

Choose a reason for hiding this comment

darincomeau commented Nov 13, 2023

darincomeau commented Nov 13, 2023

wlin7 commented Nov 13, 2023 • edited Loading

darincomeau commented Nov 13, 2023

rljacob commented Dec 5, 2023

rljacob commented Dec 7, 2023 • edited Loading

rljacob commented Jan 4, 2024

keziming commented Jan 4, 2024 • edited Loading

rljacob commented Jan 19, 2024

hwangacme left a comment

Choose a reason for hiding this comment

rljacob commented Jan 29, 2024

rljacob commented Feb 1, 2024

wlin7 commented Feb 22, 2024

wlin7 commented Feb 22, 2024

keziming commented Feb 22, 2024

zzzeamon commented May 14, 2024

keziming commented May 15, 2024

zzzeamon commented May 16, 2024

wlin7 commented Nov 6, 2023 •

edited

Loading

github-actions bot commented Nov 6, 2023 •

edited

Loading

wlin7 Nov 13, 2023 •

edited

Loading

wlin7 commented Nov 13, 2023 •

edited

Loading

rljacob commented Dec 7, 2023 •

edited

Loading

keziming commented Jan 4, 2024 •

edited

Loading