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

[WIP] Add capability to run GEOS-Chem at 0.125x0.15625 resolution #1980

Draft
wants to merge 7 commits into
base: main
Choose a base branch
from

Commits on Oct 6, 2023

  1. Add initial modifications for running GEOS-Chem at 0.125x0.15625 reso…

    …lution
    
    NASA GMAO has provided 1-hourly mass flux fields for GEOS-FP meteorology
    at c720 resolution to drive GCHP. The mass flux fields have been converted
    to derived winds fields and further regridded from c720 to 0.125x0.15625
    lat-lon resolution. The new files include the following collections and met
    fields:
    
    - GEOS.fp.asm.I1dyn : PS, SPHU
    - GEOS.fp.asm.A1dyn : U, V
    
    Signed-off-by: Melissa Sulprizio <mpayer@seas.harvard.edu>
    msulprizio committed Oct 6, 2023
    Configuration menu
    Copy the full SHA
    a093dfc View commit details
    Browse the repository at this point in the history

Commits on Nov 21, 2023

  1. Expand IF statements for GridRes to include 0.125x0.15625 resolution

    Several IF statements have been expanded to include State_Grid%GridRes =
    0.125x0.15625 resolution so that grid parameters and scale factors are
    properly determined.
    
    Also updated createRunDir.sh for GCClassic to correctly set transport and
    chemistry timesteps for 0.125x0.15625 resolution. For now the same default
    timesteps as 0.25x0.3125 and 0.5x0.625 will be used (transport=5 min,
    chemistry=10 min). This recommendation may later changed based on proper
    evaluation of the 0.125x0.15625 simulations.
    
    Signed-off-by: Melissa Sulprizio <mpayer@seas.harvard.edu>
    msulprizio committed Nov 21, 2023
    Configuration menu
    Copy the full SHA
    4d77fcb View commit details
    Browse the repository at this point in the history
  2. Modify HEMCO_Config.rc.gmao_metfields_0125 to use 0.25x0.3125 NA file…

    …s for non-dynamic fields
    
    The 0.125x0.15625 files are only available for dynamic collections A1dyn
    and I1dyn. The remaining met fields when running at this resolution should
    be read at 0.25x0.3125. This file is currently hardcoded to use the North
    America (NA) files at 0.25x0.3125 resolution for testing, but this should
    be made more flexible in the future.
    
    The 0.125x0.15625 files have also been updated here so that daily files are
    used instead of hourly files. The use of daily 24-hour met field files has
    already been implemented in GEOS-Chem, but it does not currently know how
    to handle hourly files.
    
    Signed-off-by: Melissa Sulprizio <mpayer@seas.harvard.edu>
    msulprizio committed Nov 21, 2023
    Configuration menu
    Copy the full SHA
    d22b361 View commit details
    Browse the repository at this point in the history

Commits on Mar 26, 2024

  1. Merge remote-tracking branch 'main' into feature/0.125x0.15625_resolu…

    …tion
    
    Signed-off-by: Melissa Sulprizio <mpayer@seas.harvard.edu>
    msulprizio committed Mar 26, 2024
    Configuration menu
    Copy the full SHA
    b5a8e4e View commit details
    Browse the repository at this point in the history

Commits on Apr 22, 2024

  1. Add fixes for 0.125x0.15625 grid from Xiaolin Wang

    Xiaolin Wang identified the following bugs in the 12km simulations:
    
    1. In GeosCore/pjc_pfix_window_mod.F90 for 12km code, we need to change the
       parameter dp in correspond with 12km resolution.
    2. In the input_mod.F90, we need to improve floating point precision,
       otherwise the longitude resolution is 0.1563 rather than 0.15625.
    
    Signed-off-by: Melissa Sulprizio <mpayer@seas.harvard.edu>
    msulprizio committed Apr 22, 2024
    Configuration menu
    Copy the full SHA
    cb71dfa View commit details
    Browse the repository at this point in the history

Commits on Jul 17, 2024

  1. Add fixes to read 0.125x0.15625 met fields at proper frequency from X…

    …iaolin Wang
    
    The PS and SPHU fields at 0.125x0.15625 resolution are in the I1dyn
    collection, while the U and V fields are read from the A1dyn collection.
    All of these fields should be read hourly, not 3-hourly as is the case
    for 0.25x0.3125 resolution and coarser. New subroutines FlexGrid_Read_A1dyn,
    FlexGrid_Read_I1dyn_1, FlexGrid_Read_I1dyn_2, and Copy_I1dyn_Fields have
    been added to flexgrid_mod.F90 following work by Xiaolin Wang. The code
    in hco_interface_gc_mod.F90 and main.F90 that calls these routines has
    also been updated. New I1dyn routines have also been added to time_mod.F90.
    
    Signed-off-by: Melissa Sulprizio <mpayer@seas.harvard.edu>
    msulprizio committed Jul 17, 2024
    Configuration menu
    Copy the full SHA
    c482454 View commit details
    Browse the repository at this point in the history
  2. Merge branch 'origin/main' (14.4.1) into feature/0.125x0.15625_resolu…

    …tion
    
    Signed-off-by: Melissa Sulprizio <mpayer@seas.harvard.edu>
    msulprizio committed Jul 17, 2024
    Configuration menu
    Copy the full SHA
    f23487f View commit details
    Browse the repository at this point in the history