NCAR / amwg_dev

Repo to store model sandboxes and cases used for CAM development
9 stars 2 forks source link

b.cesm3_cam058_mom_e.B1850WscMOM.ne30_L58_t061.camdev_cice5.027a #182

Open gustavo-marques opened 2 years ago

gustavo-marques commented 2 years ago

Description:

Same as 026a with a hybrid vertical coordinate system in the ocean model

In user_nl_mom we set:

INPUTDIR = "INPUT"
NK = 75
REGRIDDING_COORDINATE_MODE = "HYCOM1"
INTERPOLATION_SCHEME = "P1M_H2"
REMAPPING_SCHEME = "PPM_H4"
BOUNDARY_EXTRAPOLATION = True
ALE_COORDINATE_CONFIG = "HYBRID:hybrid_75layer_zstar2.50m-2020-11-23.nc,sigma2,dz"
REGRID_COMPRESSIBILITY_FRACTION = 0.01
MAXIMUM_INT_DEPTH_CONFIG = "FNC1:5,8000.0,1.0,.01"
MAX_LAYER_THICKNESS_CONFIG = "FNC1:400,31000.0,0.1,.01"
COORD_CONFIG = "ALE"

Directory INPUT must be copied to the $RUNDIR with appropriate symbolic links in it.

Case directory: Locally (if still available): /glade/p/cesmdata/cseg/runs/cesm2_0/b.cesm3_cam058_mom_e.B1850WscMOM.ne30_L58_t061.camdev_cice5.027a

On github:

Sandbox: Locally (if still available): /glade/work/hannay/cesm_tags/cesm3_cam6_3_058_MOM_e

On github: hash: 8f70c08

Diagnostics: AMWG diags (if available) https://webext.cgd.ucar.edu/BWsc1850MOM/$CASE/atm/

Contacts: @gustavo-marques, @cecilehannay, @JulioTBacmeister, @fobryan3, @adamrher

justin-richling commented 2 years ago

Quick ADF vs 026a for years 17-27 can be found here.

justin-richling commented 2 years ago

Timeseries:

TimeSeries_ANN

swrneale commented 2 years ago

Yikes, this is depressing. We're just delaying icing almost no matter what. I wonder if we run out 23, and even though it's warm, use that as an initial condition for subsequent runs.

gustavo-marques commented 2 years ago

Interesting how RESTOM oscillates in 27a. No other run shows such behavior.

justin-richling commented 2 years ago

27a is being plotted with only the yearly average while the others are 5-yr rolling, would that explain the appearance of the oscillating behavior?

gustavo-marques commented 2 years ago

Yes, that explains! Thanks, @justin-richling.

JulioTBacmeister commented 2 years ago

We decided to extend 25a which has somewhat more acceptable RESTOM than

  1. Let's see where that goes. I think we need to reassess our strategy this Friday. We also need to look more closely at the 100 years of 26a and 26b despite the Lab Sea freeze. What does ENSO look like ... etc. ?

On Tue, Sep 27, 2022 at 9:37 AM Rich Neale @.***> wrote:

Yikes, this is depressing. We're just delaying icing almost no matter what. I wonder if we run out 23, and even though it's warm, use that as an initial condition for subsequent runs.

— Reply to this email directly, view it on GitHub https://github.com/NCAR/amwg_dev/issues/182#issuecomment-1259685968, or unsubscribe https://github.com/notifications/unsubscribe-auth/ACGLMTVEKQHP7EPNEJUCNNDWAMICFANCNFSM6AAAAAAQUHLUZE . You are receiving this because you were mentioned.Message ID: @.***>

gustavo-marques commented 2 years ago

Sounds good. I am going to stop 27a.

JulioTBacmeister commented 2 years ago

Hmmm is convection totally dead everywhere???

adamrher commented 2 years ago

I think we need to reassess our strategy this Friday. We also need to look more closely at the 100 years of 26a and 26b despite the Lab Sea freeze. What does ENSO look like ... etc. ?

Has anyone run any of the variability packages before? I know Adam Phillip's newly pythonized variability packages requires them to be in cmip / timeseries format first. What about MDF?

gustavo-marques commented 2 years ago

