Colleen Petrik (Mar 28 2022 at 22:37): Colleen Petrik (Mar 28 2022 at 22:37):@Kristen Krumhardt , there is some interest from our NSF project in running FEISTY on the high res model. We could potentially just do the California Current System if that helps with reducing post-processing the high res outputs and storage space.
Kristen Krumhardt (Mar 29 2022 at 13:18):Hi @Colleen Petrik , Sure, this sounds great. I think I could try preparing FEISTY fields globally, but if it gets to be too much, I could subset the california current region. What hi res simulation were you thinking of using? The JRA-forced one (currently 34 years, 1958 to 1991)?
Matt Long (Mar 29 2022 at 15:44):yes, the JRA run
Colleen Petrik (Apr 19 2022 at 23:16):@Matt Long and @Kristen Krumhardt. Piers called me today wanting the output of FEISTY forced by the high res model for the NSF convergence accelerator CATCH website and the prototype demonstration. I think the deadline may be Apr 28th? Is there any chance of getting the outputs post-processed by the end of this week? Subsetting the California Current region might be easiest. If not, don't hesitate to let me know and I will tell Piers that it is not possible.
Ideally it would be good to have the high res FEISTY runs completed in time for creating figures for the proposal, which is due May 25th.
Matt Long (Apr 19 2022 at 23:17):This should be relatively straightforward. I think you and Kristen were already working on it.
Matt Long (Apr 19 2022 at 23:41):*thought
Kristen Krumhardt (Apr 20 2022 at 12:57):Hi @Colleen Petrik Sorry I didn't realize there was a close deadline for this. I will get going on this today and hopefully have the post-processed output by the end of the week. I'll keep you updated.
Matt Long (Apr 20 2022 at 12:58):@Kristen Krumhardt, I suggest running the process on a (generous) CCS subset first
Colleen Petrik (Apr 20 2022 at 17:33):No worries, @Kristen Krumhardt . The deadline was not clear to me before either.
Kristen Krumhardt (Apr 21 2022 at 16:00):Hi @Colleen Petrik , I've processed hi res output for FEISTY for the California current region. It's here: /glade/work/kristenk/fish-offline/g.e22.G1850ECO_JRA_HR.TL319_t13.004.FIESTY-forcing_hiresJRA_CAcurr.nc
Two things to note:
1) all means and depth integrals are over the top 100m (rather than the top 150m). This is because the zoo_loss variable was only output for this hi-res model run as a depth integral over the top 100m, so I just did the same with the other variables to be consistent.
2) there are explicit coccolithophores in this simulation so I added this as an additional biomass variable to the output file. You will need to consider this in your calculation of "mesozoo" biomass.
Please let me know if you have any questions. I'm also happy to make this type of FEISTY input file for other regions, just let me know :)
Colleen Petrik (Apr 21 2022 at 20:00):Thanks so much, @Kristen Krumhardt ! In your opinion (or in your 4P-4Z version), would coccos be a prey item for mesozooplankton?
Kristen Krumhardt (Apr 21 2022 at 21:21):Hi @Colleen Petrik , I think of coccos as in the transition between small phytoplankton and diatoms - some are quite large, some are small. They can be eaten by copepods (e.g., see https://www.nature.com/articles/s41598-018-28073-x) and also by microzooplankton (e.g., see https://www.sciencedirect.com/science/article/pii/S0079661117301593). Perhaps you can count around half or maybe 75% of cocco biomass going towards mesozooplankton partitioning?
Colleen Petrik (Apr 28 2022 at 00:31):@Kristen Krumhardt , I also need the bottom depth of each grid cell (HT?). The file you created has KMT, which might have been accidentally saved instead of depth.
Kristen Krumhardt (Apr 28 2022 at 14:02):Good morning, @Colleen Petrik , oh I didn't know you needed that. I found a function for calculating HT in one of Matt's notebooks and added it to my hi-res processing notebook. Now the netcdf should also contain HT: /glade/work/kristenk/fish-offline/g.e22.G1850ECO_JRA_HR.TL319_t13.004.FIESTY-forcing_hiresJRA_CAcurr.nc
Kristen Krumhardt (Aug 23 2022 at 20:23):Hi @Colleen Petrik , @Matt Long , I transferred 10 years (1964 to 1973) of FEISTY input variables from the hi-res run on Frontera to Cheyenne here: /glade/scratch/kristenk/hi-res-feisty-inputs/
Each monthly timeseries file is 1 year long.
There are the variables:
• mesozooC_zint_100m_2 (just ignore the 2 in the name, as it's to differentiate it from the daily variable)
• mesozoo_loss_zint_100m
• pocToFloor
• TEMP_BOTTOM
• TEMP_mean_100m
Please let me know if I'm missing anything. I'll keep putting more files in this directory as the run progresses... thanks!!
Matt Long (Aug 31 2022 at 19:27):@Kristen Krumhardt, we should perhaps move to campaign.
Also, can you point @Michael Levy to output from the 1° integration
Matt Long (Aug 31 2022 at 19:28):or maybe @Fred Castruccio has already done this?
Kristen Krumhardt (Aug 31 2022 at 19:32):Yes I think that @Fred Castruccio already put the first ~18 years if tseries files on campaign, which includes FEISTY vars: /glade/campaign/collections/cmip/CMIP6/iHESP/FOSI_BGC/HR/g.e22.TL319_t13.G1850ECOIAF_JRA_HR.4p2z.001
Kristen Krumhardt (Oct 12 2022 at 15:03):Hello, I have transferred a 1° companion to the hi-res run to campaign. This simulation has 2 zooplankton and was forced by JRAv1.5 over the period 1958 to 2021. There are both daily and monthly variables for inputs into FEISTY. The run is located here:
/glade/campaign/cesm/development/bgcwg/projects/CESM2-OMIP2-like-4p2z-run/g.e22.GOMIPECOIAF_JRA-1p4-2018.TL319_g17.4p2z.001branch
The above run was branched off a 5 IAF cycle spin up (forcing is JRAv1.3, 1958 to 2018), output is here:
/glade/campaign/cesm/development/bgcwg/projects/CESM2-OMIP2-like-4p2z-run/g.e22.GOMIPECOIAF_JRA-1p4-2018.TL319_g17.4p2z.001 The spin up does not have FEISTY input variables, as I did this one before Mike and I had those code changes incorporated.
Colleen Petrik (Jun 07 2023 at 17:48):@Kristen Krumhardt , are the paths above the correct place to find the high-res output and the 1 degree companion run?
Kristen Krumhardt (Jun 07 2023 at 17:52):Hey Colleen, actually the hi res run was moved here: /glade/campaign/cgd/oce/projects/FOSI_BGC/HR/g.e22.TL319_t13.G1850ECOIAF_JRA_HR.4p2z.001
... and I recently re-ran the 1 degree companion because to better match the hi res run (with some better sea ice parameters and to match the riverine BGC runoff that was implemented in the hi res). The new run is still in my scratch space, here: /glade/scratch/kristenk/archive/g.e22.GOMIPECOIAF_JRA-1p4-2018.TL319_g17.4p2z.002branch
Colleen Petrik (Jun 13 2023 at 20:05):Thanks, @Kristen Krumhardt . It looks like the FEISTY forcing variables are integrated 0-150 m in the HR, but 0-100 m in the LR.
Kristen Krumhardt (Jun 13 2023 at 20:29):Hi Colleen, yes I was just seeing that for some reason I didn't get the daily feisty forcing for this new run. I will have to re-do the simulation... but in the meantime, you can use this one /glade/campaign/cesm/development/bgcwg/projects/CESM2-OMIP2-like-4p2z-run/g.e22.GOMIPECOIAF_JRA-1p4-2018.TL319_g17.4p2z.001branch
Colleen Petrik (Jun 13 2023 at 21:39):I appear to not have permission for that one, @Kristen Krumhardt
Kristen Krumhardt (Jun 13 2023 at 21:48):Ok, yep permissions on campaign are rather strict. I'll look into getting you permissions, but in the meantime I put the files you need to force feisty here: /glade/scratch/kristenk/for_Colleen
Colleen Petrik (Jun 14 2023 at 00:02):Thank you! Please let me know when the new version with vint-150m is ready.
Kristen Krumhardt (Jun 26 2023 at 18:59):Hi @Colleen Petrik , the new run is done. Daily output (150m means/zint) that is ready for FEISTY can be found here:
/glade/scratch/kristenk/archive/g.e22.GOMIPECOIAF_JRA-1p4-2018.TL319_g17.4p2z.002branch/ocn/proc/tseries/day_1
files:
g.e22.GOMIPECOIAF_JRA-1p4-2018.TL319_g17.4p2z.002branch.pop.h.ecosys.nday1.mesozooC_zint_150m.19580101-20211231.nc
g.e22.GOMIPECOIAF_JRA-1p4-2018.TL319_g17.4p2z.002branch.pop.h.ecosys.nday1.mesozoo_loss_zint_150m.19580101-20211231.nc
g.e22.GOMIPECOIAF_JRA-1p4-2018.TL319_g17.4p2z.002branch.pop.h.ecosys.nday1.TEMP_mean_150m.19580101-20211231.nc
g.e22.GOMIPECOIAF_JRA-1p4-2018.TL319_g17.4p2z.002branch.pop.h.ecosys.nday1.TEMP_BOTTOM_2.19580101-20211231.nc
g.e22.GOMIPECOIAF_JRA-1p4-2018.TL319_g17.4p2z.002branch.pop.h.ecosys.nday1.pocToFloor_2.19580101-20211231.nc
Colleen Petrik (Jul 06 2023 at 20:39):Thanks, @Kristen Krumhardt . I got those files. We had also talked about saving the mesozooplankton quadratic mortality loss in the one degree run to see if we needed that from the high res also. Did you do that?
Kristen Krumhardt (Jul 06 2023 at 20:44):No unfortunately I don't think we got the code changes done that would be necessary to output mesozoo quadratic mortality ... I can speak with Mike about it , though and perhaps we could rerun the 1° then...
Michael Levy (Jul 10 2023 at 16:43):@Kristen Krumhardt I'm around this week and then gone for the next two -- happy to chat about adding the quadratic mortality term either in the next few days or in August
Kristen Krumhardt (Jul 10 2023 at 16:45):Hey Mike, thanks! I'm pretty swamped this week so let's plan on doing this at the beginning of August ...
Last updated: May 16 2025 at 17:14 UTC