{ "cells": [ { "cell_type": "code", "execution_count": 1, "id": "rolled-registration", "metadata": {}, "outputs": [], "source": [ "from netCDF4 import Dataset, num2date, date2num\n", "import numpy as np\n", "from scipy import stats\n", "import matplotlib.pyplot as plt\n", "import math\n", "import xarray as xr\n", "from pandas import Index" ] }, { "cell_type": "markdown", "id": "therapeutic-grammar", "metadata": {}, "source": [ "## Read FOCI-MOZ SAM index" ] }, { "cell_type": "code", "execution_count": 2, "id": "visible-ranch", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "/gxfs_work1/fs1/work-geomar/smomw418/models/foci/experiments/SAM_index\n" ] } ], "source": [ "%cd /sfs/fs1/work-geomar/smomw418/models/foci/experiments/SAM_index" ] }, { "cell_type": "code", "execution_count": 3, "id": "mental-money", "metadata": {}, "outputs": [], "source": [ "sw054 = xr.open_dataset('FOCI1.4-SW054_SAM_index_1958-2013.nc', decode_times=True)\n", "sw064 = xr.open_dataset('FOCI1.4-SW064_SAM_index_1958-2013.nc', decode_times=True)\n", "sw065 = xr.open_dataset('FOCI1.4-SW065_SAM_index_1958-2013.nc', decode_times=True)" ] }, { "cell_type": "markdown", "id": "swiss-island", "metadata": {}, "source": [ "## Read FOCI_HIST SAM index" ] }, { "cell_type": "code", "execution_count": 4, "id": "reported-english", "metadata": {}, "outputs": [], "source": [ "sw060 = xr.open_dataset('FOCI1.7-SW060_SAM_index_1958-2013.nc', decode_times=True)\n", "sw061 = xr.open_dataset('FOCI1.7-SW061_SAM_index_1958-2013.nc', decode_times=True)\n", "jh019 = xr.open_dataset('FOCI1.9.4-JH019_SAM_index_1958-2013.nc', decode_times=True)" ] }, { "cell_type": "code", "execution_count": 5, "id": "afraid-glucose", "metadata": {}, "outputs": [ { "data": { "text/html": [ "
<xarray.Dataset>\n", "Dimensions: (lat: 96, lat_SH: 48, plev: 26, time: 20454)\n", "Coordinates:\n", " * plev (plev) float64 1e+03 925.0 850.0 775.0 700.0 ... 1.0 0.5 0.2 0.1\n", " * lat (lat) float32 88.57 86.72 84.86 83.0 ... -84.86 -86.72 -88.57\n", " * lat_SH (lat_SH) float32 -0.9326 -2.798 -4.663 ... -84.86 -86.72 -88.57\n", " * time (time) datetime64[ns] 1958-01-01T23:52:00 ... 2013-12-31T23:5...\n", "Data variables:\n", " PC (time, plev) float64 ...\n", " PC_trend (time, plev) float64 ...\n", " EOF (lat_SH, plev) float64 2.097 2.017 1.992 ... -251.2 -219.3\n", " EOF_global (plev, lat) float64 0.6788 0.6984 0.7663 ... -218.7 -219.3\n", " expvar (plev) float64 69.99 70.54 70.76 70.27 ... 69.05 67.9 64.37\n", "Attributes:\n", " Description: SAM index and pattern from FOCI-MOZ experiment SW054 for th...
array([1.00e+03, 9.25e+02, 8.50e+02, 7.75e+02, 7.00e+02, 6.00e+02, 5.00e+02,\n", " 4.00e+02, 3.00e+02, 2.50e+02, 2.00e+02, 1.50e+02, 1.00e+02, 7.00e+01,\n", " 5.00e+01, 3.00e+01, 2.00e+01, 1.00e+01, 7.00e+00, 5.00e+00, 3.00e+00,\n", " 2.00e+00, 1.00e+00, 5.00e-01, 2.00e-01, 1.00e-01])
array([ 88.57217 , 86.722534, 84.86197 , 82.99894 , 81.13498 , 79.27056 ,\n", " 77.40589 , 75.54106 , 73.67613 , 71.811134, 69.94608 , 68.08099 ,\n", " 66.21587 , 64.35073 , 62.48557 , 60.620396, 58.75521 , 56.89001 ,\n", " 55.024807, 53.159595, 51.294376, 49.429153, 47.563927, 45.698692,\n", " 43.833458, 41.96822 , 40.102978, 38.237736, 36.37249 , 34.507244,\n", " 32.641994, 30.776745, 28.911493, 27.04624 , 25.180986, 23.315731,\n", " 21.450476, 19.585218, 17.719961, 15.854704, 13.989446, 12.124187,\n", " 10.258928, 8.393669, 6.528409, 4.66315 , 2.79789 , 0.93263 ,\n", " -0.93263 , -2.79789 , -4.66315 , -6.528409, -8.393669, -10.258928,\n", " -12.124187, -13.989446, -15.854704, -17.719961, -19.585218, -21.450476,\n", " -23.315731, -25.180986, -27.04624 , -28.911493, -30.776745, -32.641994,\n", " -34.507244, -36.37249 , -38.237736, -40.102978, -41.96822 , -43.833458,\n", " -45.698692, -47.563927, -49.429153, -51.294376, -53.159595, -55.024807,\n", " -56.89001 , -58.75521 , -60.620396, -62.48557 , -64.35073 , -66.21587 ,\n", " -68.08099 , -69.94608 , -71.811134, -73.67613 , -75.54106 , -77.40589 ,\n", " -79.27056 , -81.13498 , -82.99894 , -84.86197 , -86.722534, -88.57217 ],\n", " dtype=float32)
array([ -0.93263 , -2.79789 , -4.66315 , -6.528409, -8.393669, -10.258928,\n", " -12.124187, -13.989446, -15.854704, -17.719961, -19.585218, -21.450476,\n", " -23.315731, -25.180986, -27.04624 , -28.911493, -30.776745, -32.641994,\n", " -34.507244, -36.37249 , -38.237736, -40.102978, -41.96822 , -43.833458,\n", " -45.698692, -47.563927, -49.429153, -51.294376, -53.159595, -55.024807,\n", " -56.89001 , -58.75521 , -60.620396, -62.48557 , -64.35073 , -66.21587 ,\n", " -68.08099 , -69.94608 , -71.811134, -73.67613 , -75.54106 , -77.40589 ,\n", " -79.27056 , -81.13498 , -82.99894 , -84.86197 , -86.722534, -88.57217 ],\n", " dtype=float32)
array(['1958-01-01T23:52:00.000000000', '1958-01-02T23:52:00.000000000',\n", " '1958-01-03T23:51:59.999999999', ..., '2013-12-29T23:52:00.000000000',\n", " '2013-12-30T23:52:00.000000000', '2013-12-31T23:52:00.000000000'],\n", " dtype='datetime64[ns]')
[531804 values with dtype=float64]
[531804 values with dtype=float64]
array([[ 2.096768, 2.017261, 1.991643, ..., 79.021808, 77.097795,\n", " 72.737797],\n", " [ 2.114851, 2.038655, 2.011834, ..., 78.97254 , 76.850743,\n", " 72.222034],\n", " [ 2.133693, 2.060468, 2.037769, ..., 78.831859, 76.342968,\n", " 71.19606 ],\n", " ...,\n", " [ -68.94688 , -69.088987, -69.205023, ..., -273.920493, -250.937666,\n", " -217.887154],\n", " [ -69.69914 , -69.843726, -69.854963, ..., -271.870397, -251.031326,\n", " -218.747355],\n", " [ -70.450618, -70.569233, -70.539186, ..., -270.777917, -251.194593,\n", " -219.265093]])
array([[ 0.678754, 0.698426, 0.766338, ..., -68.94688 , -69.69914 ,\n", " -70.450618],\n", " [ 0.8405 , 0.829877, 0.882347, ..., -69.088987, -69.843726,\n", " -70.569233],\n", " [ 0.988098, 0.99794 , 1.055583, ..., -69.205023, -69.854963,\n", " -70.539186],\n", " ...,\n", " [ -10.606145, -10.267255, -9.733816, ..., -273.920493, -271.870397,\n", " -270.777917],\n", " [ 4.521696, 4.820326, 5.245298, ..., -250.937666, -251.031326,\n", " -251.194593],\n", " [ 19.504302, 19.613677, 19.71583 , ..., -217.887154, -218.747355,\n", " -219.265093]])
array([69.988649, 70.542177, 70.761973, 70.274043, 69.448014, 67.628528,\n", " 65.41385 , 62.317314, 59.459793, 58.735534, 59.124596, 60.828395,\n", " 65.737688, 71.231819, 73.848744, 73.602726, 71.811005, 70.677136,\n", " 70.919314, 70.814975, 69.577706, 68.630039, 68.778226, 69.052241,\n", " 67.901556, 64.373603])