SII NRT Daily CSV Comparisons

{
 "metadata": {
  "name": "daily_diffs"
 },
 "nbformat": 3,
 "nbformat_minor": 0,
 "worksheets": [
  {
   "cells": [
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "## Comparison of daily extents using SMMR (N07) and SSMI (F13, F17) pole hole files in Sea Ice Index and NSIDC-0081"
     ]
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "There are 47 days where the openwater grid cells creep into the SMMR mask in the NRT data with the max difference of 0.016 M km^2\n"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "ne_indices = nrt_pre.extent != nrt_f17.extent\n",
      "pre_ne = nrt_pre[ne_indices]\n",
      "f17_ne = nrt_f17[ne_indices]\n",
      "f13_ne = nrt_f13[ne_indices]\n",
      "print \"Number of extents different:\", len(f13_ne)\n",
      "print \"\\nYYY\\tMM\\tDD\\t0081 F13\\t\\t0081 F13\\t\\t0081 F17\\t\\tDifference\"\n",
      "print \"\\t\\t\\tSII N07\\t\\tSII None\\t\\tSII None\\t\\tSII N07 - SII None\"\n",
      "for year, month, day, pre, f13, f17 in zip(pre_ne.year, pre_ne.mm, pre_ne.dd, pre_ne.extent, f13_ne.extent, f17_ne.extent):\n",
      "    print year, \"\\t\", month, \"\\t\", day, \"\\t\", pre, \"\\t\\t\", f13, \"\\t\\t\", f17, \"\\t\\t\", pre - f17"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Number of extents different: 47\n",
        "\n",
        "YYY\tMM\tDD\t0081 F13\t\t0081 F13\t\t0081 F17\t\tDifference\n",
        "\t\t\tSII N07\t\tSII None\t\tSII None\t\tSII N07 - SII None\n",
        "2014 \t8 \t18 \t5.893 \t\t5.892 \t\t5.892 \t\t0.000999999999999\n",
        "2014 \t8 \t19 \t5.904 \t\t5.902 \t\t5.902 \t\t0.002\n",
        "2014 \t8 \t20 \t5.82 \t\t5.818 \t\t5.818 \t\t0.002\n",
        "2014 \t8 \t21 \t5.75 \t\t5.744 \t\t5.744 \t\t0.006\n",
        "2014 \t8 \t22 \t5.782 \t\t5.773 \t\t5.773 \t\t0.009\n",
        "2014 \t8 \t23 \t5.793 \t\t5.781 \t\t5.781 \t\t0.012\n",
        "2014 \t8 \t24 \t5.636 \t\t5.624 \t\t5.624 \t\t0.012\n",
        "2014 \t8 \t25 \t5.626 \t\t5.613 \t\t5.613 \t\t0.013\n",
        "2014 \t8 \t26 \t5.617 \t\t5.601 \t\t5.601 \t\t0.016\n",
        "2014 \t8 \t27 \t5.572 \t\t5.558 \t\t5.558 \t\t0.014\n",
        "2014 \t8 \t28 \t5.551 \t\t5.539 \t\t5.539 \t\t0.012\n",
        "2014 \t8 \t29 \t5.543 \t\t5.531 \t\t5.531 \t\t0.012\n",
        "2014 \t8 \t30 \t5.482 \t\t5.471 \t\t5.471 \t\t0.011\n",
        "2014 \t8 \t31 \t5.432 \t\t5.421 \t\t5.421 \t\t0.011\n",
        "2014 \t9 \t1 \t5.388 \t\t5.378 \t\t5.378 \t\t0.01\n",
        "2014 \t9 \t2 \t5.365 \t\t5.356 \t\t5.356 \t\t0.009\n",
        "2014 \t9 \t3 \t5.38 \t\t5.371 \t\t5.371 \t\t0.009\n",
        "2014 \t9 \t4 \t5.327 \t\t5.317 \t\t5.317 \t\t0.01\n",
        "2014 \t9 \t5 \t5.249 \t\t5.238 \t\t5.238 \t\t0.011\n",
        "2014 \t9 \t6 \t5.23 \t\t5.217 \t\t5.217 \t\t0.013\n",
        "2014 \t9 \t7 \t5.242 \t\t5.23 \t\t5.23 \t\t0.012\n",
        "2014 \t9 \t8 \t5.189 \t\t5.18 \t\t5.18 \t\t0.009\n",
        "2014 \t9 \t9 \t5.141 \t\t5.134 \t\t5.134 \t\t0.007\n",
        "2014 \t9 \t10 \t5.119 \t\t5.111 \t\t5.111 \t\t0.008\n",
        "2014 \t9 \t11 \t5.135 \t\t5.124 \t\t5.124 \t\t0.011\n",
        "2014 \t9 \t12 \t5.115 \t\t5.109 \t\t5.109 \t\t0.006\n",
        "2014 \t9 \t13 \t5.03 \t\t5.025 \t\t5.025 \t\t0.005\n",
        "2014 \t9 \t14 \t5.031 \t\t5.027 \t\t5.027 \t\t0.004\n",
        "2014 \t9 \t15 \t5.024 \t\t5.018 \t\t5.018 \t\t0.006"
       ]
      },
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "\n",
        "2014 \t9 \t16 \t4.983 \t\t4.979 \t\t4.979 \t\t0.004\n",
        "2014 \t9 \t17 \t5.01 \t\t5.005 \t\t5.005 \t\t0.005\n",
        "2014 \t9 \t18 \t5.063 \t\t5.056 \t\t5.056 \t\t0.007\n",
        "2014 \t9 \t19 \t5.097 \t\t5.088 \t\t5.088 \t\t0.009\n",
        "2014 \t9 \t20 \t5.115 \t\t5.104 \t\t5.104 \t\t0.011\n",
        "2014 \t9 \t21 \t5.151 \t\t5.139 \t\t5.139 \t\t0.012\n",
        "2014 \t9 \t22 \t5.141 \t\t5.131 \t\t5.131 \t\t0.01\n",
        "2014 \t9 \t23 \t5.184 \t\t5.179 \t\t5.179 \t\t0.005\n",
        "2014 \t9 \t24 \t5.237 \t\t5.232 \t\t5.232 \t\t0.005\n",
        "2014 \t9 \t25 \t5.297 \t\t5.294 \t\t5.294 \t\t0.003\n",
        "2014 \t9 \t26 \t5.368 \t\t5.366 \t\t5.366 \t\t0.002\n",
        "2014 \t9 \t27 \t5.332 \t\t5.328 \t\t5.328 \t\t0.004\n",
        "2014 \t9 \t28 \t5.394 \t\t5.39 \t\t5.39 \t\t0.004\n",
        "2014 \t9 \t29 \t5.445 \t\t5.439 \t\t5.439 \t\t0.006\n",
        "2014 \t9 \t30 \t5.53 \t\t5.517 \t\t5.517 \t\t0.013\n",
        "2014 \t10 \t1 \t5.623 \t\t5.61 \t\t5.61 \t\t0.013\n",
        "2014 \t10 \t2 \t5.662 \t\t5.654 \t\t5.654 \t\t0.008\n",
        "2014 \t10 \t4 \t5.793 \t\t5.792 \t\t5.792 \t\t0.000999999999999\n"
       ]
      }
     ],
     "prompt_number": 30
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import pandas as pd\n",
      "import numpy as np\n",
      "import matplotlib.pyplot as plt\n",
      "import matplotlib as mpl\n",
      "import os"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 15
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "base = '/sharehome/hwilcox/PM-488/'\n",
      "pre_dir = 'sii-pole-hole'\n",
      "f13_dir = 'f13-pole-hole'\n",
      "f17_dir = 'f17-pole-hole'"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 27
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "nrt_csv_pre = os.path.join(base, pre_dir, 'NH_seaice_extent_nrt.csv')\n",
      "nrt_csv_f13 = os.path.join(base, f13_dir, 'NH_seaice_extent_nrt.csv')\n",
      "nrt_csv_f17 = os.path.join(base, f17_dir, 'NH_seaice_extent_nrt.csv')"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 28
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "nrt_pre = pd.read_csv(nrt_csv_pre, header = None, skiprows=2, names=[\"year\", \"mm\", \"dd\", \"extent\", \"missing\", \"source\"])\n",
      "nrt_f13 = pd.read_csv(nrt_csv_f13, header = None, skiprows=2, names=[\"year\", \"mm\", \"dd\", \"extent\", \"missing\", \"source\"])\n",
      "nrt_f17 = pd.read_csv(nrt_csv_f17, header = None, skiprows=2, names=[\"year\", \"mm\", \"dd\", \"extent\", \"missing\", \"source\"])"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 29
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [],
     "language": "python",
     "metadata": {},
     "outputs": []
    }
   ],
   "metadata": {}
  }
 ]
}