Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Implement O3v2 in Linoz #3047

Merged
merged 7 commits into from
Aug 29, 2019
Merged

Implement O3v2 in Linoz #3047

merged 7 commits into from
Aug 29, 2019

Conversation

tangq
Copy link
Contributor

@tangq tangq commented Jul 7, 2019

Previous PR #3042 has some issues that I couldn't push new commits to the remote branch. The issue was fixed by deleting and recreating the remote branch.

Implement the newly designed Linoz O3v2

  • Add ozone sink in the boundary layer with a fix e-folding
    decay to a typical ozone concentration in the boundary layer
  • Disable the tropospheric ozone rewrite to mozart climatology
  • The parameters controlling this new ozone sink are name list variables.
  • Turn off ozone dry deposition (Currently turned off in the
    test runs via user_nl_cam, need to do it by changing compset settings
    in future commit)

[Non-BFB] - Non Bit-For-Bit

tangq added 2 commits July 6, 2019 15:57
Implement the newly designed Linoz O3v2
* Add ozone sink in the boundary layer with a fix e-folding
decay to a typical ozone concentration in the boundary layer
* Disable the tropospheric ozone rewrite to mozart climatology
* Turn off ozone dry deposition (Currently turned off in the
test runs via user_nl_cam, need to do it by changing compset settings
in future commit)
* Need to change hard-coded parameters to name list variables
with future commit.

[Non-BFB] - Non Bit-For-Bit
Change the hard-coded Linoz tuning parameters to name list variables
* The tuning parameters are for the boundary layer ozone loss

[BFB] - Bit-For-Bit
[NML] - Namelist Changing
@tangq tangq added Atmosphere non-BFB PR makes roundoff changes to answers. labels Jul 7, 2019
@tangq tangq added this to the v2.0alpha milestone Jul 7, 2019
@tangq tangq requested review from wlin7 and golaz July 7, 2019 23:51
The new Linoz O3v2 shares the tracer index with real O3. This requires
turning off the dry deposition for the redesigned Linoz O3v2's loss
in the boundary layer to work.

When the interactive tropospheric chemistry is implemented, the ozone
dry deposition should be turned on again.

[BFB] - Bit-For-Bit
[NML] - Namelist Changing
The new Linoz O3v2 shares the tracer index with real O3. This requires
turning off the dry deposition for the redesigned Linoz O3v2's loss
in the boundary layer to work.

When the interactive tropospheric chemistry is implemented, the ozone
dry deposition should be turned on again.

[BFB] - Bit-For-Bit
[NML] - Namelist Changing
@tangq
Copy link
Contributor Author

tangq commented Jul 12, 2019

After the initial implementation, the following up commits 1) change the hard-coded parameters to name list variables enabling more flexible tuning in the future, and 2) turn off the ozone dry deposition by compset settings for all CMIP6 compsets on the master.

@rljacob
Copy link
Member

rljacob commented Jul 25, 2019

Where is the code review document for this feature?

@tangq
Copy link
Contributor Author

tangq commented Jul 31, 2019

@rljacob rljacob assigned singhbalwinder and unassigned wlin7 Aug 14, 2019
@@ -5454,6 +5454,24 @@ if <varname>linoz_data_type</varname> is 'FIXED'.
Default: 0 seconds
</entry>

<entry id="linoz_lbl" type="integer" category="linoz"
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

It will be helpful to list here the actual default values -- as specified in build-namelist for linoz_lbl, linoz_tau and linoz_sfc.

Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Good catch @wlin7 . I will make these changes when cori is back.

Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

The default O3v2 parameter values are now added in the name list definitions.

@wlin7
Copy link
Contributor

wlin7 commented Aug 26, 2019

The changes look good. Special attention paid to the scoping of the dummy arrays involved in the new subroutine "lin_strat_sfcsink" and all is good.

One more comment to the drydep_list nml. Ideally, the change to the drydep list should be realized via build-namelist (which would invoke set_dep_lists in ChemNamelist). That said, the current format is ok for the handful of compsets in use.

@wlin7
Copy link
Contributor

wlin7 commented Aug 26, 2019

Hi @tangq , another note concerning the drydep_list as eventually specified in drv_flds_in. What is the consequence if "O3" remains in the nml file after the new O3v2 is in used? Can the model still run properly? It may happen to compsets not covered in this PR.

