Matt Long (Sep 03 2020 at 18:08): Matt Long (Sep 03 2020 at 18:08):In yesterday's call, we discussed extending the JRA-forced hindcasts that @Kristen Krumhardt has been running, possibly with changes in the ice albedo settings.
@David Bailey and @Marika Holland, let us know when you have these available.
I think we are otherwise poised to continue these experiments.
Stephen Yeager (Sep 03 2020 at 18:25):I'm going to have a chat (either this afternoon or tomorrow morning) with Yu-Chiao Liang (WHOI) who successfully implemented sea-ice nudging in fully-coupled CESM2. So, if we could hold off a bit...
Matt Long (Sep 03 2020 at 18:32):Sounds good. Just closing the loop on our conversation so we can move forward when ready.
David Bailey (Sep 03 2020 at 19:39):Do we have sea ice diagnostics for this run? If not, can you point me to the case so I can run them? I need to assess where the sea ice is at before suggesting tuning.
Matt Long (Sep 03 2020 at 19:46):@Kristen Krumhardt can point you to the data.
Kristen Krumhardt (Sep 03 2020 at 19:50):There are two JRA-forced hindcast runs that are possibilities. Here are the paths to the output data:
One with the standard 3 phytoplankton functional types: /glade/scratch/kristenk/archive/g.e22b05.G1850ECOIAF_JRA.TL319_g17.005
And another with 4 phytoplankton functional types: /glade/scratch/kristenk/archive/g.e22b05.G1850ECOIAF_JRA.TL319_g17.cocco.001
Kristen Krumhardt (Sep 03 2020 at 19:51):These both currently have 2 JRA cycles.
David Bailey (Sep 03 2020 at 21:02):Ok. Here are the sea ice diagnostics for the end of the second cycle. The ice is definitely too thin here. I will confer with Marika.
David Bailey (Sep 03 2020 at 21:05):Where are the case directories for these? I want to see where the sea ice tunings are set currently.
Kristen Krumhardt (Sep 03 2020 at 21:07):Here is the case directory for the 005 case:
/glade/work/kristenk/cesm_work/cesm2.2_tunings_v2_cases/g.e22b05.G1850ECOIAF_JRA.TL319_g17.005
Kristen Krumhardt (Sep 03 2020 at 21:07):and for the other case:
/glade/work/kristenk/cesm_work/cesm2.2_tunings_v2_cases/g.e22b05.G1850ECOIAF_JRA.TL319_g17.cocco.001
David Bailey (Sep 03 2020 at 21:08):Great. Thanks.
Stephen Yeager (Sep 04 2020 at 19:51):Just sent an invite for a zoom next week to discuss possible sensitivity tests branched from Kristen's runs to improve sea ice thickness. Candidate mods are: high latitude forcing tweaks, strong SST/SSS restoring under sea ice, ice albedo changes, combinations of above.
Stephen Yeager (Sep 09 2020 at 20:30):Is one of Kristen's two cases preferred for branching from for sea-ice sensitivity experiments?
Kristen Krumhardt (Sep 09 2020 at 20:44):Perhaps just use the 3 PFT one since there are less tracers and I don't think that the sea ice would differ much. So the g.e22b05.G1850ECOIAF_JRA.TL319_g17.005 case.
Stephen Yeager (Sep 10 2020 at 15:34):@David Bailey what should albedo settings be for albedo sensitivity run?
David Bailey (Sep 10 2020 at 15:42):Marika and I are discussing this in about 20 minutes.
David Bailey (Sep 10 2020 at 17:12):r_snw = 1.6
dt_mlt = 0.5
rsnw_mlt = 1000.
We are still working on forcing modifications. I am willing to do that run.
Stephen Yeager (Oct 13 2020 at 21:18):To update everyone, Who Kim is going to run 2 cycles of JRA55-do from state of rest using what we think will be our final modified configuration to improve sea ice: ice albedo changes + strong SST restoring under ice. He'll run it with the retuned kappa. Hopefully we can set up a meeting at end of the week to discuss next steps (moving forward with BGC FOSI for SMYLE).
Who Kim (Oct 13 2020 at 22:40):@Stephen Yeager, The results from this run will not be available until early next week since I cannot start this run until Thursday due to the maintenance of Cheyenne.
Marika Holland (Oct 14 2020 at 18:51):@Stephen Yeager @Who Kim I'd be happy to chat about Who's runs next week although my schedule is pretty full
Keith Lindsay (Nov 04 2020 at 18:19):I'm on target to start the 6th cycle of the JRA FOSI experiment fro SMYLE tonight. I'd like to confirm that I'm setup to generate all ocn/ice restart files needed for SMYLE. It looks like the CSL proposal calls out initializing SMYLE over years 1970-2019, though there is also mention of extensions to pre-1970 start dates. In what forcing year should I start generating monthly restart files?
Stephen Yeager (Nov 04 2020 at 19:18):Thanks @Keith Lindsay. I would save monthly restarts for the entire 6th cycle (1958-2019). We are also thinking of introducing ocean spread by starting SMYLE hindcasts from lagged ocean states (say, first 5 days of a month, rather than 0Z at 1st of the month). This is probably too complicated to accomplish in your 6th cycle on the first go, but we may want to go back to rerun 5 or 10 day segments of the 6th cycle to save these additional restarts.
Keith Lindsay (Nov 17 2020 at 04:32):Diagnostics from the FOSI experiment for generating SMYLE ICs are now online:
popdiag
popdiag for CESM2 OMIP2 (JRA)
pop adhoc_ts (005 is SMYLE FOSI run, 004 is like 005, except for 1850 aerosol dep)
icediag vs obs
icediag vs CESM2 OMIP2 (JRA)
Stephen Yeager (Nov 17 2020 at 17:59):Excellent! Thanks Keith. Shall we all have a look at these diagnostics and meet to discuss?
Stephen Yeager (Nov 17 2020 at 22:43):Where can we find the monthly POP/CICE output from SMYLE FOSI as well as CESM2 OMIP to perform some other diagnostic analysis?
Keith Lindsay (Nov 18 2020 at 00:01):SMYLE FOSI info
CASE=g.e22.GOMIPECOIAF_JRA-1p4-2018.TL319_g17.SMYLE.005
CASEROOT=/glade/work/klindsay/cesm22_cases/SMYLE/$CASE
DOUT_S_ROOT=/glade/scratch/klindsay/archive/$CASE
I don't think the monthly history files for the CESM2 OMIP2 epxeriment exist any more.
Single variable timeseries files from the CESM2 OMIP2 experiment are in the directory
/glade/collections/cdg/timeseries-cmip6/$CASE/ocn/proc/tseries/month_1
CASE=g.e21.GOMIPECOIAF_JRA.TL319_g17.CMIP6-omip2.001
Kristen Krumhardt (Nov 20 2020 at 19:42):Hi all, here are some BGC diagnostics for the SMYLE-FOSI (mean over the 6th IAF cycle):
https://nbviewer.jupyter.org/github/kristenkrumhardt/CESM2_oceanBGC_diag/blob/master/smyle_omip_diagnostics/MARBL_diagnostic_hist_3PFT.SMYLE-fosi.ipynb
...and a similar notebook for the CESM2 OMIP run for comparison:
https://nbviewer.jupyter.org/github/kristenkrumhardt/CESM2_oceanBGC_diag/blob/master/Biological_diagnostics_OMIP-CESM2.ipynb
Stephen Yeager (Dec 14 2020 at 21:53):SMYLE-FOSI data have been archived on campaign at:
/glade/campaign/cesm/development/espwg/SMYLE/SMYLE-FOSI/
This includes: spinup restart files, POP/CICE history outputs from cycle 6
Cycle 6 restarts to be used for SMYLE are on glade here:
/glade/p/cesm/espwg/CESM2-SMYLE/initial_conditions/pop_cice
Stephen Yeager (Apr 08 2021 at 13:50):@Keith Lindsay SMYLE_FOSI cycle 6 currently extends through 2018 (sim-year 0366). We'd like to extend this through 2020 (and even into 2021) using near real-time updated JRA55-do forcing. Would be good to discuss prior to the April 23 ESP meeting.
Keith Lindsay (Apr 08 2021 at 20:43):I apologize for not getting back to you on SMYLE questions.
Extending the SMYLE FOSI runs seems doable, with some nml changes. The only code change I'm seeing that might be necessary at this point is increasing the number of years for forcing timeseries extrapolation, currently 4, for cfc's. The cfc timeseries file ends at 2015.5. So we'll need to either extend the dataset or allow the code to extrapolate more than 4 years.
I'll look into availability of an extension of the cfc dataset. Unfortunately, the DOE site, CDIAC, that hosted the data was cancelled by DOE in 2017 (that was sad), and the PMEL scientist, John Bullister, that collated the cfc obs passed away in 2018 (that was really sad). But maybe it lives on elsewhere.
How recent do we have JRA forcing?
Stephen Yeager (Apr 09 2021 at 01:08):Thanks Keith. We have official JRA55-do through 2020, and we have the capacity to push JRA55-do_beta up to real time (within a week or so). So far, SMYLE is looking promising, and we'll probably want to extend our May 1st start hindcast set up through 5/1/2021 to generate an actual forecast of 2022 conditions (in time for a reveal at the CESM workshop?).
Last updated: May 16 2025 at 17:14 UTC