Adding Mixed Layer Depth Diagnostic#111
Open
breichl wants to merge 1 commit into
Open
Conversation
- Mixed Layer depth diagnostic code added for analyzing model output - Takes MLD_003, MLD_EN1, MLD_EN2, OR MLD_EN3 as inputs. - Creates map of input (either min or max of monthly means) - Compares to similar field from obs (usually Argo based).
raphaeldussin
requested changes
Aug 8, 2023
| metadata: | ||
| origin_url: '' | ||
|
|
||
| Argo_MLD_003: |
Owner
There was a problem hiding this comment.
let's make a note to copy these files into a central location. You can also include a notebook with the preprocessing steps in the om4labs/data directory.
| import xarray as xr | ||
| import xesmf as xe | ||
| import warnings | ||
| from scipy.interpolate import griddata |
| cat = open_intake_catalog(dictArgs["platform"], "obs") | ||
| ds_obs = cat["Argo_MLD_003"].to_dask() | ||
|
|
||
| ds_model = ds_input[mldvar].groupby("time.month").mean("time") |
Owner
There was a problem hiding this comment.
consider moving groupby into calculate
| LatMin = dims[dictArgs["grid"]][2] | ||
| LatMax = dims[dictArgs["grid"]][3] | ||
|
|
||
| # Extract copy of obs domain |
Owner
There was a problem hiding this comment.
is this block's purpose to select a subdomain?
|
|
||
| regridder_mod = xe.Regridder(model, ds_plot, "bilinear", periodic=True) | ||
|
|
||
| ds_plot["model"] = (("lat", "lon"), regridder_mod(model).MLD.values) |
Owner
There was a problem hiding this comment.
there might be a more concise way to do this
| ), | ||
| ) | ||
|
|
||
| # Want to also compute the metrics here (bias, RMS, r2) |
Owner
There was a problem hiding this comment.
bilinear interpolation does not conserve statistics beyond mean
| ) | ||
|
|
||
| # Note that we are area weighting our global metrics by assuming a spherical Earth. This is a decent approximation and better than not area weighting the metrics. | ||
| ds_plot["bias"] = np.nansum( |
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
2 participants
Add this suggestion to a batch that can be applied as a single commit.This suggestion is invalid because no changes were made to the code.Suggestions cannot be applied while the pull request is closed.Suggestions cannot be applied while viewing a subset of changes.Only one suggestion per line can be applied in a batch.Add this suggestion to a batch that can be applied as a single commit.Applying suggestions on deleted lines is not supported.You must change the existing code in this line in order to create a valid suggestion.Outdated suggestions cannot be applied.This suggestion has been applied or marked resolved.Suggestions cannot be applied from pending reviews.Suggestions cannot be applied on multi-line comments.Suggestions cannot be applied while the pull request is queued to merge.Suggestion cannot be applied right now. Please check back later.
This PR is a replacement for #87
In addition to updating w/ the latest om4labs, it includes two updates requested in that PR:
Potential Remaining Issue:
Command to test:
om4labs mld --grid global --method max --mldvar MLD_EN1 -s /archive/bgr/FMS2019.01.03_mom6_20201020/OM4p25_JRA55do1.5_ePBLRL19_BBL01_VS_JHLtol/gfdl.ncrc4-intel18-prod/pp/ocean_monthly/ocean_monthly.static.nc /archive/bgr/FMS2019.01.03_mom6_20201020/OM4p25_JRA55do1.5_ePBLRL19_BBL01_VS_JHLtol/gfdl.ncrc4-intel18-prod/pp/ocean_monthly/ts/monthly/5yr/ocean_monthly.199801-200212.MLD_EN1.nc /archive/bgr/FMS2019.01.03_mom6_20201020/OM4p25_JRA55do1.5_ePBLRL19_BBL01_VS_JHLtol/gfdl.ncrc4-intel18-prod/pp/ocean_monthly/ts/monthly/5yr/ocean_monthly.200301-200712.MLD_EN1.nc /archive/bgr/FMS2019.01.03_mom6_20201020/OM4p25_JRA55do1.5_ePBLRL19_BBL01_VS_JHLtol/gfdl.ncrc4-intel18-prod/pp/ocean_monthly/ts/monthly/5yr/ocean_monthly.200801-201212.MLD_EN1.nc