Skip to content

2020 MARBL Dev team meetings

Jump to bottom
Michael Levy edited this page Jun 29, 2020 · 56 revisions

June 30, 2020

POP Software Updates

  1. Upcoming CESM 2.2 freeze: need to figure out order of POP tags
    1. New tunings
    2. JRA / BGC high-res run
    3. Qing's entrainment update
  2. Added complication, the entrainment update may be more than just round-off level changes
  3. My preferred path forward
    1. Qing's entrainment update but keep default scheme in place (langmuir_opt = 'vr12-ma')
    2. New tunings (need to verify that above PR doesn't require re-tuning: another cycle? Shorter run?)
    3. JRA / BGC high-res run
      • Hoping to run a couple of 1-month simulations (one using new compset out of the box, other configured for our experiment)
      • Do we need to figure out langmuir_opt first, or is that only going to affect the 1 degree?
    4. If langmuir_opt = 'lf17' is back on the table, it should be re-tested after the tuning update
      • Listed last only because I'm uncertain if it's necessary; can be done ahead of high-res compset (that may actually be preferable)
  4. CESM 2.1.4 release still needs a few things from POP
    1. Update namelist defaults to use cdf5 files rather than netcdf-4 files
    2. Update dt_count for SSP extension compsets
    3. I'm hoping to avoid thinking about these until after my MOM6 webinar talk in August (the first is actually ready to be merged, but the second may need a little more testing); I think the 2.1.4 code freeze won't happen until after the 2.2.0 release, but I'm not 100% certain about that.

MARBL Software Updates

  1. POP PR for updated tunings is waiting on corresponding MARBL PR
  2. Will also need to update the stable branch

MOM Software Updates

  1. Can run a full month with reasonable surface forcings (includes using T and S from the model physics) with correct surface values of tracers
  2. Starting to put together talk for MOM6 webinar, but mostly focused on CESM 2.2 release

General discussion

No Meeting

  • June 16, 2020: CESM Workshop

June 2, 2020

POP Software Updates

  1. Status for CESM 2.2 release

    • Plans call for several new tags

      1. WW3 entrainment update (Alper's responsibility)
      2. iron flux forcing bug (see below)
      3. new tunings for BGC (Kristen is waiting on iron flux forcing bug fix)
      4. new compset for high-res w/ BGC (see below; will include new tunings from Kristen but I have other aspects to attend to as well)
      5. bug in selecting dt_count default (POP issue #28)

    • Open PR: update iron flux forcing (tied to PR in marbl-forcing)

    • high-res compset

      • Current definition:

        <compset>
  <!-- latest JRA forcing, ecosys, high-res -->
  <alias>GIAFECO_JRA_HR</alias>
  <lname>2000_DATM%JRA-1p4-2018_SLND_CICE%CICE4_POP2%ECO_DROF%JRA-1p4-2018_SGLC_SWAV</lname>
</compset>

      • The existing eco + interannual forcing compset is

        <compset>
  <!-- latest JRA forcing -->
  <alias>G1850ECOIAF_JRA</alias>
  <lname>1850_DATM%JRA-1p4-2018_SLND_CICE_POP2%ECO_DROF%JRA-1p4-2018_SGLC_WW3</lname>
</compset>

      • I think that means we really want our compset to be

        <compset>
  <!-- latest JRA forcing, ecosys, high-res -->
  <alias>G1850ECOIAF_JRA_HR</alias>
  <lname>1850_DATM%JRA-1p4-2018_SLND_CICE%CICE4_POP2%ECO_DROF%JRA-1p4-2018_SGLC_SWAV</lname>
</compset>

  2. JRA_HR w/ MARBL

    • using three autotrophs (no coccolithophores yet), seeing 0.77 SYPD (260k pe-hrs / simulated_year) on largest task count Alper provided:

      <decomp nproc="7507" res="tx0.1v3" >
  <maxblocks >1</maxblocks>
  <bsize_x   >25</bsize_x>
  <bsize_y   >32</bsize_y>
  <decomptype>spacecurve</decomptype>
</decomp>

    • Number above is a single day run with no output (I also ran for two days to verify initialization isn't included).

    • Just got text file from him outlining other task counts to try

    • Do we have a target SYPD? CICE is running at 1.7 SYPD (23 nodes); would need to increase task count there as well to get any faster

MOM Software Updates

  1. Still struggling with two issues in surface flux forcing

    1. With surface tracer values set to 0, I tried to set forcing fields to the following (unlisted forcings set to 0):

      u10_sqr = 2.5e5
atm_press = 1
xco2 = 284.7
xco2_alt_co2 = 284.7
sss = 35

      But the run crashes during day 5; setting sss = 0 instead lets me run for a full month.

    2. Setting surface tracer values to "true" values (CS%tr(i,j,1,m)) causes run to crash in day 7 (assuming sss=0`)

General discussion

MARBL Software Updates

No Meeting

  • May 19, 2020: Matt unavailable

May 5, 2020

General discussion

  1. JRA high-res
    • Trying to track progress in real-time on Zulip
    • CICE and CIME pull requests handle a few minor issues in those components
    • Run is throwing errors in POP
      1. NaN in tracer tendencies was traced to bad copy of DZT with partial bottom cells
      2. something in PFTs? NaN in diazChl but really small values in all PFT fields... (see #50) -- few possibilities for fixing this to discuss:
        • Keith's suggestion (if so, can we remove PAR_threshold?)
        • Move PAR_threshold to settings file, and make it resolution-dependent?

MARBL Software Updates

POP Software Updates

MOM Software Updates

  1. Cleaned up configuration of MARBL per call with GFDL folks

  2. Pulled MARBL out of submodules

    • This caused Travis CI failures (building without access to MARBL), so I added the _USE_MARBL_TRACERS cpp
    • Currently, CESM interface always builds with -D_USE_MARBL_TRACERS; will put in logic to be smarter about that (and to only build MARBL itself) towards end of project

  3. Have registered diagnostics with the model, though I may need to be smarter about pointing out 2D variables?

    "ocean_model", "ECOSYS_IFRAC"  [Unused]
    ! long_name: Ice Fraction for ecosys fluxes
    ! units: fraction
    ! cell_methods: xh:mean yh:mean zl:mean area:mean
"ocean_model", "ECOSYS_IFRAC_xyave"  [Unused]
    ! long_name: Ice Fraction for ecosys fluxes
    ! units: fraction
    ! cell_methods: zl:mean
"ocean_model_z", "ECOSYS_IFRAC"  [Unused]
    ! long_name: Ice Fraction for ecosys fluxes
    ! units: fraction
    ! cell_methods: xh:mean yh:mean z_l:mean area:mean
"ocean_model_z", "ECOSYS_IFRAC_xyave"  [Unused]
    ! long_name: Ice Fraction for ecosys fluxes
    ! units: fraction
    ! cell_methods: z_l:mean
"ocean_model_rho2", "ECOSYS_IFRAC"  [Unused]
    ! long_name: Ice Fraction for ecosys fluxes
    ! units: fraction
    ! cell_methods: xh:mean yh:mean rho2_l:mean area:mean
"ocean_model_rho2", "ECOSYS_IFRAC_xyave"  [Unused]
    ! long_name: Ice Fraction for ecosys fluxes
    ! units: fraction
    ! cell_methods: rho2_l:mean
"ocean_model_d2", "ECOSYS_IFRAC"  [Unused]
    ! long_name: Ice Fraction for ecosys fluxes
    ! units: fraction
    ! cell_methods: xh:mean yh:mean zl:mean area:mean
"ocean_model_d2", "ECOSYS_IFRAC_xyave"  [Unused]
    ! long_name: Ice Fraction for ecosys fluxes
    ! units: fraction
    ! cell_methods: zl:mean
"ocean_model_z_d2", "ECOSYS_IFRAC"  [Unused]
    ! long_name: Ice Fraction for ecosys fluxes
    ! units: fraction
    ! cell_methods: xh:mean yh:mean z_l:mean area:mean
"ocean_model_z_d2", "ECOSYS_IFRAC_xyave"  [Unused]
    ! long_name: Ice Fraction for ecosys fluxes
    ! units: fraction
    ! cell_methods: z_l:mean
"ocean_model_rho2_d2", "ECOSYS_IFRAC"  [Unused]
    ! long_name: Ice Fraction for ecosys fluxes
    ! units: fraction
    ! cell_methods: xh:mean yh:mean rho2_l:mean area:mean
"ocean_model_rho2_d2", "ECOSYS_IFRAC_xyave"  [Unused]
    ! long_name: Ice Fraction for ecosys fluxes
    ! units: fraction

April 21, 2020

General discussion

  1. JRA high-res
    • Division of labor for forcing / initial condition files? I'm happy to take some pre-existing scripts for generating x1 files and modify them for use w/ 0.1 degree, but don't want to duplicate labor if others are already on it
  2. More glade corruption?
    • Mike Mills was running into file-system troubles that seemed reminiscent of errors I saw last month when testing merge of single-column MARBL branch in CESM
    • CISL has re-opened my original ticket (if that link doesn't work, perhaps this one will)

MARBL Software Updates

POP Software Updates

MOM Software Updates

  1. Tracking progress via github project board
    • To-do: break down these big tasks into many smaller tasks. E.g. instead of add MARBL output to history file, create issues for
      1. running MARBL python script to generate diag list
      2. adding MARBL list to diag_table
      3. modifying fortran in driver to accumulate desired diagnostics correctly
  2. Emailed Alistair and Andrew about bringing in MARBL as git submodule rather than using generic tracer

April 7, 2020

General discussion

  1. Hi-res + BGC
    • #25 has been merged (brings Keith's CESM 2.1 updates for JRA to POP master) and the tag is in the plans for cesm2_2_alpha04g; this will let us start out of cesm2_2_beta04 (need to add _HR compset, put together emissions dataset, etc)
    • Keith asked to discuss output for the project

MARBL Software Updates

  1. #338 has been merged.

    • Yay!

    • Three issues that were waiting for single-column test:

      1. #53: migrate k loop (mostly done, three or four more function calls)
      2. #176: loop to kmt instead of k
      3. #336: clean up stand-alone timer results

      I could see spending some time on #53 and / or #176, though my feeling is that it's a low priority. #336 is a wishlist issue item, not something that needs attention right now

POP Software Updates

  1. Not BGC related: I was asked to help create a new compset for extending SSP runs

    • Not too much work to put together PR #27
    • Testing uncovered issue #28

    I don't want to fall down this rabbit hole, though I'm probably now in a better spot to clean this up than Alper

MOM Software Updates

  1. I did generate initial conditions for tracers on the MOM grid last summer, and updated the slides from previous meeting accordingly
  2. Not much progress to report, but I really want to stop putting out small fires and attack the big fire

No Meetings

  • March 24, 2020
  • March 10, 2020 (No urgent need to meet, other projects taking precedence)
  • February 25, 2020 (CGD Town Hall)

February 11, 2020

General discussion

MARBL Software Updates

  1. #338 is the stand-alone test of the compute() functions, just needs more documentation
    • Read through user guide, make sure it is all up-to-date with examples from the stand-alone driver (I think I finished this section last fall)
    • Link to a page with details of POP's saved state implementation from the general saved state page?
    • unit testing: need more detail on what the tests are doing
    • Write up regression testing page

POP Software Updates

MOM Software Updates

  1. Putting together slides on the process
    • Still waiting on #338 before making more progress

January 28, 2020

General discussion

  1. More talk about CESM2 papers
    • Updated nutrient plots (using different region mask for zonal means)
    • xpersist: caching data in /glade/p/cgd/oce/projects/cesm2-marbl
  2. Geocat hack-a-thon

MARBL Software Updates

POP Software Updates

MOM Software Updates

January 14, 2020

General discussion

  1. Progress on tables / plots for CESM2 papers
    • cesm2-marbl
      1. Using xpersist to store time series of global averages in flux table (will also be used for time series plots)
      2. Nutrient plots need to use the pop-tools region mask
    • Keith update?
  2. Working with Precious to get him using intake-esm for LENS study
  3. Been talking to Matt about how to turn the cesm2-marbl repository into a more general analysis package (or packages)

MARBL Software Updates

POP Software Updates

MOM Software Updates

  1. Design Document
  2. How to Submit Code Modifications
  3. Dev meeting agenda / rough notes
  4. Processes Still Being Developed
Clone this wiki locally