{ "cells": [ { "cell_type": "markdown", "id": "8d022ccc-4918-4194-8fdb-0cbd1912c716", "metadata": {}, "source": [ "## Compare equatorial ocean metrics between obs and MOM6 run\n", "1) Thermocline depth, gradient, strength \n", "2) Compare cold tongue strength \n", "3) EUC MOM6 vs obs, strength, depth, gradient\n", "\n", "\n", "Obs data: # Use something more current WOA18\n", "\n", "- /glade/p/cesm/omwg/obs_data/phc/PHC2_TEMP_tx0.66v1_34lev_ann_avg.nc \n", "- /glade/p/cesm/omwg/obs_data/phc/PHC2_SALT_tx0.66v1_34lev_ann_avg.nc \n", "- /glade/p/cesm/omwg/obs_data/johnson_pmel/meanfit_m.nc\n", "\n", "#### Problems with this notebook: \n", "1) Actually the main problem is the use of zl and not e\n", "1) the directory is hardcoded, should do something like Gustavo did with yml file\n", "2) The units for U say they are different, but they are not different by 100. What is the issue with the EUC? \n", "3) the unified plot is just a draft, should think about that some more, how to summarize the equatorial region? " ] }, { "cell_type": "code", "execution_count": 3, "id": "8d3f24df-58af-4fc8-bb52-ff6d1a6efc23", "metadata": {}, "outputs": [], "source": [ "# Frank's notebooks in MOM6-modeloutputanalysis/EquatorialPacific\n", "# use the h files \n", "# plot the model coordinates e" ] }, { "cell_type": "code", "execution_count": 1, "id": "cbc947a4-b1d6-4ef0-8b64-6cd664e190e2", "metadata": {}, "outputs": [], "source": [ "# Load required modules\n", "import warnings\n", "warnings.filterwarnings(\"ignore\") # I don't want any warnings (: \n", "\n", "# the usual suspects\n", "import numpy as np\n", "from datetime import date\n", "from matplotlib import pyplot as plt\n", "import cartopy.crs as ccrs\n", "import xarray as xr\n", "import glob\n", "import nc_time_axis # it says I need this to plot.. not sure \n", "\n", "# dask helpers \n", "from distributed import Client\n", "from ncar_jobqueue import NCARCluster\n", "\n", "# get mom6-tools \n", "import mom6_tools" ] }, { "cell_type": "code", "execution_count": 2, "id": "ae48bd1b-fc77-4c52-9c62-70c9d4b50250", "metadata": {}, "outputs": [ { "data": { "text/html": [ "
Client-5658fd75-23e7-11ed-ae24-3cecef1b12e0
\n", "| Connection method: Cluster object | \n", "Cluster type: dask_jobqueue.PBSCluster | \n", " \n", "
| \n", " Dashboard: https://jupyterhub.hpc.ucar.edu/stable/user/deppenme/proxy/8787/status\n", " | \n", "\n", " |
dcdf9042
\n", "| \n", " Dashboard: https://jupyterhub.hpc.ucar.edu/stable/user/deppenme/proxy/8787/status\n", " | \n", "\n", " Workers: 0\n", " | \n", "
| \n", " Total threads: 0\n", " | \n", "\n", " Total memory: 0 B\n", " | \n", "
Scheduler-d693d26b-f1b4-4dca-affa-3135ccca1685
\n", "| \n", " Comm: tcp://10.12.206.60:33813\n", " | \n", "\n", " Workers: 0\n", " | \n", "
| \n", " Dashboard: https://jupyterhub.hpc.ucar.edu/stable/user/deppenme/proxy/8787/status\n", " | \n", "\n", " Total threads: 0\n", " | \n", "
| \n", " Started: Just now\n", " | \n", "\n", " Total memory: 0 B\n", " | \n", "