@JulioTBacmeister: the convection is the Lab Sea died after ~ 7 years and did not recover (it's at year 30 now). I am happy to continue the run if you think it's worth it.

cecilehannay commented 2 years ago

@adamrher: you can run directly the CVDP from the ADF. I will ask @justin-richling to give it a try.

JulioTBacmeister commented 2 years ago

@gustavo-marques the feeling in the Ocean section has always been that once dead there is no way for the N Atlantic convection to recover. Is this worth revisiting with MOM? I leave it to you experts.
Another question: Is convection dead everywhere in 27a, i.e., Irminger Sea ... further east ...?

JulioTBacmeister commented 2 years ago

To start things off ... Time series of SST in Nino34 for 26a and 26b. image image

gustavo-marques commented 2 years ago

@JulioTBacmeister: no, convection still occurs elsewhere. The plot below shows winter mean MLDs averaged between years 17-27. b cesm3_cam058_mom_e B1850WscMOM ne30_L58_t061 camdev_cice5 027a_MLD_winter

JulioTBacmeister commented 2 years ago

So, is 27a still "viable"? If AMOC is being forced in different regions, but is still there is that OK? Perhaps with the biases we see in 1 degree ocean models we should accept that convection may happen in different locations than in reality. It's interesting that the mixed layer depth in 27a is bigger than estimated in several places.

On Tue, Sep 27, 2022 at 2:03 PM Gustavo Marques @.***> wrote:

@JulioTBacmeister https://github.com/JulioTBacmeister: no, convection still occurs elsewhere. The plot below shows winter mean MLDs averaged between years 17-27. [image: b cesm3_cam058_mom_e B1850WscMOM ne30_L58_t061 camdev_cice5 027a_MLD_winter] https://user-images.githubusercontent.com/11339137/192624272-80bb1cfd-01e1-493d-bf46-1d18520f9d15.png

— Reply to this email directly, view it on GitHub https://github.com/NCAR/amwg_dev/issues/182#issuecomment-1259987990, or unsubscribe https://github.com/notifications/unsubscribe-auth/ACGLMTSPY2NBAZD5JYMTVR3WANHJVANCNFSM6AAAAAAQUHLUZE . You are receiving this because you were mentioned.Message ID: @.***>

swrneale commented 2 years ago

Not to repeat exactly what you showed Julio, but ENSO looks pretty darn good, especially early on in the simulations.

image image

JulioTBacmeister commented 2 years ago

Thanks Rich, This is interesting. Is it worth looking at ENSO in some of our earlier runs even though they are shorter?

swrneale commented 2 years ago

Sure. I think if we have 30 years, we can at least have a sense of the period and maybe autocorrelation. With magnitude you're never quite sure. Which runs? @cecilehannay @JulioTBacmeister is it worth starting a 'discussion' of non standard diagnostics for some key runs (MJO, blocking, monsoon etc.). I have a few scripts for these, same as CESM2 development.

JulioTBacmeister commented 2 years ago

Why don't we look at 23, 25a and 25b. And I like the idea of setting up a new discussion on github. We can call it something like "Atmospheric means and variability in coupled runs" ...? Note also: @Justin Richling @.***> is trying to get the CVDP analyses up for the 26's.

Anyone know how to move our recent comments and figures from the 27a issues (a bad place) to this new discussion?

On Wed, Sep 28, 2022 at 10:15 AM Rich Neale @.***> wrote:

Sure. I think if we have 30 years, we can at least have a sense of the period and maybe autocorrelation. With magnitude you're never quite sure. Is it worth starting a 'discussion' of non-standard. Which runs? @cecilehannay https://github.com/cecilehannay @JulioTBacmeister https://github.com/JulioTBacmeister is it worth starting a 'discussion' of non standard diagnostics for some key runs (MJO, blocking, monsoon etc.). I have a few scripts for these, same as CESM2 development.

— Reply to this email directly, view it on GitHub https://github.com/NCAR/amwg_dev/issues/182#issuecomment-1261143119, or unsubscribe https://github.com/notifications/unsubscribe-auth/ACGLMTV2MLKJGCGHYVRWANLWARVKVANCNFSM6AAAAAAQUHLUZE . You are receiving this because you were mentioned.Message ID: @.***>

gustavo-marques commented 2 years ago

@JulioTBacmeister: sorry for the delay in responding to your question. I think 27a is not "viable". AMOC time series at both 26 N and 45 N are declining very rapidly, see figures below (please ignore the years). b cesm3_cam058_mom_e B1850WscMOM ne30_L58_t061 camdev_cice5 027a_MOC_26N_time_series b cesm3_cam058_mom_e B1850WscMOM ne30_L58_t061 camdev_cice5 027a_MOC_45N_time_series