search

NCAR / cesm-lens-aws

Examples of analysis of CESM LENS data publicly available on Amazon S3 (us-west-2 region) using xarray and dask
https://doi.org/10.26024/wt24-5j82
BSD 3-Clause "New" or "Revised" License
43 stars 23 forks source link

PD does not agree completely with sigma_0 derived from SALT and THETA #60

Open cspencerjones opened 3 years ago

cspencerjones commented 3 years ago

I've been trying to compare PD ("Potential Density Ref to Surface") with the potential density referenced to p=0 calculated using https://github.com/xgcm/fastjmd95 from SALT and THETA. These match relatively well, but not perfectly (and not enough for me to feel confident just using them). The size of the difference is correlated with depth, suggesting that maybe I am using a different equation of state.

difference 2plots_1100m zoomed_in (color is depth in the final plot and each point is a grid cell)

As far as I know, CESM uses the JMD95 equation of state. Does anyone have suggestions about what might be going wrong here?

dcherian commented 3 years ago

i think pop2 uses MDJWF 2003 though (https://github.com/NCAR/pop-tools/blob/master/pop_tools/eos.py)

cspencerjones commented 3 years ago

Thanks @dcherian , this sounded very promising, but sadly in the end it didn't help! The error has the same size and structure even if I use the pop-tools eos.

dcherian commented 3 years ago

cc @mnlevy1981