Modeled - Precip - 20 year - monthly.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "cd9e6c27",
   "metadata": {},
   "source": [
    "## WICCI Downscaled Data: Precip - 20 year - monthly\n",
    "\n",
    "Objectives\n",
    "* aggregate data for prcp, e.g. mean, standard deviation, days over 1\",2\",3\", for RCP4.5 and RCP8.5\n",
    "* run for every 20 year model window, e.g. 2021-2040 for 2030 average\n",
    "* calculate averages by whole month increments\n",
    "* creat new netcdf file(s) for aggregate data for each 20-year timeframe\n",
    "* calculate intermodel standard deviation differently -- find average of each model within the current time window and\n",
    "  then calculate standard deviation across models\n",
    "\n",
    "Eric Compas, compase@uww.edu 11/17/2021, 1/19/2022, 1/25/2022, 2/28/2022, 6/6/2022"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "b488720d",
   "metadata": {},
   "outputs": [],
   "source": [
    "import netCDF4\n",
    "import numpy as np\n",
    "import os\n",
    "import datetime\n",
    "import gc\n",
    "from netCDF4 import Dataset,num2date,date2num"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "3b7aaae1",
   "metadata": {},
   "outputs": [],
   "source": [
    "import warnings\n",
    "warnings.filterwarnings(\"ignore\", category=DeprecationWarning)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "ff913133",
   "metadata": {},
   "outputs": [],
   "source": [
    "#base_folder = \"Z:/Climate_Data\"\n",
    "base_folder = \"F:/Climate_Data\"\n",
    "#base_folder = \"/Users/ericcompas/Climate_Data\"\n",
    "if not os.path.isdir(base_folder):\n",
    "    print(\"Base folder not valid\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "c3853f08",
   "metadata": {},
   "outputs": [],
   "source": [
    "#out_folder = \"Z:/Climate_Data/aggregate_files\"\n",
    "out_folder = \"C:/Users/Eric/Dropbox/Climate_Data/aggregate_files\"\n",
    "#out_folder  = \"/Users/ericcompas/Climate_Data/output\"\n",
    "if not os.path.isdir(out_folder):\n",
    "    print(\"Out folder not valid\")"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "e76f0f3e",
   "metadata": {},
   "source": [
    "The BIG LOOP.\n",
    "\n",
    "Look across all years, models, and create aggregate file for each year."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "f6d6f533",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Processing climate scenario rcp45\n",
      "Removed GFDL-CM3\n",
      "  Processing base year 2030\n",
      "    Processing GCM ACCESS1-0\n",
      "      Processing year 2021\n",
      "      Processing year 2022\n",
      "      Processing year 2023\n",
      "      Processing year 2024\n",
      "      Processing year 2025\n",
      "      Processing year 2026\n",
      "      Processing year 2027\n",
      "      Processing year 2028\n",
      "      Processing year 2029\n",
      "      Processing year 2030\n",
      "      Processing year 2031\n",
      "      Processing year 2032\n",
      "      Processing year 2033\n",
      "      Processing year 2034\n",
      "      Processing year 2035\n",
      "      Processing year 2036\n",
      "      Processing year 2037\n",
      "      Processing year 2038\n",
      "      Processing year 2039\n",
      "      Processing year 2040\n",
      "    Processing GCM ACCESS1-3\n",
      "      Processing year 2021\n",
      "      Processing year 2022\n",
      "      Processing year 2023\n",
      "      Processing year 2024\n",
      "      Processing year 2025\n",
      "      Processing year 2026\n",
      "      Processing year 2027\n",
      "      Processing year 2028\n",
      "      Processing year 2029\n",
      "      Processing year 2030\n",
      "      Processing year 2031\n",
      "      Processing year 2032\n",
      "      Processing year 2033\n",
      "      Processing year 2034\n",
      "      Processing year 2035\n",
      "      Processing year 2036\n",
      "      Processing year 2037\n",
      "      Processing year 2038\n",
      "      Processing year 2039\n",
      "      Processing year 2040\n",
      "    Processing GCM CanESM2\n",
      "      Processing year 2021\n",
      "      Processing year 2022\n",
      "      Processing year 2023\n",
      "      Processing year 2024\n",
      "      Processing year 2025\n",
      "      Processing year 2026\n",
      "      Processing year 2027\n",
      "      Processing year 2028\n",
      "      Processing year 2029\n",
      "      Processing year 2030\n",
      "      Processing year 2031\n",
      "      Processing year 2032\n",
      "      Processing year 2033\n",
      "      Processing year 2034\n",
      "      Processing year 2035\n",
      "      Processing year 2036\n",
      "      Processing year 2037\n",
      "      Processing year 2038\n",
      "      Processing year 2039\n",
      "      Processing year 2040\n",
      "    Processing GCM CMCC-CM\n",
      "      Processing year 2021\n",
      "      Processing year 2022\n",
      "      Processing year 2023\n",
      "      Processing year 2024\n",
      "      Processing year 2025\n",
      "      Processing year 2026\n",
      "      Processing year 2027\n",
      "      Processing year 2028\n",
      "      Processing year 2029\n",
      "      Processing year 2030\n",
      "      Processing year 2031\n",
      "      Processing year 2032\n",
      "      Processing year 2033\n",
      "      Processing year 2034\n",
      "      Processing year 2035\n",
      "      Processing year 2036\n",
      "      Processing year 2037\n",
      "      Processing year 2038\n",
      "      Processing year 2039\n",
      "      Processing year 2040\n",
      "    Processing GCM CMCC-CMS\n",
      "      Processing year 2021\n",
      "      Processing year 2022\n",
      "      Processing year 2023\n",
      "      Processing year 2024\n",
      "      Processing year 2025\n",
      "      Processing year 2026\n",
      "      Processing year 2027\n",
      "      Processing year 2028\n",
      "      Processing year 2029\n",
      "      Processing year 2030\n",
      "      Processing year 2031\n",
      "      Processing year 2032\n",
      "      Processing year 2033\n",
      "      Processing year 2034\n",
      "      Processing year 2035\n",
      "      Processing year 2036\n",
      "      Processing year 2037\n",
      "      Processing year 2038\n",
      "      Processing year 2039\n",
      "      Processing year 2040\n",
      "    Processing GCM CNRM-CM5\n",
      "      Processing year 2021\n",
      "      Processing year 2022\n",
      "      Processing year 2023\n",
      "      Processing year 2024\n",
      "      Processing year 2025\n",
      "      Processing year 2026\n",
      "      Processing year 2027\n",
      "      Processing year 2028\n",
      "      Processing year 2029\n",
      "      Processing year 2030\n",
      "      Processing year 2031\n",
      "      Processing year 2032\n",
      "      Processing year 2033\n",
      "      Processing year 2034\n",
      "      Processing year 2035\n",
      "      Processing year 2036\n",
      "      Processing year 2037\n",
      "      Processing year 2038\n",
      "      Processing year 2039\n",
      "      Processing year 2040\n",
      "    Processing GCM CSIRO-Mk3-6-0\n",
      "      Processing year 2021\n",
      "      Processing year 2022\n",
      "      Processing year 2023\n",
      "      Processing year 2024\n",
      "      Processing year 2025\n",
      "      Processing year 2026\n",
      "      Processing year 2027\n",
      "      Processing year 2028\n",
      "      Processing year 2029\n",
      "      Processing year 2030\n",
      "      Processing year 2031\n",
      "      Processing year 2032\n",
      "      Processing year 2033\n",
      "      Processing year 2034\n",
      "      Processing year 2035\n",
      "      Processing year 2036\n",
      "      Processing year 2037\n",
      "      Processing year 2038\n",
      "      Processing year 2039\n",
      "      Processing year 2040\n",
      "    Processing GCM GFDL-ESM2G\n",
      "      Processing year 2021\n",
      "      Processing year 2022\n",
      "      Processing year 2023\n",
      "      Processing year 2024\n",
      "      Processing year 2025\n",
      "      Processing year 2026\n",
      "      Processing year 2027\n",
      "      Processing year 2028\n",
      "      Processing year 2029\n",
      "      Processing year 2030\n",
      "      Processing year 2031\n",
      "      Processing year 2032\n",
      "      Processing year 2033\n",
      "      Processing year 2034\n",
      "      Processing year 2035\n",
      "      Processing year 2036\n",
      "      Processing year 2037\n",
      "      Processing year 2038\n",
      "      Processing year 2039\n",
      "      Processing year 2040\n",
      "    Processing GCM GFDL-ESM2M\n",
      "      Processing year 2021\n",
      "      Processing year 2022\n",
      "      Processing year 2023\n",
      "      Processing year 2024\n",
      "      Processing year 2025\n",
      "      Processing year 2026\n",
      "      Processing year 2027\n",
      "      Processing year 2028\n",
      "      Processing year 2029\n",
      "      Processing year 2030\n",
      "      Processing year 2031\n",
      "      Processing year 2032\n",
      "      Processing year 2033\n",
      "      Processing year 2034\n",
      "      Processing year 2035\n",
      "      Processing year 2036\n",
      "      Processing year 2037\n",
      "      Processing year 2038\n",
      "      Processing year 2039\n",
      "      Processing year 2040\n",
      "    Processing GCM HadGEM2-CC\n",
      "      Processing year 2021\n",
      "      Processing year 2022\n",
      "      Processing year 2023\n",
      "      Processing year 2024\n",
      "      Processing year 2025\n",
      "      Processing year 2026\n",
      "      Processing year 2027\n",
      "      Processing year 2028\n",
      "      Processing year 2029\n",
      "      Processing year 2030\n",
      "      Processing year 2031\n",
      "      Processing year 2032\n",
      "      Processing year 2033\n",
      "      Processing year 2034\n",
      "      Processing year 2035\n",
      "      Processing year 2036\n",
      "      Processing year 2037\n",
      "      Processing year 2038\n",
      "      Processing year 2039\n",
      "      Processing year 2040\n",
      "    Processing GCM inmcm4\n",
      "      Processing year 2021\n",
      "      Processing year 2022\n",
      "      Processing year 2023\n",
      "      Processing year 2024\n",
      "      Processing year 2025\n",
      "      Processing year 2026\n",
      "      Processing year 2027\n",
      "      Processing year 2028\n",
      "      Processing year 2029\n",
      "      Processing year 2030\n",
      "      Processing year 2031\n",
      "      Processing year 2032\n",
      "      Processing year 2033\n",
      "      Processing year 2034\n",
      "      Processing year 2035\n",
      "      Processing year 2036\n",
      "      Processing year 2037\n",
      "      Processing year 2038\n",
      "      Processing year 2039\n",
      "      Processing year 2040\n",
      "    Processing GCM IPSL-CM5A-LR\n",
      "      Processing year 2021\n",
      "      Processing year 2022\n",
      "      Processing year 2023\n",
      "      Processing year 2024\n",
      "      Processing year 2025\n",
      "      Processing year 2026\n",
      "      Processing year 2027\n",
      "      Processing year 2028\n",
      "      Processing year 2029\n",
      "      Processing year 2030\n",
      "      Processing year 2031\n",
      "      Processing year 2032\n",
      "      Processing year 2033\n",
      "      Processing year 2034\n",
      "      Processing year 2035\n",
      "      Processing year 2036\n",
      "      Processing year 2037\n",
      "      Processing year 2038\n",
      "      Processing year 2039\n",
      "      Processing year 2040\n",
      "    Processing GCM IPSL-CM5A-MR\n",
      "      Processing year 2021\n",
      "      Processing year 2022\n",
      "      Processing year 2023\n",
      "      Processing year 2024\n",
      "      Processing year 2025\n",
      "      Processing year 2026\n",
      "      Processing year 2027\n",
      "      Processing year 2028\n",
      "      Processing year 2029\n",
      "      Processing year 2030\n",
      "      Processing year 2031\n",
      "      Processing year 2032\n",
      "      Processing year 2033\n",
      "      Processing year 2034\n",
      "      Processing year 2035\n",
      "      Processing year 2036\n",
      "      Processing year 2037\n",
      "      Processing year 2038\n",
      "      Processing year 2039\n",
      "      Processing year 2040\n",
      "    Processing GCM IPSL-CM5B-LR\n",
      "      Processing year 2021\n",
      "      Processing year 2022\n",
      "      Processing year 2023\n",
      "      Processing year 2024\n",
      "      Processing year 2025\n",
      "      Processing year 2026\n",
      "      Processing year 2027\n",
      "      Processing year 2028\n",
      "      Processing year 2029\n",
      "      Processing year 2030\n",
      "      Processing year 2031\n",
      "      Processing year 2032\n",
      "      Processing year 2033\n",
      "      Processing year 2034\n",
      "      Processing year 2035\n",
      "      Processing year 2036\n",
      "      Processing year 2037\n",
      "      Processing year 2038\n",
      "      Processing year 2039\n",
      "      Processing year 2040\n",
      "    Processing GCM MIROC-ESM\n",
      "      Processing year 2021\n",
      "      Processing year 2022\n",
      "      Processing year 2023\n",
      "      Processing year 2024\n",
      "      Processing year 2025\n",
      "      Processing year 2026\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "      Processing year 2027\n",
      "      Processing year 2028\n",
      "      Processing year 2029\n",
      "      Processing year 2030\n",
      "      Processing year 2031\n",
      "      Processing year 2032\n",
      "      Processing year 2033\n",
      "      Processing year 2034\n",
      "      Processing year 2035\n",
      "      Processing year 2036\n",
      "      Processing year 2037\n",
      "      Processing year 2038\n",
      "      Processing year 2039\n",
      "      Processing year 2040\n",
      "    Processing GCM MIROC-ESM-CHEM\n",
      "      Processing year 2021\n",
      "      Processing year 2022\n",
      "      Processing year 2023\n",
      "      Processing year 2024\n",
      "      Processing year 2025\n",
      "      Processing year 2026\n",
      "      Processing year 2027\n",
      "      Processing year 2028\n",
      "      Processing year 2029\n",
      "      Processing year 2030\n",
      "      Processing year 2031\n",
      "      Processing year 2032\n",
      "      Processing year 2033\n",
      "      Processing year 2034\n",
      "      Processing year 2035\n",
      "      Processing year 2036\n",
      "      Processing year 2037\n",
      "      Processing year 2038\n",
      "      Processing year 2039\n",
      "      Processing year 2040\n",
      "    Processing GCM MIROC5\n",
      "      Processing year 2021\n",
      "      Processing year 2022\n",
      "      Processing year 2023\n",
      "      Processing year 2024\n",
      "      Processing year 2025\n",
      "      Processing year 2026\n",
      "      Processing year 2027\n",
      "      Processing year 2028\n",
      "      Processing year 2029\n",
      "      Processing year 2030\n",
      "      Processing year 2031\n",
      "      Processing year 2032\n",
      "      Processing year 2033\n",
      "      Processing year 2034\n",
      "      Processing year 2035\n",
      "      Processing year 2036\n",
      "      Processing year 2037\n",
      "      Processing year 2038\n",
      "      Processing year 2039\n",
      "      Processing year 2040\n",
      "    Processing GCM MPI-ESM-LR\n",
      "      Processing year 2021\n",
      "      Processing year 2022\n",
      "      Processing year 2023\n",
      "      Processing year 2024\n",
      "      Processing year 2025\n",
      "      Processing year 2026\n",
      "      Processing year 2027\n",
      "      Processing year 2028\n",
      "      Processing year 2029\n",
      "      Processing year 2030\n",
      "      Processing year 2031\n",
      "      Processing year 2032\n",
      "      Processing year 2033\n",
      "      Processing year 2034\n",
      "      Processing year 2035\n",
      "      Processing year 2036\n",
      "      Processing year 2037\n",
      "      Processing year 2038\n",
      "      Processing year 2039\n",
      "      Processing year 2040\n",
      "    Processing GCM MPI-ESM-MR\n",
      "      Processing year 2021\n",
      "      Processing year 2022\n",
      "      Processing year 2023\n",
      "      Processing year 2024\n",
      "      Processing year 2025\n",
      "      Processing year 2026\n",
      "      Processing year 2027\n",
      "      Processing year 2028\n",
      "      Processing year 2029\n",
      "      Processing year 2030\n",
      "      Processing year 2031\n",
      "      Processing year 2032\n",
      "      Processing year 2033\n",
      "      Processing year 2034\n",
      "      Processing year 2035\n",
      "      Processing year 2036\n",
      "      Processing year 2037\n",
      "      Processing year 2038\n",
      "      Processing year 2039\n",
      "      Processing year 2040\n",
      "    Processing GCM MRI-CGCM3\n",
      "      Processing year 2021\n",
      "      Processing year 2022\n",
      "      Processing year 2023\n",
      "      Processing year 2024\n",
      "      Processing year 2025\n",
      "      Processing year 2026\n",
      "      Processing year 2027\n",
      "      Processing year 2028\n",
      "      Processing year 2029\n",
      "      Processing year 2030\n",
      "      Processing year 2031\n",
      "      Processing year 2032\n",
      "      Processing year 2033\n",
      "      Processing year 2034\n",
      "      Processing year 2035\n",
      "      Processing year 2036\n",
      "      Processing year 2037\n",
      "      Processing year 2038\n",
      "      Processing year 2039\n",
      "      Processing year 2040\n",
      "    Processing GCM NorESM1-M\n",
      "      Processing year 2021\n",
      "      Processing year 2022\n",
      "      Processing year 2023\n",
      "      Processing year 2024\n",
      "      Processing year 2025\n",
      "      Processing year 2026\n",
      "      Processing year 2027\n",
      "      Processing year 2028\n",
      "      Processing year 2029\n",
      "      Processing year 2030\n",
      "      Processing year 2031\n",
      "      Processing year 2032\n",
      "      Processing year 2033\n",
      "      Processing year 2034\n",
      "      Processing year 2035\n",
      "      Processing year 2036\n",
      "      Processing year 2037\n",
      "      Processing year 2038\n",
      "      Processing year 2039\n",
      "      Processing year 2040\n",
      "Writing time. Year = 2030\n",
      "Wrote: prcp_rcp45_2030_20yr_monthly.nc\n",
      "Cleared some memory up...restarting loop\n",
      "  Processing base year 2050\n",
      "    Processing GCM ACCESS1-0\n",
      "      Processing year 2041\n",
      "      Processing year 2042\n",
      "      Processing year 2043\n",
      "      Processing year 2044\n",
      "      Processing year 2045\n",
      "      Processing year 2046\n",
      "      Processing year 2047\n",
      "      Processing year 2048\n",
      "      Processing year 2049\n",
      "      Processing year 2050\n",
      "      Processing year 2051\n",
      "      Processing year 2052\n",
      "      Processing year 2053\n",
      "      Processing year 2054\n",
      "      Processing year 2055\n",
      "      Processing year 2056\n",
      "      Processing year 2057\n",
      "      Processing year 2058\n",
      "      Processing year 2059\n",
      "      Processing year 2060\n",
      "    Processing GCM ACCESS1-3\n",
      "      Processing year 2041\n",
      "      Processing year 2042\n",
      "      Processing year 2043\n",
      "      Processing year 2044\n",
      "      Processing year 2045\n",
      "      Processing year 2046\n",
      "      Processing year 2047\n",
      "      Processing year 2048\n",
      "      Processing year 2049\n",
      "      Processing year 2050\n",
      "      Processing year 2051\n",
      "      Processing year 2052\n",
      "      Processing year 2053\n",
      "      Processing year 2054\n",
      "      Processing year 2055\n",
      "      Processing year 2056\n",
      "      Processing year 2057\n",
      "      Processing year 2058\n",
      "      Processing year 2059\n",
      "      Processing year 2060\n",
      "    Processing GCM CanESM2\n",
      "      Processing year 2041\n",
      "      Processing year 2042\n",
      "      Processing year 2043\n",
      "      Processing year 2044\n",
      "      Processing year 2045\n",
      "      Processing year 2046\n",
      "      Processing year 2047\n",
      "      Processing year 2048\n",
      "      Processing year 2049\n",
      "      Processing year 2050\n",
      "      Processing year 2051\n",
      "      Processing year 2052\n",
      "      Processing year 2053\n",
      "      Processing year 2054\n",
      "      Processing year 2055\n",
      "      Processing year 2056\n",
      "      Processing year 2057\n",
      "      Processing year 2058\n",
      "      Processing year 2059\n",
      "      Processing year 2060\n",
      "    Processing GCM CMCC-CM\n",
      "      Processing year 2041\n",
      "      Processing year 2042\n",
      "      Processing year 2043\n",
      "      Processing year 2044\n",
      "      Processing year 2045\n",
      "      Processing year 2046\n",
      "      Processing year 2047\n",
      "      Processing year 2048\n",
      "      Processing year 2049\n",
      "      Processing year 2050\n",
      "      Processing year 2051\n",
      "      Processing year 2052\n",
      "      Processing year 2053\n",
      "      Processing year 2054\n",
      "      Processing year 2055\n",
      "      Processing year 2056\n",
      "      Processing year 2057\n",
      "      Processing year 2058\n",
      "      Processing year 2059\n",
      "      Processing year 2060\n",
      "    Processing GCM CMCC-CMS\n",
      "      Processing year 2041\n",
      "      Processing year 2042\n",
      "      Processing year 2043\n",
      "      Processing year 2044\n",
      "      Processing year 2045\n",
      "      Processing year 2046\n",
      "      Processing year 2047\n",
      "      Processing year 2048\n",
      "      Processing year 2049\n",
      "      Processing year 2050\n",
      "      Processing year 2051\n",
      "      Processing year 2052\n",
      "      Processing year 2053\n",
      "      Processing year 2054\n",
      "      Processing year 2055\n",
      "      Processing year 2056\n",
      "      Processing year 2057\n",
      "      Processing year 2058\n",
      "      Processing year 2059\n",
      "      Processing year 2060\n",
      "    Processing GCM CNRM-CM5\n",
      "      Processing year 2041\n",
      "      Processing year 2042\n",
      "      Processing year 2043\n",
      "      Processing year 2044\n",
      "      Processing year 2045\n",
      "      Processing year 2046\n",
      "      Processing year 2047\n",
      "      Processing year 2048\n",
      "      Processing year 2049\n",
      "      Processing year 2050\n",
      "      Processing year 2051\n",
      "      Processing year 2052\n",
      "      Processing year 2053\n",
      "      Processing year 2054\n",
      "      Processing year 2055\n",
      "      Processing year 2056\n",
      "      Processing year 2057\n",
      "      Processing year 2058\n",
      "      Processing year 2059\n",
      "      Processing year 2060\n",
      "    Processing GCM CSIRO-Mk3-6-0\n",
      "      Processing year 2041\n",
      "      Processing year 2042\n",
      "      Processing year 2043\n",
      "      Processing year 2044\n",
      "      Processing year 2045\n",
      "      Processing year 2046\n",
      "      Processing year 2047\n",
      "      Processing year 2048\n",
      "      Processing year 2049\n",
      "      Processing year 2050\n",
      "      Processing year 2051\n",
      "      Processing year 2052\n",
      "      Processing year 2053\n",
      "      Processing year 2054\n",
      "      Processing year 2055\n",
      "      Processing year 2056\n",
      "      Processing year 2057\n",
      "      Processing year 2058\n",
      "      Processing year 2059\n",
      "      Processing year 2060\n",
      "    Processing GCM GFDL-ESM2G\n",
      "      Processing year 2041\n",
      "      Processing year 2042\n",
      "      Processing year 2043\n",
      "      Processing year 2044\n",
      "      Processing year 2045\n",
      "      Processing year 2046\n",
      "      Processing year 2047\n",
      "      Processing year 2048\n",
      "      Processing year 2049\n",
      "      Processing year 2050\n",
      "      Processing year 2051\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "      Processing year 2052\n",
      "      Processing year 2053\n",
      "      Processing year 2054\n",
      "      Processing year 2055\n",
      "      Processing year 2056\n",
      "      Processing year 2057\n",
      "      Processing year 2058\n",
      "      Processing year 2059\n",
      "      Processing year 2060\n",
      "    Processing GCM GFDL-ESM2M\n",
      "      Processing year 2041\n",
      "      Processing year 2042\n",
      "      Processing year 2043\n",
      "      Processing year 2044\n",
      "      Processing year 2045\n",
      "      Processing year 2046\n",
      "      Processing year 2047\n",
      "      Processing year 2048\n",
      "      Processing year 2049\n",
      "      Processing year 2050\n",
      "      Processing year 2051\n",
      "      Processing year 2052\n",
      "      Processing year 2053\n",
      "      Processing year 2054\n",
      "      Processing year 2055\n",
      "      Processing year 2056\n",
      "      Processing year 2057\n",
      "      Processing year 2058\n",
      "      Processing year 2059\n",
      "      Processing year 2060\n",
      "    Processing GCM HadGEM2-CC\n",
      "      Processing year 2041\n",
      "      Processing year 2042\n",
      "      Processing year 2043\n",
      "      Processing year 2044\n",
      "      Processing year 2045\n",
      "      Processing year 2046\n",
      "      Processing year 2047\n",
      "      Processing year 2048\n",
      "      Processing year 2049\n",
      "      Processing year 2050\n",
      "      Processing year 2051\n",
      "      Processing year 2052\n",
      "      Processing year 2053\n",
      "      Processing year 2054\n",
      "      Processing year 2055\n",
      "      Processing year 2056\n",
      "      Processing year 2057\n",
      "      Processing year 2058\n",
      "      Processing year 2059\n",
      "      Processing year 2060\n",
      "    Processing GCM inmcm4\n",
      "      Processing year 2041\n",
      "      Processing year 2042\n",
      "      Processing year 2043\n",
      "      Processing year 2044\n",
      "      Processing year 2045\n",
      "      Processing year 2046\n",
      "      Processing year 2047\n",
      "      Processing year 2048\n",
      "      Processing year 2049\n",
      "      Processing year 2050\n",
      "      Processing year 2051\n",
      "      Processing year 2052\n",
      "      Processing year 2053\n",
      "      Processing year 2054\n",
      "      Processing year 2055\n",
      "      Processing year 2056\n",
      "      Processing year 2057\n",
      "      Processing year 2058\n",
      "      Processing year 2059\n",
      "      Processing year 2060\n",
      "    Processing GCM IPSL-CM5A-LR\n",
      "      Processing year 2041\n",
      "      Processing year 2042\n",
      "      Processing year 2043\n",
      "      Processing year 2044\n",
      "      Processing year 2045\n",
      "      Processing year 2046\n",
      "      Processing year 2047\n",
      "      Processing year 2048\n",
      "      Processing year 2049\n",
      "      Processing year 2050\n",
      "      Processing year 2051\n",
      "      Processing year 2052\n",
      "      Processing year 2053\n",
      "      Processing year 2054\n",
      "      Processing year 2055\n",
      "      Processing year 2056\n",
      "      Processing year 2057\n",
      "      Processing year 2058\n",
      "      Processing year 2059\n",
      "      Processing year 2060\n",
      "    Processing GCM IPSL-CM5A-MR\n",
      "      Processing year 2041\n",
      "      Processing year 2042\n",
      "      Processing year 2043\n",
      "      Processing year 2044\n",
      "      Processing year 2045\n",
      "      Processing year 2046\n",
      "      Processing year 2047\n",
      "      Processing year 2048\n",
      "      Processing year 2049\n",
      "      Processing year 2050\n",
      "      Processing year 2051\n",
      "      Processing year 2052\n",
      "      Processing year 2053\n",
      "      Processing year 2054\n",
      "      Processing year 2055\n",
      "      Processing year 2056\n",
      "      Processing year 2057\n",
      "      Processing year 2058\n",
      "      Processing year 2059\n",
      "      Processing year 2060\n",
      "    Processing GCM IPSL-CM5B-LR\n",
      "      Processing year 2041\n",
      "      Processing year 2042\n",
      "      Processing year 2043\n",
      "      Processing year 2044\n",
      "      Processing year 2045\n",
      "      Processing year 2046\n",
      "      Processing year 2047\n",
      "      Processing year 2048\n",
      "      Processing year 2049\n",
      "      Processing year 2050\n",
      "      Processing year 2051\n",
      "      Processing year 2052\n",
      "      Processing year 2053\n",
      "      Processing year 2054\n",
      "      Processing year 2055\n",
      "      Processing year 2056\n",
      "      Processing year 2057\n",
      "      Processing year 2058\n",
      "      Processing year 2059\n",
      "      Processing year 2060\n",
      "    Processing GCM MIROC-ESM\n",
      "      Processing year 2041\n",
      "      Processing year 2042\n",
      "      Processing year 2043\n",
      "      Processing year 2044\n",
      "      Processing year 2045\n",
      "      Processing year 2046\n",
      "      Processing year 2047\n",
      "      Processing year 2048\n",
      "      Processing year 2049\n",
      "      Processing year 2050\n",
      "      Processing year 2051\n",
      "      Processing year 2052\n",
      "      Processing year 2053\n",
      "      Processing year 2054\n",
      "      Processing year 2055\n",
      "      Processing year 2056\n",
      "      Processing year 2057\n",
      "      Processing year 2058\n",
      "      Processing year 2059\n",
      "      Processing year 2060\n",
      "    Processing GCM MIROC-ESM-CHEM\n",
      "      Processing year 2041\n",
      "      Processing year 2042\n",
      "      Processing year 2043\n",
      "      Processing year 2044\n",
      "      Processing year 2045\n",
      "      Processing year 2046\n",
      "      Processing year 2047\n",
      "      Processing year 2048\n",
      "      Processing year 2049\n",
      "      Processing year 2050\n",
      "      Processing year 2051\n",
      "      Processing year 2052\n",
      "      Processing year 2053\n",
      "      Processing year 2054\n",
      "      Processing year 2055\n",
      "      Processing year 2056\n",
      "      Processing year 2057\n",
      "      Processing year 2058\n",
      "      Processing year 2059\n",
      "      Processing year 2060\n",
      "    Processing GCM MIROC5\n",
      "      Processing year 2041\n",
      "      Processing year 2042\n",
      "      Processing year 2043\n",
      "      Processing year 2044\n",
      "      Processing year 2045\n",
      "      Processing year 2046\n",
      "      Processing year 2047\n",
      "      Processing year 2048\n",
      "      Processing year 2049\n",
      "      Processing year 2050\n",
      "      Processing year 2051\n",
      "      Processing year 2052\n",
      "      Processing year 2053\n",
      "      Processing year 2054\n",
      "      Processing year 2055\n",
      "      Processing year 2056\n",
      "      Processing year 2057\n",
      "      Processing year 2058\n",
      "      Processing year 2059\n",
      "      Processing year 2060\n",
      "    Processing GCM MPI-ESM-LR\n",
      "      Processing year 2041\n",
      "      Processing year 2042\n",
      "      Processing year 2043\n",
      "      Processing year 2044\n",
      "      Processing year 2045\n",
      "      Processing year 2046\n",
      "      Processing year 2047\n",
      "      Processing year 2048\n",
      "      Processing year 2049\n",
      "      Processing year 2050\n",
      "      Processing year 2051\n",
      "      Processing year 2052\n",
      "      Processing year 2053\n",
      "      Processing year 2054\n",
      "      Processing year 2055\n",
      "      Processing year 2056\n",
      "      Processing year 2057\n",
      "      Processing year 2058\n",
      "      Processing year 2059\n",
      "      Processing year 2060\n",
      "    Processing GCM MPI-ESM-MR\n",
      "      Processing year 2041\n",
      "      Processing year 2042\n",
      "      Processing year 2043\n",
      "      Processing year 2044\n",
      "      Processing year 2045\n",
      "      Processing year 2046\n",
      "      Processing year 2047\n",
      "      Processing year 2048\n",
      "      Processing year 2049\n",
      "      Processing year 2050\n",
      "      Processing year 2051\n",
      "      Processing year 2052\n",
      "      Processing year 2053\n",
      "      Processing year 2054\n",
      "      Processing year 2055\n",
      "      Processing year 2056\n",
      "      Processing year 2057\n",
      "      Processing year 2058\n",
      "      Processing year 2059\n",
      "      Processing year 2060\n",
      "    Processing GCM MRI-CGCM3\n",
      "      Processing year 2041\n",
      "      Processing year 2042\n",
      "      Processing year 2043\n",
      "      Processing year 2044\n",
      "      Processing year 2045\n",
      "      Processing year 2046\n",
      "      Processing year 2047\n",
      "      Processing year 2048\n",
      "      Processing year 2049\n",
      "      Processing year 2050\n",
      "      Processing year 2051\n",
      "      Processing year 2052\n",
      "      Processing year 2053\n",
      "      Processing year 2054\n",
      "      Processing year 2055\n",
      "      Processing year 2056\n",
      "      Processing year 2057\n",
      "      Processing year 2058\n",
      "      Processing year 2059\n",
      "      Processing year 2060\n",
      "    Processing GCM NorESM1-M\n",
      "      Processing year 2041\n",
      "      Processing year 2042\n",
      "      Processing year 2043\n",
      "      Processing year 2044\n",
      "      Processing year 2045\n",
      "      Processing year 2046\n",
      "      Processing year 2047\n",
      "      Processing year 2048\n",
      "      Processing year 2049\n",
      "      Processing year 2050\n",
      "      Processing year 2051\n",
      "      Processing year 2052\n",
      "      Processing year 2053\n",
      "      Processing year 2054\n",
      "      Processing year 2055\n",
      "      Processing year 2056\n",
      "      Processing year 2057\n",
      "      Processing year 2058\n",
      "      Processing year 2059\n",
      "      Processing year 2060\n",
      "Writing time. Year = 2050\n",
      "Wrote: prcp_rcp45_2050_20yr_monthly.nc\n",
      "Cleared some memory up...restarting loop\n",
      "  Processing base year 2070\n",
      "    Processing GCM ACCESS1-0\n",
      "      Processing year 2061\n",
      "      Processing year 2062\n",
      "      Processing year 2063\n",
      "      Processing year 2064\n",
      "      Processing year 2065\n",
      "      Processing year 2066\n",
      "      Processing year 2067\n",
      "      Processing year 2068\n",
      "      Processing year 2069\n",
      "      Processing year 2070\n",
      "      Processing year 2071\n",
      "      Processing year 2072\n",
      "      Processing year 2073\n",
      "      Processing year 2074\n",
      "      Processing year 2075\n",
      "      Processing year 2076\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "      Processing year 2077\n",
      "      Processing year 2078\n",
      "      Processing year 2079\n",
      "      Processing year 2080\n",
      "    Processing GCM ACCESS1-3\n",
      "      Processing year 2061\n",
      "      Processing year 2062\n",
      "      Processing year 2063\n",
      "      Processing year 2064\n",
      "      Processing year 2065\n",
      "      Processing year 2066\n",
      "      Processing year 2067\n",
      "      Processing year 2068\n",
      "      Processing year 2069\n",
      "      Processing year 2070\n",
      "      Processing year 2071\n",
      "      Processing year 2072\n",
      "      Processing year 2073\n",
      "      Processing year 2074\n",
      "      Processing year 2075\n",
      "      Processing year 2076\n",
      "      Processing year 2077\n",
      "      Processing year 2078\n",
      "      Processing year 2079\n",
      "      Processing year 2080\n",
      "    Processing GCM CanESM2\n",
      "      Processing year 2061\n",
      "      Processing year 2062\n",
      "      Processing year 2063\n",
      "      Processing year 2064\n",
      "      Processing year 2065\n",
      "      Processing year 2066\n",
      "      Processing year 2067\n",
      "      Processing year 2068\n",
      "      Processing year 2069\n",
      "      Processing year 2070\n",
      "      Processing year 2071\n",
      "      Processing year 2072\n",
      "      Processing year 2073\n",
      "      Processing year 2074\n",
      "      Processing year 2075\n",
      "      Processing year 2076\n",
      "      Processing year 2077\n",
      "      Processing year 2078\n",
      "      Processing year 2079\n",
      "      Processing year 2080\n",
      "    Processing GCM CMCC-CM\n",
      "      Processing year 2061\n",
      "      Processing year 2062\n",
      "      Processing year 2063\n",
      "      Processing year 2064\n",
      "      Processing year 2065\n",
      "      Processing year 2066\n",
      "      Processing year 2067\n",
      "      Processing year 2068\n",
      "      Processing year 2069\n",
      "      Processing year 2070\n",
      "      Processing year 2071\n",
      "      Processing year 2072\n",
      "      Processing year 2073\n",
      "      Processing year 2074\n",
      "      Processing year 2075\n",
      "      Processing year 2076\n",
      "      Processing year 2077\n",
      "      Processing year 2078\n",
      "      Processing year 2079\n",
      "      Processing year 2080\n",
      "    Processing GCM CMCC-CMS\n",
      "      Processing year 2061\n",
      "      Processing year 2062\n",
      "      Processing year 2063\n",
      "      Processing year 2064\n",
      "      Processing year 2065\n",
      "      Processing year 2066\n",
      "      Processing year 2067\n",
      "      Processing year 2068\n",
      "      Processing year 2069\n",
      "      Processing year 2070\n",
      "      Processing year 2071\n",
      "      Processing year 2072\n",
      "      Processing year 2073\n",
      "      Processing year 2074\n",
      "      Processing year 2075\n",
      "      Processing year 2076\n",
      "      Processing year 2077\n",
      "      Processing year 2078\n",
      "      Processing year 2079\n",
      "      Processing year 2080\n",
      "    Processing GCM CNRM-CM5\n",
      "      Processing year 2061\n",
      "      Processing year 2062\n",
      "      Processing year 2063\n",
      "      Processing year 2064\n",
      "      Processing year 2065\n",
      "      Processing year 2066\n",
      "      Processing year 2067\n",
      "      Processing year 2068\n",
      "      Processing year 2069\n",
      "      Processing year 2070\n",
      "      Processing year 2071\n",
      "      Processing year 2072\n",
      "      Processing year 2073\n",
      "      Processing year 2074\n",
      "      Processing year 2075\n",
      "      Processing year 2076\n",
      "      Processing year 2077\n",
      "      Processing year 2078\n",
      "      Processing year 2079\n",
      "      Processing year 2080\n",
      "    Processing GCM CSIRO-Mk3-6-0\n",
      "      Processing year 2061\n",
      "      Processing year 2062\n",
      "      Processing year 2063\n",
      "      Processing year 2064\n",
      "      Processing year 2065\n",
      "      Processing year 2066\n",
      "      Processing year 2067\n",
      "      Processing year 2068\n",
      "      Processing year 2069\n",
      "      Processing year 2070\n",
      "      Processing year 2071\n",
      "      Processing year 2072\n",
      "      Processing year 2073\n",
      "      Processing year 2074\n",
      "      Processing year 2075\n",
      "      Processing year 2076\n",
      "      Processing year 2077\n",
      "      Processing year 2078\n",
      "      Processing year 2079\n",
      "      Processing year 2080\n",
      "    Processing GCM GFDL-ESM2G\n",
      "      Processing year 2061\n",
      "      Processing year 2062\n",
      "      Processing year 2063\n",
      "      Processing year 2064\n",
      "      Processing year 2065\n",
      "      Processing year 2066\n",
      "      Processing year 2067\n",
      "      Processing year 2068\n",
      "      Processing year 2069\n",
      "      Processing year 2070\n",
      "      Processing year 2071\n",
      "      Processing year 2072\n",
      "      Processing year 2073\n",
      "      Processing year 2074\n",
      "      Processing year 2075\n",
      "      Processing year 2076\n",
      "      Processing year 2077\n",
      "      Processing year 2078\n",
      "      Processing year 2079\n",
      "      Processing year 2080\n",
      "    Processing GCM GFDL-ESM2M\n",
      "      Processing year 2061\n",
      "      Processing year 2062\n",
      "      Processing year 2063\n",
      "      Processing year 2064\n",
      "      Processing year 2065\n",
      "      Processing year 2066\n",
      "      Processing year 2067\n",
      "      Processing year 2068\n",
      "      Processing year 2069\n",
      "      Processing year 2070\n",
      "      Processing year 2071\n",
      "      Processing year 2072\n",
      "      Processing year 2073\n",
      "      Processing year 2074\n",
      "      Processing year 2075\n",
      "      Processing year 2076\n",
      "      Processing year 2077\n",
      "      Processing year 2078\n",
      "      Processing year 2079\n",
      "      Processing year 2080\n",
      "    Processing GCM HadGEM2-CC\n",
      "      Processing year 2061\n",
      "      Processing year 2062\n",
      "      Processing year 2063\n",
      "      Processing year 2064\n",
      "      Processing year 2065\n",
      "      Processing year 2066\n",
      "      Processing year 2067\n",
      "      Processing year 2068\n",
      "      Processing year 2069\n",
      "      Processing year 2070\n",
      "      Processing year 2071\n",
      "      Processing year 2072\n",
      "      Processing year 2073\n",
      "      Processing year 2074\n",
      "      Processing year 2075\n",
      "      Processing year 2076\n",
      "      Processing year 2077\n",
      "      Processing year 2078\n",
      "      Processing year 2079\n",
      "      Processing year 2080\n",
      "    Processing GCM inmcm4\n",
      "      Processing year 2061\n",
      "      Processing year 2062\n",
      "      Processing year 2063\n",
      "      Processing year 2064\n",
      "      Processing year 2065\n",
      "      Processing year 2066\n",
      "      Processing year 2067\n",
      "      Processing year 2068\n",
      "      Processing year 2069\n",
      "      Processing year 2070\n",
      "      Processing year 2071\n",
      "      Processing year 2072\n",
      "      Processing year 2073\n",
      "      Processing year 2074\n",
      "      Processing year 2075\n",
      "      Processing year 2076\n",
      "      Processing year 2077\n",
      "      Processing year 2078\n",
      "      Processing year 2079\n",
      "      Processing year 2080\n",
      "    Processing GCM IPSL-CM5A-LR\n",
      "      Processing year 2061\n",
      "      Processing year 2062\n",
      "      Processing year 2063\n",
      "      Processing year 2064\n",
      "      Processing year 2065\n",
      "      Processing year 2066\n",
      "      Processing year 2067\n",
      "      Processing year 2068\n",
      "      Processing year 2069\n",
      "      Processing year 2070\n",
      "      Processing year 2071\n",
      "      Processing year 2072\n",
      "      Processing year 2073\n",
      "      Processing year 2074\n",
      "      Processing year 2075\n",
      "      Processing year 2076\n",
      "      Processing year 2077\n",
      "      Processing year 2078\n",
      "      Processing year 2079\n",
      "      Processing year 2080\n",
      "    Processing GCM IPSL-CM5A-MR\n",
      "      Processing year 2061\n",
      "      Processing year 2062\n",
      "      Processing year 2063\n",
      "      Processing year 2064\n",
      "      Processing year 2065\n",
      "      Processing year 2066\n",
      "      Processing year 2067\n",
      "      Processing year 2068\n",
      "      Processing year 2069\n",
      "      Processing year 2070\n",
      "      Processing year 2071\n",
      "      Processing year 2072\n",
      "      Processing year 2073\n",
      "      Processing year 2074\n",
      "      Processing year 2075\n",
      "      Processing year 2076\n",
      "      Processing year 2077\n",
      "      Processing year 2078\n",
      "      Processing year 2079\n",
      "      Processing year 2080\n",
      "    Processing GCM IPSL-CM5B-LR\n",
      "      Processing year 2061\n",
      "      Processing year 2062\n",
      "      Processing year 2063\n",
      "      Processing year 2064\n",
      "      Processing year 2065\n",
      "      Processing year 2066\n",
      "      Processing year 2067\n",
      "      Processing year 2068\n",
      "      Processing year 2069\n",
      "      Processing year 2070\n",
      "      Processing year 2071\n",
      "      Processing year 2072\n",
      "      Processing year 2073\n",
      "      Processing year 2074\n",
      "      Processing year 2075\n",
      "      Processing year 2076\n",
      "      Processing year 2077\n",
      "      Processing year 2078\n",
      "      Processing year 2079\n",
      "      Processing year 2080\n",
      "    Processing GCM MIROC-ESM\n",
      "      Processing year 2061\n",
      "      Processing year 2062\n",
      "      Processing year 2063\n",
      "      Processing year 2064\n",
      "      Processing year 2065\n",
      "      Processing year 2066\n",
      "      Processing year 2067\n",
      "      Processing year 2068\n",
      "      Processing year 2069\n",
      "      Processing year 2070\n",
      "      Processing year 2071\n",
      "      Processing year 2072\n",
      "      Processing year 2073\n",
      "      Processing year 2074\n",
      "      Processing year 2075\n",
      "      Processing year 2076\n",
      "      Processing year 2077\n",
      "      Processing year 2078\n",
      "      Processing year 2079\n",
      "      Processing year 2080\n",
      "    Processing GCM MIROC-ESM-CHEM\n",
      "      Processing year 2061\n",
      "      Processing year 2062\n",
      "      Processing year 2063\n",
      "      Processing year 2064\n",
      "      Processing year 2065\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "      Processing year 2066\n",
      "      Processing year 2067\n",
      "      Processing year 2068\n",
      "      Processing year 2069\n",
      "      Processing year 2070\n",
      "      Processing year 2071\n",
      "      Processing year 2072\n",
      "      Processing year 2073\n",
      "      Processing year 2074\n",
      "      Processing year 2075\n",
      "      Processing year 2076\n",
      "      Processing year 2077\n",
      "      Processing year 2078\n",
      "      Processing year 2079\n",
      "      Processing year 2080\n",
      "    Processing GCM MIROC5\n",
      "      Processing year 2061\n",
      "      Processing year 2062\n",
      "      Processing year 2063\n",
      "      Processing year 2064\n",
      "      Processing year 2065\n",
      "      Processing year 2066\n",
      "      Processing year 2067\n",
      "      Processing year 2068\n",
      "      Processing year 2069\n",
      "      Processing year 2070\n",
      "      Processing year 2071\n",
      "      Processing year 2072\n",
      "      Processing year 2073\n",
      "      Processing year 2074\n",
      "      Processing year 2075\n",
      "      Processing year 2076\n",
      "      Processing year 2077\n",
      "      Processing year 2078\n",
      "      Processing year 2079\n",
      "      Processing year 2080\n",
      "    Processing GCM MPI-ESM-LR\n",
      "      Processing year 2061\n",
      "      Processing year 2062\n",
      "      Processing year 2063\n",
      "      Processing year 2064\n",
      "      Processing year 2065\n",
      "      Processing year 2066\n",
      "      Processing year 2067\n",
      "      Processing year 2068\n",
      "      Processing year 2069\n",
      "      Processing year 2070\n",
      "      Processing year 2071\n",
      "      Processing year 2072\n",
      "      Processing year 2073\n",
      "      Processing year 2074\n",
      "      Processing year 2075\n",
      "      Processing year 2076\n",
      "      Processing year 2077\n",
      "      Processing year 2078\n",
      "      Processing year 2079\n",
      "      Processing year 2080\n",
      "    Processing GCM MPI-ESM-MR\n",
      "      Processing year 2061\n",
      "      Processing year 2062\n",
      "      Processing year 2063\n",
      "      Processing year 2064\n",
      "      Processing year 2065\n",
      "      Processing year 2066\n",
      "      Processing year 2067\n",
      "      Processing year 2068\n",
      "      Processing year 2069\n",
      "      Processing year 2070\n",
      "      Processing year 2071\n",
      "      Processing year 2072\n",
      "      Processing year 2073\n",
      "      Processing year 2074\n",
      "      Processing year 2075\n",
      "      Processing year 2076\n",
      "      Processing year 2077\n",
      "      Processing year 2078\n",
      "      Processing year 2079\n",
      "      Processing year 2080\n",
      "    Processing GCM MRI-CGCM3\n",
      "      Processing year 2061\n",
      "      Processing year 2062\n",
      "      Processing year 2063\n",
      "      Processing year 2064\n",
      "      Processing year 2065\n",
      "      Processing year 2066\n",
      "      Processing year 2067\n",
      "      Processing year 2068\n",
      "      Processing year 2069\n",
      "      Processing year 2070\n",
      "      Processing year 2071\n",
      "      Processing year 2072\n",
      "      Processing year 2073\n",
      "      Processing year 2074\n",
      "      Processing year 2075\n",
      "      Processing year 2076\n",
      "      Processing year 2077\n",
      "      Processing year 2078\n",
      "      Processing year 2079\n",
      "      Processing year 2080\n",
      "    Processing GCM NorESM1-M\n",
      "      Processing year 2061\n",
      "      Processing year 2062\n",
      "      Processing year 2063\n",
      "      Processing year 2064\n",
      "      Processing year 2065\n",
      "      Processing year 2066\n",
      "      Processing year 2067\n",
      "      Processing year 2068\n",
      "      Processing year 2069\n",
      "      Processing year 2070\n",
      "      Processing year 2071\n",
      "      Processing year 2072\n",
      "      Processing year 2073\n",
      "      Processing year 2074\n",
      "      Processing year 2075\n",
      "      Processing year 2076\n",
      "      Processing year 2077\n",
      "      Processing year 2078\n",
      "      Processing year 2079\n",
      "      Processing year 2080\n",
      "Writing time. Year = 2070\n",
      "Wrote: prcp_rcp45_2070_20yr_monthly.nc\n",
      "Cleared some memory up...restarting loop\n",
      "  Processing base year 2090\n",
      "    Processing GCM ACCESS1-0\n",
      "      Processing year 2081\n",
      "      Processing year 2082\n",
      "      Processing year 2083\n",
      "      Processing year 2084\n",
      "      Processing year 2085\n",
      "      Processing year 2086\n",
      "      Processing year 2087\n",
      "      Processing year 2088\n",
      "      Processing year 2089\n",
      "      Processing year 2090\n",
      "      Processing year 2091\n",
      "      Processing year 2092\n",
      "      Processing year 2093\n",
      "      Processing year 2094\n",
      "      Processing year 2095\n",
      "      Processing year 2096\n",
      "      Processing year 2097\n",
      "      Processing year 2098\n",
      "      Processing year 2099\n",
      "      Processing year 2100\n",
      "    Processing GCM ACCESS1-3\n",
      "      Processing year 2081\n",
      "      Processing year 2082\n",
      "      Processing year 2083\n",
      "      Processing year 2084\n",
      "      Processing year 2085\n",
      "      Processing year 2086\n",
      "      Processing year 2087\n",
      "      Processing year 2088\n",
      "      Processing year 2089\n",
      "      Processing year 2090\n",
      "      Processing year 2091\n",
      "      Processing year 2092\n",
      "      Processing year 2093\n",
      "      Processing year 2094\n",
      "      Processing year 2095\n",
      "      Processing year 2096\n",
      "      Processing year 2097\n",
      "      Processing year 2098\n",
      "      Processing year 2099\n",
      "      Processing year 2100\n",
      "    Processing GCM CanESM2\n",
      "      Processing year 2081\n",
      "      Processing year 2082\n",
      "      Processing year 2083\n",
      "      Processing year 2084\n",
      "      Processing year 2085\n",
      "      Processing year 2086\n",
      "      Processing year 2087\n",
      "      Processing year 2088\n",
      "      Processing year 2089\n",
      "      Processing year 2090\n",
      "      Processing year 2091\n",
      "      Processing year 2092\n",
      "      Processing year 2093\n",
      "      Processing year 2094\n",
      "      Processing year 2095\n",
      "      Processing year 2096\n",
      "      Processing year 2097\n",
      "      Processing year 2098\n",
      "      Processing year 2099\n",
      "      Processing year 2100\n",
      "    Processing GCM CMCC-CM\n",
      "      Processing year 2081\n",
      "      Processing year 2082\n",
      "      Processing year 2083\n",
      "      Processing year 2084\n",
      "      Processing year 2085\n",
      "      Processing year 2086\n",
      "      Processing year 2087\n",
      "      Processing year 2088\n",
      "      Processing year 2089\n",
      "      Processing year 2090\n",
      "      Processing year 2091\n",
      "      Processing year 2092\n",
      "      Processing year 2093\n",
      "      Processing year 2094\n",
      "      Processing year 2095\n",
      "      Processing year 2096\n",
      "      Processing year 2097\n",
      "      Processing year 2098\n",
      "      Processing year 2099\n",
      "      Processing year 2100\n",
      "    Processing GCM CMCC-CMS\n",
      "      Processing year 2081\n",
      "      Processing year 2082\n",
      "      Processing year 2083\n",
      "      Processing year 2084\n",
      "      Processing year 2085\n",
      "      Processing year 2086\n",
      "      Processing year 2087\n",
      "      Processing year 2088\n",
      "      Processing year 2089\n",
      "      Processing year 2090\n",
      "      Processing year 2091\n",
      "      Processing year 2092\n",
      "      Processing year 2093\n",
      "      Processing year 2094\n",
      "      Processing year 2095\n",
      "      Processing year 2096\n",
      "      Processing year 2097\n",
      "      Processing year 2098\n",
      "      Processing year 2099\n",
      "      Processing year 2100\n",
      "    Processing GCM CNRM-CM5\n",
      "      Processing year 2081\n",
      "      Processing year 2082\n",
      "      Processing year 2083\n",
      "      Processing year 2084\n",
      "      Processing year 2085\n",
      "      Processing year 2086\n",
      "      Processing year 2087\n",
      "      Processing year 2088\n",
      "      Processing year 2089\n",
      "      Processing year 2090\n",
      "      Processing year 2091\n",
      "      Processing year 2092\n",
      "      Processing year 2093\n",
      "      Processing year 2094\n",
      "      Processing year 2095\n",
      "      Processing year 2096\n",
      "      Processing year 2097\n",
      "      Processing year 2098\n",
      "      Processing year 2099\n",
      "      Processing year 2100\n",
      "    Processing GCM CSIRO-Mk3-6-0\n",
      "      Processing year 2081\n",
      "      Processing year 2082\n",
      "      Processing year 2083\n",
      "      Processing year 2084\n",
      "      Processing year 2085\n",
      "      Processing year 2086\n",
      "      Processing year 2087\n",
      "      Processing year 2088\n",
      "      Processing year 2089\n",
      "      Processing year 2090\n",
      "      Processing year 2091\n",
      "      Processing year 2092\n",
      "      Processing year 2093\n",
      "      Processing year 2094\n",
      "      Processing year 2095\n",
      "      Processing year 2096\n",
      "      Processing year 2097\n",
      "      Processing year 2098\n",
      "      Processing year 2099\n",
      "      Processing year 2100\n",
      "    Processing GCM GFDL-ESM2G\n",
      "      Processing year 2081\n",
      "      Processing year 2082\n",
      "      Processing year 2083\n",
      "      Processing year 2084\n",
      "      Processing year 2085\n",
      "      Processing year 2086\n",
      "      Processing year 2087\n",
      "      Processing year 2088\n",
      "      Processing year 2089\n",
      "      Processing year 2090\n",
      "      Processing year 2091\n",
      "      Processing year 2092\n",
      "      Processing year 2093\n",
      "      Processing year 2094\n",
      "      Processing year 2095\n",
      "      Processing year 2096\n",
      "      Processing year 2097\n",
      "      Processing year 2098\n",
      "      Processing year 2099\n",
      "      Processing year 2100\n",
      "    Processing GCM GFDL-ESM2M\n",
      "      Processing year 2081\n",
      "      Processing year 2082\n",
      "      Processing year 2083\n",
      "      Processing year 2084\n",
      "      Processing year 2085\n",
      "      Processing year 2086\n",
      "      Processing year 2087\n",
      "      Processing year 2088\n",
      "      Processing year 2089\n",
      "      Processing year 2090\n",
      "      Processing year 2091\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "      Processing year 2092\n",
      "      Processing year 2093\n",
      "      Processing year 2094\n",
      "      Processing year 2095\n",
      "      Processing year 2096\n",
      "      Processing year 2097\n",
      "      Processing year 2098\n",
      "      Processing year 2099\n",
      "      Processing year 2100\n",
      "    Processing GCM HadGEM2-CC\n",
      "      Processing year 2081\n",
      "      Processing year 2082\n",
      "      Processing year 2083\n",
      "      Processing year 2084\n",
      "      Processing year 2085\n",
      "      Processing year 2086\n",
      "      Processing year 2087\n",
      "      Processing year 2088\n",
      "      Processing year 2089\n",
      "      Processing year 2090\n",
      "      Processing year 2091\n",
      "      Processing year 2092\n",
      "      Processing year 2093\n",
      "      Processing year 2094\n",
      "      Processing year 2095\n",
      "      Processing year 2096\n",
      "      Processing year 2097\n",
      "      Processing year 2098\n",
      "      Processing year 2099\n",
      "      Processing year 2100\n",
      "    Processing GCM inmcm4\n",
      "      Processing year 2081\n",
      "      Processing year 2082\n",
      "      Processing year 2083\n",
      "      Processing year 2084\n",
      "      Processing year 2085\n",
      "      Processing year 2086\n",
      "      Processing year 2087\n",
      "      Processing year 2088\n",
      "      Processing year 2089\n",
      "      Processing year 2090\n",
      "      Processing year 2091\n",
      "      Processing year 2092\n",
      "      Processing year 2093\n",
      "      Processing year 2094\n",
      "      Processing year 2095\n",
      "      Processing year 2096\n",
      "      Processing year 2097\n",
      "      Processing year 2098\n",
      "      Processing year 2099\n",
      "      Processing year 2100\n",
      "    Processing GCM IPSL-CM5A-LR\n",
      "      Processing year 2081\n",
      "      Processing year 2082\n",
      "      Processing year 2083\n",
      "      Processing year 2084\n",
      "      Processing year 2085\n",
      "      Processing year 2086\n",
      "      Processing year 2087\n",
      "      Processing year 2088\n",
      "      Processing year 2089\n",
      "      Processing year 2090\n",
      "      Processing year 2091\n",
      "      Processing year 2092\n",
      "      Processing year 2093\n",
      "      Processing year 2094\n",
      "      Processing year 2095\n",
      "      Processing year 2096\n",
      "      Processing year 2097\n",
      "      Processing year 2098\n",
      "      Processing year 2099\n",
      "      Processing year 2100\n",
      "    Processing GCM IPSL-CM5A-MR\n",
      "      Processing year 2081\n",
      "      Processing year 2082\n",
      "      Processing year 2083\n",
      "      Processing year 2084\n",
      "      Processing year 2085\n",
      "      Processing year 2086\n",
      "      Processing year 2087\n",
      "      Processing year 2088\n",
      "      Processing year 2089\n",
      "      Processing year 2090\n",
      "      Processing year 2091\n",
      "      Processing year 2092\n",
      "      Processing year 2093\n",
      "      Processing year 2094\n",
      "      Processing year 2095\n",
      "      Processing year 2096\n",
      "      Processing year 2097\n",
      "      Processing year 2098\n",
      "      Processing year 2099\n",
      "      Processing year 2100\n",
      "    Processing GCM IPSL-CM5B-LR\n",
      "      Processing year 2081\n",
      "      Processing year 2082\n",
      "      Processing year 2083\n",
      "      Processing year 2084\n",
      "      Processing year 2085\n",
      "      Processing year 2086\n",
      "      Processing year 2087\n",
      "      Processing year 2088\n",
      "      Processing year 2089\n",
      "      Processing year 2090\n",
      "      Processing year 2091\n",
      "      Processing year 2092\n",
      "      Processing year 2093\n",
      "      Processing year 2094\n",
      "      Processing year 2095\n",
      "      Processing year 2096\n",
      "      Processing year 2097\n",
      "      Processing year 2098\n",
      "      Processing year 2099\n",
      "      Processing year 2100\n",
      "    Processing GCM MIROC-ESM\n",
      "      Processing year 2081\n",
      "      Processing year 2082\n",
      "      Processing year 2083\n",
      "      Processing year 2084\n",
      "      Processing year 2085\n",
      "      Processing year 2086\n",
      "      Processing year 2087\n",
      "      Processing year 2088\n",
      "      Processing year 2089\n",
      "      Processing year 2090\n",
      "      Processing year 2091\n",
      "      Processing year 2092\n",
      "      Processing year 2093\n",
      "      Processing year 2094\n",
      "      Processing year 2095\n",
      "      Processing year 2096\n",
      "      Processing year 2097\n",
      "      Processing year 2098\n",
      "      Processing year 2099\n",
      "      Processing year 2100\n",
      "    Processing GCM MIROC-ESM-CHEM\n",
      "      Processing year 2081\n",
      "      Processing year 2082\n",
      "      Processing year 2083\n",
      "      Processing year 2084\n",
      "      Processing year 2085\n",
      "      Processing year 2086\n",
      "      Processing year 2087\n",
      "      Processing year 2088\n",
      "      Processing year 2089\n",
      "      Processing year 2090\n",
      "      Processing year 2091\n",
      "      Processing year 2092\n",
      "      Processing year 2093\n",
      "      Processing year 2094\n",
      "      Processing year 2095\n",
      "      Processing year 2096\n",
      "      Processing year 2097\n",
      "      Processing year 2098\n",
      "      Processing year 2099\n",
      "      Processing year 2100\n",
      "    Processing GCM MIROC5\n",
      "      Processing year 2081\n",
      "      Processing year 2082\n",
      "      Processing year 2083\n",
      "      Processing year 2084\n",
      "      Processing year 2085\n",
      "      Processing year 2086\n",
      "      Processing year 2087\n",
      "      Processing year 2088\n",
      "      Processing year 2089\n",
      "      Processing year 2090\n",
      "      Processing year 2091\n",
      "      Processing year 2092\n",
      "      Processing year 2093\n",
      "      Processing year 2094\n",
      "      Processing year 2095\n",
      "      Processing year 2096\n",
      "      Processing year 2097\n",
      "      Processing year 2098\n",
      "      Processing year 2099\n",
      "      Processing year 2100\n",
      "    Processing GCM MPI-ESM-LR\n",
      "      Processing year 2081\n",
      "      Processing year 2082\n",
      "      Processing year 2083\n",
      "      Processing year 2084\n",
      "      Processing year 2085\n",
      "      Processing year 2086\n",
      "      Processing year 2087\n",
      "      Processing year 2088\n",
      "      Processing year 2089\n",
      "      Processing year 2090\n",
      "      Processing year 2091\n",
      "      Processing year 2092\n",
      "      Processing year 2093\n",
      "      Processing year 2094\n",
      "      Processing year 2095\n",
      "      Processing year 2096\n",
      "      Processing year 2097\n",
      "      Processing year 2098\n",
      "      Processing year 2099\n",
      "      Processing year 2100\n",
      "    Processing GCM MPI-ESM-MR\n",
      "      Processing year 2081\n",
      "      Processing year 2082\n",
      "      Processing year 2083\n",
      "      Processing year 2084\n",
      "      Processing year 2085\n",
      "      Processing year 2086\n",
      "      Processing year 2087\n",
      "      Processing year 2088\n",
      "      Processing year 2089\n",
      "      Processing year 2090\n",
      "      Processing year 2091\n",
      "      Processing year 2092\n",
      "      Processing year 2093\n",
      "      Processing year 2094\n",
      "      Processing year 2095\n",
      "      Processing year 2096\n",
      "      Processing year 2097\n",
      "      Processing year 2098\n",
      "      Processing year 2099\n",
      "      Processing year 2100\n",
      "    Processing GCM MRI-CGCM3\n",
      "      Processing year 2081\n",
      "      Processing year 2082\n",
      "      Processing year 2083\n",
      "      Processing year 2084\n",
      "      Processing year 2085\n",
      "      Processing year 2086\n",
      "      Processing year 2087\n",
      "      Processing year 2088\n",
      "      Processing year 2089\n",
      "      Processing year 2090\n",
      "      Processing year 2091\n",
      "      Processing year 2092\n",
      "      Processing year 2093\n",
      "      Processing year 2094\n",
      "      Processing year 2095\n",
      "      Processing year 2096\n",
      "      Processing year 2097\n",
      "      Processing year 2098\n",
      "      Processing year 2099\n",
      "      Processing year 2100\n",
      "    Processing GCM NorESM1-M\n",
      "      Processing year 2081\n",
      "      Processing year 2082\n",
      "      Processing year 2083\n",
      "      Processing year 2084\n",
      "      Processing year 2085\n",
      "      Processing year 2086\n",
      "      Processing year 2087\n",
      "      Processing year 2088\n",
      "      Processing year 2089\n",
      "      Processing year 2090\n",
      "      Processing year 2091\n",
      "      Processing year 2092\n",
      "      Processing year 2093\n",
      "      Processing year 2094\n",
      "      Processing year 2095\n",
      "      Processing year 2096\n",
      "      Processing year 2097\n",
      "      Processing year 2098\n",
      "      Processing year 2099\n",
      "      Processing year 2100\n",
      "Writing time. Year = 2090\n",
      "Wrote: prcp_rcp45_2090_20yr_monthly.nc\n",
      "Cleared some memory up...restarting loop\n",
      "Processing climate scenario rcp85\n",
      "  Processing base year 2030\n",
      "    Processing GCM ACCESS1-0\n",
      "      Processing year 2021\n",
      "      Processing year 2022\n",
      "      Processing year 2023\n",
      "      Processing year 2024\n",
      "      Processing year 2025\n",
      "      Processing year 2026\n",
      "      Processing year 2027\n",
      "      Processing year 2028\n",
      "      Processing year 2029\n",
      "      Processing year 2030\n",
      "      Processing year 2031\n",
      "      Processing year 2032\n",
      "      Processing year 2033\n",
      "      Processing year 2034\n",
      "      Processing year 2035\n",
      "      Processing year 2036\n",
      "      Processing year 2037\n",
      "      Processing year 2038\n",
      "      Processing year 2039\n",
      "      Processing year 2040\n",
      "    Processing GCM ACCESS1-3\n",
      "      Processing year 2021\n",
      "      Processing year 2022\n",
      "      Processing year 2023\n",
      "      Processing year 2024\n",
      "      Processing year 2025\n",
      "      Processing year 2026\n",
      "      Processing year 2027\n",
      "      Processing year 2028\n",
      "      Processing year 2029\n",
      "      Processing year 2030\n",
      "      Processing year 2031\n",
      "      Processing year 2032\n",
      "      Processing year 2033\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "      Processing year 2034\n",
      "      Processing year 2035\n",
      "      Processing year 2036\n",
      "      Processing year 2037\n",
      "      Processing year 2038\n",
      "      Processing year 2039\n",
      "      Processing year 2040\n",
      "    Processing GCM CanESM2\n",
      "      Processing year 2021\n",
      "      Processing year 2022\n",
      "      Processing year 2023\n",
      "      Processing year 2024\n",
      "      Processing year 2025\n",
      "      Processing year 2026\n",
      "      Processing year 2027\n",
      "      Processing year 2028\n",
      "      Processing year 2029\n",
      "      Processing year 2030\n",
      "      Processing year 2031\n",
      "      Processing year 2032\n",
      "      Processing year 2033\n",
      "      Processing year 2034\n",
      "      Processing year 2035\n",
      "      Processing year 2036\n",
      "      Processing year 2037\n",
      "      Processing year 2038\n",
      "      Processing year 2039\n",
      "      Processing year 2040\n",
      "    Processing GCM CMCC-CESM\n",
      "      Processing year 2021\n",
      "      Processing year 2022\n",
      "      Processing year 2023\n",
      "      Processing year 2024\n",
      "      Processing year 2025\n",
      "      Processing year 2026\n",
      "      Processing year 2027\n",
      "      Processing year 2028\n",
      "      Processing year 2029\n",
      "      Processing year 2030\n",
      "      Processing year 2031\n",
      "      Processing year 2032\n",
      "      Processing year 2033\n",
      "      Processing year 2034\n",
      "      Processing year 2035\n",
      "      Processing year 2036\n",
      "      Processing year 2037\n",
      "      Processing year 2038\n",
      "      Processing year 2039\n",
      "      Processing year 2040\n",
      "    Processing GCM CMCC-CM\n",
      "      Processing year 2021\n",
      "      Processing year 2022\n",
      "      Processing year 2023\n",
      "      Processing year 2024\n",
      "      Processing year 2025\n",
      "      Processing year 2026\n",
      "      Processing year 2027\n",
      "      Processing year 2028\n",
      "      Processing year 2029\n",
      "      Processing year 2030\n",
      "      Processing year 2031\n",
      "      Processing year 2032\n",
      "      Processing year 2033\n",
      "      Processing year 2034\n",
      "      Processing year 2035\n",
      "      Processing year 2036\n",
      "      Processing year 2037\n",
      "      Processing year 2038\n",
      "      Processing year 2039\n",
      "      Processing year 2040\n",
      "    Processing GCM CMCC-CMS\n",
      "      Processing year 2021\n",
      "      Processing year 2022\n",
      "      Processing year 2023\n",
      "      Processing year 2024\n",
      "      Processing year 2025\n",
      "      Processing year 2026\n",
      "      Processing year 2027\n",
      "      Processing year 2028\n",
      "      Processing year 2029\n",
      "      Processing year 2030\n",
      "      Processing year 2031\n",
      "      Processing year 2032\n",
      "      Processing year 2033\n",
      "      Processing year 2034\n",
      "      Processing year 2035\n",
      "      Processing year 2036\n",
      "      Processing year 2037\n",
      "      Processing year 2038\n",
      "      Processing year 2039\n",
      "      Processing year 2040\n",
      "    Processing GCM CNRM-CM5\n",
      "      Processing year 2021\n",
      "      Processing year 2022\n",
      "      Processing year 2023\n",
      "      Processing year 2024\n",
      "      Processing year 2025\n",
      "      Processing year 2026\n",
      "      Processing year 2027\n",
      "      Processing year 2028\n",
      "      Processing year 2029\n",
      "      Processing year 2030\n",
      "      Processing year 2031\n",
      "      Processing year 2032\n",
      "      Processing year 2033\n",
      "      Processing year 2034\n",
      "      Processing year 2035\n",
      "      Processing year 2036\n",
      "      Processing year 2037\n",
      "      Processing year 2038\n",
      "      Processing year 2039\n",
      "      Processing year 2040\n",
      "    Processing GCM CSIRO-Mk3-6-0\n",
      "      Processing year 2021\n",
      "      Processing year 2022\n",
      "      Processing year 2023\n",
      "      Processing year 2024\n",
      "      Processing year 2025\n",
      "      Processing year 2026\n",
      "      Processing year 2027\n",
      "      Processing year 2028\n",
      "      Processing year 2029\n",
      "      Processing year 2030\n",
      "      Processing year 2031\n",
      "      Processing year 2032\n",
      "      Processing year 2033\n",
      "      Processing year 2034\n",
      "      Processing year 2035\n",
      "      Processing year 2036\n",
      "      Processing year 2037\n",
      "      Processing year 2038\n",
      "      Processing year 2039\n",
      "      Processing year 2040\n",
      "    Processing GCM GFDL-CM3\n",
      "      Processing year 2021\n",
      "      Processing year 2022\n",
      "      Processing year 2023\n",
      "      Processing year 2024\n",
      "      Processing year 2025\n",
      "      Processing year 2026\n",
      "      Processing year 2027\n",
      "      Processing year 2028\n",
      "      Processing year 2029\n",
      "      Processing year 2030\n",
      "      Processing year 2031\n",
      "      Processing year 2032\n",
      "      Processing year 2033\n",
      "      Processing year 2034\n",
      "      Processing year 2035\n",
      "      Processing year 2036\n",
      "      Processing year 2037\n",
      "      Processing year 2038\n",
      "      Processing year 2039\n",
      "      Processing year 2040\n",
      "    Processing GCM GFDL-ESM2G\n",
      "      Processing year 2021\n",
      "      Processing year 2022\n",
      "      Processing year 2023\n",
      "      Processing year 2024\n",
      "      Processing year 2025\n",
      "      Processing year 2026\n",
      "      Processing year 2027\n",
      "      Processing year 2028\n",
      "      Processing year 2029\n",
      "      Processing year 2030\n",
      "      Processing year 2031\n",
      "      Processing year 2032\n",
      "      Processing year 2033\n",
      "      Processing year 2034\n",
      "      Processing year 2035\n",
      "      Processing year 2036\n",
      "      Processing year 2037\n",
      "      Processing year 2038\n",
      "      Processing year 2039\n",
      "      Processing year 2040\n",
      "    Processing GCM GFDL-ESM2M\n",
      "      Processing year 2021\n",
      "      Processing year 2022\n",
      "      Processing year 2023\n",
      "      Processing year 2024\n",
      "      Processing year 2025\n",
      "      Processing year 2026\n",
      "      Processing year 2027\n",
      "      Processing year 2028\n",
      "      Processing year 2029\n",
      "      Processing year 2030\n",
      "      Processing year 2031\n",
      "      Processing year 2032\n",
      "      Processing year 2033\n",
      "      Processing year 2034\n",
      "      Processing year 2035\n",
      "      Processing year 2036\n",
      "      Processing year 2037\n",
      "      Processing year 2038\n",
      "      Processing year 2039\n",
      "      Processing year 2040\n",
      "    Processing GCM HadGEM2-CC\n",
      "      Processing year 2021\n",
      "      Processing year 2022\n",
      "      Processing year 2023\n",
      "      Processing year 2024\n",
      "      Processing year 2025\n",
      "      Processing year 2026\n",
      "      Processing year 2027\n",
      "      Processing year 2028\n",
      "      Processing year 2029\n",
      "      Processing year 2030\n",
      "      Processing year 2031\n",
      "      Processing year 2032\n",
      "      Processing year 2033\n",
      "      Processing year 2034\n",
      "      Processing year 2035\n",
      "      Processing year 2036\n",
      "      Processing year 2037\n",
      "      Processing year 2038\n",
      "      Processing year 2039\n",
      "      Processing year 2040\n",
      "    Processing GCM inmcm4\n",
      "      Processing year 2021\n",
      "      Processing year 2022\n",
      "      Processing year 2023\n",
      "      Processing year 2024\n",
      "      Processing year 2025\n",
      "      Processing year 2026\n",
      "      Processing year 2027\n",
      "      Processing year 2028\n",
      "      Processing year 2029\n",
      "      Processing year 2030\n",
      "      Processing year 2031\n",
      "      Processing year 2032\n",
      "      Processing year 2033\n",
      "      Processing year 2034\n",
      "      Processing year 2035\n",
      "      Processing year 2036\n",
      "      Processing year 2037\n",
      "      Processing year 2038\n",
      "      Processing year 2039\n",
      "      Processing year 2040\n",
      "    Processing GCM IPSL-CM5A-LR\n",
      "      Processing year 2021\n",
      "      Processing year 2022\n",
      "      Processing year 2023\n",
      "      Processing year 2024\n",
      "      Processing year 2025\n",
      "      Processing year 2026\n",
      "      Processing year 2027\n",
      "      Processing year 2028\n",
      "      Processing year 2029\n",
      "      Processing year 2030\n",
      "      Processing year 2031\n",
      "      Processing year 2032\n",
      "      Processing year 2033\n",
      "      Processing year 2034\n",
      "      Processing year 2035\n",
      "      Processing year 2036\n",
      "      Processing year 2037\n",
      "      Processing year 2038\n",
      "      Processing year 2039\n",
      "      Processing year 2040\n",
      "    Processing GCM IPSL-CM5A-MR\n",
      "      Processing year 2021\n",
      "      Processing year 2022\n",
      "      Processing year 2023\n",
      "      Processing year 2024\n",
      "      Processing year 2025\n",
      "      Processing year 2026\n",
      "      Processing year 2027\n",
      "      Processing year 2028\n",
      "      Processing year 2029\n",
      "      Processing year 2030\n",
      "      Processing year 2031\n",
      "      Processing year 2032\n",
      "      Processing year 2033\n",
      "      Processing year 2034\n",
      "      Processing year 2035\n",
      "      Processing year 2036\n",
      "      Processing year 2037\n",
      "      Processing year 2038\n",
      "      Processing year 2039\n",
      "      Processing year 2040\n",
      "    Processing GCM IPSL-CM5B-LR\n",
      "      Processing year 2021\n",
      "      Processing year 2022\n",
      "      Processing year 2023\n",
      "      Processing year 2024\n",
      "      Processing year 2025\n",
      "      Processing year 2026\n",
      "      Processing year 2027\n",
      "      Processing year 2028\n",
      "      Processing year 2029\n",
      "      Processing year 2030\n",
      "      Processing year 2031\n",
      "      Processing year 2032\n",
      "      Processing year 2033\n",
      "      Processing year 2034\n",
      "      Processing year 2035\n",
      "      Processing year 2036\n",
      "      Processing year 2037\n",
      "      Processing year 2038\n",
      "      Processing year 2039\n",
      "      Processing year 2040\n",
      "    Processing GCM MIROC-ESM\n",
      "      Processing year 2021\n",
      "      Processing year 2022\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "      Processing year 2023\n",
      "      Processing year 2024\n",
      "      Processing year 2025\n",
      "      Processing year 2026\n",
      "      Processing year 2027\n",
      "      Processing year 2028\n",
      "      Processing year 2029\n",
      "      Processing year 2030\n",
      "      Processing year 2031\n",
      "      Processing year 2032\n",
      "      Processing year 2033\n",
      "      Processing year 2034\n",
      "      Processing year 2035\n",
      "      Processing year 2036\n",
      "      Processing year 2037\n",
      "      Processing year 2038\n",
      "      Processing year 2039\n",
      "      Processing year 2040\n",
      "    Processing GCM MIROC-ESM-CHEM\n",
      "      Processing year 2021\n",
      "      Processing year 2022\n",
      "      Processing year 2023\n",
      "      Processing year 2024\n",
      "      Processing year 2025\n",
      "      Processing year 2026\n",
      "      Processing year 2027\n",
      "      Processing year 2028\n",
      "      Processing year 2029\n",
      "      Processing year 2030\n",
      "      Processing year 2031\n",
      "      Processing year 2032\n",
      "      Processing year 2033\n",
      "      Processing year 2034\n",
      "      Processing year 2035\n",
      "      Processing year 2036\n",
      "      Processing year 2037\n",
      "      Processing year 2038\n",
      "      Processing year 2039\n",
      "      Processing year 2040\n",
      "    Processing GCM MIROC5\n",
      "      Processing year 2021\n",
      "      Processing year 2022\n",
      "      Processing year 2023\n",
      "      Processing year 2024\n",
      "      Processing year 2025\n",
      "      Processing year 2026\n",
      "      Processing year 2027\n",
      "      Processing year 2028\n",
      "      Processing year 2029\n",
      "      Processing year 2030\n",
      "      Processing year 2031\n",
      "      Processing year 2032\n",
      "      Processing year 2033\n",
      "      Processing year 2034\n",
      "      Processing year 2035\n",
      "      Processing year 2036\n",
      "      Processing year 2037\n",
      "      Processing year 2038\n",
      "      Processing year 2039\n",
      "      Processing year 2040\n",
      "    Processing GCM MPI-ESM-LR\n",
      "      Processing year 2021\n",
      "      Processing year 2022\n",
      "      Processing year 2023\n",
      "      Processing year 2024\n",
      "      Processing year 2025\n",
      "      Processing year 2026\n",
      "      Processing year 2027\n",
      "      Processing year 2028\n",
      "      Processing year 2029\n",
      "      Processing year 2030\n",
      "      Processing year 2031\n",
      "      Processing year 2032\n",
      "      Processing year 2033\n",
      "      Processing year 2034\n",
      "      Processing year 2035\n",
      "      Processing year 2036\n",
      "      Processing year 2037\n",
      "      Processing year 2038\n",
      "      Processing year 2039\n",
      "      Processing year 2040\n",
      "    Processing GCM MPI-ESM-MR\n",
      "      Processing year 2021\n",
      "      Processing year 2022\n",
      "      Processing year 2023\n",
      "      Processing year 2024\n",
      "      Processing year 2025\n",
      "      Processing year 2026\n",
      "      Processing year 2027\n",
      "      Processing year 2028\n",
      "      Processing year 2029\n",
      "      Processing year 2030\n",
      "      Processing year 2031\n",
      "      Processing year 2032\n",
      "      Processing year 2033\n",
      "      Processing year 2034\n",
      "      Processing year 2035\n",
      "      Processing year 2036\n",
      "      Processing year 2037\n",
      "      Processing year 2038\n",
      "      Processing year 2039\n",
      "      Processing year 2040\n",
      "    Processing GCM MRI-CGCM3\n",
      "      Processing year 2021\n",
      "      Processing year 2022\n",
      "      Processing year 2023\n",
      "      Processing year 2024\n",
      "      Processing year 2025\n",
      "      Processing year 2026\n",
      "      Processing year 2027\n",
      "      Processing year 2028\n",
      "      Processing year 2029\n",
      "      Processing year 2030\n",
      "      Processing year 2031\n",
      "      Processing year 2032\n",
      "      Processing year 2033\n",
      "      Processing year 2034\n",
      "      Processing year 2035\n",
      "      Processing year 2036\n",
      "      Processing year 2037\n",
      "      Processing year 2038\n",
      "      Processing year 2039\n",
      "      Processing year 2040\n",
      "    Processing GCM MRI-ESM1\n",
      "      Processing year 2021\n",
      "      Processing year 2022\n",
      "      Processing year 2023\n",
      "      Processing year 2024\n",
      "      Processing year 2025\n",
      "      Processing year 2026\n",
      "      Processing year 2027\n",
      "      Processing year 2028\n",
      "      Processing year 2029\n",
      "      Processing year 2030\n",
      "      Processing year 2031\n",
      "      Processing year 2032\n",
      "      Processing year 2033\n",
      "      Processing year 2034\n",
      "      Processing year 2035\n",
      "      Processing year 2036\n",
      "      Processing year 2037\n",
      "      Processing year 2038\n",
      "      Processing year 2039\n",
      "      Processing year 2040\n",
      "    Processing GCM NorESM1-M\n",
      "      Processing year 2021\n",
      "      Processing year 2022\n",
      "      Processing year 2023\n",
      "      Processing year 2024\n",
      "      Processing year 2025\n",
      "      Processing year 2026\n",
      "      Processing year 2027\n",
      "      Processing year 2028\n",
      "      Processing year 2029\n",
      "      Processing year 2030\n",
      "      Processing year 2031\n",
      "      Processing year 2032\n",
      "      Processing year 2033\n",
      "      Processing year 2034\n",
      "      Processing year 2035\n",
      "      Processing year 2036\n",
      "      Processing year 2037\n",
      "      Processing year 2038\n",
      "      Processing year 2039\n",
      "      Processing year 2040\n",
      "Writing time. Year = 2030\n",
      "Wrote: prcp_rcp85_2030_20yr_monthly.nc\n",
      "Cleared some memory up...restarting loop\n",
      "  Processing base year 2050\n",
      "    Processing GCM ACCESS1-0\n",
      "      Processing year 2041\n",
      "      Processing year 2042\n",
      "      Processing year 2043\n",
      "      Processing year 2044\n",
      "      Processing year 2045\n",
      "      Processing year 2046\n",
      "      Processing year 2047\n",
      "      Processing year 2048\n",
      "      Processing year 2049\n",
      "      Processing year 2050\n",
      "      Processing year 2051\n",
      "      Processing year 2052\n",
      "      Processing year 2053\n",
      "      Processing year 2054\n",
      "      Processing year 2055\n",
      "      Processing year 2056\n",
      "      Processing year 2057\n",
      "      Processing year 2058\n",
      "      Processing year 2059\n",
      "      Processing year 2060\n",
      "    Processing GCM ACCESS1-3\n",
      "      Processing year 2041\n",
      "      Processing year 2042\n",
      "      Processing year 2043\n",
      "      Processing year 2044\n",
      "      Processing year 2045\n",
      "      Processing year 2046\n",
      "      Processing year 2047\n",
      "      Processing year 2048\n",
      "      Processing year 2049\n",
      "      Processing year 2050\n",
      "      Processing year 2051\n",
      "      Processing year 2052\n",
      "      Processing year 2053\n",
      "      Processing year 2054\n",
      "      Processing year 2055\n",
      "      Processing year 2056\n",
      "      Processing year 2057\n",
      "      Processing year 2058\n",
      "      Processing year 2059\n",
      "      Processing year 2060\n",
      "    Processing GCM CanESM2\n",
      "      Processing year 2041\n",
      "      Processing year 2042\n",
      "      Processing year 2043\n",
      "      Processing year 2044\n",
      "      Processing year 2045\n",
      "      Processing year 2046\n",
      "      Processing year 2047\n",
      "      Processing year 2048\n",
      "      Processing year 2049\n",
      "      Processing year 2050\n",
      "      Processing year 2051\n",
      "      Processing year 2052\n",
      "      Processing year 2053\n",
      "      Processing year 2054\n",
      "      Processing year 2055\n",
      "      Processing year 2056\n",
      "      Processing year 2057\n",
      "      Processing year 2058\n",
      "      Processing year 2059\n",
      "      Processing year 2060\n",
      "    Processing GCM CMCC-CESM\n",
      "      Processing year 2041\n",
      "      Processing year 2042\n",
      "      Processing year 2043\n",
      "      Processing year 2044\n",
      "      Processing year 2045\n",
      "      Processing year 2046\n",
      "      Processing year 2047\n",
      "      Processing year 2048\n",
      "      Processing year 2049\n",
      "      Processing year 2050\n",
      "      Processing year 2051\n",
      "      Processing year 2052\n",
      "      Processing year 2053\n",
      "      Processing year 2054\n",
      "      Processing year 2055\n",
      "      Processing year 2056\n",
      "      Processing year 2057\n",
      "      Processing year 2058\n",
      "      Processing year 2059\n",
      "      Processing year 2060\n",
      "    Processing GCM CMCC-CM\n",
      "      Processing year 2041\n",
      "      Processing year 2042\n",
      "      Processing year 2043\n",
      "      Processing year 2044\n",
      "      Processing year 2045\n",
      "      Processing year 2046\n",
      "      Processing year 2047\n",
      "      Processing year 2048\n",
      "      Processing year 2049\n",
      "      Processing year 2050\n",
      "      Processing year 2051\n",
      "      Processing year 2052\n",
      "      Processing year 2053\n",
      "      Processing year 2054\n",
      "      Processing year 2055\n",
      "      Processing year 2056\n",
      "      Processing year 2057\n",
      "      Processing year 2058\n",
      "      Processing year 2059\n",
      "      Processing year 2060\n",
      "    Processing GCM CMCC-CMS\n",
      "      Processing year 2041\n",
      "      Processing year 2042\n",
      "      Processing year 2043\n",
      "      Processing year 2044\n",
      "      Processing year 2045\n",
      "      Processing year 2046\n",
      "      Processing year 2047\n",
      "      Processing year 2048\n",
      "      Processing year 2049\n",
      "      Processing year 2050\n",
      "      Processing year 2051\n",
      "      Processing year 2052\n",
      "      Processing year 2053\n",
      "      Processing year 2054\n",
      "      Processing year 2055\n",
      "      Processing year 2056\n",
      "      Processing year 2057\n",
      "      Processing year 2058\n",
      "      Processing year 2059\n",
      "      Processing year 2060\n",
      "    Processing GCM CNRM-CM5\n",
      "      Processing year 2041\n",
      "      Processing year 2042\n",
      "      Processing year 2043\n",
      "      Processing year 2044\n",
      "      Processing year 2045\n",
      "      Processing year 2046\n",
      "      Processing year 2047\n",
      "      Processing year 2048\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "      Processing year 2049\n",
      "      Processing year 2050\n",
      "      Processing year 2051\n",
      "      Processing year 2052\n",
      "      Processing year 2053\n",
      "      Processing year 2054\n",
      "      Processing year 2055\n",
      "      Processing year 2056\n",
      "      Processing year 2057\n",
      "      Processing year 2058\n",
      "      Processing year 2059\n",
      "      Processing year 2060\n",
      "    Processing GCM CSIRO-Mk3-6-0\n",
      "      Processing year 2041\n",
      "      Processing year 2042\n",
      "      Processing year 2043\n",
      "      Processing year 2044\n",
      "      Processing year 2045\n",
      "      Processing year 2046\n",
      "      Processing year 2047\n",
      "      Processing year 2048\n",
      "      Processing year 2049\n",
      "      Processing year 2050\n",
      "      Processing year 2051\n",
      "      Processing year 2052\n",
      "      Processing year 2053\n",
      "      Processing year 2054\n",
      "      Processing year 2055\n",
      "      Processing year 2056\n",
      "      Processing year 2057\n",
      "      Processing year 2058\n",
      "      Processing year 2059\n",
      "      Processing year 2060\n",
      "    Processing GCM GFDL-CM3\n",
      "      Processing year 2041\n",
      "      Processing year 2042\n",
      "      Processing year 2043\n",
      "      Processing year 2044\n",
      "      Processing year 2045\n",
      "      Processing year 2046\n",
      "      Processing year 2047\n",
      "      Processing year 2048\n",
      "      Processing year 2049\n",
      "      Processing year 2050\n",
      "      Processing year 2051\n",
      "      Processing year 2052\n",
      "      Processing year 2053\n",
      "      Processing year 2054\n",
      "      Processing year 2055\n",
      "      Processing year 2056\n",
      "      Processing year 2057\n",
      "      Processing year 2058\n",
      "      Processing year 2059\n",
      "      Processing year 2060\n",
      "    Processing GCM GFDL-ESM2G\n",
      "      Processing year 2041\n",
      "      Processing year 2042\n",
      "      Processing year 2043\n",
      "      Processing year 2044\n",
      "      Processing year 2045\n",
      "      Processing year 2046\n",
      "      Processing year 2047\n",
      "      Processing year 2048\n",
      "      Processing year 2049\n",
      "      Processing year 2050\n",
      "      Processing year 2051\n",
      "      Processing year 2052\n",
      "      Processing year 2053\n",
      "      Processing year 2054\n",
      "      Processing year 2055\n",
      "      Processing year 2056\n",
      "      Processing year 2057\n",
      "      Processing year 2058\n",
      "      Processing year 2059\n",
      "      Processing year 2060\n",
      "    Processing GCM GFDL-ESM2M\n",
      "      Processing year 2041\n",
      "      Processing year 2042\n",
      "      Processing year 2043\n",
      "      Processing year 2044\n",
      "      Processing year 2045\n",
      "      Processing year 2046\n",
      "      Processing year 2047\n",
      "      Processing year 2048\n",
      "      Processing year 2049\n",
      "      Processing year 2050\n",
      "      Processing year 2051\n",
      "      Processing year 2052\n",
      "      Processing year 2053\n",
      "      Processing year 2054\n",
      "      Processing year 2055\n",
      "      Processing year 2056\n",
      "      Processing year 2057\n",
      "      Processing year 2058\n",
      "      Processing year 2059\n",
      "      Processing year 2060\n",
      "    Processing GCM HadGEM2-CC\n",
      "      Processing year 2041\n",
      "      Processing year 2042\n",
      "      Processing year 2043\n",
      "      Processing year 2044\n",
      "      Processing year 2045\n",
      "      Processing year 2046\n",
      "      Processing year 2047\n",
      "      Processing year 2048\n",
      "      Processing year 2049\n",
      "      Processing year 2050\n",
      "      Processing year 2051\n",
      "      Processing year 2052\n",
      "      Processing year 2053\n",
      "      Processing year 2054\n",
      "      Processing year 2055\n",
      "      Processing year 2056\n",
      "      Processing year 2057\n",
      "      Processing year 2058\n",
      "      Processing year 2059\n",
      "      Processing year 2060\n",
      "    Processing GCM inmcm4\n",
      "      Processing year 2041\n",
      "      Processing year 2042\n",
      "      Processing year 2043\n",
      "      Processing year 2044\n",
      "      Processing year 2045\n",
      "      Processing year 2046\n",
      "      Processing year 2047\n",
      "      Processing year 2048\n",
      "      Processing year 2049\n",
      "      Processing year 2050\n",
      "      Processing year 2051\n",
      "      Processing year 2052\n",
      "      Processing year 2053\n",
      "      Processing year 2054\n",
      "      Processing year 2055\n",
      "      Processing year 2056\n",
      "      Processing year 2057\n",
      "      Processing year 2058\n",
      "      Processing year 2059\n",
      "      Processing year 2060\n",
      "    Processing GCM IPSL-CM5A-LR\n",
      "      Processing year 2041\n",
      "      Processing year 2042\n",
      "      Processing year 2043\n",
      "      Processing year 2044\n",
      "      Processing year 2045\n",
      "      Processing year 2046\n",
      "      Processing year 2047\n",
      "      Processing year 2048\n",
      "      Processing year 2049\n",
      "      Processing year 2050\n",
      "      Processing year 2051\n",
      "      Processing year 2052\n",
      "      Processing year 2053\n",
      "      Processing year 2054\n",
      "      Processing year 2055\n",
      "      Processing year 2056\n",
      "      Processing year 2057\n",
      "      Processing year 2058\n",
      "      Processing year 2059\n",
      "      Processing year 2060\n",
      "    Processing GCM IPSL-CM5A-MR\n",
      "      Processing year 2041\n",
      "      Processing year 2042\n",
      "      Processing year 2043\n",
      "      Processing year 2044\n",
      "      Processing year 2045\n",
      "      Processing year 2046\n",
      "      Processing year 2047\n",
      "      Processing year 2048\n",
      "      Processing year 2049\n",
      "      Processing year 2050\n",
      "      Processing year 2051\n",
      "      Processing year 2052\n",
      "      Processing year 2053\n",
      "      Processing year 2054\n",
      "      Processing year 2055\n",
      "      Processing year 2056\n",
      "      Processing year 2057\n",
      "      Processing year 2058\n",
      "      Processing year 2059\n",
      "      Processing year 2060\n",
      "    Processing GCM IPSL-CM5B-LR\n",
      "      Processing year 2041\n",
      "      Processing year 2042\n",
      "      Processing year 2043\n",
      "      Processing year 2044\n",
      "      Processing year 2045\n",
      "      Processing year 2046\n",
      "      Processing year 2047\n",
      "      Processing year 2048\n",
      "      Processing year 2049\n",
      "      Processing year 2050\n",
      "      Processing year 2051\n",
      "      Processing year 2052\n",
      "      Processing year 2053\n",
      "      Processing year 2054\n",
      "      Processing year 2055\n",
      "      Processing year 2056\n",
      "      Processing year 2057\n",
      "      Processing year 2058\n",
      "      Processing year 2059\n",
      "      Processing year 2060\n",
      "    Processing GCM MIROC-ESM\n",
      "      Processing year 2041\n",
      "      Processing year 2042\n",
      "      Processing year 2043\n",
      "      Processing year 2044\n",
      "      Processing year 2045\n",
      "      Processing year 2046\n",
      "      Processing year 2047\n",
      "      Processing year 2048\n",
      "      Processing year 2049\n",
      "      Processing year 2050\n",
      "      Processing year 2051\n",
      "      Processing year 2052\n",
      "      Processing year 2053\n",
      "      Processing year 2054\n",
      "      Processing year 2055\n",
      "      Processing year 2056\n",
      "      Processing year 2057\n",
      "      Processing year 2058\n",
      "      Processing year 2059\n",
      "      Processing year 2060\n",
      "    Processing GCM MIROC-ESM-CHEM\n",
      "      Processing year 2041\n",
      "      Processing year 2042\n",
      "      Processing year 2043\n",
      "      Processing year 2044\n",
      "      Processing year 2045\n",
      "      Processing year 2046\n",
      "      Processing year 2047\n",
      "      Processing year 2048\n",
      "      Processing year 2049\n",
      "      Processing year 2050\n",
      "      Processing year 2051\n",
      "      Processing year 2052\n",
      "      Processing year 2053\n",
      "      Processing year 2054\n",
      "      Processing year 2055\n",
      "      Processing year 2056\n",
      "      Processing year 2057\n",
      "      Processing year 2058\n",
      "      Processing year 2059\n",
      "      Processing year 2060\n",
      "    Processing GCM MIROC5\n",
      "      Processing year 2041\n",
      "      Processing year 2042\n",
      "      Processing year 2043\n",
      "      Processing year 2044\n",
      "      Processing year 2045\n",
      "      Processing year 2046\n",
      "      Processing year 2047\n",
      "      Processing year 2048\n",
      "      Processing year 2049\n",
      "      Processing year 2050\n",
      "      Processing year 2051\n",
      "      Processing year 2052\n",
      "      Processing year 2053\n",
      "      Processing year 2054\n",
      "      Processing year 2055\n",
      "      Processing year 2056\n",
      "      Processing year 2057\n",
      "      Processing year 2058\n",
      "      Processing year 2059\n",
      "      Processing year 2060\n",
      "    Processing GCM MPI-ESM-LR\n",
      "      Processing year 2041\n",
      "      Processing year 2042\n",
      "      Processing year 2043\n",
      "      Processing year 2044\n",
      "      Processing year 2045\n",
      "      Processing year 2046\n",
      "      Processing year 2047\n",
      "      Processing year 2048\n",
      "      Processing year 2049\n",
      "      Processing year 2050\n",
      "      Processing year 2051\n",
      "      Processing year 2052\n",
      "      Processing year 2053\n",
      "      Processing year 2054\n",
      "      Processing year 2055\n",
      "      Processing year 2056\n",
      "      Processing year 2057\n",
      "      Processing year 2058\n",
      "      Processing year 2059\n",
      "      Processing year 2060\n",
      "    Processing GCM MPI-ESM-MR\n",
      "      Processing year 2041\n",
      "      Processing year 2042\n",
      "      Processing year 2043\n",
      "      Processing year 2044\n",
      "      Processing year 2045\n",
      "      Processing year 2046\n",
      "      Processing year 2047\n",
      "      Processing year 2048\n",
      "      Processing year 2049\n",
      "      Processing year 2050\n",
      "      Processing year 2051\n",
      "      Processing year 2052\n",
      "      Processing year 2053\n",
      "      Processing year 2054\n",
      "      Processing year 2055\n",
      "      Processing year 2056\n",
      "      Processing year 2057\n",
      "      Processing year 2058\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "      Processing year 2059\n",
      "      Processing year 2060\n",
      "    Processing GCM MRI-CGCM3\n",
      "      Processing year 2041\n",
      "      Processing year 2042\n",
      "      Processing year 2043\n",
      "      Processing year 2044\n",
      "      Processing year 2045\n",
      "      Processing year 2046\n",
      "      Processing year 2047\n",
      "      Processing year 2048\n",
      "      Processing year 2049\n",
      "      Processing year 2050\n",
      "      Processing year 2051\n",
      "      Processing year 2052\n",
      "      Processing year 2053\n",
      "      Processing year 2054\n",
      "      Processing year 2055\n",
      "      Processing year 2056\n",
      "      Processing year 2057\n",
      "      Processing year 2058\n",
      "      Processing year 2059\n",
      "      Processing year 2060\n",
      "    Processing GCM MRI-ESM1\n",
      "      Processing year 2041\n",
      "      Processing year 2042\n",
      "      Processing year 2043\n",
      "      Processing year 2044\n",
      "      Processing year 2045\n",
      "      Processing year 2046\n",
      "      Processing year 2047\n",
      "      Processing year 2048\n",
      "      Processing year 2049\n",
      "      Processing year 2050\n",
      "      Processing year 2051\n",
      "      Processing year 2052\n",
      "      Processing year 2053\n",
      "      Processing year 2054\n",
      "      Processing year 2055\n",
      "      Processing year 2056\n",
      "      Processing year 2057\n",
      "      Processing year 2058\n",
      "      Processing year 2059\n",
      "      Processing year 2060\n",
      "    Processing GCM NorESM1-M\n",
      "      Processing year 2041\n",
      "      Processing year 2042\n",
      "      Processing year 2043\n",
      "      Processing year 2044\n",
      "      Processing year 2045\n",
      "      Processing year 2046\n",
      "      Processing year 2047\n",
      "      Processing year 2048\n",
      "      Processing year 2049\n",
      "      Processing year 2050\n",
      "      Processing year 2051\n",
      "      Processing year 2052\n",
      "      Processing year 2053\n",
      "      Processing year 2054\n",
      "      Processing year 2055\n",
      "      Processing year 2056\n",
      "      Processing year 2057\n",
      "      Processing year 2058\n",
      "      Processing year 2059\n",
      "      Processing year 2060\n",
      "Writing time. Year = 2050\n",
      "Wrote: prcp_rcp85_2050_20yr_monthly.nc\n",
      "Cleared some memory up...restarting loop\n",
      "  Processing base year 2070\n",
      "    Processing GCM ACCESS1-0\n",
      "      Processing year 2061\n",
      "      Processing year 2062\n",
      "      Processing year 2063\n",
      "      Processing year 2064\n",
      "      Processing year 2065\n",
      "      Processing year 2066\n",
      "      Processing year 2067\n",
      "      Processing year 2068\n",
      "      Processing year 2069\n",
      "      Processing year 2070\n",
      "      Processing year 2071\n",
      "      Processing year 2072\n",
      "      Processing year 2073\n",
      "      Processing year 2074\n",
      "      Processing year 2075\n",
      "      Processing year 2076\n",
      "      Processing year 2077\n",
      "      Processing year 2078\n",
      "      Processing year 2079\n",
      "      Processing year 2080\n",
      "    Processing GCM ACCESS1-3\n",
      "      Processing year 2061\n",
      "      Processing year 2062\n",
      "      Processing year 2063\n",
      "      Processing year 2064\n",
      "      Processing year 2065\n",
      "      Processing year 2066\n",
      "      Processing year 2067\n",
      "      Processing year 2068\n",
      "      Processing year 2069\n",
      "      Processing year 2070\n",
      "      Processing year 2071\n",
      "      Processing year 2072\n",
      "      Processing year 2073\n",
      "      Processing year 2074\n",
      "      Processing year 2075\n",
      "      Processing year 2076\n",
      "      Processing year 2077\n",
      "      Processing year 2078\n",
      "      Processing year 2079\n",
      "      Processing year 2080\n",
      "    Processing GCM CanESM2\n",
      "      Processing year 2061\n",
      "      Processing year 2062\n",
      "      Processing year 2063\n",
      "      Processing year 2064\n",
      "      Processing year 2065\n",
      "      Processing year 2066\n",
      "      Processing year 2067\n",
      "      Processing year 2068\n",
      "      Processing year 2069\n",
      "      Processing year 2070\n",
      "      Processing year 2071\n",
      "      Processing year 2072\n",
      "      Processing year 2073\n",
      "      Processing year 2074\n",
      "      Processing year 2075\n",
      "      Processing year 2076\n",
      "      Processing year 2077\n",
      "      Processing year 2078\n",
      "      Processing year 2079\n",
      "      Processing year 2080\n",
      "    Processing GCM CMCC-CESM\n",
      "      Processing year 2061\n",
      "      Processing year 2062\n",
      "      Processing year 2063\n",
      "      Processing year 2064\n",
      "      Processing year 2065\n",
      "      Processing year 2066\n",
      "      Processing year 2067\n",
      "      Processing year 2068\n",
      "      Processing year 2069\n",
      "      Processing year 2070\n",
      "      Processing year 2071\n",
      "      Processing year 2072\n",
      "      Processing year 2073\n",
      "      Processing year 2074\n",
      "      Processing year 2075\n",
      "      Processing year 2076\n",
      "      Processing year 2077\n",
      "      Processing year 2078\n",
      "      Processing year 2079\n",
      "      Processing year 2080\n",
      "    Processing GCM CMCC-CM\n",
      "      Processing year 2061\n",
      "      Processing year 2062\n",
      "      Processing year 2063\n",
      "      Processing year 2064\n",
      "      Processing year 2065\n",
      "      Processing year 2066\n",
      "      Processing year 2067\n",
      "      Processing year 2068\n",
      "      Processing year 2069\n",
      "      Processing year 2070\n",
      "      Processing year 2071\n",
      "      Processing year 2072\n",
      "      Processing year 2073\n",
      "      Processing year 2074\n",
      "      Processing year 2075\n",
      "      Processing year 2076\n",
      "      Processing year 2077\n",
      "      Processing year 2078\n",
      "      Processing year 2079\n",
      "      Processing year 2080\n",
      "    Processing GCM CMCC-CMS\n",
      "      Processing year 2061\n",
      "      Processing year 2062\n",
      "      Processing year 2063\n",
      "      Processing year 2064\n",
      "      Processing year 2065\n",
      "      Processing year 2066\n",
      "      Processing year 2067\n",
      "      Processing year 2068\n",
      "      Processing year 2069\n",
      "      Processing year 2070\n",
      "      Processing year 2071\n",
      "      Processing year 2072\n",
      "      Processing year 2073\n",
      "      Processing year 2074\n",
      "      Processing year 2075\n",
      "      Processing year 2076\n",
      "      Processing year 2077\n",
      "      Processing year 2078\n",
      "      Processing year 2079\n",
      "      Processing year 2080\n",
      "    Processing GCM CNRM-CM5\n",
      "      Processing year 2061\n",
      "      Processing year 2062\n",
      "      Processing year 2063\n",
      "      Processing year 2064\n",
      "      Processing year 2065\n",
      "      Processing year 2066\n",
      "      Processing year 2067\n",
      "      Processing year 2068\n",
      "      Processing year 2069\n",
      "      Processing year 2070\n",
      "      Processing year 2071\n",
      "      Processing year 2072\n",
      "      Processing year 2073\n",
      "      Processing year 2074\n",
      "      Processing year 2075\n",
      "      Processing year 2076\n",
      "      Processing year 2077\n",
      "      Processing year 2078\n",
      "      Processing year 2079\n",
      "      Processing year 2080\n",
      "    Processing GCM CSIRO-Mk3-6-0\n",
      "      Processing year 2061\n",
      "      Processing year 2062\n",
      "      Processing year 2063\n",
      "      Processing year 2064\n",
      "      Processing year 2065\n",
      "      Processing year 2066\n",
      "      Processing year 2067\n",
      "      Processing year 2068\n",
      "      Processing year 2069\n",
      "      Processing year 2070\n",
      "      Processing year 2071\n",
      "      Processing year 2072\n",
      "      Processing year 2073\n",
      "      Processing year 2074\n",
      "      Processing year 2075\n",
      "      Processing year 2076\n",
      "      Processing year 2077\n",
      "      Processing year 2078\n",
      "      Processing year 2079\n",
      "      Processing year 2080\n",
      "    Processing GCM GFDL-CM3\n",
      "      Processing year 2061\n",
      "      Processing year 2062\n",
      "      Processing year 2063\n",
      "      Processing year 2064\n",
      "      Processing year 2065\n",
      "      Processing year 2066\n",
      "      Processing year 2067\n",
      "      Processing year 2068\n",
      "      Processing year 2069\n",
      "      Processing year 2070\n",
      "      Processing year 2071\n",
      "      Processing year 2072\n",
      "      Processing year 2073\n",
      "      Processing year 2074\n",
      "      Processing year 2075\n",
      "      Processing year 2076\n",
      "      Processing year 2077\n",
      "      Processing year 2078\n",
      "      Processing year 2079\n",
      "      Processing year 2080\n",
      "    Processing GCM GFDL-ESM2G\n",
      "      Processing year 2061\n",
      "      Processing year 2062\n",
      "      Processing year 2063\n",
      "      Processing year 2064\n",
      "      Processing year 2065\n",
      "      Processing year 2066\n",
      "      Processing year 2067\n",
      "      Processing year 2068\n",
      "      Processing year 2069\n",
      "      Processing year 2070\n",
      "      Processing year 2071\n",
      "      Processing year 2072\n",
      "      Processing year 2073\n",
      "      Processing year 2074\n",
      "      Processing year 2075\n",
      "      Processing year 2076\n",
      "      Processing year 2077\n",
      "      Processing year 2078\n",
      "      Processing year 2079\n",
      "      Processing year 2080\n",
      "    Processing GCM GFDL-ESM2M\n",
      "      Processing year 2061\n",
      "      Processing year 2062\n",
      "      Processing year 2063\n",
      "      Processing year 2064\n",
      "      Processing year 2065\n",
      "      Processing year 2066\n",
      "      Processing year 2067\n",
      "      Processing year 2068\n",
      "      Processing year 2069\n",
      "      Processing year 2070\n",
      "      Processing year 2071\n",
      "      Processing year 2072\n",
      "      Processing year 2073\n",
      "      Processing year 2074\n",
      "      Processing year 2075\n",
      "      Processing year 2076\n",
      "      Processing year 2077\n",
      "      Processing year 2078\n",
      "      Processing year 2079\n",
      "      Processing year 2080\n",
      "    Processing GCM HadGEM2-CC\n",
      "      Processing year 2061\n",
      "      Processing year 2062\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "      Processing year 2063\n",
      "      Processing year 2064\n",
      "      Processing year 2065\n",
      "      Processing year 2066\n",
      "      Processing year 2067\n",
      "      Processing year 2068\n",
      "      Processing year 2069\n",
      "      Processing year 2070\n",
      "      Processing year 2071\n",
      "      Processing year 2072\n",
      "      Processing year 2073\n",
      "      Processing year 2074\n",
      "      Processing year 2075\n",
      "      Processing year 2076\n",
      "      Processing year 2077\n",
      "      Processing year 2078\n",
      "      Processing year 2079\n",
      "      Processing year 2080\n",
      "    Processing GCM inmcm4\n",
      "      Processing year 2061\n",
      "      Processing year 2062\n",
      "      Processing year 2063\n",
      "      Processing year 2064\n",
      "      Processing year 2065\n",
      "      Processing year 2066\n",
      "      Processing year 2067\n",
      "      Processing year 2068\n",
      "      Processing year 2069\n",
      "      Processing year 2070\n",
      "      Processing year 2071\n",
      "      Processing year 2072\n",
      "      Processing year 2073\n",
      "      Processing year 2074\n",
      "      Processing year 2075\n",
      "      Processing year 2076\n",
      "      Processing year 2077\n",
      "      Processing year 2078\n",
      "      Processing year 2079\n",
      "      Processing year 2080\n",
      "    Processing GCM IPSL-CM5A-LR\n",
      "      Processing year 2061\n",
      "      Processing year 2062\n",
      "      Processing year 2063\n",
      "      Processing year 2064\n",
      "      Processing year 2065\n",
      "      Processing year 2066\n",
      "      Processing year 2067\n",
      "      Processing year 2068\n",
      "      Processing year 2069\n",
      "      Processing year 2070\n",
      "      Processing year 2071\n",
      "      Processing year 2072\n",
      "      Processing year 2073\n",
      "      Processing year 2074\n",
      "      Processing year 2075\n",
      "      Processing year 2076\n",
      "      Processing year 2077\n",
      "      Processing year 2078\n",
      "      Processing year 2079\n",
      "      Processing year 2080\n",
      "    Processing GCM IPSL-CM5A-MR\n",
      "      Processing year 2061\n",
      "      Processing year 2062\n",
      "      Processing year 2063\n",
      "      Processing year 2064\n",
      "      Processing year 2065\n",
      "      Processing year 2066\n",
      "      Processing year 2067\n",
      "      Processing year 2068\n",
      "      Processing year 2069\n",
      "      Processing year 2070\n",
      "      Processing year 2071\n",
      "      Processing year 2072\n",
      "      Processing year 2073\n",
      "      Processing year 2074\n",
      "      Processing year 2075\n",
      "      Processing year 2076\n",
      "      Processing year 2077\n",
      "      Processing year 2078\n",
      "      Processing year 2079\n",
      "      Processing year 2080\n",
      "    Processing GCM IPSL-CM5B-LR\n",
      "      Processing year 2061\n",
      "      Processing year 2062\n",
      "      Processing year 2063\n",
      "      Processing year 2064\n",
      "      Processing year 2065\n",
      "      Processing year 2066\n",
      "      Processing year 2067\n",
      "      Processing year 2068\n",
      "      Processing year 2069\n",
      "      Processing year 2070\n",
      "      Processing year 2071\n",
      "      Processing year 2072\n",
      "      Processing year 2073\n",
      "      Processing year 2074\n",
      "      Processing year 2075\n",
      "      Processing year 2076\n",
      "      Processing year 2077\n",
      "      Processing year 2078\n",
      "      Processing year 2079\n",
      "      Processing year 2080\n",
      "    Processing GCM MIROC-ESM\n",
      "      Processing year 2061\n",
      "      Processing year 2062\n",
      "      Processing year 2063\n",
      "      Processing year 2064\n",
      "      Processing year 2065\n",
      "      Processing year 2066\n",
      "      Processing year 2067\n",
      "      Processing year 2068\n",
      "      Processing year 2069\n",
      "      Processing year 2070\n",
      "      Processing year 2071\n",
      "      Processing year 2072\n",
      "      Processing year 2073\n",
      "      Processing year 2074\n",
      "      Processing year 2075\n",
      "      Processing year 2076\n",
      "      Processing year 2077\n",
      "      Processing year 2078\n",
      "      Processing year 2079\n",
      "      Processing year 2080\n",
      "    Processing GCM MIROC-ESM-CHEM\n",
      "      Processing year 2061\n",
      "      Processing year 2062\n",
      "      Processing year 2063\n",
      "      Processing year 2064\n",
      "      Processing year 2065\n",
      "      Processing year 2066\n",
      "      Processing year 2067\n",
      "      Processing year 2068\n",
      "      Processing year 2069\n",
      "      Processing year 2070\n",
      "      Processing year 2071\n",
      "      Processing year 2072\n",
      "      Processing year 2073\n",
      "      Processing year 2074\n",
      "      Processing year 2075\n",
      "      Processing year 2076\n",
      "      Processing year 2077\n",
      "      Processing year 2078\n",
      "      Processing year 2079\n",
      "      Processing year 2080\n",
      "    Processing GCM MIROC5\n",
      "      Processing year 2061\n",
      "      Processing year 2062\n",
      "      Processing year 2063\n",
      "      Processing year 2064\n",
      "      Processing year 2065\n",
      "      Processing year 2066\n",
      "      Processing year 2067\n",
      "      Processing year 2068\n",
      "      Processing year 2069\n",
      "      Processing year 2070\n",
      "      Processing year 2071\n",
      "      Processing year 2072\n",
      "      Processing year 2073\n",
      "      Processing year 2074\n",
      "      Processing year 2075\n",
      "      Processing year 2076\n",
      "      Processing year 2077\n",
      "      Processing year 2078\n",
      "      Processing year 2079\n",
      "      Processing year 2080\n",
      "    Processing GCM MPI-ESM-LR\n",
      "      Processing year 2061\n",
      "      Processing year 2062\n",
      "      Processing year 2063\n",
      "      Processing year 2064\n",
      "      Processing year 2065\n",
      "      Processing year 2066\n",
      "      Processing year 2067\n",
      "      Processing year 2068\n",
      "      Processing year 2069\n",
      "      Processing year 2070\n",
      "      Processing year 2071\n",
      "      Processing year 2072\n",
      "      Processing year 2073\n",
      "      Processing year 2074\n",
      "      Processing year 2075\n",
      "      Processing year 2076\n",
      "      Processing year 2077\n",
      "      Processing year 2078\n",
      "      Processing year 2079\n",
      "      Processing year 2080\n",
      "    Processing GCM MPI-ESM-MR\n",
      "      Processing year 2061\n",
      "      Processing year 2062\n",
      "      Processing year 2063\n",
      "      Processing year 2064\n",
      "      Processing year 2065\n",
      "      Processing year 2066\n",
      "      Processing year 2067\n",
      "      Processing year 2068\n",
      "      Processing year 2069\n",
      "      Processing year 2070\n",
      "      Processing year 2071\n",
      "      Processing year 2072\n",
      "      Processing year 2073\n",
      "      Processing year 2074\n",
      "      Processing year 2075\n",
      "      Processing year 2076\n",
      "      Processing year 2077\n",
      "      Processing year 2078\n",
      "      Processing year 2079\n",
      "      Processing year 2080\n",
      "    Processing GCM MRI-CGCM3\n",
      "      Processing year 2061\n",
      "      Processing year 2062\n",
      "      Processing year 2063\n",
      "      Processing year 2064\n",
      "      Processing year 2065\n",
      "      Processing year 2066\n",
      "      Processing year 2067\n",
      "      Processing year 2068\n",
      "      Processing year 2069\n",
      "      Processing year 2070\n",
      "      Processing year 2071\n",
      "      Processing year 2072\n",
      "      Processing year 2073\n",
      "      Processing year 2074\n",
      "      Processing year 2075\n",
      "      Processing year 2076\n",
      "      Processing year 2077\n",
      "      Processing year 2078\n",
      "      Processing year 2079\n",
      "      Processing year 2080\n",
      "    Processing GCM MRI-ESM1\n",
      "      Processing year 2061\n",
      "      Processing year 2062\n",
      "      Processing year 2063\n",
      "      Processing year 2064\n",
      "      Processing year 2065\n",
      "      Processing year 2066\n",
      "      Processing year 2067\n",
      "      Processing year 2068\n",
      "      Processing year 2069\n",
      "      Processing year 2070\n",
      "      Processing year 2071\n",
      "      Processing year 2072\n",
      "      Processing year 2073\n",
      "      Processing year 2074\n",
      "      Processing year 2075\n",
      "      Processing year 2076\n",
      "      Processing year 2077\n",
      "      Processing year 2078\n",
      "      Processing year 2079\n",
      "      Processing year 2080\n",
      "    Processing GCM NorESM1-M\n",
      "      Processing year 2061\n",
      "      Processing year 2062\n",
      "      Processing year 2063\n",
      "      Processing year 2064\n",
      "      Processing year 2065\n",
      "      Processing year 2066\n",
      "      Processing year 2067\n",
      "      Processing year 2068\n",
      "      Processing year 2069\n",
      "      Processing year 2070\n",
      "      Processing year 2071\n",
      "      Processing year 2072\n",
      "      Processing year 2073\n",
      "      Processing year 2074\n",
      "      Processing year 2075\n",
      "      Processing year 2076\n",
      "      Processing year 2077\n",
      "      Processing year 2078\n",
      "      Processing year 2079\n",
      "      Processing year 2080\n",
      "Writing time. Year = 2070\n",
      "Wrote: prcp_rcp85_2070_20yr_monthly.nc\n",
      "Cleared some memory up...restarting loop\n",
      "  Processing base year 2090\n",
      "    Processing GCM ACCESS1-0\n",
      "      Processing year 2081\n",
      "      Processing year 2082\n",
      "      Processing year 2083\n",
      "      Processing year 2084\n",
      "      Processing year 2085\n",
      "      Processing year 2086\n",
      "      Processing year 2087\n",
      "      Processing year 2088\n",
      "      Processing year 2089\n",
      "      Processing year 2090\n",
      "      Processing year 2091\n",
      "      Processing year 2092\n",
      "      Processing year 2093\n",
      "      Processing year 2094\n",
      "      Processing year 2095\n",
      "      Processing year 2096\n",
      "      Processing year 2097\n",
      "      Processing year 2098\n",
      "      Processing year 2099\n",
      "      Processing year 2100\n",
      "    Processing GCM ACCESS1-3\n",
      "      Processing year 2081\n",
      "      Processing year 2082\n",
      "      Processing year 2083\n",
      "      Processing year 2084\n",
      "      Processing year 2085\n",
      "      Processing year 2086\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "      Processing year 2087\n",
      "      Processing year 2088\n",
      "      Processing year 2089\n",
      "      Processing year 2090\n",
      "      Processing year 2091\n",
      "      Processing year 2092\n",
      "      Processing year 2093\n",
      "      Processing year 2094\n",
      "      Processing year 2095\n",
      "      Processing year 2096\n",
      "      Processing year 2097\n",
      "      Processing year 2098\n",
      "      Processing year 2099\n",
      "      Processing year 2100\n",
      "    Processing GCM CanESM2\n",
      "      Processing year 2081\n",
      "      Processing year 2082\n",
      "      Processing year 2083\n",
      "      Processing year 2084\n",
      "      Processing year 2085\n",
      "      Processing year 2086\n",
      "      Processing year 2087\n",
      "      Processing year 2088\n",
      "      Processing year 2089\n",
      "      Processing year 2090\n",
      "      Processing year 2091\n",
      "      Processing year 2092\n",
      "      Processing year 2093\n",
      "      Processing year 2094\n",
      "      Processing year 2095\n",
      "      Processing year 2096\n",
      "      Processing year 2097\n",
      "      Processing year 2098\n",
      "      Processing year 2099\n",
      "      Processing year 2100\n",
      "    Processing GCM CMCC-CESM\n",
      "      Processing year 2081\n",
      "      Processing year 2082\n",
      "      Processing year 2083\n",
      "      Processing year 2084\n",
      "      Processing year 2085\n",
      "      Processing year 2086\n",
      "      Processing year 2087\n",
      "      Processing year 2088\n",
      "      Processing year 2089\n",
      "      Processing year 2090\n",
      "      Processing year 2091\n",
      "      Processing year 2092\n",
      "      Processing year 2093\n",
      "      Processing year 2094\n",
      "      Processing year 2095\n",
      "      Processing year 2096\n",
      "      Processing year 2097\n",
      "      Processing year 2098\n",
      "      Processing year 2099\n",
      "      Processing year 2100\n",
      "    Processing GCM CMCC-CM\n",
      "      Processing year 2081\n",
      "      Processing year 2082\n",
      "      Processing year 2083\n",
      "      Processing year 2084\n",
      "      Processing year 2085\n",
      "      Processing year 2086\n",
      "      Processing year 2087\n",
      "      Processing year 2088\n",
      "      Processing year 2089\n",
      "      Processing year 2090\n",
      "      Processing year 2091\n",
      "      Processing year 2092\n",
      "      Processing year 2093\n",
      "      Processing year 2094\n",
      "      Processing year 2095\n",
      "      Processing year 2096\n",
      "      Processing year 2097\n",
      "      Processing year 2098\n",
      "      Processing year 2099\n",
      "      Processing year 2100\n",
      "    Processing GCM CMCC-CMS\n",
      "      Processing year 2081\n",
      "      Processing year 2082\n",
      "      Processing year 2083\n",
      "      Processing year 2084\n",
      "      Processing year 2085\n",
      "      Processing year 2086\n",
      "      Processing year 2087\n",
      "      Processing year 2088\n",
      "      Processing year 2089\n",
      "      Processing year 2090\n",
      "      Processing year 2091\n",
      "      Processing year 2092\n",
      "      Processing year 2093\n",
      "      Processing year 2094\n",
      "      Processing year 2095\n",
      "      Processing year 2096\n",
      "      Processing year 2097\n",
      "      Processing year 2098\n",
      "      Processing year 2099\n",
      "      Processing year 2100\n",
      "    Processing GCM CNRM-CM5\n",
      "      Processing year 2081\n",
      "      Processing year 2082\n",
      "      Processing year 2083\n",
      "      Processing year 2084\n",
      "      Processing year 2085\n",
      "      Processing year 2086\n",
      "      Processing year 2087\n",
      "      Processing year 2088\n",
      "      Processing year 2089\n",
      "      Processing year 2090\n",
      "      Processing year 2091\n",
      "      Processing year 2092\n",
      "      Processing year 2093\n",
      "      Processing year 2094\n",
      "      Processing year 2095\n",
      "      Processing year 2096\n",
      "      Processing year 2097\n",
      "      Processing year 2098\n",
      "      Processing year 2099\n",
      "      Processing year 2100\n",
      "    Processing GCM CSIRO-Mk3-6-0\n",
      "      Processing year 2081\n",
      "      Processing year 2082\n",
      "      Processing year 2083\n",
      "      Processing year 2084\n",
      "      Processing year 2085\n",
      "      Processing year 2086\n",
      "      Processing year 2087\n",
      "      Processing year 2088\n",
      "      Processing year 2089\n",
      "      Processing year 2090\n",
      "      Processing year 2091\n",
      "      Processing year 2092\n",
      "      Processing year 2093\n",
      "      Processing year 2094\n",
      "      Processing year 2095\n",
      "      Processing year 2096\n",
      "      Processing year 2097\n",
      "      Processing year 2098\n",
      "      Processing year 2099\n",
      "      Processing year 2100\n",
      "    Processing GCM GFDL-CM3\n",
      "      Processing year 2081\n",
      "      Processing year 2082\n",
      "      Processing year 2083\n",
      "      Processing year 2084\n",
      "      Processing year 2085\n",
      "      Processing year 2086\n",
      "      Processing year 2087\n",
      "      Processing year 2088\n",
      "      Processing year 2089\n",
      "      Processing year 2090\n",
      "      Processing year 2091\n",
      "      Processing year 2092\n",
      "      Processing year 2093\n",
      "      Processing year 2094\n",
      "      Processing year 2095\n",
      "      Processing year 2096\n",
      "      Processing year 2097\n",
      "      Processing year 2098\n",
      "      Processing year 2099\n",
      "      Processing year 2100\n",
      "    Processing GCM GFDL-ESM2G\n",
      "      Processing year 2081\n",
      "      Processing year 2082\n",
      "      Processing year 2083\n",
      "      Processing year 2084\n",
      "      Processing year 2085\n",
      "      Processing year 2086\n",
      "      Processing year 2087\n",
      "      Processing year 2088\n",
      "      Processing year 2089\n",
      "      Processing year 2090\n",
      "      Processing year 2091\n",
      "      Processing year 2092\n",
      "      Processing year 2093\n",
      "      Processing year 2094\n",
      "      Processing year 2095\n",
      "      Processing year 2096\n",
      "      Processing year 2097\n",
      "      Processing year 2098\n",
      "      Processing year 2099\n",
      "      Processing year 2100\n",
      "    Processing GCM GFDL-ESM2M\n",
      "      Processing year 2081\n",
      "      Processing year 2082\n",
      "      Processing year 2083\n",
      "      Processing year 2084\n",
      "      Processing year 2085\n",
      "      Processing year 2086\n",
      "      Processing year 2087\n",
      "      Processing year 2088\n",
      "      Processing year 2089\n",
      "      Processing year 2090\n",
      "      Processing year 2091\n",
      "      Processing year 2092\n",
      "      Processing year 2093\n",
      "      Processing year 2094\n",
      "      Processing year 2095\n",
      "      Processing year 2096\n",
      "      Processing year 2097\n",
      "      Processing year 2098\n",
      "      Processing year 2099\n",
      "      Processing year 2100\n",
      "    Processing GCM HadGEM2-CC\n",
      "      Processing year 2081\n",
      "      Processing year 2082\n",
      "      Processing year 2083\n",
      "      Processing year 2084\n",
      "      Processing year 2085\n",
      "      Processing year 2086\n",
      "      Processing year 2087\n",
      "      Processing year 2088\n",
      "      Processing year 2089\n",
      "      Processing year 2090\n",
      "      Processing year 2091\n",
      "      Processing year 2092\n",
      "      Processing year 2093\n",
      "      Processing year 2094\n",
      "      Processing year 2095\n",
      "      Processing year 2096\n",
      "      Processing year 2097\n",
      "      Processing year 2098\n",
      "      Processing year 2099\n",
      "      Processing year 2100\n",
      "    Processing GCM inmcm4\n",
      "      Processing year 2081\n",
      "      Processing year 2082\n",
      "      Processing year 2083\n",
      "      Processing year 2084\n",
      "      Processing year 2085\n",
      "      Processing year 2086\n",
      "      Processing year 2087\n",
      "      Processing year 2088\n",
      "      Processing year 2089\n",
      "      Processing year 2090\n",
      "      Processing year 2091\n",
      "      Processing year 2092\n",
      "      Processing year 2093\n",
      "      Processing year 2094\n",
      "      Processing year 2095\n",
      "      Processing year 2096\n",
      "      Processing year 2097\n",
      "      Processing year 2098\n",
      "      Processing year 2099\n",
      "      Processing year 2100\n",
      "    Processing GCM IPSL-CM5A-LR\n",
      "      Processing year 2081\n",
      "      Processing year 2082\n",
      "      Processing year 2083\n",
      "      Processing year 2084\n",
      "      Processing year 2085\n",
      "      Processing year 2086\n",
      "      Processing year 2087\n",
      "      Processing year 2088\n",
      "      Processing year 2089\n",
      "      Processing year 2090\n",
      "      Processing year 2091\n",
      "      Processing year 2092\n",
      "      Processing year 2093\n",
      "      Processing year 2094\n",
      "      Processing year 2095\n",
      "      Processing year 2096\n",
      "      Processing year 2097\n",
      "      Processing year 2098\n",
      "      Processing year 2099\n",
      "      Processing year 2100\n",
      "    Processing GCM IPSL-CM5A-MR\n",
      "      Processing year 2081\n",
      "      Processing year 2082\n",
      "      Processing year 2083\n",
      "      Processing year 2084\n",
      "      Processing year 2085\n",
      "      Processing year 2086\n",
      "      Processing year 2087\n",
      "      Processing year 2088\n",
      "      Processing year 2089\n",
      "      Processing year 2090\n",
      "      Processing year 2091\n",
      "      Processing year 2092\n",
      "      Processing year 2093\n",
      "      Processing year 2094\n",
      "      Processing year 2095\n",
      "      Processing year 2096\n",
      "      Processing year 2097\n",
      "      Processing year 2098\n",
      "      Processing year 2099\n",
      "      Processing year 2100\n",
      "    Processing GCM IPSL-CM5B-LR\n",
      "      Processing year 2081\n",
      "      Processing year 2082\n",
      "      Processing year 2083\n",
      "      Processing year 2084\n",
      "      Processing year 2085\n",
      "      Processing year 2086\n",
      "      Processing year 2087\n",
      "      Processing year 2088\n",
      "      Processing year 2089\n",
      "      Processing year 2090\n",
      "      Processing year 2091\n",
      "      Processing year 2092\n",
      "      Processing year 2093\n",
      "      Processing year 2094\n",
      "      Processing year 2095\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "      Processing year 2096\n",
      "      Processing year 2097\n",
      "      Processing year 2098\n",
      "      Processing year 2099\n",
      "      Processing year 2100\n",
      "    Processing GCM MIROC-ESM\n",
      "      Processing year 2081\n",
      "      Processing year 2082\n",
      "      Processing year 2083\n",
      "      Processing year 2084\n",
      "      Processing year 2085\n",
      "      Processing year 2086\n",
      "      Processing year 2087\n",
      "      Processing year 2088\n",
      "      Processing year 2089\n",
      "      Processing year 2090\n",
      "      Processing year 2091\n",
      "      Processing year 2092\n",
      "      Processing year 2093\n",
      "      Processing year 2094\n",
      "      Processing year 2095\n",
      "      Processing year 2096\n",
      "      Processing year 2097\n",
      "      Processing year 2098\n",
      "      Processing year 2099\n",
      "      Processing year 2100\n",
      "    Processing GCM MIROC-ESM-CHEM\n",
      "      Processing year 2081\n",
      "      Processing year 2082\n",
      "      Processing year 2083\n",
      "      Processing year 2084\n",
      "      Processing year 2085\n",
      "      Processing year 2086\n",
      "      Processing year 2087\n",
      "      Processing year 2088\n",
      "      Processing year 2089\n",
      "      Processing year 2090\n",
      "      Processing year 2091\n",
      "      Processing year 2092\n",
      "      Processing year 2093\n",
      "      Processing year 2094\n",
      "      Processing year 2095\n",
      "      Processing year 2096\n",
      "      Processing year 2097\n",
      "      Processing year 2098\n",
      "      Processing year 2099\n",
      "      Processing year 2100\n",
      "    Processing GCM MIROC5\n",
      "      Processing year 2081\n",
      "      Processing year 2082\n",
      "      Processing year 2083\n",
      "      Processing year 2084\n",
      "      Processing year 2085\n",
      "      Processing year 2086\n",
      "      Processing year 2087\n",
      "      Processing year 2088\n",
      "      Processing year 2089\n",
      "      Processing year 2090\n",
      "      Processing year 2091\n",
      "      Processing year 2092\n",
      "      Processing year 2093\n",
      "      Processing year 2094\n",
      "      Processing year 2095\n",
      "      Processing year 2096\n",
      "      Processing year 2097\n",
      "      Processing year 2098\n",
      "      Processing year 2099\n",
      "      Processing year 2100\n",
      "    Processing GCM MPI-ESM-LR\n",
      "      Processing year 2081\n",
      "      Processing year 2082\n",
      "      Processing year 2083\n",
      "      Processing year 2084\n",
      "      Processing year 2085\n",
      "      Processing year 2086\n",
      "      Processing year 2087\n",
      "      Processing year 2088\n",
      "      Processing year 2089\n",
      "      Processing year 2090\n",
      "      Processing year 2091\n",
      "      Processing year 2092\n",
      "      Processing year 2093\n",
      "      Processing year 2094\n",
      "      Processing year 2095\n",
      "      Processing year 2096\n",
      "      Processing year 2097\n",
      "      Processing year 2098\n",
      "      Processing year 2099\n",
      "      Processing year 2100\n",
      "    Processing GCM MPI-ESM-MR\n",
      "      Processing year 2081\n",
      "      Processing year 2082\n",
      "      Processing year 2083\n",
      "      Processing year 2084\n",
      "      Processing year 2085\n",
      "      Processing year 2086\n",
      "      Processing year 2087\n",
      "      Processing year 2088\n",
      "      Processing year 2089\n",
      "      Processing year 2090\n",
      "      Processing year 2091\n",
      "      Processing year 2092\n",
      "      Processing year 2093\n",
      "      Processing year 2094\n",
      "      Processing year 2095\n",
      "      Processing year 2096\n",
      "      Processing year 2097\n",
      "      Processing year 2098\n",
      "      Processing year 2099\n",
      "      Processing year 2100\n",
      "    Processing GCM MRI-CGCM3\n",
      "      Processing year 2081\n",
      "      Processing year 2082\n",
      "      Processing year 2083\n",
      "      Processing year 2084\n",
      "      Processing year 2085\n",
      "      Processing year 2086\n",
      "      Processing year 2087\n",
      "      Processing year 2088\n",
      "      Processing year 2089\n",
      "      Processing year 2090\n",
      "      Processing year 2091\n",
      "      Processing year 2092\n",
      "      Processing year 2093\n",
      "      Processing year 2094\n",
      "      Processing year 2095\n",
      "      Processing year 2096\n",
      "      Processing year 2097\n",
      "      Processing year 2098\n",
      "      Processing year 2099\n",
      "      Processing year 2100\n",
      "    Processing GCM MRI-ESM1\n",
      "      Processing year 2081\n",
      "      Processing year 2082\n",
      "      Processing year 2083\n",
      "      Processing year 2084\n",
      "      Processing year 2085\n",
      "      Processing year 2086\n",
      "      Processing year 2087\n",
      "      Processing year 2088\n",
      "      Processing year 2089\n",
      "      Processing year 2090\n",
      "      Processing year 2091\n",
      "      Processing year 2092\n",
      "      Processing year 2093\n",
      "      Processing year 2094\n",
      "      Processing year 2095\n",
      "      Processing year 2096\n",
      "      Processing year 2097\n",
      "      Processing year 2098\n",
      "      Processing year 2099\n",
      "      Processing year 2100\n",
      "    Processing GCM NorESM1-M\n",
      "      Processing year 2081\n",
      "      Processing year 2082\n",
      "      Processing year 2083\n",
      "      Processing year 2084\n",
      "      Processing year 2085\n",
      "      Processing year 2086\n",
      "      Processing year 2087\n",
      "      Processing year 2088\n",
      "      Processing year 2089\n",
      "      Processing year 2090\n",
      "      Processing year 2091\n",
      "      Processing year 2092\n",
      "      Processing year 2093\n",
      "      Processing year 2094\n",
      "      Processing year 2095\n",
      "      Processing year 2096\n",
      "      Processing year 2097\n",
      "      Processing year 2098\n",
      "      Processing year 2099\n",
      "      Processing year 2100\n",
      "Writing time. Year = 2090\n",
      "Wrote: prcp_rcp85_2090_20yr_monthly.nc\n",
      "Cleared some memory up...restarting loop\n"
     ]
    }
   ],
   "source": [
    "# loop across models\n",
    "models = [\"rcp45\",\"rcp85\"]\n",
    "for m in models:\n",
    "    print(\"Processing climate scenario \"+m)\n",
    "\n",
    "    # get subfolders/GCMs for current model\n",
    "    GCMs = os.listdir(os.path.join(base_folder,m))\n",
    "    \n",
    "    # remove 'GFDL-CM3' from RCP45 scenario -- incorrect precip via WICCI Climate Working Group\n",
    "    if m == \"rcp45\":\n",
    "        GCMs.remove('GFDL-CM3')\n",
    "        print(\"Removed GFDL-CM3\")\n",
    "\n",
    "    # loop across years in 20-year increments\n",
    "    for y in range(2030,2091,20):    # 2030,2091,20\n",
    "        print(\"  Processing base year \"+str(y))\n",
    "\n",
    "        # get mean and std for prcp for current 20-year window\n",
    "        # loop through GCMS first and then years in current window\n",
    "        # (create large blank masked arrays more memory efficient?)\n",
    "        num_models = len(GCMs)\n",
    "        years = 20\n",
    "        realizations = 3\n",
    "        prcp_all_array = np.ma.empty([num_models*years*realizations, 365, 48, 62])\n",
    "        prcp_means_array = np.ma.empty([num_models, 365, 48, 62])\n",
    "        \n",
    "         # indices for adding to blank arrays\n",
    "        i = 0\n",
    "        j = 0\n",
    "        \n",
    "        for gcm in GCMs:      \n",
    "            print(\"    Processing GCM \"+gcm)\n",
    "            prcp_gcm_array = []\n",
    "            for yr in range(y-9,y+11):\n",
    "                print(\"      Processing year \"+str(yr))\n",
    "                realizations = [\"01\",\"02\",\"03\"]\n",
    "                for r in realizations:\n",
    "                    nf = os.path.join(base_folder,m,gcm,\"r1i1p1\",\"prcp_\"+r+\"_\"+str(yr)+\".nc\")\n",
    "                    try:\n",
    "                        n = netCDF4.Dataset(nf)\n",
    "                    except:\n",
    "                        #print(\"File not found: \"+nf)\n",
    "                        q = 0\n",
    "                    else:\n",
    "                        prcp = n.variables['prcp']\n",
    "                        \n",
    "                        # convert to inches\n",
    "                        prcp = prcp[:,:,:] * 0.0393701\n",
    "                        \n",
    "                        # get mask - to reapply below\n",
    "                        mask = np.ma.getmask(prcp[0])\n",
    "                        \n",
    "                        # handle leap years -- average Feb 28 and 29\n",
    "                        if prcp.shape[0] == 366:\n",
    "                            # get and average Feb 28 and 29 for prcp\n",
    "                            prcp_leap = prcp[58:60,:,:]\n",
    "                            prcp_mean_leap = np.ma.mean(prcp_leap,axis=0)\n",
    "                            prcp_part1 = np.ma.append(prcp[0:58,:,:],np.ma.expand_dims(prcp_mean_leap,axis=0),axis=0)\n",
    "                            prcp = np.ma.append(prcp_part1,prcp[60:367,:,:],axis=0)\n",
    "                        \n",
    "                        # add prcp array for gcm, year, realization to total array\n",
    "                        prcp_all_array[i] = prcp[:,:,:]\n",
    "                        \n",
    "                        # add prcp array for gcm, year, realization to gcm-specific array\n",
    "                        prcp_gcm_array.append(prcp[:,:,:])\n",
    "                        \n",
    "                        i += 1\n",
    "\n",
    "            # calculate mean precip for gcm over time window (for std calculation)\n",
    "            prcp_model_mean = np.ma.mean(prcp_gcm_array,axis=0)\n",
    "\n",
    "            # add this to an array of means (one for each gcm)\n",
    "            prcp_means_array[j] = prcp_model_mean[:,:,:]\n",
    "            \n",
    "            j += 1\n",
    "\n",
    "        ##################################################\n",
    "        ## calculate metrics for current 20-year window ##\n",
    "        ##################################################\n",
    "        \n",
    "        # calculate daily mean from total array    \n",
    "        prcp_mean = np.ma.mean(prcp_all_array,axis=0)\n",
    "\n",
    "        # calculate monthly mean from total array\n",
    "        prcp_months_total = np.ma.empty([12, 48, 62])\n",
    "        prcp_months_std = np.ma.empty([12, 48, 62])\n",
    "        prcp_months_gt1 = np.ma.empty([12, 48, 62])\n",
    "        prcp_months_gt2 = np.ma.empty([12, 48, 62])\n",
    "        prcp_months_gt3 = np.ma.empty([12, 48, 62])\n",
    "        prcp_months_norain = np.ma.empty([12, 48, 62])\n",
    "        \n",
    "        # index for julian/ordinal calendar (start/end of each month)\n",
    "        julian = [0,31,59,90,120,151,181,212,243,273,304,334,366]\n",
    "        \n",
    "        # loop through months\n",
    "        for month in range(1,13):\n",
    "            \n",
    "            # calculate month precip total for each month\n",
    "            prcp_month_total = np.ma.sum(prcp_mean[(julian[month-1]):(julian[month])],axis=0)\n",
    "            prcp_months_total[month-1] = prcp_month_total\n",
    "            \n",
    "            # calculate inter-model standard deviation\n",
    "            all_std_nparray = np.array(prcp_means_array)\n",
    "            month_std_nparray = all_std_nparray[:,(julian[month-1]):(julian[month]),:,:]\n",
    "            # loop over each model -- add to total monthly sample array\n",
    "            for i in range(0, (np.shape(all_std_nparray)[0])):\n",
    "                model_array = month_std_nparray[i]\n",
    "                if i==0:\n",
    "                    concat_array = model_array\n",
    "                else:\n",
    "                    concat_array = np.ma.concatenate((concat_array,model_array),axis=0)\n",
    "            std_month = np.ma.std(concat_array,axis=0)\n",
    "            prcp_months_std[month-1] = np.ma.masked_array(std_month,mask)\n",
    "            \n",
    "            # calc extreme estimates - prcp\n",
    "            all_nparray = np.array(prcp_all_array)\n",
    "            month_nparray = all_nparray[:,(julian[month-1]):(julian[month]),:,:]\n",
    "            # get number of days in month and number of \"samples\"(needed for estimate formula)\n",
    "            num_days = np.shape(month_nparray)[1]\n",
    "            num_samples = np.shape(month_nparray)[0]*np.shape(month_nparray)[1]     \n",
    "            \n",
    "            # prcp extreme - prcp gt 1\" (25.4 mm)\n",
    "            # get number of days meeting threshold (by iteration)\n",
    "            extreme_sum_days = (month_nparray >= 1).sum(axis=0)\n",
    "            # add number of days (by month)\n",
    "            extreme_sum = extreme_sum_days.sum(axis=0)\n",
    "            # estimate of number of monthly extreme days\n",
    "            extreme_days = (extreme_sum / num_samples) * num_days\n",
    "            # add to array\n",
    "            prcp_months_gt1[month-1] = np.ma.masked_array(extreme_days,mask)\n",
    "            \n",
    "            # prcp extreme - prcp gt 2\" (50.8 mm)     \n",
    "            # get number of days meeting threshold (by iteration)\n",
    "            extreme_sum_days = (month_nparray >= 2).sum(axis=0)\n",
    "            # add number of days (by month)\n",
    "            extreme_sum = extreme_sum_days.sum(axis=0)\n",
    "            # estimate of number of monthly extreme days\n",
    "            extreme_days = (extreme_sum / num_samples) * num_days\n",
    "            # add to array\n",
    "            prcp_months_gt2[month-1] = np.ma.masked_array(extreme_days,mask)\n",
    "            \n",
    "            # prcp extreme - prcp gt 3\" (76.2 mm)     \n",
    "            # get number of days meeting threshold (by iteration)\n",
    "            extreme_sum_days = (month_nparray >= 3).sum(axis=0)\n",
    "            # add number of days (by month)\n",
    "            extreme_sum = extreme_sum_days.sum(axis=0)\n",
    "            # estimate of number of monthly extreme days\n",
    "            extreme_days = (extreme_sum / num_samples) * num_days\n",
    "            # add to array\n",
    "            prcp_months_gt3[month-1] = np.ma.masked_array(extreme_days,mask)\n",
    "            \n",
    "            # prcp extreme(?) - prcp is zero    \n",
    "            # get number of days meeting threshold (by iteration)\n",
    "            extreme_sum_days = (month_nparray == 0).sum(axis=0)\n",
    "            # add number of days (by month)\n",
    "            extreme_sum = extreme_sum_days.sum(axis=0)\n",
    "            # estimate of number of monthly extreme days\n",
    "            extreme_days = (extreme_sum / num_samples) * num_days\n",
    "            # add to array\n",
    "            prcp_months_norain[month-1] = np.ma.masked_array(extreme_days,mask)\n",
    "\n",
    "        # write netcdf files with results\n",
    "        filename = \"prcp_\"+m+\"_\"+str(y)+\"_20yr_monthly.nc\"\n",
    "        newfile = os.path.join(out_folder,filename)\n",
    "        ncfile = netCDF4.Dataset(newfile,mode='w',format='NETCDF4_CLASSIC')\n",
    "        lat_dim = ncfile.createDimension('lat', 48)     # latitude axis\n",
    "        lon_dim = ncfile.createDimension('lon', 62)    # longitude axis\n",
    "        time_dim = ncfile.createDimension('time', None) # unlimited axis (can be appended to).\n",
    "\n",
    "        ncfile.title='Aggregate monthly prcp values for WICCI downscaled climate data for all GCMs for '+m+' and 20-year window around year '+str(y)\n",
    "        ncfile.subtitle=\"Data source: UW-Madison WICCI; Data aggregation: Eric Compas, compase@uww.edu\"\n",
    "        lat = ncfile.createVariable('lat', np.float64, ('lat',))\n",
    "        lat.units = 'degrees_north'\n",
    "        lat.long_name = 'latitude'\n",
    "        lon = ncfile.createVariable('lon', np.float64, ('lon',))\n",
    "        lon.units = 'degrees_east'\n",
    "        lon.long_name = 'longitude'\n",
    "        time = ncfile.createVariable('time', np.float64, ('time',))\n",
    "        timeunits = 'days since '+str(y)+'-01-01'\n",
    "        time.units = timeunits\n",
    "        time.long_name = 'time'\n",
    "\n",
    "        prcp_total = ncfile.createVariable('prcp_total',np.float32,('time','lat','lon')) # note: unlimited dimension is leftmost\n",
    "        prcp_total.units = 'inches' # inches\n",
    "        prcp_total.standard_name = 'total of mean of daily precipitation (inches) per month across 20-year window'\n",
    "        prcp_total.missing_value = -32768\n",
    "        \n",
    "        prcp_std = ncfile.createVariable('prcp_std',np.float32,('time','lat','lon')) # note: unlimited dimension is leftmost\n",
    "        prcp_std.units = 'inches' # inches\n",
    "        prcp_std.standard_name = 'standard deviation of daily precipitation (inches) per month across 20-year window'\n",
    "        prcp_std.missing_value = -32768\n",
    "        \n",
    "        prcp_gt1 = ncfile.createVariable('prcp_gt1',np.float32,('time','lat','lon')) # note: unlimited dimension is leftmost\n",
    "        prcp_gt1.units = 'days' # number of days\n",
    "        prcp_gt1.standard_name = 'estimated number of days 1\" or more of rain (25.4 mm) per month across 20-year window'\n",
    "        prcp_gt1.missing_value = -32768\n",
    "        \n",
    "        prcp_gt2 = ncfile.createVariable('prcp_gt2',np.float32,('time','lat','lon')) # note: unlimited dimension is leftmost\n",
    "        prcp_gt2.units = 'days' # number of days\n",
    "        prcp_gt2.standard_name = 'estimated number of days 2\" or more of rain (50.8 mm) per month across 20-year window'\n",
    "        prcp_gt2.missing_value = -32768\n",
    "        \n",
    "        prcp_gt3 = ncfile.createVariable('prcp_gt3',np.float32,('time','lat','lon')) # note: unlimited dimension is leftmost\n",
    "        prcp_gt3.units = 'days' # number of days\n",
    "        prcp_gt3.standard_name = 'estimated number of days 3\" or more of rain (76.2 mm) per month across 20-year window'\n",
    "        prcp_gt3.missing_value = -32768\n",
    "        \n",
    "        prcp_norain = ncfile.createVariable('prcp_norain',np.float32,('time','lat','lon')) # note: unlimited dimension is leftmost\n",
    "        prcp_norain.units = 'days' # number of days\n",
    "        prcp_norain.standard_name = 'estimated number of days with no rain (0 mm) per month across 20-year window'\n",
    "        prcp_norain.missing_value = -32768 \n",
    "\n",
    "        # Write latitudes, longitudes\n",
    "        # Note: the \":\" is necessary in these \"write\" statements\n",
    "        n_lat = n.variables['lat']\n",
    "        n_lon = n.variables['lon']\n",
    "        lat[:] = n_lat[:]\n",
    "        lon[:] = n_lon[:]\n",
    "\n",
    "        # write temp variables\n",
    "        prcp_total[:,:,:] = prcp_months_total\n",
    "        prcp_std[:,:,:] = prcp_months_std\n",
    "        prcp_gt1[:,:,:] = prcp_months_gt1\n",
    "        prcp_gt2[:,:,:] = prcp_months_gt2\n",
    "        prcp_gt3[:,:,:] = prcp_months_gt3\n",
    "        prcp_norain[:,:,:] = prcp_months_norain\n",
    "\n",
    "        # write time\n",
    "        yystart = y\n",
    "        nyears = 1\n",
    "        ntime = 12\n",
    "        print(\"Writing time. Year = \"+str(y))\n",
    "        datesout = [datetime.datetime(y,mm,15,0) for yy in range(yystart,yystart+nyears) for mm in range(1,13)]\n",
    "        time[:] = date2num(datesout,timeunits)\n",
    "\n",
    "        # close file\n",
    "        ncfile.close()\n",
    "\n",
    "        # report progress\n",
    "        print(\"Wrote: \"+filename)\n",
    "        \n",
    "        \n",
    "        # clear some memory, invoke Python garbage collector\n",
    "        del prcp_all_array\n",
    "        del prcp_means_array\n",
    "        del prcp_months_total\n",
    "        del prcp_months_std\n",
    "        del prcp_months_gt1\n",
    "        del prcp_months_gt2\n",
    "        del prcp_months_gt3\n",
    "        del prcp_months_norain\n",
    "        del prcp_mean\n",
    "        gc.collect()\n",
    "        print(\"Cleared some memory up...restarting loop\")\n",
    "        "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "7cf18b66",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "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.8.12"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}