CS_Figures_PosterOSM16_13Feb16.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "## generate colors of time dep phase and adv phase of CS transport and flux \n",
    "\n",
    "\n",
    "from math import *\n",
    "\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "import matplotlib as mpl\n",
    "\n",
    "import matplotlib.colors as mcolors\n",
    "\n",
    "from netCDF4 import Dataset\n",
    "\n",
    "import numpy as np\n",
    "\n",
    "import os \n",
    "\n",
    "import pandas as pd\n",
    "\n",
    "import pylab as pl\n",
    "\n",
    "import scipy.io\n",
    "\n",
    "import scipy as spy\n",
    "\n",
    "import seaborn as sns\n",
    "\n",
    "import sys\n",
    "\n",
    "\n",
    "lib_path = os.path.abspath('../../Building_canyon/BuildCanyon/PythonModulesMITgcm') # Add absolute path to my python scripts\n",
    "#lib_path = os.path.abspath('../BuildCanyon/PythonModulesMITgcm') # Add absolute path to my python scripts\n",
    "\n",
    "sys.path.append(lib_path)\n",
    "\n",
    "import ReadOutTools_MITgcm as rout \n",
    "import MetricsPythonTools as mpt\n",
    "\n",
    "\n",
    "\n",
    "# Files\n",
    "CGrid = '/ocean/kramosmu/MITgcm/TracerExperiments/CNTDIFF/run03/gridGlob.nc' # Smallest volume grid, closed bdy, no canyon.\n",
    "CGridOut = Dataset(CGrid)\n",
    "\n",
    "FluxTR01 = ('/ocean/kramosmu/MITgcm/TracerExperiments/CNTDIFF/run03/FluxTR01Glob.nc' )\n",
    "FluxTR013D = ('/ocean/kramosmu/MITgcm/TracerExperiments/3DDIFF/run06/FluxTR01Glob.nc' )\n",
    "\n",
    "# General input\n",
    "nx = 360\n",
    "ny = 360\n",
    "nz = 90\n",
    "nt = 19 # t dimension size \n",
    "\n",
    "rc = CGridOut.variables['RC']\n",
    "xc = rout.getField(CGrid, 'XC') # x coords tracer cells\n",
    "yc = rout.getField(CGrid, 'YC') # y coords tracer cells\n",
    "depth = rout.getField(CGrid, 'Depth') # y coords tracer cells\n",
    "\n",
    "drF = CGridOut.variables['drF'] # vertical distance between faces\n",
    "drC = CGridOut.variables['drC'] # vertical distance between centers\n",
    "\n",
    "times = np.arange(0,nt,1)\n",
    "MaskCan = rout.getMask(CGrid,'HFacC') \n",
    "hFacCCan = rout.getField(CGrid,'HFacC') \n",
    "\n",
    "\n",
    "#Transect definitions (indices x,y,z,t)\n",
    "  \n",
    "CS1 = [0,40,227,227,0,30,0,18]\n",
    "CS2 = [40,120,227,227,0,30,0,18]\n",
    "CS3 = [120,240,267,267,0,30,0,18]\n",
    "CS3sb = [120,240,227,227,0,30,0,18]\n",
    "CS4 = [240,320,227,227,0,30,0,18]\n",
    "CS5 = [320,359,227,227,0,30,0,18]\n",
    "AS1 = [120,120,227,267,0,30,0,18]\n",
    "AS2 = [240,240,227,267,0,30,0,18]\n",
    "LID1 = [120,180,227,267,30,30,0,18]\n",
    "LID2 = [180,240,227,267,30,30,0,18]\n",
    "  \n",
    "day = [0.5, 1., 1.5, 2., 2.5, 3., 3.5, 4., 4.5, 5., 5.5,  6., 6.5,  7., 7.5,  8., 8.5,  9.] # Fluxes are calculated between two outputs\n",
    "\n",
    "TracerList = ['Tr1']  \n",
    "fluxfile = [FluxTR01]\n",
    "fluxtr = ['1']\n",
    "  \n",
    "for f,tr,trstr in zip (fluxfile,fluxtr,TracerList):\n",
    "    \n",
    "    keyw = ('DFrITr0%s' %tr)\n",
    "    keywe = ('DFrETr0%s' %tr)\n",
    "    keyv = ('DFyETr0%s' %tr)\n",
    "    keyu = ('DFxETr0%s' %tr)\n",
    "    \n",
    "    Wnm,Vnm,Unm = mpt.get_TRAC(f, keyw ,keyv, keyu)\n",
    "    Wnme,Vnm,Unm = mpt.get_TRAC(f, keywe ,keyv, keyu)\n",
    "    \n",
    "    MaskExp = mpt.maskExpand(MaskCan,Unm)\n",
    "    U = np.ma.MaskedArray(Unm,mask=MaskExp)\n",
    "    V = np.ma.MaskedArray(Vnm,mask=MaskExp)\n",
    "    W = np.ma.MaskedArray(Wnm,mask=MaskExp)\n",
    "    We = np.ma.MaskedArray(Wnme,mask=MaskExp)\n",
    "\n",
    "\n",
    "    #Get slices\n",
    "    V_CS1a = mpt.slice_TRAC(V,CS1[0],CS1[1],CS1[2],CS1[3],CS1[4],CS1[5],CS1[6],CS1[7])\n",
    "    V_CS2a = mpt.slice_TRAC(V,CS2[0],CS2[1],CS2[2],CS2[3],CS2[4],CS2[5],CS2[6],CS2[7])\n",
    "    V_CS3a = mpt.slice_TRAC(V,CS3[0],CS3[1],CS3[2],CS3[3],CS3[4],CS3[5],CS3[6],CS3[7])\n",
    "    V_CS4a = mpt.slice_TRAC(V,CS4[0],CS4[1],CS4[2],CS4[3],CS4[4],CS4[5],CS4[6],CS4[7])\n",
    "    V_CS5a = mpt.slice_TRAC(V,CS5[0],CS5[1],CS5[2],CS5[3],CS5[4],CS5[5],CS5[6],CS5[7])\n",
    "    V_CS3sba = mpt.slice_TRAC(V,CS3sb[0],CS3sb[1],CS3sb[2],CS3sb[3],CS3sb[4],CS3sb[5],CS3sb[6],CS3sb[7])\n",
    "    U_AS1a = mpt.slice_TRAC(U,AS1[0],AS1[1],AS1[2],AS1[3],AS1[4],AS1[5],AS1[6],AS1[7])\n",
    "    U_AS2a = mpt.slice_TRAC(U,AS2[0],AS2[1],AS2[2],AS2[3],AS2[4],AS2[5],AS2[6],AS2[7])\n",
    "    W_LID1a = mpt.slice_TRAC(W,LID1[0],LID1[1],LID1[2],LID1[3],LID1[4],LID1[5],LID1[6],LID1[7])\n",
    "    W_LID2a = mpt.slice_TRAC(W,LID2[0],LID2[1],LID2[2],LID2[3],LID2[4],LID2[5],LID2[6],LID2[7])\n",
    "    W_LID1ae = mpt.slice_TRAC(We,LID1[0],LID1[1],LID1[2],LID1[3],LID1[4],LID1[5],LID1[6],LID1[7])\n",
    "    W_LID2ae = mpt.slice_TRAC(We,LID2[0],LID2[1],LID2[2],LID2[3],LID2[4],LID2[5],LID2[6],LID2[7])\n",
    "\n",
    "    \n",
    "# Concatenate arrays to plot\n",
    "Up = (np.concatenate((V_CS1a,V_CS2a),axis = 2))*1000.0\n",
    "ASup = (-U_AS1a)*1000.0\n",
    "Head = (V_CS3a)*1000.0\n",
    "ASdown = (U_AS2a[:,:,::-1])*1000.0\n",
    "Down= (np.concatenate((V_CS4a,V_CS5a),axis = 2))*1000.0\n",
    "Vert = (np.concatenate((W_LID1a,W_LID2a),axis = 2))*1000.0\n",
    "Verte = (np.concatenate((W_LID1ae,W_LID2ae),axis = 2))*1000.0\n",
    "\n",
    "# 3D case\n",
    "TracerList = ['Tr1']  \n",
    "fluxfile = [FluxTR013D]\n",
    "fluxtr = ['1']\n",
    "\n",
    "for f,tr,trstr in zip (fluxfile,fluxtr,TracerList):\n",
    "    \n",
    "    keyw = ('DFrITr0%s' %tr)\n",
    "    keywe = ('DFrETr0%s' %tr)\n",
    "    keyv = ('DFyETr0%s' %tr)\n",
    "    keyu = ('DFxETr0%s' %tr)\n",
    "    \n",
    "    Wnm,Vnm,Unm = mpt.get_TRAC(f, keyw ,keyv, keyu)\n",
    "    Wnme,Vnm,Unm = mpt.get_TRAC(f, keywe ,keyv, keyu)\n",
    "    \n",
    "    MaskExp = mpt.maskExpand(MaskCan,Unm)\n",
    "    U = np.ma.MaskedArray(Unm,mask=MaskExp)\n",
    "    V = np.ma.MaskedArray(Vnm,mask=MaskExp)\n",
    "    W = np.ma.MaskedArray(Wnm,mask=MaskExp)\n",
    "    We = np.ma.MaskedArray(Wnme,mask=MaskExp)\n",
    "\n",
    "\n",
    "    #Get slices\n",
    "    V_CS1b = mpt.slice_TRAC(V,CS1[0],CS1[1],CS1[2],CS1[3],CS1[4],CS1[5],CS1[6],CS1[7])\n",
    "    V_CS2b = mpt.slice_TRAC(V,CS2[0],CS2[1],CS2[2],CS2[3],CS2[4],CS2[5],CS2[6],CS2[7])\n",
    "    V_CS3b = mpt.slice_TRAC(V,CS3[0],CS3[1],CS3[2],CS3[3],CS3[4],CS3[5],CS3[6],CS3[7])\n",
    "    V_CS4b = mpt.slice_TRAC(V,CS4[0],CS4[1],CS4[2],CS4[3],CS4[4],CS4[5],CS4[6],CS4[7])\n",
    "    V_CS5b = mpt.slice_TRAC(V,CS5[0],CS5[1],CS5[2],CS5[3],CS5[4],CS5[5],CS5[6],CS5[7])\n",
    "    V_CS3sbb = mpt.slice_TRAC(V,CS3sb[0],CS3sb[1],CS3sb[2],CS3sb[3],CS3sb[4],CS3sb[5],CS3sb[6],CS3sb[7])\n",
    "    U_AS1b = mpt.slice_TRAC(U,AS1[0],AS1[1],AS1[2],AS1[3],AS1[4],AS1[5],AS1[6],AS1[7])\n",
    "    U_AS2b = mpt.slice_TRAC(U,AS2[0],AS2[1],AS2[2],AS2[3],AS2[4],AS2[5],AS2[6],AS2[7])\n",
    "    W_LID1b = mpt.slice_TRAC(W,LID1[0],LID1[1],LID1[2],LID1[3],LID1[4],LID1[5],LID1[6],LID1[7])\n",
    "    W_LID2b = mpt.slice_TRAC(W,LID2[0],LID2[1],LID2[2],LID2[3],LID2[4],LID2[5],LID2[6],LID2[7])\n",
    "    W_LID1be = mpt.slice_TRAC(We,LID1[0],LID1[1],LID1[2],LID1[3],LID1[4],LID1[5],LID1[6],LID1[7])\n",
    "    W_LID2be = mpt.slice_TRAC(We,LID2[0],LID2[1],LID2[2],LID2[3],LID2[4],LID2[5],LID2[6],LID2[7])\n",
    "\n",
    "    \n",
    "# Concatenate arrays to plot\n",
    "Upb = (np.concatenate((V_CS1b,V_CS2b),axis = 2))*1000.0\n",
    "ASupb = (-U_AS1b)*1000.0\n",
    "Headb = (V_CS3b)*1000.0\n",
    "ASdownb = (U_AS2b[:,:,::-1])*1000.0\n",
    "Downb= (np.concatenate((V_CS4b,V_CS5b),axis = 2))*1000.0\n",
    "Vertb = (np.concatenate((W_LID1b,W_LID2b),axis = 2))*1000.0\n",
    "Vertbe = (np.concatenate((W_LID1be,W_LID2be),axis = 2))*1000.0\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "\n",
    "## FIGURE\n",
    "sns.set_context(\"talk\", font_scale=0.9, rc={\"lines.linewidth\": 2.5})\n",
    "sns.set_palette( sns.hls_palette(11, l=.4, s=.8)) \n",
    "sns.set_style('dark')\n",
    "sns.set_context('talk')\n",
    "\n",
    "\n",
    "fig45=plt.figure(figsize=(25,12))\n",
    "\n",
    "tt = 12\n",
    "#ax1 = plt.subplot(2,6,1)\n",
    "#vmini = np.min(np.mean(Upb[10:,:,:],axis=0))\n",
    "#vmaxi = -np.min(np.mean(Upb[10:,:,:],axis=0))\n",
    "#cn = ax1.contourf(xc[227,0:120]/1000.0,rc[0:30],np.mean(Up[10:,:,:],axis=0)/(1E1),15, vmax = vmaxi/(1E1),\n",
    "#                  vmin = vmini/(1E1), cmap = 'RdYlBu_r')\n",
    "#plt.ylabel('Depth (m) ')\n",
    "\n",
    "\n",
    "#ax2 = plt.subplot(2,6,2)\n",
    "#vmaxi =  -np.min(np.mean(ASupb[10:,:,:],axis=0))\n",
    "#vmini = np.min(np.mean(ASupb[10:,:,:],axis=0))\n",
    "#cn = ax2.contourf(yc[227:267,120]/1000.0,rc[0:30],np.mean(ASup[10:,:,:],axis=0)/(1E1),15, vmax = vmaxi/(1E1),\n",
    "#                  vmin = vmini/(1E1),cmap = 'RdYlBu_r')\n",
    "\n",
    "ax3 = plt.subplot(2,6,3)\n",
    "vmaxi =  np.max(np.mean(Headb[10:,:,:],axis=0))\n",
    "vmini = -np.max(np.mean(Headb[10:,:,:],axis=0))\n",
    "cn = ax3.contourf(xc[267,120:240]/1000.0,rc[0:30],np.mean(Head[10:,:,:],axis=0)/(1E1),15, vmax = vmaxi/(1E1),\n",
    "                  vmin = vmini/(1E1),cmap = 'RdYlBu_r')\n",
    "plt.ylabel('Depth (m) ')\n",
    "\n",
    "ax4 = plt.subplot(2,6,4)\n",
    "vmaxi =  np.max(np.mean(ASdownb[10:,:,:],axis=0))\n",
    "vmini = -np.max(np.mean(ASdownb[10:,:,:],axis=0))\n",
    "plt.title(' $K_{v}=10^{-5}m^2s^{-1}$, Mean cross-shelf diffusive flux (advective phase)')\n",
    "cn = ax4.contourf(yc[227:267,120]/1000.0,rc[0:30],np.mean(ASdown[10:,:,:],axis=0)/(1E1),15, vmax = vmaxi/(1E1),\n",
    "                  vmin = vmini/(1E1),cmap = 'RdYlBu_r')\n",
    "\n",
    "ax5 = plt.subplot(2,6,5)\n",
    "vmaxi = -np.min(np.mean(Downb[10:,:,:],axis=0))\n",
    "vmini = np.min(np.mean(Downb[10:,:,:],axis=0))\n",
    "cn = ax5.contourf(xc[227,240:-1]/1000.0,rc[0:30],np.mean(Down[10:,:,:],axis=0)/(1E1),15, vmax = vmaxi/(1E1),\n",
    "                  vmin = vmini/(1E1),cmap = 'RdYlBu_r')\n",
    "\n",
    "ax6 = plt.subplot(2,6,6)\n",
    "vmaxi = np.max(np.mean(Vertb[10:,:,:],axis=0))\n",
    "vmini = -np.max(np.mean(Vertb[10:,:,:],axis=0))\n",
    "cn = ax6.contourf(xc[227,120:240]/1000.0,yc[227:267,120]/1000.0,np.mean(Vert[10:,:,:]+Verte[10:,:,:],axis=0)/(1E1),15, vmax = vmaxi/(1E1),\n",
    "                  vmin = vmini/(1E1),cmap = 'RdYlBu_r')\n",
    "shelfbreakline = ax6.contour(xc[227,120:240]/1000.0,yc[227:267,120]/1000.0,depth[227:267,120:240],[152.0],colors='k') \n",
    "plt.ylabel('Cross-shore distance (km) ')\n",
    "plt.title('Vertical (shelf-break depth)')\n",
    "\n",
    "\n",
    "### Second row\n",
    "#ax7 = plt.subplot(2,6,7)\n",
    "#vmini = np.min(np.mean(Upb[10:,:,:],axis=0))\n",
    "#vmaxi = -np.min(np.mean(Upb[10:,:,:],axis=0))\n",
    "#plt.ylabel('Depth (m) ')\n",
    "#plt.xlabel('Alongshore distance (km) ')\n",
    "#cn = ax7.contourf(xc[227,0:120]/1000.0,rc[0:30],np.mean(Upb[10:,:,:],axis=0)/(1E1),15, vmax = vmaxi/(1E1),\n",
    "#                  vmin = vmini/(1E1), cmap = 'RdYlBu_r')\n",
    "#cax,kw = mpl.colorbar.make_axes([ax1,ax7],location='bottom')\n",
    "#cb = plt.colorbar(cn, cax=cax,ticks=[np.linspace(vmini/(1E1), vmaxi/(1E1),10) ],format='%.1f',**kw)\n",
    "#cb.set_label('$10^1$ Mol/s')\n",
    "\n",
    "\n",
    "#ax8 = plt.subplot(2,6,8)\n",
    "#vmaxi =  -np.min(np.mean(ASupb[10:,:,:],axis=0))\n",
    "#vmini = np.min(np.mean(ASupb[10:,:,:],axis=0))\n",
    "#cn = ax8.contourf(yc[227:267,120]/1000.0,rc[0:30],np.mean(ASupb[10:,:,:],axis=0)/(1E1),15, vmax= vmaxi/(1E1),\n",
    "#                  vmin = vmini/(1E1),cmap = 'RdYlBu_r')\n",
    "#plt.xlabel('Cross-shore distance (km) ')\n",
    "#cax,kw = mpl.colorbar.make_axes([ax2,ax8],location='bottom')\n",
    "#cb = plt.colorbar(cn, cax=cax,ticks=[np.linspace(vmini/(1E1), vmaxi/(1E1),10) ],format='%.1f',**kw)\n",
    "#cb.set_label('$10^1$ Mol/s')\n",
    "\n",
    "ax9 = plt.subplot(2,6,9)\n",
    "vmaxi = np.max(np.mean(Headb[10:,:,:],axis=0))\n",
    "vmini = -np.max(np.mean(Headb[10:,:,:],axis=0))\n",
    "cn = ax9.contourf(xc[267,120:240]/1000.0,rc[0:30],np.mean(Headb[10:,:,:],axis=0)/(1E1),15, vmax = vmaxi/(1E1),\n",
    "                  vmin = vmini/(1E1),cmap = 'RdYlBu_r')\n",
    "plt.xlabel('Alongshore distance (km) ')\n",
    "plt.title(' $K_{v} = 10^{-3} m^2s^{-1}$')\n",
    "plt.ylabel('Depth (m) ')\n",
    "cax,kw = mpl.colorbar.make_axes([ax3,ax9],location='bottom')\n",
    "cb = plt.colorbar(cn, cax=cax,ticks=[np.linspace(vmini/(1E1), vmaxi/(1E1),10) ],format='%.1f',**kw)\n",
    "cb.set_label('$10^1$ Mol/s')\n",
    "\n",
    "ax10 = plt.subplot(2,6,10)\n",
    "vmaxi = np.max(np.mean(ASdownb[10:,:,:],axis=0))\n",
    "vmini = -np.max(np.mean(ASdownb[10:,:,:],axis=0))\n",
    "cn = ax10.contourf(yc[227:267,120]/1000.0,rc[0:30],np.mean(ASdownb[10:,:,:],axis=0)/(1E1),15, vmax = vmaxi/(1E1),\n",
    "                  vmin = vmini/(1E1),cmap = 'RdYlBu_r')\n",
    "plt.xlabel('Cross-shore distance (km) ')\n",
    "cax,kw = mpl.colorbar.make_axes([ax4,ax10],location='bottom')\n",
    "cb = plt.colorbar(cn, cax=cax,ticks=[np.linspace(vmini/(1E1), vmaxi/(1E1),10) ],format='%.1f',**kw)\n",
    "cb.set_label('$10^1$ Mol/s')\n",
    "\n",
    "ax11 = plt.subplot(2,6,11)\n",
    "vmaxi = -np.min(np.mean(Downb[10:,:,:],axis=0))\n",
    "vmini = np.min(np.mean(Downb[10:,:,:],axis=0))\n",
    "cn = ax11.contourf(xc[227,240:-1]/1000.0,rc[0:30],np.mean(Downb[10:,:,:],axis=0)/(1E1),15, vmax = vmaxi/(1E1),\n",
    "                  vmin = vmini/(1E1),cmap = 'RdYlBu_r')\n",
    "plt.xlabel('Alongshore distance (km) ')\n",
    "cax,kw = mpl.colorbar.make_axes([ax5,ax11],location='bottom')\n",
    "cb = plt.colorbar(cn, cax=cax,ticks=[np.linspace(vmini/(1E1), vmaxi/(1E1),10) ],format='%.1f',**kw)\n",
    "cb.set_label('$10^1$ Mol/s')\n",
    "\n",
    "ax12 = plt.subplot(2,6,12)\n",
    "vmaxi = np.max(np.mean(Vertb[10:,:,:],axis=0))\n",
    "vmini = -np.max(np.mean(Vertb[10:,:,:],axis=0))\n",
    "cn = ax12.contourf(xc[227,120:240]/1000.0,yc[227:267,120]/1000.0,np.mean(Vertb[10:,:,:]+Vertbe[10:,:,:],axis=0)/(1E1),15, vmax = vmaxi/(1E1),\n",
    "                  vmin = vmini/(1E1),cmap = 'RdYlBu_r')\n",
    "shelfbreakline = ax12.contour(xc[227,120:240]/1000.0,yc[227:267,120]/1000.0,depth[227:267,120:240],[152.0],colors='k') \n",
    "\n",
    "plt.xlabel('Alongshore distance (km) ')\n",
    "plt.ylabel('Cross-shore distance (km) ')\n",
    "cax,kw = mpl.colorbar.make_axes([ax6,ax12],location='bottom')\n",
    "cb = plt.colorbar(cn,cax=cax,ticks=np.linspace(vmini/(1E1), vmaxi/(1E1),11),format='%.1f',**kw)\n",
    "cb.update_ticks()\n",
    "cb.set_label('$10^1$ Mol/s')\n",
    "\n",
    "plt.show()\n",
    "fig45.savefig('results/figures/PosterOSM16/CS_DIFF_AdvPh_CNTTr13Drun06.eps', format='eps', dpi=1000, bbox_inches='tight')\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "\n",
    "\n",
    "## FIGURE\n",
    "sns.set_context(\"talk\", font_scale=0.9, rc={\"lines.linewidth\": 2.5})\n",
    "sns.set_palette( sns.hls_palette(11, l=.4, s=.8)) \n",
    "sns.set_style('dark')\n",
    "sns.set_context('talk')\n",
    "\n",
    "\n",
    "fig45=plt.figure(figsize=(25,12))\n",
    "\n",
    "tt = 12\n",
    "#ax1 = plt.subplot(2,6,1)\n",
    "#vmini =  -np.max(np.mean(Upb[:10,:,:],axis=0))\n",
    "#vmaxi = np.max(np.mean(Upb[:10,:,:],axis=0))\n",
    "#cn = ax1.contourf(xc[227,0:120]/1000.0,rc[0:30],np.mean(Up[:10,:,:],axis=0),15, vmax = vmaxi,\n",
    "                  vmin = vmini, cmap = 'RdYlBu_r')\n",
    "#plt.ylabel('Depth (m) ')\n",
    "\n",
    "\n",
    "#ax2 = plt.subplot(2,6,2)\n",
    "#vmaxi = -np.min(np.mean(ASupb[:10,:,:],axis=0))\n",
    "#vmini = np.min(np.mean(ASupb[:10,:,:],axis=0))\n",
    "#cn = ax2.contourf(yc[227:267,120]/1000.0,rc[0:30],np.mean(ASup[:10,:,:],axis=0),15, vmax = vmaxi,\n",
    "                  vmin = vmini,cmap = 'RdYlBu_r')\n",
    "\n",
    "ax3 = plt.subplot(2,6,3)\n",
    "vmaxi = np.max(np.mean(Headb[:10,:,:],axis=0))\n",
    "vmini = -np.max(np.mean(Headb[:10,:,:],axis=0))\n",
    "cn = ax3.contourf(xc[267,120:240]/1000.0,rc[0:30],np.mean(Head[:10,:,:],axis=0),15, vmax = vmaxi,\n",
    "                  vmin = vmini,cmap = 'RdYlBu_r')\n",
    "plt.title(' $K_{v}=10^{-5}m^2s^{-1}$, Mean cross-shelf advective flux time-dependent phase')\n",
    "\n",
    "ax4 = plt.subplot(2,6,4)\n",
    "vmaxi = np.max(np.mean(ASdownb[:10,:,:],axis=0))\n",
    "vmini = -np.max(np.mean(ASdownb[:10,:,:],axis=0))\n",
    "cn = ax4.contourf(yc[227:267,120]/1000.0,rc[0:30],np.mean(ASdown[:10,:,:],axis=0),15, vmax = vmaxi,\n",
    "                  vmin = vmini,cmap = 'RdYlBu_r')\n",
    "\n",
    "ax5 = plt.subplot(2,6,5)\n",
    "vmaxi = np.max(np.mean(Downb[:10,:,:],axis=0))\n",
    "vmini = -np.max(np.mean(Downb[:10,:,:],axis=0))\n",
    "cn = ax5.contourf(xc[227,240:-1]/1000.0,rc[0:30],np.mean(Down[:10,:,:],axis=0),15, vmax = vmaxi,\n",
    "                  vmin = vmini,cmap = 'RdYlBu_r')\n",
    "\n",
    "ax6 = plt.subplot(2,6,6)\n",
    "vmaxi = np.max(np.mean(Vertb[:10,:,:],axis=0))\n",
    "vmini = -np.max(np.mean(Vertb[:10,:,:],axis=0))\n",
    "cn = ax6.contourf(xc[227,120:240]/1000.0,yc[227:267,120]/1000.0,np.mean(Vert[:10,:,:],axis=0),15, vmax = vmaxi,\n",
    "                  vmin = vmini,cmap = 'RdYlBu_r')\n",
    "shelfbreakline = ax6.contour(xc[227,120:240]/1000.0,yc[227:267,120]/1000.0,depth[227:267,120:240],[152.0],colors='k') \n",
    "plt.ylabel('Cross-shore distance (km) ')\n",
    "plt.title('Vertical (shelf-break depth)')\n",
    "\n",
    "\n",
    "### Second row\n",
    "#ax7 = plt.subplot(2,6,7)\n",
    "#vmini = -np.max(np.mean(Upb[:10,:,:],axis=0))\n",
    "#vmaxi = np.max(np.mean(Upb[:10,:,:],axis=0))\n",
    "#plt.ylabel('Depth (m) ')\n",
    "#plt.xlabel('Alongshore distance (km) ')\n",
    "#cn = ax7.contourf(xc[227,0:120]/1000.0,rc[0:30],np.mean(Upb[:10,:,:],axis=0),15, vmax = vmaxi,\n",
    "#                  vmin = vmini, cmap = 'RdYlBu_r')\n",
    "#cax,kw = mpl.colorbar.make_axes([ax1,ax7],location='bottom')\n",
    "#cb = plt.colorbar(cn, cax=cax,ticks=[np.linspace(vmini, vmaxi,10) ],format='%.2e',**kw)\n",
    "#cb.set_label('Mol/s')\n",
    "\n",
    "\n",
    "#ax8 = plt.subplot(2,6,8)\n",
    "#vmaxi = -np.min(np.mean(ASupb[:10,:,:],axis=0))\n",
    "#vmini = np.min(np.mean(ASupb[:10,:,:],axis=0))\n",
    "#cn = ax8.contourf(yc[227:267,120]/1000.0,rc[0:30],np.mean(ASupb[:10,:,:],axis=0),15, vmax = vmaxi,\n",
    "#                  vmin = vmini,cmap = 'RdYlBu_r')\n",
    "#plt.xlabel('Cross-shore distance (km) ')\n",
    "#cax,kw = mpl.colorbar.make_axes([ax2,ax8],location='bottom')\n",
    "#cb = plt.colorbar(cn, cax=cax,ticks=[np.linspace(vmini, vmaxi,10) ],format='%.2e',**kw)\n",
    "#cb.set_label('Mol/s')\n",
    "\n",
    "ax9 = plt.subplot(2,6,9)\n",
    "vmaxi = np.max(np.mean(Headb[:10,:,:],axis=0))\n",
    "vmini = -np.max(np.mean(Headb[:10,:,:],axis=0))\n",
    "cn = ax9.contourf(xc[267,120:240]/1000.0,rc[0:30],np.mean(Headb[:10,:,:],axis=0),15, vmax = vmaxi,\n",
    "                  vmin = vmini,cmap = 'RdYlBu_r')\n",
    "plt.xlabel('Alongshore distance (km) ')\n",
    "plt.title(' $K_{v} = 10^{-3} m^2s^{-1}$')\n",
    "cax,kw = mpl.colorbar.make_axes([ax3,ax9],location='bottom')\n",
    "cb = plt.colorbar(cn, cax=cax,ticks=[np.linspace(vmini, vmaxi,10) ],format='%.2e',**kw)\n",
    "cb.set_label('Mol/s')\n",
    "\n",
    "ax10 = plt.subplot(2,6,10)\n",
    "vmaxi = np.max(np.mean(ASdownb[:10,:,:],axis=0))\n",
    "vmini = -np.max(np.mean(ASdownb[:10,:,:],axis=0))\n",
    "cn = ax10.contourf(yc[227:267,120]/1000.0,rc[0:30],np.mean(ASdownb[:10,:,:],axis=0),15, vmax = vmaxi,\n",
    "                  vmin = vmini,cmap = 'RdYlBu_r')\n",
    "plt.xlabel('Cross-shore distance (km) ')\n",
    "cax,kw = mpl.colorbar.make_axes([ax4,ax10],location='bottom')\n",
    "cb = plt.colorbar(cn, cax=cax,ticks=[np.linspace(vmini, vmaxi,10) ],format='%.2e',**kw)\n",
    "cb.set_label('Mol/s')\n",
    "\n",
    "ax11 = plt.subplot(2,6,11)\n",
    "vmaxi = np.max(np.mean(Downb[:10,:,:],axis=0))\n",
    "vmini = -np.max(np.mean(Downb[:10,:,:],axis=0))\n",
    "cn = ax11.contourf(xc[227,240:-1]/1000.0,rc[0:30],np.mean(Downb[:10,:,:],axis=0),15, vmax = vmaxi,\n",
    "                  vmin = vmini,cmap = 'RdYlBu_r')\n",
    "plt.xlabel('Alongshore distance (km) ')\n",
    "cax,kw = mpl.colorbar.make_axes([ax5,ax11],location='bottom')\n",
    "cb = plt.colorbar(cn, cax=cax,ticks=[np.linspace(vmini, vmaxi,10) ],format='%.2e',**kw)\n",
    "cb.set_label('Mol/s')\n",
    "\n",
    "ax12 = plt.subplot(2,6,12)\n",
    "vmaxi = np.max(np.mean(Vertb[:10,:,:],axis=0))\n",
    "vmini = -np.max(np.mean(Vertb[:10,:,:],axis=0))\n",
    "cn = ax12.contourf(xc[227,120:240]/1000.0,yc[227:267,120]/1000.0,np.mean(Vertb[:10,:,:],axis=0),15, vmax = vmaxi,\n",
    "                  vmin = vmini,cmap = 'RdYlBu_r')\n",
    "shelfbreakline = ax12.contour(xc[227,120:240]/1000.0,yc[227:267,120]/1000.0,depth[227:267,120:240],[152.0],colors='k') \n",
    "\n",
    "plt.xlabel('Alongshore distance (km) ')\n",
    "plt.ylabel('Cross-shore distance (km) ')\n",
    "cax,kw = mpl.colorbar.make_axes([ax6,ax12],location='bottom')\n",
    "cb = plt.colorbar(cn, cax=cax,ticks=[np.linspace(vmini, vmaxi,4) ],format='%0.2e',**kw)\n",
    "cb.set_label('Mol/s')\n",
    "\n",
    "\n",
    "plt.show()\n",
    "fig45.savefig('results/figures/PosterOSM16/CS_ADV_TimeDep_CNTTr13Drun06.eps', format='eps', dpi=1000, bbox_inches='tight')\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "anaconda-cloud": {},
  "kernelspec": {
   "display_name": "Python [conda root]",
   "language": "python",
   "name": "conda-root-py"
  },
  "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.5.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}