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

V10.19.3 ctm jk #40

Open
wants to merge 5 commits into
base: main
Choose a base branch
from
Open

V10.19.3 ctm jk #40

wants to merge 5 commits into from

Conversation

JulesKouatchou
Copy link
Contributor

Made the following changes in CTM Env to align the Icarus version of the code with the Jason one:

  • Added an option to produce the courant numbers and mass fluxes needed for Advection.
    - The flag output_forcingData needs to be set to true in CTM_GridComp.rc.
    - The fields MFX, MFY, CX and CY should be included in HISTORY.rc
  • Added an option to read the courant number and mass fluxes.
    - The flag read_advCoreFields needs to be set to true in CTM_GridComp.rc
  • CTM Convection was refactored.
  • Added an option to do RAS calculations.
    • The flag enable_rasCalculations needs to be set to true in CTM_GridComp.rc
  • Added an option to do Advection or not.
    • If choosing not to do Advection, the flag do_ctmAdvection needs to be set to false in CTM_GridComp.rc
  • CTM Env was refactored to add an internal state and to include calculations to derive fields into subroutines.
  • CTM GridComp was edited to include an internal state.
  • The atmospheric pressure is now calculated using the surface pressure and AK and BK. This is the most significant change among all the modifications.
  • Included the setting Ext_AllowExtrap: .TRUE. (in the ExtData rc file templates) to allow reading data files that have one record only.

The following tests were done to validate the changes:

  • Ran the Idealized Passive Tracers configuration using MERRA2 forcings and showed that there is total mass conservation.
  • Compared the regridding outputs of this version of the code and that of Icarus and showed that both versions of the code produce the same total budget.
  • Did a one-month TR+GOCART experiment.
  • Ran the Idealized Passive Tracers configuration using FP collections where the courant numbers (CX, CY), mass fluxes (MFX, MFY), surface pressure (PS) and specific humidity (Q) are read on the cubed-sphere grid and the remaining fields on the lat-lon grid. The results showed mass conservation.

The following modifications were made:

- CTM Env has the option to produce the courant numbers and mass fluxes needed for Advection.
     - The flag output_forcingData needs to be set to true in CTM_GridComp.rc.
     - The fields MFX, MFY, CX and CY should be included in HISTORY.rc
- CTM Env has the option to read the courant number and mass fluxes.
     - The flag read_advCoreFields needs to be set to true in CTM_GridComp.rc
- CTM Convection was refactored.
- CTM Env has the option to do RAS calculations.
     - The flag enable_rasCalculations needs to be set to true in CTM_GridComp.rc
- CTM has the option to do Advection or not.
     - If choosing not to do Advection, the flag do_ctmAdvection needs to be set to false in CTM_GridComp.rc
- CTM Env was refactored to add an internal state and to include calculations to derive fields into
  subroutines.
- CTM GridComp was edited to include an internal state.
- In CTM Env, the atmospheric pressure is now calculate using the surface pressure and AK and BK.
  This is the most significant change among all the modifications.
@JulesKouatchou JulesKouatchou requested review from a team as code owners December 2, 2021 13:36
Copy link
Collaborator

@tclune tclune left a comment

Choose a reason for hiding this comment

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

cmake changes ok

USE MAPL
USE GmiConvectionMethod_mod ! GMI Convection component
USE GenericConvectionMethod_mod ! Generic Convection component
USE MAPL_Mod
Copy link
Member

Choose a reason for hiding this comment

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

@JulesKouatchou Is there a reason you went back to MAPL_Mod? @tclune and the SI Team moved GEOS to use MAPL as the default in GEOS.

Copy link
Member

Choose a reason for hiding this comment

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

Note there are other moves to use MAPL_Mod in other files in this PR as well

Copy link
Contributor Author

Choose a reason for hiding this comment

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

  • MAPL_Mod can be replaced by MAPL.
  • The changes I do not produce zero-diff because the atmospheric pressure is now calculated instead of being read from file.

- Fully implmented the internal state in CTM_ConvectionGridCompMod.F90
- Removed GenericConvectionMethod_mod.F90 and GmiConvectionMethod_mod.F90
  from CMakeLists.txt. The two files do not need to be compiled at this time.
  The file GenericConvectionMethod_mod.F90 might be completely removed.
  GmiConvectionMethod_mod.F90 needs to be kept and used if we decide to
  implement GMI specific calculations.
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

3 participants