Change needs to be made to one other file at least: 2000_cam5_av1c-04p2.xml. This is the use case file for FC5AV1C-04P2. I know AV1C- compsets are deprecated and no need to maintain. But FC5AV1C-L is still used for many of the nightly tests. It is an alias to FC5AV1C-04P2. Until another compset is selected for the tests, the change to 2000_cam5_av1c-04p2.xml is needed.

@tangq
Copy link
Contributor Author

tangq commented Aug 26, 2019

It will give the wrong results for ozone. I will add the change to 2000_cam5_av1c-04p2.xml when cori's back.

This drydep change won't be needed till the interactive tropospheric chemistry is implemented later.

@tangq
Copy link
Contributor Author

tangq commented Aug 27, 2019

@wlin7 , I added the changes you suggested yesterday - name list default values and turning off ozone dry deposition in 2000_cam5_av1c-04p2.xml. It should be ready to go now.

@@ -162,6 +162,9 @@
<linoz_data_path >atm/cam/chem/trop_mozart/ub</linoz_data_path>
<linoz_data_type >CYCLICAL</linoz_data_type>

<!-- Turn off ozone dry deposition, as Linoz O3v2 and ozone are not separated for now. Need to turn on ozone dry deposition when interactive tropospheric chemistry is implemented -->
<drydep_list >'H2O2', 'H2SO4', 'SO2'</drydep_list>
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Reminder: this change will also be needed for the several F-2010 compsets, either via this PR or a future one, depending on the order of merging them to master.

@wlin7
Copy link
Contributor

wlin7 commented Aug 27, 2019

Thanks @tangq . All is good now. F2010SC5-CMIP6 is in the master now. I can add the change to drydep_list for it after merging this PR to master.

@golaz , can we have your approval to get this PR in? Thanks.

@rljacob
Copy link
Member

rljacob commented Aug 27, 2019

wlin7 added a commit that referenced this pull request Aug 28, 2019
Implement O3v2 in Linoz

The following elements are added:

 . Add ozone sink in the boundary layer with a fix e-folding
decay to a typical ozone concentration in the boundary layer

 . Disable the tropospheric ozone rewrite to mozart climatology

 . The parameters controlling this new ozone sink are name list variables.

 . Turn off ozone dry deposition (Currently turned off in the
test runs via user_nl_cam, need to do it by changing compset settings
in future commit)

[Non-BFB] - Non Bit-For-Bit
[NML] drv_flds_in changed on the list for drydep_list.
@wlin7
Copy link
Contributor

wlin7 commented Aug 28, 2019

Merged to next.

wlin7 added a commit that referenced this pull request Aug 29, 2019
Implement O3v2 in Linoz

The following elements are implemented in the newly designed Linoz O3v2:

. Add ozone sink in the boundary layer with a fix e-folding
decay to a typical ozone concentration in the boundary layer

. Disable the tropospheric ozone rewrite to mozart climatology

. The parameters controlling this new ozone sink are name list variables.

. Turn off ozone dry deposition (Currently turned off in the
test runs via user_nl_cam, need to do it by changing compset settings
in future commit)

[Non-BFB] - Non Bit-For-Bit
[NML] drydep_list in drv_flds_in changed
@wlin7 wlin7 merged commit 031b28f into master Aug 29, 2019
@wlin7
Copy link
Contributor

wlin7 commented Aug 29, 2019

Merged to master. Requests submitted to bless DIFFs: sandistoss3 (38), melvin(7), cori-knl(2). Two FAILs on anvil will need to be blessed after inputdata issue is resolved.

  • SMS_Ld2.ne30_oECv3_ICG.A_WCYCL1850S_CMIP6.anvil_intel.allactive-v1cmip6
  • SMS_Ln5.ne30_ne30.FC5AV1C-L.anvil_intel.cam-cosplite

@wlin7
Copy link
Contributor

wlin7 commented Sep 5, 2019

Request submitted to bless the anvil DIFFs, which were FAILs previously due to machine issue. Note that PR #3057 could also cause DIFF for SMS_Ld2.ne30_oECv3_ICG.A_WCYCL1850S_CMIP6.anvil_intel.allactive-v1cmip6.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Atmosphere non-BFB PR makes roundoff changes to answers.
Projects
None yet
Development

Successfully merging this pull request may close these issues.

None yet

4 participants