diff --git a/Bottom_property_bias_WOA18.ipynb b/Bottom_property_bias_WOA18.ipynb index 5886eb8..2b7d266 100644 --- a/Bottom_property_bias_WOA18.ipynb +++ b/Bottom_property_bias_WOA18.ipynb @@ -20,7 +20,7 @@ }, { "cell_type": "code", - "execution_count": 1, + "execution_count": null, "id": "68078c46-4b7b-4b2f-b8c1-2d2538f46ebc", "metadata": {}, "outputs": [], @@ -31,7 +31,7 @@ }, { "cell_type": "code", - "execution_count": 2, + "execution_count": null, "id": "50093a3b-2193-4041-8409-de3dc594cf7d", "metadata": {}, "outputs": [], @@ -63,7 +63,7 @@ }, { "cell_type": "code", - "execution_count": 3, + "execution_count": null, "id": "3ed173d4-cc6d-4693-934a-ff02dd96e632", "metadata": {}, "outputs": [], @@ -86,7 +86,7 @@ }, { "cell_type": "code", - "execution_count": 4, + "execution_count": null, "id": "dfcdc828-ad45-4ef0-92bc-05df9c036270", "metadata": {}, "outputs": [], @@ -106,7 +106,7 @@ }, { "cell_type": "code", - "execution_count": 5, + "execution_count": null, "id": "806241bd-ebea-4640-ab69-cea6edba62de", "metadata": {}, "outputs": [], @@ -126,7 +126,7 @@ }, { "cell_type": "code", - "execution_count": 6, + "execution_count": null, "id": "04229c22-3781-400f-a390-8c5eaaef5189", "metadata": {}, "outputs": [], @@ -153,35 +153,10 @@ }, { "cell_type": "code", - "execution_count": 7, + "execution_count": null, "id": "382835f3-a3d7-419f-a92f-e65808230c59", "metadata": {}, - "outputs": [ - { - "name": "stdout", - "output_type": "stream", - "text": [ - "/archive/Olga.Sergienko/fre_om5/FMS2023.01_om5_20231122/om5_b03_standard_GEBCO_BMv3Ant/gfdl.ncrc5-intel23-prod/pp/ocean_monthly/ts/monthly/5yr/ocean_monthly.197301-197712.tob.nc\n", - "/archive/Olga.Sergienko/fre_om5/FMS2023.01_om5_20231122/om5_b03_standard_GEBCO_BMv3Ant/gfdl.ncrc5-intel23-prod/pp/ocean_monthly/ts/monthly/5yr/ocean_monthly.197801-198212.tob.nc\n", - "/archive/Olga.Sergienko/fre_om5/FMS2023.01_om5_20231122/om5_b03_standard_GEBCO_BMv3Ant/gfdl.ncrc5-intel23-prod/pp/ocean_monthly/ts/monthly/5yr/ocean_monthly.198301-198712.tob.nc\n", - "/archive/Olga.Sergienko/fre_om5/FMS2023.01_om5_20231122/om5_b03_standard_GEBCO_BMv3Ant/gfdl.ncrc5-intel23-prod/pp/ocean_monthly/ts/monthly/5yr/ocean_monthly.198801-199212.tob.nc\n", - "/archive/Olga.Sergienko/fre_om5/FMS2023.01_om5_20231122/om5_b03_standard_GEBCO_BMv3Ant/gfdl.ncrc5-intel23-prod/pp/ocean_monthly/ts/monthly/5yr/ocean_monthly.199301-199712.tob.nc\n", - "/archive/Olga.Sergienko/fre_om5/FMS2023.01_om5_20231122/om5_b03_standard_GEBCO_BMv3Ant/gfdl.ncrc5-intel23-prod/pp/ocean_monthly/ts/monthly/5yr/ocean_monthly.199801-200212.tob.nc\n", - "/archive/Olga.Sergienko/fre_om5/FMS2023.01_om5_20231122/om5_b03_standard_GEBCO_BMv3Ant/gfdl.ncrc5-intel23-prod/pp/ocean_monthly/ts/monthly/5yr/ocean_monthly.200301-200712.tob.nc\n", - "/archive/Olga.Sergienko/fre_om5/FMS2023.01_om5_20231122/om5_b03_standard_GEBCO_BMv3Ant/gfdl.ncrc5-intel23-prod/pp/ocean_monthly/ts/monthly/5yr/ocean_monthly.200801-201212.tob.nc\n", - "/archive/Olga.Sergienko/fre_om5/FMS2023.01_om5_20231122/om5_b03_standard_GEBCO_BMv3Ant/gfdl.ncrc5-intel23-prod/pp/ocean_monthly/ts/monthly/5yr/ocean_monthly.201301-201712.tob.nc\n", - "/archive/Olga.Sergienko/fre_om5/FMS2023.01_om5_20231122/om5_b03_standard_GEBCO_BMv3Ant/gfdl.ncrc5-intel23-prod/pp/ocean_monthly/ts/monthly/5yr/ocean_monthly.197301-197712.sob.nc\n", - "/archive/Olga.Sergienko/fre_om5/FMS2023.01_om5_20231122/om5_b03_standard_GEBCO_BMv3Ant/gfdl.ncrc5-intel23-prod/pp/ocean_monthly/ts/monthly/5yr/ocean_monthly.197801-198212.sob.nc\n", - "/archive/Olga.Sergienko/fre_om5/FMS2023.01_om5_20231122/om5_b03_standard_GEBCO_BMv3Ant/gfdl.ncrc5-intel23-prod/pp/ocean_monthly/ts/monthly/5yr/ocean_monthly.198301-198712.sob.nc\n", - "/archive/Olga.Sergienko/fre_om5/FMS2023.01_om5_20231122/om5_b03_standard_GEBCO_BMv3Ant/gfdl.ncrc5-intel23-prod/pp/ocean_monthly/ts/monthly/5yr/ocean_monthly.198801-199212.sob.nc\n", - "/archive/Olga.Sergienko/fre_om5/FMS2023.01_om5_20231122/om5_b03_standard_GEBCO_BMv3Ant/gfdl.ncrc5-intel23-prod/pp/ocean_monthly/ts/monthly/5yr/ocean_monthly.199301-199712.sob.nc\n", - "/archive/Olga.Sergienko/fre_om5/FMS2023.01_om5_20231122/om5_b03_standard_GEBCO_BMv3Ant/gfdl.ncrc5-intel23-prod/pp/ocean_monthly/ts/monthly/5yr/ocean_monthly.199801-200212.sob.nc\n", - "/archive/Olga.Sergienko/fre_om5/FMS2023.01_om5_20231122/om5_b03_standard_GEBCO_BMv3Ant/gfdl.ncrc5-intel23-prod/pp/ocean_monthly/ts/monthly/5yr/ocean_monthly.200301-200712.sob.nc\n", - "/archive/Olga.Sergienko/fre_om5/FMS2023.01_om5_20231122/om5_b03_standard_GEBCO_BMv3Ant/gfdl.ncrc5-intel23-prod/pp/ocean_monthly/ts/monthly/5yr/ocean_monthly.200801-201212.sob.nc\n", - "/archive/Olga.Sergienko/fre_om5/FMS2023.01_om5_20231122/om5_b03_standard_GEBCO_BMv3Ant/gfdl.ncrc5-intel23-prod/pp/ocean_monthly/ts/monthly/5yr/ocean_monthly.201301-201712.sob.nc\n" - ] - } - ], + "outputs": [], "source": [ "# Determine what files are needed (leave this up to the developer for flexibility)\n", "# This is an example of what someone might do:\n", @@ -232,18 +207,10 @@ }, { "cell_type": "code", - "execution_count": 8, + "execution_count": null, "id": "7096dc5b-559b-4221-a189-0836d77afd12", "metadata": {}, - "outputs": [ - { - "name": "stdout", - "output_type": "stream", - "text": [ - "Calling dmget on files ...\n" - ] - } - ], + "outputs": [], "source": [ "if not \"DORA_EXECUTE\" in os.environ.keys():\n", " print(\"Calling dmget on files ...\")\n", @@ -261,7 +228,7 @@ }, { "cell_type": "code", - "execution_count": 9, + "execution_count": null, "id": "70b1b444-20e4-4e63-b8a6-4d076d977eb2", "metadata": {}, "outputs": [], @@ -289,19 +256,20 @@ }, { "cell_type": "code", - "execution_count": 10, + "execution_count": null, "id": "3f322d2e-1c88-479f-86c7-5910d71cf3f9", "metadata": {}, "outputs": [], "source": [ "if \"om4\" in model_type:\n", - " dsobs = \"/archive/jpk/datasets/OM5/obs/WOA18/WOA18_z_35level_OM4_1080x1440_annual_v20240224.nc\"\n", + " dsobs = \"/archive/jpk/datasets/OM5/obs/WOA18/WOA18_z_35level_OM4_1080x1440_annual_v20240602.nc\"\n", "elif \"om5\" in model_type:\n", - " dsobs = \"/archive/jpk/datasets/OM5/obs/WOA18/WOA18_z_35level_OM5_1161x1440_annual_v20240224.nc\"\n", + " dsobs = \"/archive/jpk/datasets/OM5/obs/WOA18/WOA18_z_35level_OM5_1161x1440_annual_v20240602.nc\"\n", "else:\n", " raise ValueError(f\"Model type {model_type} not supported\")\n", "\n", "dsobs = xr.open_dataset(dsobs).isel(time=0)\n", + "dsobs = dsobs.assign_coords({\"xh\": ds.xh, \"yh\": ds.yh})\n", "\n", "# Get select bottom most numerical value\n", "dsobs['tob'] = dsobs['thetao'].isel(z_l=np.isnan(dsobs['thetao']).argmax(\"z_l\")-1)\n", @@ -318,7 +286,7 @@ }, { "cell_type": "code", - "execution_count": 11, + "execution_count": null, "id": "c4fa45cd-e682-4aac-a98c-a0af350de2a3", "metadata": {}, "outputs": [], @@ -338,7 +306,7 @@ }, { "cell_type": "code", - "execution_count": 12, + "execution_count": null, "id": "294f6a30-47b8-4938-8fa4-f6390e841234", "metadata": {}, "outputs": [], @@ -354,7 +322,7 @@ }, { "cell_type": "code", - "execution_count": 13, + "execution_count": null, "id": "574269bd-afb5-437f-ad30-703b70489d85", "metadata": {}, "outputs": [], @@ -377,7 +345,7 @@ }, { "cell_type": "code", - "execution_count": 14, + "execution_count": null, "id": "80be5109-0e21-450c-827f-2bed52263b00", "metadata": {}, "outputs": [], @@ -397,7 +365,7 @@ }, { "cell_type": "code", - "execution_count": 15, + "execution_count": null, "id": "aadf7139-8861-4a4a-8a2a-1a6c6244cc37", "metadata": {}, "outputs": [], @@ -409,856 +377,6 @@ " )" ] }, - { - "cell_type": "code", - "execution_count": 16, - "id": "64a9c87a-dc17-4c3c-ba71-db76f86dad3e", - "metadata": {}, - "outputs": [ - { - "data": { - "text/html": [ - "
<xarray.Dataset> Size: 2GB\n", - "Dimensions: (xh: 1440, yh: 1161, time: 41, xq: 1441, yq: 1162)\n", - "Coordinates:\n", - " * xh (xh) float64 12kB 0.5 1.5 2.5 ... 1.438e+03 1.438e+03 1.44e+03\n", - " * yh (yh) float64 9kB 0.5 1.5 2.5 3.5 ... 1.158e+03 1.16e+03 1.16e+03\n", - " geolon (yh, xh) float32 7MB -299.9 -299.6 -299.4 ... 60.0 60.0 60.0\n", - " geolat (yh, xh) float32 7MB -88.52 -88.52 -88.52 ... 64.33 64.22 64.11\n", - " areacello (yh, xh) float32 7MB 7.708e+06 7.708e+06 ... 2.384e+06 1.201e+06\n", - " * xq (xq) float64 12kB 1.0 2.0 3.0 ... 1.439e+03 1.44e+03 1.441e+03\n", - " * yq (yq) float64 9kB 1.0 2.0 3.0 4.0 ... 1.16e+03 1.161e+03 1.162e+03\n", - " * time (time) object 328B 1975-07-02 12:00:00 ... 2015-07-02 12:00:00\n", - "Data variables:\n", - " tob (time, yh, xh) float64 548MB dask.array<chunksize=(1, 1161, 1440), meta=np.ndarray>\n", - " geolat_c (time, yq, xq) float64 549MB -88.57 -88.57 -88.57 ... 64.17 64.06\n", - " sob (time, yh, xh) float64 548MB dask.array<chunksize=(1, 1161, 1440), meta=np.ndarray>\n", - " geolon_c (time, yq, xq) float64 549MB -300.0 -299.8 -299.5 ... 60.0 60.0\n", - "Attributes:\n", - " title: om5_b03_standard_GEBCO_BMv3Ant\n", - " associated_files: areacello: 19730101.ocean_static.nc\n", - " grid_type: regular\n", - " grid_tile: N/A\n", - " external_variables: areacello