{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "%matplotlib inline" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import netCDF4 as nc\n", "import matplotlib.pylab as plt\n", "import imp\n", "import csv\n", "import pandas as pd\n", "from io import StringIO" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "# Read data\n", "# AISM_VUB\n", "\n", "fname=\"../ComputeProjections4OneIceModel_SSPs/EnsembleSingleModelProjections/SL_wTd_nos_base_AISM_VUB_RCP70.nc\"\n", "ncf = nc.Dataset(fname, \"r\")\n", "\n", "Time = ncf.variables[\"Time\"][:]\n", "SL_wTd_nos_base_AISM_VUB_R0_RCP70 = ncf.variables[\"Antarctica\"][:]\n", "SL_wTd_nos_base_AISM_VUB_R1_RCP70 = ncf.variables[\"EAIS\"][:]\n", "SL_wTd_nos_base_AISM_VUB_R2_RCP70 = ncf.variables[\"Ross\"][:]\n", "SL_wTd_nos_base_AISM_VUB_R3_RCP70 = ncf.variables[\"Amundsen\"][:]\n", "SL_wTd_nos_base_AISM_VUB_R4_RCP70 = ncf.variables[\"Weddell\"][:]\n", "SL_wTd_nos_base_AISM_VUB_R5_RCP70 = ncf.variables[\"Peninsula\"][:]\n", "SL_wTd_nos_base_AISM_VUB_SU_RCP70 = SL_wTd_nos_base_AISM_VUB_R1_RCP70+SL_wTd_nos_base_AISM_VUB_R2_RCP70+SL_wTd_nos_base_AISM_VUB_R3_RCP70+SL_wTd_nos_base_AISM_VUB_R4_RCP70+SL_wTd_nos_base_AISM_VUB_R5_RCP70\n", "ncf.close()\n", "\n" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "# Read data\n", "# BISI_LBL\n", "\n", "fname=\"../ComputeProjections4OneIceModel_SSPs/EnsembleSingleModelProjections/SL_wTd_nos_base_BISI_LBL_RCP70.nc\"\n", "ncf = nc.Dataset(fname, \"r\")\n", "\n", "Time = ncf.variables[\"Time\"][:]\n", "SL_wTd_nos_base_BISI_LBL_R0_RCP70 = ncf.variables[\"Antarctica\"][:]\n", "SL_wTd_nos_base_BISI_LBL_R1_RCP70 = ncf.variables[\"EAIS\"][:]\n", "SL_wTd_nos_base_BISI_LBL_R2_RCP70 = ncf.variables[\"Ross\"][:]\n", "SL_wTd_nos_base_BISI_LBL_R3_RCP70 = ncf.variables[\"Amundsen\"][:]\n", "SL_wTd_nos_base_BISI_LBL_R4_RCP70 = ncf.variables[\"Weddell\"][:]\n", "SL_wTd_nos_base_BISI_LBL_R5_RCP70 = ncf.variables[\"Peninsula\"][:]\n", "SL_wTd_nos_base_BISI_LBL_SU_RCP70 = SL_wTd_nos_base_BISI_LBL_R1_RCP70+SL_wTd_nos_base_BISI_LBL_R2_RCP70+SL_wTd_nos_base_BISI_LBL_R3_RCP70+SL_wTd_nos_base_BISI_LBL_R4_RCP70+SL_wTd_nos_base_BISI_LBL_R5_RCP70\n", "ncf.close()\n", "\n" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "# Read data\n", "# CISM_NCA\n", "\n", "fname=\"../ComputeProjections4OneIceModel_SSPs/EnsembleSingleModelProjections/SL_wTd_nos_base_CISM_NCA_RCP70.nc\"\n", "ncf = nc.Dataset(fname, \"r\")\n", "\n", "Time = ncf.variables[\"Time\"][:]\n", "SL_wTd_nos_base_CISM_NCA_R0_RCP70 = ncf.variables[\"Antarctica\"][:]\n", "SL_wTd_nos_base_CISM_NCA_R1_RCP70 = ncf.variables[\"EAIS\"][:]\n", "SL_wTd_nos_base_CISM_NCA_R2_RCP70 = ncf.variables[\"Ross\"][:]\n", "SL_wTd_nos_base_CISM_NCA_R3_RCP70 = ncf.variables[\"Amundsen\"][:]\n", "SL_wTd_nos_base_CISM_NCA_R4_RCP70 = ncf.variables[\"Weddell\"][:]\n", "SL_wTd_nos_base_CISM_NCA_R5_RCP70 = ncf.variables[\"Peninsula\"][:]\n", "SL_wTd_nos_base_CISM_NCA_SU_RCP70 = SL_wTd_nos_base_CISM_NCA_R1_RCP70+SL_wTd_nos_base_CISM_NCA_R2_RCP70+SL_wTd_nos_base_CISM_NCA_R3_RCP70+SL_wTd_nos_base_CISM_NCA_R4_RCP70+SL_wTd_nos_base_CISM_NCA_R5_RCP70\n", "ncf.close()\n", "\n" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "# Read data\n", "# FETI_VUB\n", "\n", "fname=\"../ComputeProjections4OneIceModel_SSPs/EnsembleSingleModelProjections/SL_wTd_nos_base_FETI_VUB_RCP70.nc\"\n", "ncf = nc.Dataset(fname, \"r\")\n", "\n", "Time = ncf.variables[\"Time\"][:]\n", "SL_wTd_nos_base_FETI_VUB_R0_RCP70 = ncf.variables[\"Antarctica\"][:]\n", "SL_wTd_nos_base_FETI_VUB_R1_RCP70 = ncf.variables[\"EAIS\"][:]\n", "SL_wTd_nos_base_FETI_VUB_R2_RCP70 = ncf.variables[\"Ross\"][:]\n", "SL_wTd_nos_base_FETI_VUB_R3_RCP70 = ncf.variables[\"Amundsen\"][:]\n", "SL_wTd_nos_base_FETI_VUB_R4_RCP70 = ncf.variables[\"Weddell\"][:]\n", "SL_wTd_nos_base_FETI_VUB_R5_RCP70 = ncf.variables[\"Peninsula\"][:]\n", "SL_wTd_nos_base_FETI_VUB_SU_RCP70 = SL_wTd_nos_base_FETI_VUB_R1_RCP70+SL_wTd_nos_base_FETI_VUB_R2_RCP70+SL_wTd_nos_base_FETI_VUB_R3_RCP70+SL_wTd_nos_base_FETI_VUB_R4_RCP70+SL_wTd_nos_base_FETI_VUB_R5_RCP70\n", "ncf.close()\n", "\n" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "ename": "FileNotFoundError", "evalue": "[Errno 2] No such file or directory: b'../ComputeProjections4OneIceModel_SSPs/EnsembleSingleModelProjections/SL_wTd_nos_base_GRIS_LSC_RCP70.nc'", "output_type": "error", "traceback": [ "\u001b[1;31m---------------------------------------------------------------------------\u001b[0m", "\u001b[1;31mFileNotFoundError\u001b[0m Traceback (most recent call last)", "\u001b[1;32m\u001b[0m in \u001b[0;36m\u001b[1;34m\u001b[0m\n\u001b[0;32m 3\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 4\u001b[0m \u001b[0mfname\u001b[0m\u001b[1;33m=\u001b[0m\u001b[1;34m\"../ComputeProjections4OneIceModel_SSPs/EnsembleSingleModelProjections/SL_wTd_nos_base_GRIS_LSC_RCP70.nc\"\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m----> 5\u001b[1;33m \u001b[0mncf\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mnc\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mDataset\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mfname\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;34m\"r\"\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 6\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 7\u001b[0m \u001b[0mTime\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mncf\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mvariables\u001b[0m\u001b[1;33m[\u001b[0m\u001b[1;34m\"Time\"\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m[\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n", "\u001b[1;32mnetCDF4\\_netCDF4.pyx\u001b[0m in \u001b[0;36mnetCDF4._netCDF4.Dataset.__init__\u001b[1;34m()\u001b[0m\n", "\u001b[1;32mnetCDF4\\_netCDF4.pyx\u001b[0m in \u001b[0;36mnetCDF4._netCDF4._ensure_nc_success\u001b[1;34m()\u001b[0m\n", "\u001b[1;31mFileNotFoundError\u001b[0m: [Errno 2] No such file or directory: b'../ComputeProjections4OneIceModel_SSPs/EnsembleSingleModelProjections/SL_wTd_nos_base_GRIS_LSC_RCP70.nc'" ] } ], "source": [ "# Read data\n", "# GRIS_LSC\n", "\n", "fname=\"../ComputeProjections4OneIceModel_SSPs/EnsembleSingleModelProjections/SL_wTd_nos_base_GRIS_LSC_RCP70.nc\"\n", "ncf = nc.Dataset(fname, \"r\")\n", "\n", "Time = ncf.variables[\"Time\"][:]\n", "SL_wTd_nos_base_GRIS_LSC_R0_RCP70 = ncf.variables[\"Antarctica\"][:]\n", "SL_wTd_nos_base_GRIS_LSC_R1_RCP70 = ncf.variables[\"EAIS\"][:]\n", "SL_wTd_nos_base_GRIS_LSC_R2_RCP70 = ncf.variables[\"Ross\"][:]\n", "SL_wTd_nos_base_GRIS_LSC_R3_RCP70 = ncf.variables[\"Amundsen\"][:]\n", "SL_wTd_nos_base_GRIS_LSC_R4_RCP70 = ncf.variables[\"Weddell\"][:]\n", "SL_wTd_nos_base_GRIS_LSC_R5_RCP70 = ncf.variables[\"Peninsula\"][:]\n", "SL_wTd_nos_base_GRIS_LSC_SU_RCP70 = SL_wTd_nos_base_GRIS_LSC_R1_RCP70+SL_wTd_nos_base_GRIS_LSC_R2_RCP70+SL_wTd_nos_base_GRIS_LSC_R3_RCP70+SL_wTd_nos_base_GRIS_LSC_R4_RCP70+SL_wTd_nos_base_GRIS_LSC_R5_RCP70\n", "ncf.close()\n", "\n" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [ "# Read data\n", "# IMAU_VUB\n", "\n", "fname=\"../ComputeProjections4OneIceModel_SSPs/EnsembleSingleModelProjections/SL_wTd_nos_base_IMAU_VUB_RCP70.nc\"\n", "ncf = nc.Dataset(fname, \"r\")\n", "\n", "Time = ncf.variables[\"Time\"][:]\n", "SL_wTd_nos_base_IMAU_VUB_R0_RCP70 = ncf.variables[\"Antarctica\"][:]\n", "SL_wTd_nos_base_IMAU_VUB_R1_RCP70 = ncf.variables[\"EAIS\"][:]\n", "SL_wTd_nos_base_IMAU_VUB_R2_RCP70 = ncf.variables[\"Ross\"][:]\n", "SL_wTd_nos_base_IMAU_VUB_R3_RCP70 = ncf.variables[\"Amundsen\"][:]\n", "SL_wTd_nos_base_IMAU_VUB_R4_RCP70 = ncf.variables[\"Weddell\"][:]\n", "SL_wTd_nos_base_IMAU_VUB_R5_RCP70 = ncf.variables[\"Peninsula\"][:]\n", "SL_wTd_nos_base_IMAU_VUB_SU_RCP70 = SL_wTd_nos_base_IMAU_VUB_R1_RCP70+SL_wTd_nos_base_IMAU_VUB_R2_RCP70+SL_wTd_nos_base_IMAU_VUB_R3_RCP70+SL_wTd_nos_base_IMAU_VUB_R4_RCP70+SL_wTd_nos_base_IMAU_VUB_R5_RCP70\n", "ncf.close()\n", "\n" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [ "# Read data\n", "# ISSM_JPL\n", "\n", "fname=\"../ComputeProjections4OneIceModel_SSPs/EnsembleSingleModelProjections/SL_wTd_nos_base_ISSM_JPL_RCP70.nc\"\n", "ncf = nc.Dataset(fname, \"r\")\n", "\n", "Time = ncf.variables[\"Time\"][:]\n", "SL_wTd_nos_base_ISSM_JPL_R0_RCP70 = ncf.variables[\"Antarctica\"][:]\n", "SL_wTd_nos_base_ISSM_JPL_R1_RCP70 = ncf.variables[\"EAIS\"][:]\n", "SL_wTd_nos_base_ISSM_JPL_R2_RCP70 = ncf.variables[\"Ross\"][:]\n", "SL_wTd_nos_base_ISSM_JPL_R3_RCP70 = ncf.variables[\"Amundsen\"][:]\n", "SL_wTd_nos_base_ISSM_JPL_R4_RCP70 = ncf.variables[\"Weddell\"][:]\n", "SL_wTd_nos_base_ISSM_JPL_R5_RCP70 = ncf.variables[\"Peninsula\"][:]\n", "SL_wTd_nos_base_ISSM_JPL_SU_RCP70 = SL_wTd_nos_base_ISSM_JPL_R1_RCP70+SL_wTd_nos_base_ISSM_JPL_R2_RCP70+SL_wTd_nos_base_ISSM_JPL_R3_RCP70+SL_wTd_nos_base_ISSM_JPL_R4_RCP70+SL_wTd_nos_base_ISSM_JPL_R5_RCP70\n", "ncf.close()\n", "\n" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [ "# Read data\n", "# ISSM_UCI\n", "\n", "fname=\"../ComputeProjections4OneIceModel_SSPs/EnsembleSingleModelProjections/SL_wTd_nos_base_ISSM_UCI_RCP70.nc\"\n", "ncf = nc.Dataset(fname, \"r\")\n", "\n", "Time = ncf.variables[\"Time\"][:]\n", "SL_wTd_nos_base_ISSM_UCI_R0_RCP70 = ncf.variables[\"Antarctica\"][:]\n", "SL_wTd_nos_base_ISSM_UCI_R1_RCP70 = ncf.variables[\"EAIS\"][:]\n", "SL_wTd_nos_base_ISSM_UCI_R2_RCP70 = ncf.variables[\"Ross\"][:]\n", "SL_wTd_nos_base_ISSM_UCI_R3_RCP70 = ncf.variables[\"Amundsen\"][:]\n", "SL_wTd_nos_base_ISSM_UCI_R4_RCP70 = ncf.variables[\"Weddell\"][:]\n", "SL_wTd_nos_base_ISSM_UCI_R5_RCP70 = ncf.variables[\"Peninsula\"][:]\n", "SL_wTd_nos_base_ISSM_UCI_SU_RCP70 = SL_wTd_nos_base_ISSM_UCI_R1_RCP70+SL_wTd_nos_base_ISSM_UCI_R2_RCP70+SL_wTd_nos_base_ISSM_UCI_R3_RCP70+SL_wTd_nos_base_ISSM_UCI_R4_RCP70+SL_wTd_nos_base_ISSM_UCI_R5_RCP70\n", "ncf.close()\n", "\n" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [ "# Read data\n", "# MALI_LAN\n", "\n", "fname=\"../ComputeProjections4OneIceModel_SSPs/EnsembleSingleModelProjections/SL_wTd_nos_base_MALI_LAN_RCP70.nc\"\n", "ncf = nc.Dataset(fname, \"r\")\n", "\n", "Time = ncf.variables[\"Time\"][:]\n", "SL_wTd_nos_base_MALI_LAN_R0_RCP70 = ncf.variables[\"Antarctica\"][:]\n", "SL_wTd_nos_base_MALI_LAN_R1_RCP70 = ncf.variables[\"EAIS\"][:]\n", "SL_wTd_nos_base_MALI_LAN_R2_RCP70 = ncf.variables[\"Ross\"][:]\n", "SL_wTd_nos_base_MALI_LAN_R3_RCP70 = ncf.variables[\"Amundsen\"][:]\n", "SL_wTd_nos_base_MALI_LAN_R4_RCP70 = ncf.variables[\"Weddell\"][:]\n", "SL_wTd_nos_base_MALI_LAN_R5_RCP70 = ncf.variables[\"Peninsula\"][:]\n", "SL_wTd_nos_base_MALI_LAN_SU_RCP70 = SL_wTd_nos_base_MALI_LAN_R1_RCP70+SL_wTd_nos_base_MALI_LAN_R2_RCP70+SL_wTd_nos_base_MALI_LAN_R3_RCP70+SL_wTd_nos_base_MALI_LAN_R4_RCP70+SL_wTd_nos_base_MALI_LAN_R5_RCP70\n", "ncf.close()\n", "\n" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [ "# Read data\n", "# PISM_AWI\n", "\n", "fname=\"../ComputeProjections4OneIceModel_SSPs/EnsembleSingleModelProjections/SL_wTd_nos_base_PISM_AWI_RCP70.nc\"\n", "ncf = nc.Dataset(fname, \"r\")\n", "\n", "Time = ncf.variables[\"Time\"][:]\n", "SL_wTd_nos_base_PISM_AWI_R0_RCP70 = ncf.variables[\"Antarctica\"][:]\n", "SL_wTd_nos_base_PISM_AWI_R1_RCP70 = ncf.variables[\"EAIS\"][:]\n", "SL_wTd_nos_base_PISM_AWI_R2_RCP70 = ncf.variables[\"Ross\"][:]\n", "SL_wTd_nos_base_PISM_AWI_R3_RCP70 = ncf.variables[\"Amundsen\"][:]\n", "SL_wTd_nos_base_PISM_AWI_R4_RCP70 = ncf.variables[\"Weddell\"][:]\n", "SL_wTd_nos_base_PISM_AWI_R5_RCP70 = ncf.variables[\"Peninsula\"][:]\n", "SL_wTd_nos_base_PISM_AWI_SU_RCP70 = SL_wTd_nos_base_PISM_AWI_R1_RCP70+SL_wTd_nos_base_PISM_AWI_R2_RCP70+SL_wTd_nos_base_PISM_AWI_R3_RCP70+SL_wTd_nos_base_PISM_AWI_R4_RCP70+SL_wTd_nos_base_PISM_AWI_R5_RCP70\n", "ncf.close()\n", "\n" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [ "# Read data\n", "# PISM_DMI\n", "\n", "fname=\"../ComputeProjections4OneIceModel_SSPs/EnsembleSingleModelProjections/SL_wTd_nos_base_PISM_DMI_RCP70.nc\"\n", "ncf = nc.Dataset(fname, \"r\")\n", "\n", "Time = ncf.variables[\"Time\"][:]\n", "SL_wTd_nos_base_PISM_DMI_R0_RCP70 = ncf.variables[\"Antarctica\"][:]\n", "SL_wTd_nos_base_PISM_DMI_R1_RCP70 = ncf.variables[\"EAIS\"][:]\n", "SL_wTd_nos_base_PISM_DMI_R2_RCP70 = ncf.variables[\"Ross\"][:]\n", "SL_wTd_nos_base_PISM_DMI_R3_RCP70 = ncf.variables[\"Amundsen\"][:]\n", "SL_wTd_nos_base_PISM_DMI_R4_RCP70 = ncf.variables[\"Weddell\"][:]\n", "SL_wTd_nos_base_PISM_DMI_R5_RCP70 = ncf.variables[\"Peninsula\"][:]\n", "SL_wTd_nos_base_PISM_DMI_SU_RCP70 = SL_wTd_nos_base_PISM_DMI_R1_RCP70+SL_wTd_nos_base_PISM_DMI_R2_RCP70+SL_wTd_nos_base_PISM_DMI_R3_RCP70+SL_wTd_nos_base_PISM_DMI_R4_RCP70+SL_wTd_nos_base_PISM_DMI_R5_RCP70\n", "ncf.close()\n", "\n" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [ "# Read data\n", "# PISM_PIK\n", "\n", "fname=\"../ComputeProjections4OneIceModel_SSPs/EnsembleSingleModelProjections/SL_wTd_nos_base_PISM_PIK_RCP70.nc\"\n", "ncf = nc.Dataset(fname, \"r\")\n", "\n", "Time = ncf.variables[\"Time\"][:]\n", "SL_wTd_nos_base_PISM_PIK_R0_RCP70 = ncf.variables[\"Antarctica\"][:]\n", "SL_wTd_nos_base_PISM_PIK_R1_RCP70 = ncf.variables[\"EAIS\"][:]\n", "SL_wTd_nos_base_PISM_PIK_R2_RCP70 = ncf.variables[\"Ross\"][:]\n", "SL_wTd_nos_base_PISM_PIK_R3_RCP70 = ncf.variables[\"Amundsen\"][:]\n", "SL_wTd_nos_base_PISM_PIK_R4_RCP70 = ncf.variables[\"Weddell\"][:]\n", "SL_wTd_nos_base_PISM_PIK_R5_RCP70 = ncf.variables[\"Peninsula\"][:]\n", "SL_wTd_nos_base_PISM_PIK_SU_RCP70 = SL_wTd_nos_base_PISM_PIK_R1_RCP70+SL_wTd_nos_base_PISM_PIK_R2_RCP70+SL_wTd_nos_base_PISM_PIK_R3_RCP70+SL_wTd_nos_base_PISM_PIK_R4_RCP70+SL_wTd_nos_base_PISM_PIK_R5_RCP70\n", "ncf.close()\n", "\n" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [ "# Read data\n", "# PISM_VUW\n", "\n", "fname=\"../ComputeProjections4OneIceModel_SSPs/EnsembleSingleModelProjections/SL_wTd_nos_base_PISM_VUW_RCP70.nc\"\n", "ncf = nc.Dataset(fname, \"r\")\n", "\n", "Time = ncf.variables[\"Time\"][:]\n", "SL_wTd_nos_base_PISM_VUW_R0_RCP70 = ncf.variables[\"Antarctica\"][:]\n", "SL_wTd_nos_base_PISM_VUW_R1_RCP70 = ncf.variables[\"EAIS\"][:]\n", "SL_wTd_nos_base_PISM_VUW_R2_RCP70 = ncf.variables[\"Ross\"][:]\n", "SL_wTd_nos_base_PISM_VUW_R3_RCP70 = ncf.variables[\"Amundsen\"][:]\n", "SL_wTd_nos_base_PISM_VUW_R4_RCP70 = ncf.variables[\"Weddell\"][:]\n", "SL_wTd_nos_base_PISM_VUW_R5_RCP70 = ncf.variables[\"Peninsula\"][:]\n", "SL_wTd_nos_base_PISM_VUW_SU_RCP70 = SL_wTd_nos_base_PISM_VUW_R1_RCP70+SL_wTd_nos_base_PISM_VUW_R2_RCP70+SL_wTd_nos_base_PISM_VUW_R3_RCP70+SL_wTd_nos_base_PISM_VUW_R4_RCP70+SL_wTd_nos_base_PISM_VUW_R5_RCP70\n", "ncf.close()\n", "\n" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [ "# Read data\n", "# PS3D_PSU\n", "\n", "fname=\"../ComputeProjections4OneIceModel_SSPs/EnsembleSingleModelProjections/SL_wTd_nos_base_PS3D_PSU_RCP70.nc\"\n", "ncf = nc.Dataset(fname, \"r\")\n", "\n", "Time = ncf.variables[\"Time\"][:]\n", "SL_wTd_nos_base_PS3D_PSU_R0_RCP70 = ncf.variables[\"Antarctica\"][:]\n", "SL_wTd_nos_base_PS3D_PSU_R1_RCP70 = ncf.variables[\"EAIS\"][:]\n", "SL_wTd_nos_base_PS3D_PSU_R2_RCP70 = ncf.variables[\"Ross\"][:]\n", "SL_wTd_nos_base_PS3D_PSU_R3_RCP70 = ncf.variables[\"Amundsen\"][:]\n", "SL_wTd_nos_base_PS3D_PSU_R4_RCP70 = ncf.variables[\"Weddell\"][:]\n", "SL_wTd_nos_base_PS3D_PSU_R5_RCP70 = ncf.variables[\"Peninsula\"][:]\n", "SL_wTd_nos_base_PS3D_PSU_SU_RCP70 = SL_wTd_nos_base_PS3D_PSU_R1_RCP70+SL_wTd_nos_base_PS3D_PSU_R2_RCP70+SL_wTd_nos_base_PS3D_PSU_R3_RCP70+SL_wTd_nos_base_PS3D_PSU_R4_RCP70+SL_wTd_nos_base_PS3D_PSU_R5_RCP70\n", "ncf.close()\n", "\n" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [ "# Read data\n", "# SICO_UHO\n", "\n", "fname=\"../ComputeProjections4OneIceModel_SSPs/EnsembleSingleModelProjections/SL_wTd_nos_base_SICO_UHO_RCP70.nc\"\n", "ncf = nc.Dataset(fname, \"r\")\n", "\n", "Time = ncf.variables[\"Time\"][:]\n", "SL_wTd_nos_base_SICO_UHO_R0_RCP70 = ncf.variables[\"Antarctica\"][:]\n", "SL_wTd_nos_base_SICO_UHO_R1_RCP70 = ncf.variables[\"EAIS\"][:]\n", "SL_wTd_nos_base_SICO_UHO_R2_RCP70 = ncf.variables[\"Ross\"][:]\n", "SL_wTd_nos_base_SICO_UHO_R3_RCP70 = ncf.variables[\"Amundsen\"][:]\n", "SL_wTd_nos_base_SICO_UHO_R4_RCP70 = ncf.variables[\"Weddell\"][:]\n", "SL_wTd_nos_base_SICO_UHO_R5_RCP70 = ncf.variables[\"Peninsula\"][:]\n", "SL_wTd_nos_base_SICO_UHO_SU_RCP70 = SL_wTd_nos_base_SICO_UHO_R1_RCP70+SL_wTd_nos_base_SICO_UHO_R2_RCP70+SL_wTd_nos_base_SICO_UHO_R3_RCP70+SL_wTd_nos_base_SICO_UHO_R4_RCP70+SL_wTd_nos_base_SICO_UHO_R5_RCP70\n", "ncf.close()\n", "\n" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [ "# Read data\n", "# UA_UNN\n", "\n", "fname=\"../ComputeProjections4OneIceModel_SSPs/EnsembleSingleModelProjections/SL_wTd_nos_base_UA_UNN_RCP70.nc\"\n", "ncf = nc.Dataset(fname, \"r\")\n", "\n", "Time = ncf.variables[\"Time\"][:]\n", "SL_wTd_nos_base_UA_UNN_R0_RCP70 = ncf.variables[\"Antarctica\"][:]\n", "SL_wTd_nos_base_UA_UNN_R1_RCP70 = ncf.variables[\"EAIS\"][:]\n", "SL_wTd_nos_base_UA_UNN_R2_RCP70 = ncf.variables[\"Ross\"][:]\n", "SL_wTd_nos_base_UA_UNN_R3_RCP70 = ncf.variables[\"Amundsen\"][:]\n", "SL_wTd_nos_base_UA_UNN_R4_RCP70 = ncf.variables[\"Weddell\"][:]\n", "SL_wTd_nos_base_UA_UNN_R5_RCP70 = ncf.variables[\"Peninsula\"][:]\n", "SL_wTd_nos_base_UA_UNN_SU_RCP70 = SL_wTd_nos_base_UA_UNN_R1_RCP70+SL_wTd_nos_base_UA_UNN_R2_RCP70+SL_wTd_nos_base_UA_UNN_R3_RCP70+SL_wTd_nos_base_UA_UNN_R4_RCP70+SL_wTd_nos_base_UA_UNN_R5_RCP70\n", "ncf.close()\n", "\n" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "ename": "NameError", "evalue": "name 'SL_wTd_nos_base_GRIS_LSC_SU_RCP70' is not defined", "output_type": "error", "traceback": [ "\u001b[1;31m---------------------------------------------------------------------------\u001b[0m", "\u001b[1;31mNameError\u001b[0m Traceback (most recent call last)", "\u001b[1;32m\u001b[0m in \u001b[0;36m\u001b[1;34m\u001b[0m\n\u001b[1;32m----> 1\u001b[1;33m \u001b[0mSL_wTd_nos_base_RCP70\u001b[0m \u001b[1;33m=\u001b[0m\u001b[0mnp\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mvstack\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m[\u001b[0m\u001b[0mSL_wTd_nos_base_AISM_VUB_SU_RCP70\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mSL_wTd_nos_base_BISI_LBL_SU_RCP70\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mSL_wTd_nos_base_CISM_NCA_SU_RCP70\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mSL_wTd_nos_base_FETI_VUB_SU_RCP70\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mSL_wTd_nos_base_GRIS_LSC_SU_RCP70\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mSL_wTd_nos_base_IMAU_VUB_SU_RCP70\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mSL_wTd_nos_base_ISSM_JPL_SU_RCP70\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mSL_wTd_nos_base_ISSM_UCI_SU_RCP70\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mSL_wTd_nos_base_MALI_LAN_SU_RCP70\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mSL_wTd_nos_base_PISM_AWI_SU_RCP70\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mSL_wTd_nos_base_PISM_DMI_SU_RCP70\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mSL_wTd_nos_base_PISM_PIK_SU_RCP70\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mSL_wTd_nos_base_PISM_VUW_SU_RCP70\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mSL_wTd_nos_base_PS3D_PSU_SU_RCP70\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mSL_wTd_nos_base_SICO_UHO_SU_RCP70\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mSL_wTd_nos_base_UA_UNN_SU_RCP70\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 2\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 3\u001b[0m \u001b[0mSL_wTd_nos_base_R1_RCP70\u001b[0m \u001b[1;33m=\u001b[0m\u001b[0mnp\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mvstack\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m[\u001b[0m\u001b[0mSL_wTd_nos_base_AISM_VUB_R1_RCP70\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mSL_wTd_nos_base_BISI_LBL_R1_RCP70\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mSL_wTd_nos_base_CISM_NCA_R1_RCP70\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mSL_wTd_nos_base_FETI_VUB_R1_RCP70\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mSL_wTd_nos_base_GRIS_LSC_R1_RCP70\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mSL_wTd_nos_base_IMAU_VUB_R1_RCP70\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mSL_wTd_nos_base_ISSM_JPL_R1_RCP70\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mSL_wTd_nos_base_ISSM_UCI_R1_RCP70\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mSL_wTd_nos_base_MALI_LAN_R1_RCP70\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mSL_wTd_nos_base_PISM_AWI_R1_RCP70\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mSL_wTd_nos_base_PISM_DMI_R1_RCP70\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mSL_wTd_nos_base_PISM_PIK_R1_RCP70\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mSL_wTd_nos_base_PISM_VUW_R1_RCP70\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mSL_wTd_nos_base_PS3D_PSU_R1_RCP70\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mSL_wTd_nos_base_SICO_UHO_R1_RCP70\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mSL_wTd_nos_base_UA_UNN_R1_RCP70\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 4\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 5\u001b[0m \u001b[0mSL_wTd_nos_base_R2_RCP70\u001b[0m \u001b[1;33m=\u001b[0m\u001b[0mnp\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mvstack\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m[\u001b[0m\u001b[0mSL_wTd_nos_base_AISM_VUB_R2_RCP70\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mSL_wTd_nos_base_BISI_LBL_R2_RCP70\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mSL_wTd_nos_base_CISM_NCA_R2_RCP70\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mSL_wTd_nos_base_FETI_VUB_R2_RCP70\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mSL_wTd_nos_base_GRIS_LSC_R2_RCP70\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mSL_wTd_nos_base_IMAU_VUB_R2_RCP70\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mSL_wTd_nos_base_ISSM_JPL_R2_RCP70\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mSL_wTd_nos_base_ISSM_UCI_R2_RCP70\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mSL_wTd_nos_base_MALI_LAN_R2_RCP70\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mSL_wTd_nos_base_PISM_AWI_R2_RCP70\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mSL_wTd_nos_base_PISM_DMI_R2_RCP70\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mSL_wTd_nos_base_PISM_PIK_R2_RCP70\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mSL_wTd_nos_base_PISM_VUW_R2_RCP70\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mSL_wTd_nos_base_PS3D_PSU_R2_RCP70\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mSL_wTd_nos_base_SICO_UHO_R2_RCP70\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mSL_wTd_nos_base_UA_UNN_R2_RCP70\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n", "\u001b[1;31mNameError\u001b[0m: name 'SL_wTd_nos_base_GRIS_LSC_SU_RCP70' is not defined" ] } ], "source": [ "SL_wTd_nos_base_RCP70 =np.vstack([SL_wTd_nos_base_AISM_VUB_SU_RCP70,SL_wTd_nos_base_BISI_LBL_SU_RCP70,SL_wTd_nos_base_CISM_NCA_SU_RCP70,SL_wTd_nos_base_FETI_VUB_SU_RCP70,SL_wTd_nos_base_GRIS_LSC_SU_RCP70,SL_wTd_nos_base_IMAU_VUB_SU_RCP70,SL_wTd_nos_base_ISSM_JPL_SU_RCP70,SL_wTd_nos_base_ISSM_UCI_SU_RCP70,SL_wTd_nos_base_MALI_LAN_SU_RCP70,SL_wTd_nos_base_PISM_AWI_SU_RCP70,SL_wTd_nos_base_PISM_DMI_SU_RCP70,SL_wTd_nos_base_PISM_PIK_SU_RCP70,SL_wTd_nos_base_PISM_VUW_SU_RCP70,SL_wTd_nos_base_PS3D_PSU_SU_RCP70,SL_wTd_nos_base_SICO_UHO_SU_RCP70,SL_wTd_nos_base_UA_UNN_SU_RCP70])\n", "\n", "SL_wTd_nos_base_R1_RCP70 =np.vstack([SL_wTd_nos_base_AISM_VUB_R1_RCP70,SL_wTd_nos_base_BISI_LBL_R1_RCP70,SL_wTd_nos_base_CISM_NCA_R1_RCP70,SL_wTd_nos_base_FETI_VUB_R1_RCP70,SL_wTd_nos_base_GRIS_LSC_R1_RCP70,SL_wTd_nos_base_IMAU_VUB_R1_RCP70,SL_wTd_nos_base_ISSM_JPL_R1_RCP70,SL_wTd_nos_base_ISSM_UCI_R1_RCP70,SL_wTd_nos_base_MALI_LAN_R1_RCP70,SL_wTd_nos_base_PISM_AWI_R1_RCP70,SL_wTd_nos_base_PISM_DMI_R1_RCP70,SL_wTd_nos_base_PISM_PIK_R1_RCP70,SL_wTd_nos_base_PISM_VUW_R1_RCP70,SL_wTd_nos_base_PS3D_PSU_R1_RCP70,SL_wTd_nos_base_SICO_UHO_R1_RCP70,SL_wTd_nos_base_UA_UNN_R1_RCP70])\n", "\n", "SL_wTd_nos_base_R2_RCP70 =np.vstack([SL_wTd_nos_base_AISM_VUB_R2_RCP70,SL_wTd_nos_base_BISI_LBL_R2_RCP70,SL_wTd_nos_base_CISM_NCA_R2_RCP70,SL_wTd_nos_base_FETI_VUB_R2_RCP70,SL_wTd_nos_base_GRIS_LSC_R2_RCP70,SL_wTd_nos_base_IMAU_VUB_R2_RCP70,SL_wTd_nos_base_ISSM_JPL_R2_RCP70,SL_wTd_nos_base_ISSM_UCI_R2_RCP70,SL_wTd_nos_base_MALI_LAN_R2_RCP70,SL_wTd_nos_base_PISM_AWI_R2_RCP70,SL_wTd_nos_base_PISM_DMI_R2_RCP70,SL_wTd_nos_base_PISM_PIK_R2_RCP70,SL_wTd_nos_base_PISM_VUW_R2_RCP70,SL_wTd_nos_base_PS3D_PSU_R2_RCP70,SL_wTd_nos_base_SICO_UHO_R2_RCP70,SL_wTd_nos_base_UA_UNN_R2_RCP70])\n", "\n", "SL_wTd_nos_base_R3_RCP70 =np.vstack([SL_wTd_nos_base_AISM_VUB_R3_RCP70,SL_wTd_nos_base_BISI_LBL_R3_RCP70,SL_wTd_nos_base_CISM_NCA_R3_RCP70,SL_wTd_nos_base_FETI_VUB_R3_RCP70,SL_wTd_nos_base_GRIS_LSC_R3_RCP70,SL_wTd_nos_base_IMAU_VUB_R3_RCP70,SL_wTd_nos_base_ISSM_JPL_R3_RCP70,SL_wTd_nos_base_ISSM_UCI_R3_RCP70,SL_wTd_nos_base_MALI_LAN_R3_RCP70,SL_wTd_nos_base_PISM_AWI_R3_RCP70,SL_wTd_nos_base_PISM_DMI_R3_RCP70,SL_wTd_nos_base_PISM_PIK_R3_RCP70,SL_wTd_nos_base_PISM_VUW_R3_RCP70,SL_wTd_nos_base_PS3D_PSU_R3_RCP70,SL_wTd_nos_base_SICO_UHO_R3_RCP70,SL_wTd_nos_base_UA_UNN_R3_RCP70])\n", "\n", "SL_wTd_nos_base_R4_RCP70 =np.vstack([SL_wTd_nos_base_AISM_VUB_R4_RCP70,SL_wTd_nos_base_BISI_LBL_R4_RCP70,SL_wTd_nos_base_CISM_NCA_R4_RCP70,SL_wTd_nos_base_FETI_VUB_R4_RCP70,SL_wTd_nos_base_GRIS_LSC_R4_RCP70,SL_wTd_nos_base_IMAU_VUB_R4_RCP70,SL_wTd_nos_base_ISSM_JPL_R4_RCP70,SL_wTd_nos_base_ISSM_UCI_R4_RCP70,SL_wTd_nos_base_MALI_LAN_R4_RCP70,SL_wTd_nos_base_PISM_AWI_R4_RCP70,SL_wTd_nos_base_PISM_DMI_R4_RCP70,SL_wTd_nos_base_PISM_PIK_R4_RCP70,SL_wTd_nos_base_PISM_VUW_R4_RCP70,SL_wTd_nos_base_PS3D_PSU_R4_RCP70,SL_wTd_nos_base_SICO_UHO_R4_RCP70,SL_wTd_nos_base_UA_UNN_R4_RCP70])\n", "\n", "SL_wTd_nos_base_R5_RCP70 =np.vstack([SL_wTd_nos_base_AISM_VUB_R5_RCP70,SL_wTd_nos_base_BISI_LBL_R5_RCP70,SL_wTd_nos_base_CISM_NCA_R5_RCP70,SL_wTd_nos_base_FETI_VUB_R5_RCP70,SL_wTd_nos_base_GRIS_LSC_R5_RCP70,SL_wTd_nos_base_IMAU_VUB_R5_RCP70,SL_wTd_nos_base_ISSM_JPL_R5_RCP70,SL_wTd_nos_base_ISSM_UCI_R5_RCP70,SL_wTd_nos_base_MALI_LAN_R5_RCP70,SL_wTd_nos_base_PISM_AWI_R5_RCP70,SL_wTd_nos_base_PISM_DMI_R5_RCP70,SL_wTd_nos_base_PISM_PIK_R5_RCP70,SL_wTd_nos_base_PISM_VUW_R5_RCP70,SL_wTd_nos_base_PS3D_PSU_R5_RCP70,SL_wTd_nos_base_SICO_UHO_R5_RCP70,SL_wTd_nos_base_UA_UNN_R5_RCP70])\n" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "ename": "NameError", "evalue": "name 'SL_wTd_nos_base_RCP70' is not defined", "output_type": "error", "traceback": [ "\u001b[1;31m---------------------------------------------------------------------------\u001b[0m", "\u001b[1;31mNameError\u001b[0m Traceback (most recent call last)", "\u001b[1;32m\u001b[0m in \u001b[0;36m\u001b[1;34m\u001b[0m\n\u001b[0;32m 1\u001b[0m \u001b[1;31m# compute cumulative probability distributions\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 2\u001b[0m \u001b[0mcdfnum\u001b[0m \u001b[1;33m=\u001b[0m \u001b[1;36m1000\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m----> 3\u001b[1;33m \u001b[0mcdfstep\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mint\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mlen\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mSL_wTd_nos_base_RCP70\u001b[0m\u001b[1;33m[\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m,\u001b[0m\u001b[1;36m0\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m/\u001b[0m\u001b[0mcdfnum\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 4\u001b[0m \u001b[0mprint\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mcdfstep\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 5\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n", "\u001b[1;31mNameError\u001b[0m: name 'SL_wTd_nos_base_RCP70' is not defined" ] } ], "source": [ "# compute cumulative probability distributions\n", "cdfnum = 1000\n", "cdfstep = int(len(SL_wTd_nos_base_RCP70[:,0])/cdfnum)\n", "print(cdfstep)\n", "\n", "SL_wTd_nos_base_R0_RCP70_cdf = [0] * (cdfnum+1)\n", "for t in range(len(SL_wTd_nos_base_RCP70[1,:])):\n", " print(\"R0: \",t)\n", " sortind = np.argsort(SL_wTd_nos_base_RCP70[:,t])\n", "\n", " slcdf = []\n", " slcdf.append(0)\n", " for i in range(1,cdfnum+1):\n", " cdfval = float(i/cdfnum)\n", " slval = SL_wTd_nos_base_RCP70[sortind[i*cdfstep],t]\n", " slcdf.append(slval)\n", " SL_wTd_nos_base_R0_RCP70_cdf=np.vstack([SL_wTd_nos_base_R0_RCP70_cdf, slcdf])\n", "\n", "\n", "SL_wTd_nos_base_R1_RCP70_cdf = [0] * (cdfnum+1)\n", "for t in range(len(SL_wTd_nos_base_R1_RCP70[1,:])):\n", " print(\"R1: \",t)\n", " sortind = np.argsort(SL_wTd_nos_base_R1_RCP70[:,t])\n", "\n", " slcdf = []\n", " slcdf.append(0)\n", " for i in range(1,cdfnum+1):\n", " cdfval = float(i/cdfnum)\n", " slval = SL_wTd_nos_base_R1_RCP70[sortind[i*cdfstep],t]\n", " slcdf.append(slval)\n", " SL_wTd_nos_base_R1_RCP70_cdf=np.vstack([SL_wTd_nos_base_R1_RCP70_cdf, slcdf])\n", "\n", "SL_wTd_nos_base_R2_RCP70_cdf = [0] * (cdfnum+1)\n", "for t in range(len(SL_wTd_nos_base_R2_RCP70[1,:])):\n", " print(\"R2: \",t)\n", " sortind = np.argsort(SL_wTd_nos_base_R2_RCP70[:,t])\n", "\n", " slcdf = []\n", " slcdf.append(0)\n", " for i in range(1,cdfnum+1):\n", " cdfval = float(i/cdfnum)\n", " slval = SL_wTd_nos_base_R2_RCP70[sortind[i*cdfstep],t]\n", " slcdf.append(slval)\n", " SL_wTd_nos_base_R2_RCP70_cdf=np.vstack([SL_wTd_nos_base_R2_RCP70_cdf, slcdf])\n", "\n", "SL_wTd_nos_base_R3_RCP70_cdf = [0] * (cdfnum+1)\n", "for t in range(len(SL_wTd_nos_base_R3_RCP70[1,:])):\n", " print(\"R3: \",t)\n", " sortind = np.argsort(SL_wTd_nos_base_R3_RCP70[:,t])\n", "\n", " slcdf = []\n", " slcdf.append(0)\n", " for i in range(1,cdfnum+1):\n", " cdfval = float(i/cdfnum)\n", " slval = SL_wTd_nos_base_R3_RCP70[sortind[i*cdfstep],t]\n", " slcdf.append(slval)\n", " SL_wTd_nos_base_R3_RCP70_cdf=np.vstack([SL_wTd_nos_base_R3_RCP70_cdf, slcdf])\n", "\n", "SL_wTd_nos_base_R4_RCP70_cdf = [0] * (cdfnum+1)\n", "for t in range(len(SL_wTd_nos_base_R4_RCP70[1,:])):\n", " print(\"R4: \",t)\n", " sortind = np.argsort(SL_wTd_nos_base_R4_RCP70[:,t])\n", "\n", " slcdf = []\n", " slcdf.append(0)\n", " for i in range(1,cdfnum+1):\n", " cdfval = float(i/cdfnum)\n", " slval = SL_wTd_nos_base_R4_RCP70[sortind[i*cdfstep],t]\n", " slcdf.append(slval)\n", " SL_wTd_nos_base_R4_RCP70_cdf=np.vstack([SL_wTd_nos_base_R4_RCP70_cdf, slcdf])\n", "\n", "SL_wTd_nos_base_R5_RCP70_cdf = [0] * (cdfnum+1)\n", "for t in range(len(SL_wTd_nos_base_R5_RCP70[1,:])):\n", " print(\"R5: \",t)\n", " sortind = np.argsort(SL_wTd_nos_base_R5_RCP70[:,t])\n", "\n", " slcdf = []\n", " slcdf.append(0)\n", " for i in range(1,cdfnum+1):\n", " cdfval = float(i/cdfnum)\n", " slval = SL_wTd_nos_base_R5_RCP70[sortind[i*cdfstep],t]\n", " slcdf.append(slval)\n", " SL_wTd_nos_base_R5_RCP70_cdf=np.vstack([SL_wTd_nos_base_R5_RCP70_cdf, slcdf])\n", "\n", "\n", "Percentile = np.arange(0,float((cdfnum+1)/cdfnum),float(1/cdfnum))\n", "\n", "# write cdfs\n", "ncfile = nc.Dataset('Cdfs/SL_wTd_nos_base_RCP70_cdf.nc','w', format='NETCDF4')\n", "ncfile.createDimension('Time', None)\n", "ncfile.createDimension('Percentile', None)\n", "\n", "SL_wTd_weighted_base_R0 = ncfile.createVariable('Antarctica', 'f4', ('Time','Percentile'))\n", "SL_wTd_weighted_base_R1 = ncfile.createVariable('EAIS', 'f4', ('Time','Percentile'))\n", "SL_wTd_weighted_base_R2 = ncfile.createVariable('Ross', 'f4', ('Time','Percentile'))\n", "SL_wTd_weighted_base_R3 = ncfile.createVariable('Amundsen', 'f4', ('Time','Percentile'))\n", "SL_wTd_weighted_base_R4 = ncfile.createVariable('Weddell', 'f4', ('Time','Percentile'))\n", "SL_wTd_weighted_base_R5 = ncfile.createVariable('Peninsula', 'f4', ('Time','Percentile'))\n", "p = ncfile.createVariable('Percentile', 'f4', 'Percentile')\n", "t = ncfile.createVariable('Time', 'f4', 'Time')\n", "\n", "t[:] = Time\n", "t.units = 'years'\n", "SL_wTd_weighted_base_R0[:,:] = SL_wTd_nos_base_R0_RCP70_cdf\n", "SL_wTd_weighted_base_R1[:,:] = SL_wTd_nos_base_R1_RCP70_cdf\n", "SL_wTd_weighted_base_R2[:,:] = SL_wTd_nos_base_R2_RCP70_cdf\n", "SL_wTd_weighted_base_R3[:,:] = SL_wTd_nos_base_R3_RCP70_cdf\n", "SL_wTd_weighted_base_R4[:,:] = SL_wTd_nos_base_R4_RCP70_cdf\n", "SL_wTd_weighted_base_R5[:,:] = SL_wTd_nos_base_R5_RCP70_cdf\n", "p[:] = Percentile\n", "\n", "SL_wTd_weighted_base_R0.units = 'meter'\n", "SL_wTd_weighted_base_R1.units = 'meter'\n", "SL_wTd_weighted_base_R2.units = 'meter'\n", "SL_wTd_weighted_base_R3.units = 'meter'\n", "SL_wTd_weighted_base_R4.units = 'meter'\n", "SL_wTd_weighted_base_R5.units = 'meter'\n", "\n", "p.units = 'percent'\n", "\n", "ncfile.close()" ] }, { "cell_type": "code", "execution_count": 24, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "200\n", "200\n" ] }, { "data": { "text/plain": [ "[]" ] }, "execution_count": 24, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD8CAYAAACMwORRAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJzs3Xl8XNV9///XubPPaF8tS7bkfcfYmMUQ9h0CDhACZC8ktM1WkjTfpE2b9pv+0rRJm6TZQ9YvSQMlCSFmC5BiFhPb2Abvq7zLkq1dmn3mzj2/P+bKHgnZlu1ZJPnzfDyu7zJ35n58NXrP1bl3zlVaa4QQQowvRqELEEIIkX0S7kIIMQ5JuAshxDgk4S6EEOOQhLsQQoxDEu5CCDEOSbgLIcQ4JOEuhBDjkIS7EEKMQ85Cbbiqqko3NTUVavNCCDEmrV+/vlNrXX2q9QoW7k1NTaxbt65QmxdCiDFJKXVgJOtJs4wQQoxDEu5CCDEOSbgLIcQ4JOEuhBDjkIS7EEKMQxLuQggxDkm4CyHEOCThLoQQeaItTe8ze0kcDuV8WxLuQgiRJ4mWIKHXDpNsj+R8WxLuQgiRJ9EtneBQ+GZX5HxbEu5CCJEHWmuimzvxTi/D8OW+5xcJdyGEyINka5hUTxzf/Kq8bE/CXQgh8iC6pRMM8M6tzMv2JNyFECLHBppkPFPLcARcedmmhLsQQuRYsi2M2RnFtyA/TTIg4S6EEDkX3ZRukvHNy0+TDEi4CyFETmmtiWzuwDOtDEeRO2/blXAXQogcSraGSXXF8C845Z3xskrCXQghcii6qSN9lUwem2RAwl0IIXIm3STTiWd6ed6ukhlwynBXSv1MKdWulNpygseVUurbSqlmpdQmpdTi7JcphBBjT/JwiFR3DH8er5IZMJIj918AN53k8ZuBGfbwIPCDsy9LCCHGvsimTjBUXq+SGXDKcNdavwp0n2SVZcAjOm01UKaUqstWgUIIMRalv7jUgXdGGYY/v00ykJ0293rgUMZ8i71MCCHOWcmWULovmQI0yUB2wl0Ns0wPu6JSDyql1iml1nV0dGRh00IIMTpFNneku/fNU18yQ2Uj3FuASRnzDUDrcCtqrR/WWi/RWi+prs7vNZ9CCJEvWmuim+zufQvQJAPZCfflwAftq2YuAfq01m1ZeF0hhBiTEoeCpHrj+M4r3EHsKXuMV0o9ClwFVCmlWoB/AlwAWusfAs8CtwDNQAT4i1wVK4QQY0F0U2dBm2RgBOGutb7vFI9r4ONZq0gIIcYwbdl3XJpRnpc7Lp2IfENVCCGyKHEoSKovju+8wlwlM0DCXQghsii6qbBXyQyQcBdCiCzRKU1kUwfemeUY3sI1yYCEuxBCZE18Ty9WMIl/cW2hS5FwF0KIbAm/eRTldeKbU1HoUiTchRAiG6yYSWxrF/6FVShn4aO18BUIIcQ4EN3SiU5ao6JJBiTchRAiKyJvtuOs9OKeXFzoUgAJdyGEOGtmT4z43j78i2tRari+FPNPwl0IIc5SZEM7AP5FNQWu5DgJdyGEOAtaayJvtuNuKsFZ4S10OcdIuAshxFlItoQwO6IERsmJ1AES7kIIcRbCbx4Fpyp4XzJDSbgLIcQZ0qZFdGMHvrmVBe9uYCgJdyGEOEOxnd1YEXPUXNueScJdCCHOUPjNdowiF94Z5YUu5W0k3IUQ4gxYkSSxHd34z69BOUbHte2ZJNyFEOIMRDZ1QErjXzx6rm3PJOEuhBBnIPJmO64Jflx1gUKXMiwJdyGEOE3JjgiJg8FR1d3AUBLuQghxmiJvtYMC//nVhS7lhCTchRDiNGgr3d2AZ0Y5jhJPocs5IQl3IYQ4DYn9faR64wRG6YnUARLuQghxGkJvHEF5HXjnVha6lJOScBdCiBFKhRJEN3cSWFyL4XYUupyTknAXQogRCq87CilN4OIJhS7llCTchRBiBLSlCb9xBPeUUly1o/Pa9kwS7kIIMQLx3T2kumMUXVJX6FJGZEThrpS6SSm1UynVrJT6wjCPT1ZKrVBKvaWU2qSUuiX7pQohROGEVrdhFLnwzRvdJ1IHnDLclVIO4HvAzcBc4D6l1Nwhq/0D8LjWehFwL/D9bBcqhBCFYvbGiO3oJrBkAso5Nho8RlLlRUCz1nqv1joBPAYsG7KOBkrs6VKgNXslCiFEYYXfOAJA4KLRfyJ1wEjCvR44lDHfYi/L9M/A+5VSLcCzwCeHeyGl1INKqXVKqXUdHR1nUK4QQuSXTlmE1x7BO6tiVN0A+1RGEu7D9Yqjh8zfB/xCa90A3AL8Uin1ttfWWj+stV6itV5SXT16+2QQQogB0W1dWMEkgTFyInXASMK9BZiUMd/A25tdHgAeB9BarwK8wOi6W6wQQpyB8Oo2HGUevDNH392WTmYk4b4WmKGUmqKUcpM+Ybp8yDoHgWsBlFJzSIe7tLsIIca0ZHuE+J4+AhfXoYzR2bXviZwy3LXWJvAJ4HlgO+mrYrYqpb6slLrdXu2zwEeVUhuBR4EPa62HNt0IIcSYElrVCk5F4MLRdwPsU3GOZCWt9bOkT5RmLvtSxvQ24LLsliaEEIVjxUwi69vxn1eNo8hd6HJO29i4YFMIIfIssv4oOpGi6NKJhS7ljEi4CyHEENrShFa14Z5cjLuhuNDlnBEJdyGEGCLe3IvZGaVo6dg8agcJdyGEeJvQn1vT/cgsGLtXdEu4CyFEBrMrSmxnN4GLxk4/MsMZu5ULIUQOhFa3gVJjpmvfE5FwF0IIm5VIEV57FN/8ShwlnkKXc1Yk3IUQwhbZ0I6OmWP28sdMEu5CCIF9+ePKw7jqArgbS079hFFOwl0IIYDYrh7M9ijFVzSg1NjqR2Y4Eu5CCAGEXm3BUerGd97Yvfwxk4S7EOKcl2gJEt/bR9Fl9SjH+IjF8fG/EEKIsxB87TDK4xhTt9E7FQl3IcQ5zeyJEd3cQeDiCRjeEXWUOyZIuAshzmmhlYcBRdGlQ28NPbZJuAshzlmpUILwG0fwn1+Ns2xsf2lpKAl3IcQ5K7SyFW1aFF816dQrjzES7kKIc5IVSRJa1YpvQRWuGn+hy8k6CXchxDkp9OdWdDxF8dWTC11KTki4CyHOOVbcJPh6K945FbjrAoUuJyck3IUQ55zw6jZ01KTkmvF51A4S7kKIc4yVSBF87TCeGWW4J43N+6OOhIS7EOKcEl57BCuUHNdH7SDhLoQ4h2jTIvRKC+4pJXimlBa6nJyScBdCnDPC646S6k9QMk6vkMkk4S6EOCdo0yK44hDuycV4ZpQVupyck3AXQpwTwuuOkuqLU3Jd47i4GcepjCjclVI3KaV2KqWalVJfOME671FKbVNKbVVK/Tq7ZQohxJk7147aAU7Zv6VSygF8D7geaAHWKqWWa623ZawzA/g74DKtdY9SqiZXBQshxOkKr08ftZffNeOcOGqHkR25XwQ0a633aq0TwGPAsiHrfBT4nta6B0Br3Z7dMoUQ4sxo0yL40rl11A4jC/d64FDGfIu9LNNMYKZS6nWl1Gql1E3ZKlAIIc5GeN2RdFv7tZPPmaN2GFm4D7c39JB5JzADuAq4D/iJUuptH5FKqQeVUuuUUus6OjpOt1YhhDgtViJF//8ewt1UgmdmeaHLQWvNn3/z3/S1H835tkYS7i1AZmfHDUDrMOv8QWud1FrvA3aSDvtBtNYPa62XaK2XVFdXn2nNQggxIuFVrVjBBKU3No2Ko/Zdq1ey6rePsn/j+pxvayThvhaYoZSaopRyA/cCy4es8yRwNYBSqop0M83ebBYqhBCnw4qZ9L/cgndW+aj4NmoyHuOVX/6M6qapLLj2xpxv75ThrrU2gU8AzwPbgce11luVUl9WSt1ur/Y80KWU2gasAD6nte7KVdFCCHEqwVdb0j0/3tBU6FIAeOMPvyXY1cE1H34Qw3DkfHsjutW31vpZ4Nkhy76UMa2Bz9iDEEIUVCqUILTyML7zqnDXFxW6HHqOtLJ2+e+YfdmVNMyZn5dtyjdUhRDjTnDFIXTSouT6xkKXgtaal37+IxxOJ1d+4IG8bVfCXQgxrpi9MUKr2/BfUIuruvD3Rm1+YxX7N6zn0rvfT1F5Rd62K+EuhBhX+v90EICS6wrf82MyFmPF//sx1ZObWHTTO/O6bQl3IcS4kTgcIrL+KEVLJ+Is8xa6HFY98RjBrg6ufeBjGI7cn0TNJOEuhBgXtKXp/UMzRsBFybWFP2rvOLif9U//nnlXXUf97Ll5376EuxBiXIi81U7iYJDSm6Zg+EZ0IWDOaMvixYe/g8cf4Mr331+QGiTchRBjnhUz6XtuH+7JxfgXF75T2o1/+iNtu3dy1Qc/gq+4pCA1SLgLIca8/hcPYIWTlN0+DWUUtpuBUE83r/36F0yev5A5l19dsDok3IUQY1ryaJjQqlYCF07A3VBc6HJY8YuHSZlJrvvIxwran42EuxBizNJa0/vMPpTbQckNhf/C0t4317Jr9UouufNeyuuG9oyeXxLuQogxK7azh/iuHkqunYyjyF3QWpKxGH/66fepbJjMhbffWdBaQMJdCDFG6ZRF3zN7cVb5KFo6sdDl8Oqvf0Gws4PrPvpxHE5XocuRcBdCjE2hVW2YHVFK3zkV5SxslB3auokNzz/N4ptvp2H2vILWMkDCXQgx5qTCSfr/dBDPzHK8swp7h6VELMoff/BflE2o4x33fbCgtWSScBdCjDn9Lx5AJ0zKbp1S8Dssvfqrn9Pf2c5Nf/1pXJ7Cd3kwQMJdCDGmJI+ECa9po+iSibhqAwWt5cCmDWx88VkuuGVZQboYOBkJdyHEmKEtTc/vmzF8zoL3+hiPRHj+R/9FeV09l937gYLWMhwJdyHEmBFe3UbiQD+l75yK4S/sFSmv/PInhLq6uOljn8bl9hS0luFIuAshxgSzJ0bfH/fjmVmOf1Fh+4/ZvXYVm196gSW33cHEmbMLWsuJSLgLIUY9rdPNMaApv2N6QU+ihrq7eOFH36FmyjQuu+f9BavjVCTchRCjXmRDB/FdPZTe2ISzvHBXpGjL4rnvfQMzEefWT31uVHxZ6UQk3IUQo1oqlKDvqT24G0sIFPibqOueeZKDWzZy9Yc+SsXEhoLWcioS7kKIUa33qb1Y8RTld80oaHe+R/c2s/LRR5h+4VIWXHNjweoYKQl3IcSoFdncSXRjByXXTMZV4y9YHclYjGe+/XX8paXc8JefLPgXp0ZCwl0IMSqlggl6f78bV0MRxVcVrglEa83//uyH9Bxp5eaPf7Zgd1Y6XRLuQohRR2tNz+92YyUsKt4zC+UoXFRtfukFtr7yJy65814mzz+vYHWcLgl3IcSoE157hNiObkpvbipoc8zRfXt46ec/pPG8RSx9970Fq+NMjCjclVI3KaV2KqWalVJfOMl671ZKaaXUkuyVKIQ4l5idUfqe3otnellB+2mPhUM89c2v4isp5ZZP/i2G4ShYLWfilOGulHIA3wNuBuYC9yml3tZDjlKqGPgUsCbbRQohzg06ZdH12A4wDMrfPbNgV8dorfnj979JsLOD2x76PP6S0oLUcTZGcuR+EdCstd6rtU4AjwHLhlnvX4CvAbEs1ieEOIf0v3iQZEuI8rum4ywrXH8ta5f/jj3r1nDl++9n4sw5BavjbIwk3OuBQxnzLfayY5RSi4BJWuuns1ibEOIcEtvTS/CVQwQunIB/QXXB6ji0bTMrH32EmUsvZ9HNtxesjrM1knAf7u8ifexBpQzgm8BnT/lCSj2olFqnlFrX0dEx8iqFEONaKpyk53924qz0UXrb1ILVEerp5ulv/TtldRO5cYxcz34iIwn3FmBSxnwD0JoxXwzMB15WSu0HLgGWD3dSVWv9sNZ6idZ6SXV14T6ZhRCjh7Y0PY/vJBVOUnHfbAx3YU5cmokET33jqySiUW7/9Bdw+wp3lU42jCTc1wIzlFJTlFJu4F5g+cCDWus+rXWV1rpJa90ErAZu11qvy0nFQohxJfhqC7GdPZS9cyru+qKC1KC15oWHv0Prru3c/PFPUzW5qSB1ZNMpw11rbQKfAJ4HtgOPa623KqW+rJQauw1SQoiCi+/ro/+F/fgWVBG4pK5gdbzx5G/Y/toKLnvP+5l5yTsKVkc2OUeyktb6WeDZIcu+dIJ1rzr7soQQ410qlKDr0R04y73pTsEK1L69a83rrHzsEea84youvvOegtSQC/INVSFE3mlL0/0/O7EiSSreNwfDO6LjzKw7ureZ5777DepmzOKGv/xUXj5gutvCaK1PveJZknAXQuRdcMUh4rt7KbttGu6JhWln7+/s4MmvfRlfSQnL/vYfcLrdOd/moR3dPP6VtWz830OnXvksSbgLIfIquq2L/j8dwH9+NYGLJhSkhlgoxBNf/ScSsRh3/J8vESgrz/k2j+zt49kfbKa0xsfsPJxfkHAXQuRN8kiY7sd24qovKlg7ezIR58mv/wu9R1pZ9rf/QHXjlJxvs+NQkKe+s5FAiZvb/+Z8vEW5vz2fhLsQIi9S4SSdj2xDeRxUfWAuypX/69ktK8Wz3/4PDu/cxs2f+GxeuvDtbg3z1Lc34PY6uP2h8wmU5qdbBQl3IUTO6ZRF939vJ9Ufp+qDc3HkKeAG1aA1L/38YZrXruLqD36EWUsvz/k2e46EefJbb6EMxbKHFlFS6cv5NgdIuAshcq53+R7ie/uouGsm7knFBalhzRP/w8YXnuHC2+9i8S3D9X2YXT1Hwjz5jbcAeNenF1FWm99vvEq4CyFyKrSqlfCaIxRf2YB/UU1Banjzuad4/fFfMffyq7n8vg/lfHu9RyM8+c230FrzrocWUT4hkPNtDiXhLoTImVhzL71P7cE7p4KSG5sKUsOWFS+y4hc/YvqFS7nxrx9CGbmNvb6OdLBbKc2yhxZRMTH/wQ4j/IaqEEKcrkRriK5fbcdZ7afi3lkFufHGzlUreeFH36HxvEXc+jf/B8OR25O4PUfC/OFbG0glLd71mUVUFqivHJBwF0LkQLI9QudPt2B4HFR9eB6GJ/9Rs/ettTz7na8zcdZsln32izhdub38sONQkKe+vQGUKniwg4S7ECLLkh0ROn68GRRUfWQ+znJv3ms4sHkDT/3nV6lunMIdn/8nXN7c1nBkbx9Pf3cjLo+DZQ/l/+TpcCTchRBZY3ZG08FuaaofXICrOv8ht3/Dev7wH1+hbEIdd/7d/8Xjz22b9+GdPTz9/U34S9wse+j8vF7ueDIS7kKIrDC7onT8eBOkLKo/eh6u2vyfSNz71lqW/8dXqGiYzLu/+C85v7H1vk2dPP/jLZRU+ViWxy8ojYSEuxDirJndMTp+vBmdtKj6yAJcBbj0b8/6NSz/z69SNbmRd//D/4evKLfX02997TCv/Hon1ZOLeecnF+Iryn3HY6dDwl0IcVbMnhgdP96EFU9R/ZEFBenlcfcbf+bpb32NmqYp3PX3/4K3KHc1aK1Z+/Q+1j6zn8nzKrnxo/NwF6jL4pMZfRUJIcYMszdOx483Y0XNdLAX4AqRLSte5IWHv8OEaTO46++/nNM2ditl8cqju9i2spXZl9Zx1ftm4XCMzq8LSbgLIc5Iqi+ePmIPJ9PB3pDfbgW01rzx5G9Y+dgjNJ63iNs/+/e4vbk7mZmImbz4063s39zFBTc3cvHtUwt296iRkHAXQpy2VH8ifcQeSlL1wPy89xejLYuXH/kJbz63nNmXXclNH3sIhzN317H3d0V59vub6W4Lc+V9M5l/ZUPOtpUtEu5CiNOSCibo+PEmUv0Jqh6Yj2dySX63byb54/e/xY7XX2HxLcu46gMP5LRLgfRNNjaRMjXv/MR5TJ5bmbNtZZOEuxBixJLtETp/sRUrlKDqL+bjacxvsEdDQZ76xlc5tHUTl7/3w1x4+105bRrZueYIL/1yO0XlXu74+HkF6QDsTEm4CyFGJLanl65fbkc5FVUfWZD3I/bu1sM8+bX/S39HOzd/4rPMvfzqnG3LSlmsWb6XN58/SP3MMm56cEFe7p6UTRLuQohTCq8/Ss8Tu3FWeqn68HycFfntUuDglk089Y1/RRkGd//jv1I/e27OthUNJnjhp1tp2dHDvMsncvk9M3E4R+cVMScj4S6EOCFtafpfPEBwxSE808uofN8cDF9+Y2Pzihf404+/R9mEidzx+X+irDZ3N9VuP9DPcz/cTDSY5JoPzmbOpRNztq1ck3AXQgxLJy26f7uL6MYO/EtqKb9jOiqP13SnzCQr/t9P2PjCMzSet4h3PvR5vIHcXUe/7fVWXn10F/4SN3d+bjE1eT6fkG0S7kKIt0mFEnT9cjuJA/2U3NRE8ZUNeb2mO9TdxfJvfpW2XTu44J13cMV7P5yzvtgTMZNXHt3JrjVHaZhdzg0fmTfquhI4ExLuQohBBq6ISfUnqHjfbPwLqvO6/ZZtW3jqW/9GMhbjnQ99Pqc3sm4/0M8LP9lKf2eUi26bwgU3N2EU4KYiuSDhLoQ4JvOKmOoH83tFjNaa9c88yav//XPKaidw9z9+hapJjTnb1sb/PcSq3+/BX+LmXZ9ZzMQZZTnZVqGMKNyVUjcB/wU4gJ9orf9tyOOfAT4CmEAHcL/W+kCWaxVC5Ii2NMFXW+h/YT/OKl/er4iJ9Pfxx+9/k31vrWP6hZdw08c+nbM+YsK9cV765XYObu1mysIqrvnAnDF3meNInDLclVIO4HvA9UALsFYptVxrvS1jtbeAJVrriFLqr4GvAffkomAhRHalwkl6Ht9JbGcPvvOqKL9zBkYeezk8uGUTz373P4iFglxz/19x/g235qx9f/fao7zy6E5SSYsr7p3J/CvrR3X/MGdjJD/Bi4BmrfVeAKXUY8Ay4Fi4a61XZKy/Gnh/NosUQuRG/EA/3b/eQSqUoGzZNAKX1OUt7FKmyarfPsqaJx+nvK6eO7/wz9Q0Tc3JtmKhJK88upPm9e3UTinhug/PHRW3wsulkYR7PXAoY74FuPgk6z8APHc2RQkhcktrTWhlK33P7cNR5qHmrxfmtVfHzoP7+eMPvsXRvc3Mv/p6rvnwX+bkPqdaa/Zu6ODVR3cRCye55F1TWXT9ZIxR2k1vNo0k3If7GNfDrqjU+4ElwJUnePxB4EGAyZMnj7BEIUQ2pYIJep7YTWx7N965lVTcPTNvX0yyUinWPvUEq37z37j9AW77zN8x8+LLcrKtUE+MVx/bxb6NnVRNKuK2Ty2kKs/dEhfSSH6iLcCkjPkGoHXoSkqp64AvAldqrePDvZDW+mHgYYAlS5YM+wEhhMgNrTWR9e30Pr0XbaYofedUii6bmLdmmK7Dh3j++9+irXknMy6+lOs+8vGc3OPUsjRbXmlh9ZN70Zbm0juns/DahnPiaD3TSMJ9LTBDKTUFOAzcC7w3cwWl1CLgR8BNWuv2rFcphDgrqb443b/dRXx3L+6mEsrvmoGrOj9tzmYyyRtP/oY3nnwcl9fHrZ/6HLMuvSInHyqdLSFW/GoH7fv7mTy3givfO4uSqtzdwGM0O2W4a61NpdQngOdJXwr5M631VqXUl4F1WuvlwNeBIuA39g/soNb69hzWLYQYocjmDnqeaAbTSp80vbgOlacv6hzcsok//eR79LQdZtalV3D1hz5KoKw869uJR5KsfWY/m1e04Ak4uf7+ucy4sHZUXgljxWIojyfntY2ooU1r/Szw7JBlX8qYvi7LdQkhzlIqnKTvmb1E3mzH1VBExb2zceXpKDbc28Orv/oZ215bQWntBO76+y/TtHBx1rdjWZrtr7eyZvleoqEkcy+tY+kd00fldes6maT3d7+j8/s/oPaLX6Tkxhtyuj35hqoQ44y2NJH1R+l7bh9WzKT4mkmUXDs5L51+JeMx1j/9JG8s/x2pZJKL77iHi+98Dy63J+vbat3dw2uP76bzUIi66aXc9smZVE8efSdMdTJJ3/Kn6PzhD0keOoRv8WJcE2pzvl0JdyHGkURLkN6n9pI40I+7sYTyO6bjysPdg7Rlse21Fax87BFC3V1Mv3Apl7/3w1RMrM/6tjpbgqz+w14ObO6iqNzDDR+Zx/QLakZdE4w2zeOhfvAg3nnzqP3hDyi68sq81CrhLsQ4YPbF6X9+P5E32zECLsrvmoH/gtq8tK0f3LKRl3/5Uzr276V26gxu/eTnaJg7P+vb6euIsGb5PnavO4rH52TpHdNYcHUDLndueos8U1Y8Tt8f/kDXT39K8sBBPHPn0PD971N09VV5/QCScBdiDLPiKYKvthB6tQWtNcVXNVB81aS8dB9wZM9uXn/8V+zfsJ7iqmpu+eTfMvvSK7J+s+pwX5x1z+xn28pWDIdi8Y2NLLp+Mt7A6GpXTwWD9Dz6GN2/fIRUR2f6SP1736XommsK8leFhLsQY5CVSBFe1Ubw1RascBLfwmpKb2zKS2dfHQf38+fHf0Xz2tV4i4q5/L0fZtHNt2W9XT3YHWPDiwfZtrIVK6WZe/lEltzSRKA0++33ZyN59CjdjzxC72P/gxUOE7jsMiq/9jX8l1xS0KYiCXchxhArZhJa3UbotRassIlnZjkl103OS9e87fv38saTv2Hn6pW4vT4uvft9LL5lGR5/dq+X722P8ObzB9i5+ghomHlxLUtuaaI0T9flj1R8zx66fvYz+pY/BakUJTffTOUD9+Odm7v7u54OCXchxoBUKEHo9VZCq1rRsVTeQl1rzYHNG1j31BMc2PQWLq+Pi991Nxe88w58Rdm9MqXrcIj1z+2neX07hsNg3jsmcv4NkympHD1fQtKmSXDFCnp+/Wsiq1ajvF7K776bivv/AndDQ6HLG0TCXYhRzOyLE3qlhfDaI2jTwje/iuIrG3LeyVfKNNm16jXWPv17OvbvJVBWzjvu+xALr7sZb1H27mNqWZr9mzrZ/HILLTt6cHkcnH/9ZBZeO2lUNb+YXV30/uY39PzP45htbTjr6qj+9Kcpu/vdOCsqCl3esCTchRiFUv1xgi+3EFrTBhr8i2sovrIh510GJKIRNr/0IuuffZJgZwcV9ZO44a8+xZx3XI3Tlb0TmLFwku2vt7H5lRaCXTGKyj1cvGwq86+oHzUnSrVlEVmzht7fPUH/889DMkng0qVM+OLfU3TVVSjn6I6pX71aAAAeeklEQVTP0V2dEOeY5NEwoddbCb/ZDpYmsKSW4qsn4SzP7YnS9v172fSnP7J95QoS0SgNc+Zz7f1/zdRFS7J29YvWmo6DQba+1squNUcwkxYTZ5Rx6Z3TmXp+1ajp2CvRcpi+J5+k74knSLa2YhQXU37PPZS/9z48U3PT33wuSLgLUWBWzCS6qZPwuiMkDgbBqQgsqqX4qgacOWxvTsZj7Fy1kk0vPkdb806cLjczl76D82+4lboZs7K2nVgoyc43jrD99Ta6DodwuAxmXlTLeVc3jJoueK1wmOBLL9H7xBNEVq0GpQgsXUr1Zz5D8XXXYuSgr/lck3AXogC01iT29xNed5Topg500sJZ46P0lin4F9fgKHLnbLttu3ewfeXLbH/tZeKRMBUTG7j6Qx9lzhXXZO0kqWVpWnZ0s/31NvZu7MAyNTWNxVx530xmXFiLx1/4phcrFiP08iv0P/ccoVdeQcdiuBoaqPrUJyl717twTZxY6BLPioS7EHmitSbZFia6uZPopg7MrhjK7cB/fg3+C2txTyrO2XXRnYcOsOP1V9jx+iv0tR/F6XIz/aKlLLzuZurnzMvKdi1L07a7l+Y329nzVgfR/gSegJP5V9Qz59K6UXGUbiUShFeupP/Z5wi99BJWJIKjspKyO++k5Jab8S1enPUvYRWKhLsQOWZ2RYm81U5kQwdmZxQUeKaWUnz1ZHznVWHk6Ovz/R3t7Pjzq+xY+TIdB/ejDIPGBeez9N3vZfqFS7Nyfbpladqae2lefzzQnW6DxvlVTL+ghinnVeFwFTYsU6EQ4ZUrCa1YQfClFVjBII7SUkpuvZWSW2/Bf+GFKMfo6sIgGyTchciBVDhJdEsnkTfbSRzoTwf6lFKKLq/HN68yJ80uWmu6Dh1gz/o3aF63miPNuwComzmba/7iL5l5yTuy0pe6ZWmO7OmleV060CP9CZwug8YFlUy/oJbG+ZW4PIUNy0RLC6GXVhBc8RKRdeshmcRRWkrxtddScsvNBJYuRWXx6p+R0lpzMHiQUncpZd6ynG5Lwl2ILNCWJtkaIrajm9iuHhKHgqDBWeOj5KYm/OfX4CzL/nXbKdPk8I5t7Fm/hj3rVtPXfhSACdNm8I57P8isS6+grHbCWW8nFk5ycFsXB7Z0cXBrN7FQ8ligT1tcQ9OCqoIGupVIEH1rQ/oI/eUVxHc3A+CeOpWKD36A4quvxnf++QW5fLEt1Mbao2tZe2Qta9rW0BZu4x8v+UfeM+s9Od2uhLsQZ0BbGrM9QnxfH/G96cEKJ0GBq6GY4msm45tbiWtiIKvt6Nqy6Dx0gINbNnJwy0YObdtCMhbF4XIxef5CLlp2N1MXX0hRReVZbceyNJ2Hghzc1s3BLV0c2duH1uAtctE4r5LGBZU0zq/EnYcOyoajLYv4jh2EV60ivGo1kXXr0LEYOJ34L7iAmi/cRfHVV+NubMx7bZlhvvbIWg6HDgNQ6illSe0SHpj/AJfXX57zOiTchRgBbaVPhsb39hHf10difx9WxATAUerBO7Mcz8xyvDPKst7k0td+hAObN3Bw80YObt1EtL8PgPK6icy9/CoaFyyiceEi3N4zv2xSa013W5jDO3to2dFD6+5e4vb/r6axmAtuaaJxfiU1jSUYebpF36D6LIvEnj1E1r9JeM1qIqtWk+rtBcA9fRpl7343gUuX4r/wQhzF+Ttxa1omu3t2s7FjIxs7NvJW+1tvC/MPzP0AS2qXMKN8BobK3/kHCXchhqFTFonDIRIDR+b7+9HxFACOCi/eOZV4ppbimVKKozx798PUWtPT1srhnVtp3bmdQ1s3HWtqCZRX0LRwMZPnL2Ty/PMoqao58+1Ymp6jEdqae9OBvrOHaDAJQEmVl6mLqmmYVU79rPKCdANgxWLENm8m8uZbRN98k8iGDVh96Q81Z00NRVdemQ7zS5biqj3z/XC6emO9bOrcxIb2DWzs2Mjmzs1EzSgAVb4qFlYvLFiYDyXhLs55WmtSfXGSbWGSrWHi+/tIHOhHJywAnNU+/Aur8UwtxT2lFGcWwy7S18vRfXto37eHtuZdtO7afuzI3FtUTMOceVxw67uYPP98KuobzvhDJBZO0nEoyNG9fbTt6efovr5jR+b+UjeT5lRQP6uchlnllOTpPqsDtGWR2H+A2NYtxLZsIbphI9Ft2yCZ/rBxT51K8fXX4V+0GN/iRbibmvLSlW5ntJNtXdvY3rWd7d3b2d61ndZwKwAO5WBWxSzeNf1dLKxeyPk15zMxMHFU3Q1Kwl2cU6y4SfJIhOSRcDrMj6QHHUsdW8dZ68e/uPb4kXlxdppZQt1dHN23h6N7m2nfv4ej+/YQ6uo89njZhDqmLlrCxFlzqZ81l4qJ9Wd0zXU8atJxMEj7gX57HKS/I3rs8YqJAaYtqmbCtFImTC2lrNaft1DSlkXy4EGiW7YS27qV2JYtxLZtwwqHAVAeD965c6n80AfxLV6Mb9EinOVnf4XPyaSsFIdDh9ndu5sd3TvSYd61nfZo+7F1GksaWVi9kHtm38OCqgXMq5yH33WGl5LG+kEZ4MleB2zDkXAX45K2NGZ37HiA2+NUd+zYOsrjwDUhgP/8GlwTArjqArhq/Wd9F6NkIk5P62G6Dh+i69BB2venj8zDvT32hhUVdfU0zJ5H7dTp1E6ZRnXTVLyB0/tlt1IWfR1Rug6H6W4L090aorMlRF/78SAvrvBS01jM3MvqqJlcQk1TcV6+Haq1xmxvJ75rN/HdGcOePehouj7lduOZPZvSZbfjnTcf7/x5eKZNy9kVLVpr2sJtNPc209zbzJ7ePezu2c2+vn3EUun3haEMppRM4eK6i5lTOYc5FXOYXTGbIvcZBHG4Ezp2QufO9LhjJ3Tugv7DcPt3YPEHs/w/HEzCXYxZ2tKkgglSXVHM7tjxoTOKeTSCTqabVVDgrPLhri/CdUFtOsQnBM66rTweidB9+BBdLQfpOnwoPX34ULqNXOv0pg2DyvpJNC1cTM2UaemhaeppnfyMhZL0tkfobY/Q1x6ltz1CT1uEnqNhLFMf+z+WVPmonBhg9iV11DQWU91YjC9H3RgM0KZJsrWVxP79JPYfIL5vrx3kzcfayAEc1VV4Z8yg/D1345kxA++8eXimT8/Jteb9iX4O9R/iQP8BDgQPHJve07eHcDJ8bL0aXw3TyqZx96y7mV42nWll05hRNuP0jsiTMejZB91700NXM3TsSgd6pOv4eq4AVM+EpsvT4/oLsvg/Hp7S9psw35YsWaLXrVtXkG2LscOKm6R645i9cVK9cVJ2eCc7o6S6Y8cDHECBo8yDs9KHq9Z/LMRdtX6U6/SvwU6ZJsGuTvqOHqGv40h63H6Uvo6j9B09QjTYf2xdh9NJ+cQGKuonUVk/icqG9Lisrn5EXeXGoyZ9doD3Ho3a0+nxQNs4gFJQXOmlrDZA5cQAFfUBKuoClNcFcnajaCsWI9naSvJwK8nDLSQOHLTDfD+JlpZjbeMARkkJnhkz8MyYbo/TQzabVrTWdMW6aAu1cSh46HiA2+OeeM+g9Wv9tTSWNDKtbBrTy6YfC/JST+mpN5YyIXQE+lrs4RD07IeuPdC9L30UTkaG+iqgehZUzUyPq2dB1SwoqYcsdWuglFqvtV5yqvXkyF0UjLY0ViiB2WMHd28cszeWMR1HR83BT3IonBVenJU+vDPKcValp50VXhxlHtQIu41NmUnCvT2EursJ93YT6ukm3NNDqKeL/o52+tqPEuzsQOvjHx6Gw0FJVQ2ltROYcdGllNZOsMO8gdKaCRgn+Aq7ZWkifXGC3XFCPTGC3TFCPXFC9jjYHSMWSg56TlGFh7IaP9OX1FJW46O0xk9ZjY+SKh8OZ3avwEiFwphtrekj8MOHSR4+fDzMW1tJdXYOWl95PLgnT8YzfTrF112Hu6kJd1Mj7qYmHBUVZ91+n0glOBI+Qlu4jdZQ6/HpsD0daiNhJQY9Z0JgAo3FjVzbeC2NxY1MKplEY3EjDcUNeJ0n6NFRa4j1Qt/h48Hd15IO7IEw728FnRr8PH8lVEyFpnekx5XToGJKetp38g+xSMKkIxin1OeizJ/bv6rkyF1klTYtrHCSVDiJFUqmp0NJrHDCHtvLgglS/QlIDX7/Ka8DZ1k6qAcGZ7kHx8CyYjfqBNdZa8siFgkT7e8j2t9PuK+HcM/x4D4e4t2DjrqPbVsZ+MvKKKmuobS6lrLaCZTU1FJWM4HSmgkUVVQOCnArZRENJon0J+whnh73pefDvXGCPTHCvQm0Nfj/6fY6KKrwUlTupajCQ2mVj7IaP6U1PkqrfTjP8ihca40VDmO2d2C2t2N2nHhsRSKD94PLhXNiHe76elz19bgmTkyP7Wlnbe0ZneiNmlE6I510RDvoiHbQGe2kM9pJRyQ9PbCsO9b9tudW+6qpK6qjLlDHxMBEJgQmUBeoY1LxpLcHeMqESCcEj0DoqD1uTx+BZ06H2sGMDd6Q4YLSeihpgNLMYZI9rgfP4OvoLUvTE0nQHozTEYzTHozTHowdm+7IGELx9MHKV+6Yz/suPrMvWI30yF3CXZyUTmmsyEBIJ46H9aDgPj6tY+bwL2QojIALR5ELo8iFI+DKCHAvTnt64GSm1hozEScaDBIN9qcDO2McGZgP9hPt7z82rS3r7Zt2OPCXlVNUVk6gvJKi8nICZRUEyisoKk+P/SVlGE4/iahFLJQkFk4P0VCSeDhpLzOPLw8miIaSg/4iH+D2OvCXevCXuCmuSId3Ubk3PV3uoajCi8c38j+atdboSASzp5dUby+pnh5SvT2kenrT495ezJ6Befvxnh50IvG211JeL86aGpw11Tirq3HV1OCsrsY5oQ5X/URcE+txVledMryTVpK+eB+9sV564j30xnvTgz3fF++jJ2aP4z30xHoIJUNvex2nclLpq6TKV0W1r5oqfxU1vhrqitIhXheYQK3DjzsehGg3RLrTbdnHxvYQ7UmPw53pYNdvfx/gK4ei2vRQPAGKaqBoQjqwSyeRLJpIv6OcvliKvmiSvmiS/piZHtvzfZGB5elxVyhBZyiOab39jRBwO6gp8VJd5KG6xEN1kYcae3xhUwVNVYERvwcG/QyzGe5KqZuA/wIcwE+01v825HEP8AhwAdAF3KO13n+y15RwLwxtZYb14GC2wkmsUGLQUbcVOUFYKzACruOBHXDhKHKj/E6U34HlsjCNJKZKkiRO3IyQiEWJR8IkolHi0QiJSJh4JEIiGjk2zpy2UqnhN60MvMXF+EtK8RYV4w2U4PIV4fYW4XQHcLj8OFwBDEcAw1GMVh7MmCYRM0nETJKxlD1tj6MpEkObfzIYhsJT5MIbcOENOPEGXPhL3Omh1IO/2I2/ND3vK3EPavvWWqOjUVLBEFYoiBUKHZtOBYNYwVB6WWhgOmg/HsIKBknZ4+GC2t4ZOEpLcZSV4SgvPz4uL8NZUZEO7oEAr6nBKCrC0haxVIxIMkLEjBBKhAgmgwQTwfR0IkgwmZ7uT/QTSoQIJe3liSB98T6CyeAJ95ff6afMU0aZp5QydzFlzgDlTj9VDh/VhocqHFRpRXVKU2YmMOJBiPenh9iQcbQXrOTwG1IO8Fekm0l8FaR8FSQ95UQ9VYRdlQRdlfQalXSrMjoooy9hHAvlY+GdEdzhxPDvtwEep0Gpz3VsKPG5qAy4qS72UFPsobrYeyy8q4s9BDy5afXOWrgrpRzALuB6oAVYC9yntd6Wsc7HgPO01n+llLoXuENrfc/JXlfC/cxprdEJCx0zsWImViyFFTUHzeuoPR01sUJ2M4k9DHe0CYBHoT1guXQ6nB1mOpxVnARx4laUuBkiaoaJJUIkEzGSsTjJeJRkLEYyHicZj50wlDMpw4HL48Ph8uF0eXG4vBhOD4bDi+HwoAwPynBjGH4wfCjlA3xY2oOZdGPGUySiKVLmMEdob9sYuDwO3F4nbq8Dt9eBy23gcilcbnA5we0Cj9vC40wPbsPEZSRx6ziOZAwdi6CjMaxoFB2LYmVOR6JYsRhW9Pg6ViyKtpczgv1hFBWhAgFUUQCK/OiADx3wYQW8pPxezBIfZpGPZLGXWJGbeJGbWMBFyANRHSdiRoiaUaJmlEjy+PTQ+YgZIZ6Kn3qfAUUOD0WGh2LDTbHhpBgnRSjKMCizoNzSlKZSlJtJysw4ZYkYZfEonmQYEmFIneADaQjTGSDpKibhLCLuKCLuCBA1AkRUgJAK0GeU0kcxPbqYLl1Eh1VMV8pPl+klYmqiiRT9MZPECN4LPpfjbQFdOmhwUup3UeJ1vW097xmclM+FbJ5QvQho1lrvtV/4MWAZsC1jnWXAP9vTvwW+q5RSulBtPnmitYaUTrenmhba0uiUvSxlQUpjJVNYppkeJ1JYCRMrbmIlUuiEiZWw0PZynbCwkql0cJsakhbatCCZfl1lalQKlKlQnPyklYWFSZKkTpCwYsStCFEzTNwMEUmGiJkh4qkw8VSEWCpCwoqiT5j6oAwnynChlBulXCjlBOUm/RYqRlOO1u50e63LBcqNUp70OhnTx8d28wuQtCAZB+JgKAuHsjLGKQwsHKQw6MepTbwkceoEDp3AoeI4UgmcVgxHMorTjGIkIjgTYZzxCEYihCMWQSUT6as6EknUMH9CZzLt4URSToOU24HpdmC6DUy3QdJlkHQrTLdBoggSLkXS6SXu9hL1KCJuCLstQm5NyGMRdFn0uVP0uy363SamigEx0n/4nsTAaoPPceLQ4EPh0wqvBp89+LVFlWXht1L4rRS+lEmRZeKzND6t8WuLIktTbFnHhiLLImBphosyEydx5SWqvMSUlyheonjowUuLLiasPUS0m6DlJWS5CWkPETxEtIcgfoLab499BPETwofF8E1ALofC63LgczmOj90OfC4Dn9/B1IxlxV4nJV4XxV4nRR4nxfZ0sddJscde7nXiGiX3ac2HkYR7PXAoY74FuPhE62itTaVUH1DJ296CZ+/ZT36dKe7Zx+YHh5wifaJeZSwZ/l9U5nOHRqU6vq4a7tUUhjIwVPY+yVOWiamTpLRJSicz5pOY2iRlmfZ0goQVJ3lsSJC0kiQsk6Q2Ma0USZ0ihQIcdog6QLlQygXKCbhRqhxUDbhcKJcLFw4UBgYGyjIwUBjawKEVytI40CjLxNApDCuZHlLpscMyMayEPYQwLNNePng9wzJxWAn78cSx9Y4P5qCfg2mA6QDLgJSClGFPDxlMA+IGJB1gOhVJBySdkPSBWWRPO9KvlZ42hlmWnk44B5alXyfhgrgT4m5IOTQpp0YZ4NQmTkycOh2CDq1xanCicQwZOzU4tcZtD2VArda4tMYT07ijxx/zaI1Lg0drDEvh0AqH/XNwaCP9M7EMlDZQ2oFLO3BYBuDExIGpHemxPcRwE9cu4rhIkB7HcdGv3cem47iOrRPHTRInlsOD5fCQcngxHT5STj8phw/D5cHtNHA5DDxOA7cjPe12GoOXOw1cDkWJw0GVPT2w3O00jgX2QGj73I6MZQZel+OcCuJcGEm4D3eIOPTQZyTroJR6EHgQYPLkySPY9Nul3Ba9iW7IPM60/0DQdiGZR6CDj0btOT3co0OPW/Wx0XCvYWkLCwtLW2idQqOxdOr4MrT9mMYiPU5pC4sUKZ3C0ik7yNPT6e1pUBqlQQ1s0/5wUfa20xeKaHveQimNQmMoCzcWHodGYWGpFKgUkB5rZaGVCVjp5Splr2OhlQVKo5VC2xvTCkylSNrz1qDHFNoA7GV6YJk9zbFlRvq5hkI7FMoAbRhYSqEMF5ZyoYwA2l6OocBIv3a6w6X0h6xhj9VwY22P1eB1wTg251AO+8NLYeDAUA4UDvujzIFSjhOOFU577ABloJWBxh6UQiuHvUzZy9KDUgYoA2U40tOGA2U40Q4nGC6U4QSHC204MRxOcLjBcKIc6eWmcuBwGGhDYdn7P2UoDKVwGGCo9LS2lxkKHIbCO7COUhj2eg5D2bt2YFphGPY6Kn1OwWUcD15HAXp9FNk3knBvASZlzDcArSdYp0WlDxVLgbddz6S1fhh4GNJt7mdS8G3/+fkzeZoQQpxTRvJ3z1pghlJqilLKDdwLLB+yznLgQ/b0u4GXxnt7uxBCjGanPHK329A/ATxP+lLIn2mttyqlvgys01ovB34K/FIp1Uz6iP3eXBYthBDi5EZ0IabW+lng2SHLvpQxHQPuzm5pQgghzpScjhZCiHFIwl0IIcYhCXchhBiHJNyFEGIcknAXQohxqGBd/iqlOoADZ/j0KnLQtUGWjNbapK7TI3WdvtFa23irq1FrXX2qlQoW7mdDKbVuJL2iFcJorU3qOj1S1+kbrbWdq3VJs4wQQoxDEu5CCDEOjdVwf7jQBZzEaK1N6jo9UtfpG621nZN1jck2dyGEECc3Vo/chRBCnMSoCXel1M+UUu1KqS0ZyxYqpVYppTYrpZ5SSpVkPPZ3SqlmpdROpdSNGctvspc1K6W+kM+6lFLXK6XW28vXK6WuyXjOy3ZdG+yhJo91NSmlohnb/mHGcy6w129WSn1bKXVWd2o4zbrel1HTBqWUpZQ6334s2/trklJqhVJqu1Jqq1Lqb+zlFUqpF5VSu+1xub1c2fujWSm1SSm1OOO1PmSvv1sp9aETbTOHtb3PrmmTUurPSqmFGa+1397PG5RSZ3WT4jOo6yqlVF/Gz+xLGa+Vtd/LM6jrcxk1bVFKpZRSFfZj+dhfd9vzllJqyZDn5C7HtNajYgCuABYDWzKWrQWutKfvB/7Fnp4LbAQ8wBRgD+nuiB329FTAba8zN491LQIm2tPzgcMZz3kZWFKg/dWUud6Q13kDWEr6Zk/PATfnq64hz1sA7M3h/qoDFtvTxaRv+j4X+BrwBXv5F4B/t6dvsfeHAi4B1tjLK4C99rjcni7Pc22XDmwTuHmgNnt+P1BVoH12FfD0MK+T1d/L061ryHNvI32/iXzurznArKHvaXKcY1n5xcnWwJAQAvo5fl5gErDNnv474O8y1nuedEAtBZ7PWD5ovVzXNeQ5ivTdjj32/KAfbJ7316D1hrwZd2TM3wf8qED761+Br2TMZ31/DdneH4DrgZ1AXcb+2GlP/wi4L2P9nfbjg/bR0PXyUduQdcsZfBCxnyyF1Rnss6sYPtxz8nt5hvvr18BH87m/MuYHvaeH7geynGOjplnmBLYAt9vTd3P8dn/D3bS7/iTL81VXpruAt7TW8YxlP7f//PvHs23+OIO6piil3lJKvaKUutxeVk96Hw0o5P66B3h0yLKc7C+lVBPpv7LWALVa6zYAezzQ/FOQ99gIa8v0AOm/MAZo4AWVbhZ8sAB1LVVKbVRKPaeUmmcvy9k+O539pZTyAzcBv8tYnI/9dSI5fY+N9nC/H/i4Umo96T9zEvbyE92Qe0Q36s5hXQDYb+p/B/4yY/H7tNYLgMvt4QN5rKsNmKy1XgR8Bvi1Srd7j5b9dTEQ0VpvyVick/2llCoi/cv9kNa6/2SrDrMsp++x06htYP2rSYd75o2FL9NaLybdXPNxpdQVeazrTdJfjV8IfAd4cuAlhln3rPfZ6e4v0k0yr2utM+/vXMj9ldP32KgOd631Dq31DVrrC0gf1e2xHzrRTbtHcjPvXNaFUqoB+D3wQa31noznHLbHQdJ/Gl6Ur7q01nGtdZc9vd5ePpP0/mrIeIm87y/bvQw5as/F/lJKuUj/0v231voJe/FRpVSd/Xgd0G4vz+t77DRrQyl1HvATYNnAzxZAa91qj9tJvw/Par+dTl1a636tdciefhZwKaWqyME+O939ZRvufZaP/XUiuX2P5aKt6SzaqJoY3FZbY48N4BHgfnt+HoNPROwlfRLCaU9P4fiJiHl5rKvM3uZdQ57vxG7XA1zAb4G/ymNd1YDDnp4KHAYq7Pm1pE8YDpxQvSVfdWUsawGm5nJ/2f+/R4BvDVn+dQafhPuaPX0rg0+ovmEvrwD2kW7rLrenK/Jc22SgGbh0yPoBoDhj+s/ATXmsawLHz61cBBy0XyOrv5enW5c9X0r6/s6BfO+vjMdfZnCbe05z7Kx+ibM5kP5EbQOS9i/7A8DfkD7jvAv4t4E3jr3+F0kfAe4k4woP0lc57LIf+2I+6wL+AQgDGzKGGvuNsx7YBGwF/gs7bPNU1132djeS/tP5tozXWUK6TXwP8N3MfZynn+NVwOohr5GL/fUO0n/absr42dwCVAL/C+y2xwMfegr4nr1fNg/5pbyfdLg2A3+RhffY6db2E6AnY9119vKp9s94o73fzur9fwZ1fSLjfbaajA8fsvh7ebp12c/5MPDYkNfJ1/66w/5diANHGXyyNGc5Jt9QFUL8/+3YMQ0AAAzDMP6sx6HXFNkgcoSg188dgI24AwSJO0CQuAMEiTtAkLgDBIk7QJC4AwQdiByx8riK/5QAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "\n", "print(len(Time))\n", "print(len(SL_wTd_nos_base_R0_RCP70_cdf[0:-1,500]))\n", "plt.plot(Time,SL_wTd_nos_base_R0_RCP70_cdf[0:-1,10])\n", "plt.plot(Time,SL_wTd_nos_base_R0_RCP70_cdf[0:-1,50])\n", "plt.plot(Time,SL_wTd_nos_base_R0_RCP70_cdf[0:-1,166])\n", "plt.plot(Time,SL_wTd_nos_base_R0_RCP70_cdf[0:-1,500])\n", "plt.plot(Time,SL_wTd_nos_base_R0_RCP70_cdf[0:-1,833])\n", "plt.plot(Time,SL_wTd_nos_base_R0_RCP70_cdf[0:-1,950])\n", "plt.plot(Time,SL_wTd_nos_base_R0_RCP70_cdf[0:-1,990])\n" ] }, { "cell_type": "code", "execution_count": 25, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "0.3733869194984436\n", "0.6061061024665833\n", "1.1127036809921265\n" ] } ], "source": [ "print(SL_wTd_nos_base_R0_RCP70_cdf[-1,833])\n", "print(SL_wTd_nos_base_R0_RCP70_cdf[-1,950])\n", "print(SL_wTd_nos_base_R0_RCP70_cdf[-1,990])" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.7.1" } }, "nbformat": 4, "nbformat_minor": 2 }