Skip to content

V10.19.3 ctm jk#40

Open
JulesKouatchou wants to merge 5 commits intomainfrom
v10.19.3-CTM_jk
Open

V10.19.3 ctm jk#40
JulesKouatchou wants to merge 5 commits intomainfrom
v10.19.3-CTM_jk

Conversation

@JulesKouatchou
Copy link

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.

Jules Kouatchou added 4 commits September 17, 2021 08:03
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
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