p3_1_octavekernel.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Вернемся к Matlab/Octave (Занятие 3_1 в Jupyter)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "close all\n",
    "clear all"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "% octave specific function\n",
    "pkg load netcdf\n",
    "pkg load io\n",
    "pkg load statistics"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "fname4='/home/mizyuk/DataW2/AO_MexicanGulf_by_Dop5/SST_Reinolds_7D_AO_MexicanGulf_all_auto.mat.nc'; % имя файла"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "ans =\n",
      "\n",
      "  scalar structure containing the fields:\n",
      "\n",
      "    Filename = /home/mizyuk/DataW2/AO_MexicanGulf_by_Dop5/SST_Reinolds_7D_AO_MexicanGulf_all_auto.mat.nc\n",
      "    Name = /\n",
      "    Dimensions =\n",
      "\n",
      "      1x3 struct array containing the fields:\n",
      "\n",
      "        Name\n",
      "        Length\n",
      "        Unlimited\n",
      "\n",
      "    Variables =\n",
      "\n",
      "      1x4 struct array containing the fields:\n",
      "\n",
      "        Name\n",
      "        Dimensions\n",
      "        Size\n",
      "        Datatype\n",
      "        Attributes\n",
      "        ChunkSize\n",
      "        FillValue\n",
      "        DeflateLevel\n",
      "        Shuffle\n",
      "        Checksum\n",
      "\n",
      "    Attributes =\n",
      "\n",
      "      1x5 struct array containing the fields:\n",
      "\n",
      "        Name\n",
      "        Value\n",
      "\n",
      "    Groups = [](0x0)\n",
      "    Format = netcdf4\n",
      "\n"
     ]
    }
   ],
   "source": [
    "ncinfo(fname4) % инфо к nc файлу"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "загружаем перменные"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "x=ncread(fname4,'lon'); % долгота\n",
    "y=ncread(fname4,'lat'); % широта\n",
    "d=ncread(fname4,'time'); % дата\n",
    "aa=ncread(fname4,'aa'); % массив температуры"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "ans =\n",
      "\n",
      "    56    72   476\n",
      "\n"
     ]
    }
   ],
   "source": [
    "size(aa)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "aa(aa==0)=nan; % 0all zeros to Nan (Not A Number)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "ans = 1\n"
     ]
    }
   ],
   "source": [
    "picdir='Results/SST_Reinolds_7D/'; % название директории\n",
    "mkdir(picdir) % создаем директорию"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [],
   "source": [
    "%% среднее за выбранный год\n",
    "dv=datevec(d); % дата переводится в вектор (год, месяц, число, день, час, минута)\n",
    "year=dv(:,1); % год\n",
    "\n",
    "yu1=2013; % будем строить за 2013\n",
    "\n",
    "inye=find(year==yu1); % ищем индекс даты год которой равен 2013\n",
    "aa1=aa(:,:,inye); % выбираем только эти индексы\n",
    "\n",
    "aam=nanmean(aa1,3); % среднее по времени (3-ая размерность)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [],
   "source": [
    "aam(aam<1)=nan;"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "'print' is a function from the file /usr/share/octave/4.2.2/m/plot/util/print.m\r\n",
       "\r\n",
       " -- print ()\r\n",
       " -- print (OPTIONS)\r\n",
       " -- print (FILENAME, OPTIONS)\r\n",
       " -- print (H, FILENAME, OPTIONS)\r\n",
       "     Print a plot, or save it to a file.\r\n",
       "\r\n",
       "     Both output formatted for printing (PDF and PostScript), and many\r\n",
       "     bitmapped and vector image formats are supported.\r\n",
       "\r\n",
       "     FILENAME defines the name of the output file.  If the filename has\r\n",
       "     no suffix, one is inferred from the specified device and appended\r\n",
       "     to the filename.  If no filename is specified, the output is sent\r\n",
       "     to the printer.\r\n",
       "\r\n",
       "     H specifies the handle of the figure to print.  If no handle is\r\n",
       "     specified the current figure is used.\r\n",
       "\r\n",
       "     For output to a printer, PostScript file, or PDF file, the paper\r\n",
       "     size is specified by the figure's 'papersize' property.  The\r\n",
       "     location and size of the image on the page are specified by the\r\n",
       "     figure's 'paperposition' property.  The orientation of the page is\r\n",
       "     specified by the figure's 'paperorientation' property.\r\n",
       "\r\n",
       "     The width and height of images are specified by the figure's\r\n",
       "     'paperposition(3:4)' property values.\r\n",
       "\r\n",
       "     The 'print' command supports many OPTIONS:\r\n",
       "\r\n",
       "     '-fH'\r\n",
       "          Specify the handle, H, of the figure to be printed.  The\r\n",
       "          default is the current figure.\r\n",
       "\r\n",
       "     '-PPRINTER'\r\n",
       "          Set the PRINTER name to which the plot is sent if no FILENAME\r\n",
       "          is specified.\r\n",
       "\r\n",
       "     '-GGHOSTSCRIPT_COMMAND'\r\n",
       "          Specify the command for calling Ghostscript.  For Unix and\r\n",
       "          Windows the defaults are \"gs\" and \"gswin32c\", respectively.\r\n",
       "\r\n",
       "     '-color'\r\n",
       "     '-mono'\r\n",
       "          Color or monochrome output.\r\n",
       "\r\n",
       "     '-solid'\r\n",
       "     '-dashed'\r\n",
       "          Force all lines to be solid or dashed, respectively.\r\n",
       "\r\n",
       "     '-portrait'\r\n",
       "     '-landscape'\r\n",
       "          Specify the orientation of the plot for printed output.  For\r\n",
       "          non-printed output the aspect ratio of the output corresponds\r\n",
       "          to the plot area defined by the \"paperposition\" property in\r\n",
       "          the orientation specified.  This option is equivalent to\r\n",
       "          changing the figure's \"paperorientation\" property.\r\n",
       "\r\n",
       "     '-TextAlphaBits=N'\r\n",
       "     '-GraphicsAlphaBits=N'\r\n",
       "          Octave is able to produce output for various printers,\r\n",
       "          bitmaps, and vector formats by using Ghostscript.  For bitmap\r\n",
       "          and printer output anti-aliasing is applied using\r\n",
       "          Ghostscript's TextAlphaBits and GraphicsAlphaBits options.\r\n",
       "          The default number of bits are 4 and 1 respectively.  Allowed\r\n",
       "          values for N are 1, 2, or 4.\r\n",
       "\r\n",
       "     '-dDEVICE'\r\n",
       "          The available output format is specified by the option DEVICE,\r\n",
       "          and is one of:\r\n",
       "\r\n",
       "          'ps'\r\n",
       "          'ps2'\r\n",
       "          'psc'\r\n",
       "          'psc2'\r\n",
       "               PostScript (level 1 and 2, mono and color).  The\r\n",
       "               OpenGL-based toolkits always generate PostScript level\r\n",
       "               3.0.\r\n",
       "\r\n",
       "          'eps'\r\n",
       "          'eps2'\r\n",
       "          'epsc'\r\n",
       "          'epsc2'\r\n",
       "               Encapsulated PostScript (level 1 and 2, mono and color).\r\n",
       "               The OpenGL-based toolkits always generate PostScript\r\n",
       "               level 3.0.\r\n",
       "\r\n",
       "          'pslatex'\r\n",
       "          'epslatex'\r\n",
       "          'pdflatex'\r\n",
       "          'pslatexstandalone'\r\n",
       "          'epslatexstandalone'\r\n",
       "          'pdflatexstandalone'\r\n",
       "               Generate a LaTeX file 'FILENAME.tex' for the text\r\n",
       "               portions of a plot and a file 'FILENAME.(ps|eps|pdf)' for\r\n",
       "               the remaining graphics.  The graphics file suffix\r\n",
       "               .ps|eps|pdf is determined by the specified device type.\r\n",
       "               The LaTeX file produced by the 'standalone' option can be\r\n",
       "               processed directly by LaTeX.  The file generated without\r\n",
       "               the 'standalone' option is intended to be included from\r\n",
       "               another LaTeX document.  In either case, the LaTeX file\r\n",
       "               contains an '\\includegraphics' command so that the\r\n",
       "               generated graphics file is automatically included when\r\n",
       "               the LaTeX file is processed.  The text that is written to\r\n",
       "               the LaTeX file contains the strings *exactly* as they\r\n",
       "               were specified in the plot.  If any special characters of\r\n",
       "               the TeX mode interpreter were used, the file must be\r\n",
       "               edited before LaTeX processing.  Specifically, the\r\n",
       "               special characters must be enclosed with dollar signs ('$\r\n",
       "               ... $'), and other characters that are recognized by\r\n",
       "               LaTeX may also need editing (.e.g., braces).  The\r\n",
       "               'pdflatex' device, and any of the 'standalone' formats,\r\n",
       "               are not available with the Gnuplot toolkit.\r\n",
       "\r\n",
       "          'epscairo'\r\n",
       "          'pdfcairo'\r\n",
       "          'epscairolatex'\r\n",
       "          'pdfcairolatex'\r\n",
       "          'epscairolatexstandalone'\r\n",
       "          'pdfcairolatexstandalone'\r\n",
       "               Generate Cairo based output when using the Gnuplot\r\n",
       "               graphics toolkit.  The 'epscairo' and 'pdfcairo' devices\r\n",
       "               are synonymous with the 'epsc' device.  The LaTeX\r\n",
       "               variants generate a LaTeX file, 'FILENAME.tex', for the\r\n",
       "               text portions of a plot, and an image file,\r\n",
       "               'FILENAME.(eps|pdf)', for the graph portion of the plot.\r\n",
       "               The 'standalone' variants behave as described for\r\n",
       "               'epslatexstandalone' above.\r\n",
       "\r\n",
       "          'ill'\r\n",
       "          'aifm'\r\n",
       "               Adobe Illustrator (Obsolete for Gnuplot versions > 4.2)\r\n",
       "\r\n",
       "          'canvas'\r\n",
       "               Javascript-based drawing on HTML5 canvas viewable in a\r\n",
       "               web browser (only available for the Gnuplot graphics\r\n",
       "               toolkit).\r\n",
       "\r\n",
       "          'cdr'\r\n",
       "          'corel'\r\n",
       "               CorelDraw\r\n",
       "\r\n",
       "          'dxf'\r\n",
       "               AutoCAD\r\n",
       "\r\n",
       "          'emf'\r\n",
       "          'meta'\r\n",
       "               Microsoft Enhanced Metafile\r\n",
       "\r\n",
       "          'fig'\r\n",
       "               XFig.  For the Gnuplot graphics toolkit, the additional\r\n",
       "               options '-textspecial' or '-textnormal' can be used to\r\n",
       "               control whether the special flag should be set for the\r\n",
       "               text in the figure.  (default is '-textnormal')\r\n",
       "\r\n",
       "          'gif'\r\n",
       "               GIF image (only available for the Gnuplot graphics\r\n",
       "               toolkit)\r\n",
       "\r\n",
       "          'hpgl'\r\n",
       "               HP plotter language\r\n",
       "\r\n",
       "          'jpg'\r\n",
       "          'jpeg'\r\n",
       "               JPEG image\r\n",
       "\r\n",
       "          'latex'\r\n",
       "               LaTeX picture environment (only available for the Gnuplot\r\n",
       "               graphics toolkit).\r\n",
       "\r\n",
       "          'mf'\r\n",
       "               Metafont\r\n",
       "\r\n",
       "          'png'\r\n",
       "               Portable network graphics\r\n",
       "\r\n",
       "          'pbm'\r\n",
       "               PBMplus\r\n",
       "\r\n",
       "          'pdf'\r\n",
       "               Portable document format\r\n",
       "\r\n",
       "     'svg'\r\n",
       "          Scalable vector graphics\r\n",
       "\r\n",
       "     'tikz'\r\n",
       "     'tikzstandalone'\r\n",
       "          Generate a LaTeX file using PGF/TikZ format.  The OpenGL-based\r\n",
       "          toolkits create a PGF file while Gnuplot creates a TikZ file.\r\n",
       "          The 'tikzstandalone' device produces a LaTeX document which\r\n",
       "          includes the TikZ file ('tikzstandalone' and is only available\r\n",
       "          for the Gnuplot graphics toolkit).\r\n",
       "\r\n",
       "          If the device is omitted, it is inferred from the file\r\n",
       "          extension, or if there is no filename it is sent to the\r\n",
       "          printer as PostScript.\r\n",
       "\r\n",
       "     '-dGHOSTSCRIPT_DEVICE'\r\n",
       "          Additional devices are supported by Ghostscript.  Some\r\n",
       "          examples are;\r\n",
       "\r\n",
       "          'pdfwrite'\r\n",
       "               Produces pdf output from eps\r\n",
       "\r\n",
       "          'ljet2p'\r\n",
       "               HP LaserJet IIP\r\n",
       "\r\n",
       "          'pcx24b'\r\n",
       "               24-bit color PCX file format\r\n",
       "\r\n",
       "          'ppm'\r\n",
       "               Portable Pixel Map file format\r\n",
       "\r\n",
       "          For a complete list, type 'system (\"gs -h\")' to see what\r\n",
       "          formats and devices are available.\r\n",
       "\r\n",
       "          When Ghostscript output is sent to a printer the size is\r\n",
       "          determined by the figure's \"papersize\" property.  When the\r\n",
       "          output is sent to a file the size is determined by the plot\r\n",
       "          box defined by the figure's \"paperposition\" property.\r\n",
       "\r\n",
       "     '-append'\r\n",
       "          Append PostScript or PDF output to a pre-existing file of the\r\n",
       "          same type.\r\n",
       "\r\n",
       "     '-rNUM'\r\n",
       "          Resolution of bitmaps in pixels per inch.  For both metafiles\r\n",
       "          and SVG the default is the screen resolution; for other\r\n",
       "          formats it is 150 dpi.  To specify screen resolution, use\r\n",
       "          \"-r0\".\r\n",
       "\r\n",
       "     '-loose'\r\n",
       "     '-tight'\r\n",
       "          Force a tight or loose bounding box for eps files.  The\r\n",
       "          default is loose.\r\n",
       "\r\n",
       "     '-PREVIEW'\r\n",
       "          Add a preview to eps files.  Supported formats are:\r\n",
       "\r\n",
       "          '-interchange'\r\n",
       "               Provide an interchange preview.\r\n",
       "\r\n",
       "          '-metafile'\r\n",
       "               Provide a metafile preview.\r\n",
       "\r\n",
       "          '-pict'\r\n",
       "               Provide pict preview.\r\n",
       "\r\n",
       "          '-tiff'\r\n",
       "               Provide a tiff preview.\r\n",
       "\r\n",
       "     '-SXSIZE,YSIZE'\r\n",
       "          Plot size in pixels for EMF, GIF, JPEG, PBM, PNG, and SVG.\r\n",
       "          For PS, EPS, PDF, and other vector formats the plot size is in\r\n",
       "          points.  This option is equivalent to changing the size of the\r\n",
       "          plot box associated with the \"paperposition\" property.  When\r\n",
       "          using the command form of the print function you must quote\r\n",
       "          the XSIZE,YSIZE option.  For example, by writing \"-S640,480\".\r\n",
       "\r\n",
       "     '-FFONTNAME'\r\n",
       "     '-FFONTNAME:SIZE'\r\n",
       "     '-F:SIZE'\r\n",
       "          Use FONTNAME and/or FONTSIZE for all text.  FONTNAME is\r\n",
       "          ignored for some devices: dxf, fig, hpgl, etc.\r\n",
       "\r\n",
       "     The filename and options can be given in any order.\r\n",
       "\r\n",
       "     Example: Print to a file using the pdf device.\r\n",
       "\r\n",
       "          figure (1);\r\n",
       "          clf ();\r\n",
       "          surf (peaks);\r\n",
       "          print figure1.pdf\r\n",
       "\r\n",
       "     Example: Print to a file using jpg device.\r\n",
       "\r\n",
       "          clf ();\r\n",
       "          surf (peaks);\r\n",
       "          print -djpg figure2.jpg\r\n",
       "\r\n",
       "     Example: Print to printer named PS_printer using ps format.\r\n",
       "\r\n",
       "          clf ();\r\n",
       "          surf (peaks);\r\n",
       "          print -dpswrite -PPS_printer\r\n",
       "\r\n",
       "     See also: saveas, hgsave, orient, figure.\r\n",
       "\r\n",
       "Additional help for built-in functions and operators is\r\n",
       "available in the online version of the manual.  Use the command\r\n",
       "'doc <topic>' to search the manual index.\r\n",
       "\r\n",
       "Help and information about Octave is also available on the WWW\r\n",
       "at http://www.octave.org and via the help@octave.org\r\n",
       "mailing list.\r\n"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "print?"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "c = -7.9522\n",
      "warning: print.m: fig2dev binary is not available.\n",
      "Some output formats are not available.\n",
      "warning: called from\n",
      "    __print_parse_opts__ at line 388 column 9\n",
      "    print at line 316 column 8\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjAAAAGkCAIAAACgjIjwAAAJMmlDQ1BkZWZhdWx0X3JnYi5pY2MAAEiJlZVnUJNZF8fv8zzphUASQodQQ5EqJYCUEFoo0quoQOidUEVsiLgCK4qINEWQRQEXXJUia0UUC4uCAhZ0gywCyrpxFVFBWXDfGZ33HT+8/5l7z2/+c+bec8/5cAEgiINlwct7YlK6wNvJjhkYFMwE3yiMn5bC8fR0A9/VuxEArcR7ut/P+a4IEZFp/OW4uLxy+SmCdACg7GXWzEpPWeGjy0wPj//CZ1dYsFzgMt9Y4eh/eexLzr8s+pLj681dfhUKABwp+hsO/4b/c++KVDiC9NioyGymT3JUelaYIJKZttIJHpfL9BQkR8UmRH5T8P+V/B2lR2anr0RucsomQWx0TDrzfw41MjA0BF9n8cbrS48hRv9/z2dFX73kegDYcwAg+7564ZUAdO4CQPrRV09tua+UfAA67vAzBJn/eqiVDQ0IgALoQAYoAlWgCXSBETADlsAWOAAX4AF8QRDYAPggBiQCAcgCuWAHKABFYB84CKpALWgATaAVnAad4Dy4Aq6D2+AuGAaPgRBMgpdABN6BBQiCsBAZokEykBKkDulARhAbsoYcIDfIGwqCQqFoKAnKgHKhnVARVApVQXVQE/QLdA66At2EBqGH0Dg0A/0NfYQRmATTYQVYA9aH2TAHdoV94fVwNJwK58D58F64Aq6HT8Id8BX4NjwMC+GX8BwCECLCQJQRXYSNcBEPJBiJQgTIVqQQKUfqkVakG+lD7iFCZBb5gMKgaCgmShdliXJG+aH4qFTUVlQxqgp1AtWB6kXdQ42jRKjPaDJaHq2DtkDz0IHoaHQWugBdjm5Et6OvoYfRk+h3GAyGgWFhzDDOmCBMHGYzphhzGNOGuYwZxExg5rBYrAxWB2uF9cCGYdOxBdhK7EnsJewQdhL7HkfEKeGMcI64YFwSLg9XjmvGXcQN4aZwC3hxvDreAu+Bj8BvwpfgG/Dd+Dv4SfwCQYLAIlgRfAlxhB2ECkIr4RphjPCGSCSqEM2JXsRY4nZiBfEU8QZxnPiBRCVpk7ikEFIGaS/pOOky6SHpDZlM1iDbkoPJ6eS95CbyVfJT8nsxmpieGE8sQmybWLVYh9iQ2CsKnqJO4VA2UHIo5ZQzlDuUWXG8uIY4VzxMfKt4tfg58VHxOQmahKGEh0SiRLFEs8RNiWkqlqpBdaBGUPOpx6hXqRM0hKZK49L4tJ20Bto12iQdQ2fRefQ4ehH9Z/oAXSRJlTSW9JfMlqyWvCApZCAMDQaPkcAoYZxmjDA+SilIcaQipfZItUoNSc1Ly0nbSkdKF0q3SQ9Lf5RhyjjIxMvsl+mUeSKLktWW9ZLNkj0ie012Vo4uZynHlyuUOy33SB6W15b3lt8sf0y+X35OQVHBSSFFoVLhqsKsIkPRVjFOsUzxouKMEk3JWilWqUzpktILpiSTw0xgVjB7mSJleWVn5QzlOuUB5QUVloqfSp5Km8oTVYIqWzVKtUy1R1WkpqTmrpar1qL2SB2vzlaPUT+k3qc+r8HSCNDYrdGpMc2SZvFYOawW1pgmWdNGM1WzXvO+FkaLrRWvdVjrrjasbaIdo12tfUcH1jHVidU5rDO4Cr3KfFXSqvpVo7okXY5upm6L7rgeQ89NL0+vU++Vvpp+sP5+/T79zwYmBgkGDQaPDamGLoZ5ht2GfxtpG/GNqo3uryavdly9bXXX6tfGOsaRxkeMH5jQTNxNdpv0mHwyNTMVmLaazpipmYWa1ZiNsulsT3Yx+4Y52tzOfJv5efMPFqYW6RanLf6y1LWMt2y2nF7DWhO5pmHNhJWKVZhVnZXQmmkdan3UWmijbBNmU2/zzFbVNsK20XaKo8WJ45zkvLIzsBPYtdvNcy24W7iX7RF7J/tC+wEHqoOfQ5XDU0cVx2jHFkeRk4nTZqfLzmhnV+f9zqM8BR6f18QTuZi5bHHpdSW5+rhWuT5z03YTuHW7w+4u7gfcx9aqr01a2+kBPHgeBzyeeLI8Uz1/9cJ4eXpVez33NvTO9e7zofls9Gn2eedr51vi+9hP0y/Dr8ef4h/i3+Q/H2AfUBogDNQP3BJ4O0g2KDaoKxgb7B/cGDy3zmHdwXWTISYhBSEj61nrs9ff3CC7IWHDhY2UjWEbz4SiQwNCm0MXwzzC6sPmwnnhNeEiPpd/iP8ywjaiLGIm0iqyNHIqyiqqNGo62ir6QPRMjE1MecxsLDe2KvZ1nHNcbdx8vEf88filhICEtkRcYmjiuSRqUnxSb7JicnbyYIpOSkGKMNUi9WCqSOAqaEyD0tandaXTlz/F/gzNjF0Z45nWmdWZ77P8s85kS2QnZfdv0t60Z9NUjmPOT5tRm/mbe3KVc3fkjm/hbKnbCm0N39qzTXVb/rbJ7U7bT+wg7Ijf8VueQV5p3tudATu78xXyt+dP7HLa1VIgViAoGN1tubv2B9QPsT8M7Fm9p3LP58KIwltFBkXlRYvF/OJbPxr+WPHj0t6ovQMlpiVH9mH2Je0b2W+z/0SpRGlO6cQB9wMdZcyywrK3BzcevFluXF57iHAo45Cwwq2iq1Ktcl/lYlVM1XC1XXVbjXzNnpr5wxGHh47YHmmtVagtqv14NPbogzqnuo56jfryY5hjmceeN/g39P3E/qmpUbaxqPHT8aTjwhPeJ3qbzJqamuWbS1rgloyWmZMhJ+/+bP9zV6tua10bo63oFDiVcerFL6G/jJx2Pd1zhn2m9az62Zp2WnthB9SxqUPUGdMp7ArqGjzncq6n27K7/Ve9X4+fVz5ffUHyQslFwsX8i0uXci7NXU65PHsl+spEz8aex1cDr97v9eoduOZ67cZ1x+tX+zh9l25Y3Th/0+LmuVvsW523TW939Jv0t/9m8lv7gOlAxx2zO113ze92D64ZvDhkM3Tlnv296/d5928Prx0eHPEbeTAaMip8EPFg+mHCw9ePMh8tPN4+hh4rfCL+pPyp/NP637V+bxOaCi+M24/3P/N59niCP/Hyj7Q/Fifzn5Ofl08pTTVNG02fn3Gcufti3YvJlykvF2YL/pT4s+aV5quzf9n+1S8KFE2+Frxe+rv4jcyb42+N3/bMec49fZf4bmG+8L3M+xMf2B/6PgZ8nFrIWsQuVnzS+tT92fXz2FLi0tI/QiyQvpTNDAsAAAAJcEhZcwAACxMAAAsTAQCanBgAAAAddEVYdFNvZnR3YXJlAEdQTCBHaG9zdHNjcmlwdCA5LjI2WJButwAAIABJREFUeJzt3UFsG9edP/Cf/vtHsNjDkL748IeoyAlKAjsCqo0bFKQO9VYKSv730FQbi0IulQCR3j3sygUoKwVKRZawiB0SaAQsEJgUIPVimAzizV4qFnaw6UFDbGunKTAsQBVxFA3hQw6xhocA+7/of3jNeGZIPg2pmeET+f2gKGbmzcx7HI71y5v58b2Rk5MTAgAA6Lf/1e8GAAAAECEgAQCAIBCQAABACAhIAAAgBAQkAAAQAgISDKCHDx863AgA4kBAAhc8e/ZsZGRkZGTkyZMn/W3J6urqyMjIa6+9NjIyMjMz8+zZs7YbZ2ZmRqxWV1eNkzx8+HBkZMRJAFtdXb1w4cKFCxeMw1u3EBGr8fHjx5x9AIBOAM7szp077Ha6detW35tx69atBw8e3Lp1i4hu3LjRduOjR48ePHiQSqWMos8//9w4z9dff/3gwYOvv/7aYXXstOw85i3stDdu3Lh06RJbPTk5YUV37ty5ceMGEZXLZW8vCsD5gR4SuKBcLr/yyiuvvPKKEZmuXr360ksvseWXXnrp6tWr9G1PxeiOsI4I6y7Qt90ItoPRuSkUCiMjIy+99FI6nZ6ZmWEntJ3H8PnnnxPRzMzMzMzM6urqrVu3Xn755bYbL1++PDMzw+ply0Zriejx48evvfba48ePWQvT6TSrrlAotFa3urq6uroaDAYfPHjAAhLbQkRsldVuHHXhwoVUKpVOp69du8Y+tStfAcAg6HdEhHOP/V02ugWPHj06+bb38OjRI1Z6584dVloul8vlMpn6E5cuXbpx48bJt10HY4c7d+48evSIiFKpFNsyPT3dupu5T/bo0aNgMEhEr7zySiqVunPnTqeNDOujsI6LGWuYEWDeeOONBw8evPLKK53+vbB23rp1a3p6+tKlS2zjK6+8whpsPqG5ChbCT+2HAQwP9JDgrD744AMiunr1KusGsTiRTqeDwWC5XGZvYq5evfrgwYNLly49evSIRQij93D79u3bt28T0erqKitlwezzzz9np7p9+/bVq1ffeOMNtn+n8xDR5cuXHz9+fOvWrcuXLz98+PDatWvpdLrtxq4+4LVr11gHq9PHn5mZSaVSrFfEel3mhU5eeumlTz/9FKkWAAYEJDgrFjZefvnll19+mb6NT0R09erVhw8fPnjw4I033rD9dU6n03Nzc2zZKLp9+/b3vvc9ImLPsgxsB/MjtbbnYVXfuXNndXW1UCg8efJkenq6WCy23ejSR6fV1dW5ubm33nqLPc27fPnyp59+yoo+/vjjtm1+/Pjx48ePZ2ZmPvjgg0uXLrGrBwCEgARn9Pjx408//fTGjRvsAdeNGze++OILlk722muvsR7Aa6+9RkSXL19+9uzZW2+99dZbbz158qT13Qnr+rBStoVFuNu3bz958sR4f8M5z5MnT959993bt28/e/bs8ePHT548uXTpUtuNrnz2QqHw7rvvssduLOAZIZn1e1h8tSmXy9/73vfYC6ovvviibdACGFL9fmYI5xt7DWOkqLE3Ruyd0MnJCXt5w16TfP311+w1DBGxP+K2NytGX2F6eto4ifkQ9kqm9TxGY77++mvjyR4RsSd7bTeaG3/qOyS2A9vZvBtL0jOw5rHGG6u2E7JLZEREvEMCMBs5wfQT4KOHDx9euHDh8uXLbUtZj8coffLkibHl1q1bT548Md64cM7DjrKVtt3oEdZB5FfEvw4AwwkBCcT15MmTy5cvX7hwYWZmplgs3rlzp9t8BAA4RxCQQGjPnj374IMPnj17Zv5hEwAMJAQkAAAQArLsAABACAhIAAAgBAQkAAAQAgISAAAIAQEJAACEgIAEAABCQEACAAAhICABAIAQhA5IiqLoum7eoqqqpmn9ag8AAHjnr9bX1/vdhjZ0XX/zzTdPTk7+7d/+7eLFi9/5znfYlq+++uqjjz7685//PDU11e82AgCAmwQdOohNoZZKpTRNy2azu7u7xWLx+Ph4ZWWFiKanpz/++ON+txEAANz0v/vdgPaMmWZUVQ2FQkQUj8fZFttDPAAAGAyCBiQml8tVKhU24wALS4qi5PN5J3MQRCIRz9sHAOdHvV7vbwNc/6PU90/kOkEf2Rl0XZ+Zmfn9739PRLlcrtFoZDIZFpz4IpHImwcH3jfQYp1ofdBrHJ5K/a+xL5X6X+MZK12nX/ZwVDj8ft//fLv7R+luONz3T+Q6QbPsstns3t4eEQUCAbalVCo1Go2trS0n0QgAAM4dQR/ZJZPJ5eVlVVVrtRp7QMcSvhcWFtgOu7u7fWweAAC4TtCANDExcf/+/VqtNj8/z7pEm5ub/W4UAAB4SNBHdkQUCARisdj5ekAXDocHvsbhqRQfc8AqBfEJ2kMCABg8gX43QNf1Wq1GRLIsG2/oO9E0jY2ME4vF/GgcAhIACK63tDoiIhp3sRnnHfsVTaPRMLbIspxKpRKJhJOd5+bm0um014+sxH1kBwAAZ6fr+uzs7Pb2tjnAEFGtVrt+/Xo2m3Wyc7lcnp2dVVXV06YiIAEADLJCocAe08Xj8Q8//LBer9fr9ffee0+WZSIql8ulUql157m5uYcPH7KdNzY2JElqNpuLi4uejpWDgAQAMLA0Tdve3iaiaDS6tbU1MTHBticSiZ2dndHRUSIqFArG/mzneDy+ublpPKBLJpNbW1tE1Gw2K5WKd61FQAIA8EnQvf85pCgKW8hkMraiQCAwPz9PRI1Ggz2LM3Y2xg41GHkNR0dHPXxwhxCQAAAGlhE/jL6RGXtqR0TNZtO8nZOAZ9vTXciyA4ABNfn68+VvVvvXjn4aGxuLRqOdAozthZDxjG5/f9+W6m3MjDo2NuZBM/8CPSQAgIGVTCZ3d3fZG6BW7BmdJEks/IRCIfawrlwuG4/viEjX9eXlZbZn69M8F6GHBAAwjEqlUrlcJqK5uTlj48bGBhFVKpXFxUVZlqPRaK1Wq9VqzWZTkqSdnR1Pf4qEgAQAMHRyuRxLqJNl2TzDXCAQyGQyiqI0m00WioyiTCbT9kWUixCQAAB8smha/o8uj/2JdbXn0QAVRclms+x3r7Is7+zsmN8wlUqlfD7P+kPxeFySJCJiozasra2pqurpONcISAAAPuk2CHGO7SFJQ9O0bDZbrVbZ6tLSUjqdNkcjTdNYNLIFqpWVlWKxmM/ny+WyJEkrKys9fobTICABgFjOMHjd65a1t03LQ5pk91yxWCwUCixpOx6Pt516m+3A3hXZEvNSqdTR0VG5XN7e3raFMRchyw4AYMBls1lz16fT1NtsFIZ4PN423hiJ4OYXS+5CDwkAYJDlcjkjm47/Boj1n9h7o1Ye9YrMEJAAAHzifMgftxhj2Z0ajYhodHS00WjYxvk2eDqsKoNHdgAAA4s9hZMkqXUsu1bsoZyiKMa4DG1P5d18feghAQAMrP39fSKSJMk8pHer+fn5UCiUTCbL5XKz2VxeXra9Z2JT9pH1V7SuQ0ACgEGxMG5em3z9M2N5WIey+0sCQqPRYA/uOpmamgqFQhMTE5lMJp/P12q1mZmZeDzO5qcwZo+1/YrWdQhIAAADq9vBuVOplCRJLCXPNvXR3NxcJpPxNLUBAQkAwCf+JzXU6/VuD0kmk/F4vFKpHB0dNRqN0dHRYDAYj8c9HcWOQUACAACLQCCQTCb9rxdZdgAAIAQEJAAAEAIe2QHA+XXdsrZsWfsp/cpYft+X1sAZoYcEAABCQA8JAMAn/mfZnS/oIQEAgBAQkAAAQAgISAAAIAS8QwJB2aYNDXeTJ3WGKUctNZrPs04/49Z40s25jzssU5h+1M15hpPpRYx1htjgpOViXqFPjGVk2Z0LCEgAAD5BUgMfHtkBAIAQEJAAAEAICEgAACAEBCQAABACkhpAFKemxnFy3s6QVjduXT3uVLpOX1hLrAfa3lbzX14fm4o/s++6Tn8wrR1ai37CPa9T5svVVfqiP077csefLwZtBYfm1aDl23zBjaadFZIa+NBDAgAAIQgdkBRF0XXdvEVVVdsWAAAYDIIGJF3XZ2dn9/f3Z2dn9/b22MaFhYVSqbS4uKgoSn+bBwAArhM0IJXL5UQisbKysru7WyqViGhvby8UCm1ubm5tbRUKhX43EAAAXCZoUkMqlWILqqqGQiG2MDExQUShUKharfazcQAA4IGRk5OuxuDyVS6Xq1Qq6XQ6mUzmcrmpqalYLEZECwsLu7u7/GMjkcjBwYF5Szgc9q6pACAa81+AcDhcr9f72BgiikQi9a8OTt/P4dku9v8TuU7QHhKzsrKSTqdnZmaSySQRaZrGtjvsIa3b1g9cuxU6uRsOv+l9Lf2t8VxUys0Cv2Jdney0XzgcOTgw/YMf5x5nWw10rp+IvjQtf2at9Km1Ultxr1ngnAsSDr9/cPDPnQ/kDSnbs7N9m1eeLy5YrvvkjuVyPaO/N5a/H+mqgdAfgr5DymazLJchEPjLv+yJiYmjoyMi0jRNluV+Ng4AADwgaA8pmUwuLy+rqlqr1dLpNBElEolisZjL5arVqvGGCQAABoagAWliYuL+/fu1Wm1+fp4lNRDR/fv3FUUxbwEAgIEhaEAiokAgwFIYzFq3AACcGxg7iEvQd0gAADBsxO0hQSe2kTFdma67XS0dB7jkV+pRXhZfF2l149Z8OFt2nPk/YBWiBfOB1j1ftK7aSk8ZXLXzed63Tsv9ibV91qFf1+m/nNVBtvQ852wXti9fLgwP9JAAAEAICEgAACAEPLIDAPCLi0kN37h3KmGghwQAAEJAQAIAACHgkZ2IPEqca3lecNx+r7O1oU95WbaPZkpOm7Qmqr1u3fG7nU+jEP20cw1dzVluY77wtgPfJ/qxaXXcWmpLuju0ndZ83o5zsbcrdaovX25LwuezTnse42c+5xx6SAAAIAQEJAAAEAIe2QEA+AVZdlzoIQEAgBAQkAAAQAh4ZNc37qXS2Z4CjHde7TkVjFoyupyeqmW8NWueWNB0nnHeecLfRNYnTzoW2xprruTH1iJrll0wyMs3C17pMRvtFPyLZ26hbag928B3X1pXPzOd99hax6F1z2PO9bLt3UUqpl8ZlaYmHVo+5uHxuGU1aF4dxCdcAwcBCQBgWOi6XqvViEiWZWM+bnEgIAEA+MXFpIan3e2ey+UqlUqj0TC2yLKcSqUSiYR5Hxau+BKJRDKZ7K56ZxCQAAAGma7ri4uLrZGmVqtdv35dUZTNzU1jS7VaPfWEsiy730oiQkACABhshUKBRaN4PJ5KpSYmJohob2+vWCzWarVyuTwxMcF6PIlEgh9stre3iYidwQsISAAAA0vTNBZFotHo1taWsT2RSMRisdnZ2UajUSgUWEDiP4jLZrNENDc3Z37K5y4EJP90k1bnPHHufeu0pu4l2dn2NucvPSWaNFVqy8PiJ6aNd161Fdlae7dlGDozW/qZKWssOGlp0Lg14SzYublPiSbpM2O1q3HSbDtzarEVPSWaDJoqnbRmkY2PWw7+zLJmuQi2BLxD6+pnts/CmTfXdqSNN1mIXOt06fnyZ9bEy08sa5+9bv5ciodtEpii/OWDZzIZW1EgEJifn8/n841GQ1VVfr9nb2+vXC6Pjo62nsdFCEgAAH7xffTXo6MjttA23hgP6JrNJuckuq6vra0R0ebmpqe5eQhIAAADa2xsLBqNdooiuq47OUmhUGg2m/F4PBaLudo6OwQkAICBlUwmOW+G2AM9SZI4kUZRlO3tbUmSPH1Yx2DoIACAYVQqlcrlMhHNzc1xdisUCmyfUCjkdZPQQwIAGDq5XI5l38mynE6nO+1WKpWq1aokSZx9XISA1DX++F3rtGNa488qOm4tHO9UwjtOIbrO3TnQuajn96urRMum1a4yrThtGOfueddaadC2b8dUOudpdTZPrTvbznOWyUmdp/Yd2m8Sy9rhFWupec2WgMdp7J+Ixm3nNSendTXL8H9ZV/+es7M7uAMu/tEyH7AQWXYju8+XTxZ6P5aIwuEe26AoSjabZUM2yLK8s7PDyVMwukf+jDOEgAQA4JOT66fv4/DYyK+7PoOmadls1hiLYWlpKZ1OcyJNqVRqNBq+dY8IAQkAYBgUi0WWLEdE8Xg8k8mc+k5ob2+P7ezbMKwISAAAAy6bzbL8BVmWM5mMk+xtVVVZX8qjcVTbQkACABhkuVzOyKYzxlE9FesejY6OejdyXSukfQMADCxjLLuuohERVSoVIorH4161rB30kE7HH4Nunf5gLIcpaRlwzZa8NG49kjMZqG1PToqUQvRT7s7Oc8Gc77lqHz+vZ5w0ttZB3l4PftSplDtYnKPfordSiCbpj53Lv+xc1Dulm2tiWzWPdHcctKchWpjfCPyJ6Iq11JyhdzhuKeKmKK7TCK/YA7Ya14OWoe3OkgnpFd9bxOJKtz9rVVWVpeH52T0iBCQAgAG2v79PRJIksQTuTubn5805Dux5HRF5N7B3WwhIAAADi82E1Gg02IO7TqampswBiR0VjUa9bp4N3iEBAAws/jDenbD8Ou9mhu0EPSQAgIFVr9d9O+rsEJAAAPzi0w9MzysEJCfGravW9LhJU+k3RNdNaTTfJeue3LMGzYtdDA9nmxHV+eSkPfuG6Iptbs7OnFfKTyF7ys15C3SRdGcp4idivWhKcuvq6vFPyy81pwXaUgSPuW0wt/DQPDQi0eHkOHHYCjmV8K/BIbfUB9bm2UcCBOHhHRIAAAhB6ICkqqqmabYtDqc4BACA80XQR3a6ri8uLsqyrGmaLMsrKyvGFkVR0um0n8MrAQCADwTtIZXL5Wg0urm5ubu7y35pXKlU2Jb79+8bP9oCAICBIWgPyRhAyXhAJ8tyqVRSFKVWq/mfHQ8A4ALxBjMSysjJycnpe/WJoij5fD6ZTCaTSV3Xl5eXQ6FQrVZjW/jHRiKRg4MD85ZwzzMsAsA5ZP4LEA6H+/XbGkMkEqn//OD0/Rye7Z3+fyLXCdpDIqJcLtdoNLa2ttiAFoVCIZFIsDg0PT3t5B3Sum39wOmtYB4vlYjs+doL1sIfPF8MvxM5+M3zWyQ43nFebepmeFBOlvDTyLW/rd/uVOqRP0VWY/Wfdyp1K7nclvF8N3L3n+v/t1Mtzi+mc9ciT0v1v+ntnLavjLNqK3onovy8HutUaqNbSw9NA/TaMp4/s97Dx4fPDwz/KHJw2/p3zZxdf2itkj+46kdOB1e9Gw6/6fjfo3Pr65b/vL7y9ifG8tPINderA9cJGpDY1LlbW1v9bggAAPhE0IDEEr4XFhbY6u7u7vz8/MLCwtHRUa1W83mKDgAA8IGgAal1IqlQKPTxxx8rimIbJh0A4NxAUgOXoAGpEydTwQMAwHl0zgKSR9bpv6wbuFkMthlar5iW36HJ8efTbdpeg49zh/rq7T3809NO64U/WVvb85SsNoHTroC1UqfX9gwJDi+YTxs87iapIdhjUgNxP8tpuRLPR+7kZ3kc8+cy5uAnNVhzCtbXfZ9A1lbj2+vG4l2fmwI9EfSHsQAAMGwQkAAAQAh4ZAcA4BckNXChhwQAAEJAQAIAACEgIAEAgBDwDokZt6xdsRb+2LIWvGJJfZ2k53neT62rNvz8456Tp1+kL3s70EbnPt625WS7VWlXE66bt/AHBnRYdNrO4eChadV2Gtuq9eIFg9asa+uohty07xecp33bBqwzfym2PU/5bcAkr9CCexrRvHi+mgvoIQEAgCDQQwIA8Auy7LjQQwIAACEgIAEAgBAQkAAAQAhD+g5pnX5p3cAda3LSVnjYafVpNzldnLQ6/jCjXsyO2u15eq7UeTpc656WoU67uURd7HlIHVf5p+Qm3dkrNeXg2fLxiMLOs+yc41zbP7Wkhn422Tnrzla/O7eeV/wfdxjOaEgDEgBAHyCpgQuP7AAAQAgISAAAIAQEJAAAEAICEgAACAFJDUw3OVJczoek46TSdZVC5laWXVfcyrLjzEre4gW3Bqyz4I9Qd9y5qKUFvPMcdt659V47bNli7GtNyRsPWnblXMyW1ed3aWuWnZk9487WNv6M5n+wzmj+d37PaP7bkU+er4TDPtcOPUBAAgDwC7LsuPDIDgAAhICABAAAQkBAAgAAISAgAQCAEIY0qWGdfmZdfd1SfEic1eNxp+8luxqSjlPk1rhtfD2Pk8bnTxKgrfFd5OPxs+y6aYGtGqenbS067HweeyXW+yToPMvu+epdou/SH60t6jih7eH4OFl35a1CKyQ1cKGHBAAAQkBAAgAAISAgAQCAEIb0HRIAwBDSdb1WqxGRLMuBQODU/VVVbTaboVAoFAp53zoEJAAAv3iUOuRELperVCqNRsPYIstyKpVKJBKtO+u6XigUyuVys9k0ds5kMrFYzNNGIiAx1rG8Phu3rB5a1vi3VKDzeGI2PY/M5tZYdvzMNLf+5TjPf3O+Z1fOlFbnVtIdp7T1MjuvtHNqn23Uu8ngZ9YdzaUv2GZW/ZJeNNXAuw2Og0HOqs3uyRfPq4z8yLy6MHKJcyCcka7ri4uLrGNkVqvVrl+/rijK5uamebuqqsvLy+bQxXZeXFzc2NhIJpPeNRXvkAAABlmhUGDRKB6Pf/jhh/V6vV6vv/fee7IsE1G5XC6VSsbOuq6vra01Gg1JkjY2NtjOOzs7o6OjRJTP5zVN866pCEgAAANL07Tt7W0iikajW1tbExMTbHsikTDCTKFQMPY3otfOzo7RGYrFYru7u0TUbDbv3bvnXWsRkAAABpaiKGwhk8nYigKBwPz8PBE1Gg1VVdnGSqVCRHNzc0boYkKhUDweJ6Jqtepda/EOCQBgYB0dHbEFW4Bh2FM7ImLJC6qqsldHbTMdtra2vGrltxCQAAB84n+W3djYWDQa7ZThreuWCUX39vbYgtfZdJ0gIDG9J3SdYQpU3nkcFp1ayvkHwE+r8yfprhvfOP8stip7X+2ZR6O69ZwTaL1y45OHprWwLcvukMZNp+Rm2VlL+Tubv6NvfBnk8O3rz5fv/trr2gSVTCY5eXHsgZ4kSeYIZHSbVFWtVqvHx8djY2OxWMyHnyIhIAEADKNSqVQul4lobm6ObWHpDJIkqaq6trZmyxRfWlpKp9NOfk7bM6EDkqqqgUDAHJZZxqE/vxkGABhUuVyOZd/JspxOp81FmqYtLi42m01JklhvqVarNZvN7e3tarV6//5971olaEBiv+SSZVnTNFmWV1ZWiCibzRKRpmmJRMLTH2cBAHjh0sihsfzFyXjPxxJROBzurQ2KomSzWZa8IMvyzs6OrdPDijKZTCqVYlt0Xc/n8+VyuVar5XI59gfZC4IGpHK5HI1G2ceenp5eWVlhzzo3NzfZ77YQkADg3PnDyaSx3O07NPOxRJSMfNNt7ZqmZbNZI2+b8whuaWnJiEZEFAgENjc3NU2rVqvlcnnoAhJLeCdTEsj+/v7Y2Bj7RbEP2YcAAIOkWCwWCgWW3h2PxzOZDOfdh+0hHpNIJKrVarPZVBTFozQ8QQMSu1KKouTzeePS3Lt3L51OHx0dLSwssJ8N861bVzk93DC9Y93wDmfNtmpL3nknopzaMHf9KPICt9z5f0bx93xqXlmN/MnxaV3D/U9C20W42GG5O5F3enwqchaRVb8rbbmwv+6wfIqubsRvIs8fctzt9ekT3/p7B8ayNzWcG9lsluUv8MdIZb2l0dHRtt2moc6yy+VyjUZja2vLuArxeJw9qZudnXVyhnXb+sFBu73Ynn+wrE9ausb0S2vhFcs4la/TfxrLdyN3b9f/1ljtObHV+YHJyDe/qf+/3mrpKqPX7FrkqfljesR2Ea5Fnpbqf9OplLd6aD2v7dJ2Lo28E67/vOM9c8pXpHNLzf/YrZc58k64fvugU2l3bTCznWfcVOOPwv9d/8pc+AldMZb/SJZ/C4emcVepy7Rvs6eRa39bv21qzqGlsSM/c3gePnPad2RY076JKJfLGdl0tnFUbdhIQsYI3/4TdOigUqlki0ZjY2NGqe3HXAAA0JYxlt2p0Yi+Hc2h2Wy2HUHV2ChJktvN/AtBe0iqqmqatrCwwFZ3d3eTyeTs7Gwul6vVam2fbwIACM7VkRocPY1nY9NJktQ6ll2rRCKxtrbGRlBtzVxg4ziMjo62HYXIFYIGpLaR/P79+4qizM/P43dIAABO7O/vE5EkSeYhvVsZf1fn5ua2t7e3t7fHxsbMycxGeh4bj9UjggakTvo1whIAwHnERltoNBrswV0nU1NTLCCtrKxUq9Varba2tlYqlaLRKBEZU83Ozc2Z08Fdd84CkleC1iyG162lVyxrk9bpZV+0vY/1PpFhULl5BThn6mosuy9daoDtOY3zd6BdDV7neCy7lkLLkeYUg9MyX3ocR/GptZag9Yq8eLJgXv1yZJdTC/D1kKGws7OzvLzMwpJ59KClpSXvfoHEICABAAyser3e7SGBQGB3d1dV1b29PTaAUDAYjMfjQ532DQAA/TIxMeFd8kInCEgAAD5xNcvu6em7nDeC/g4JAACGDQISAAAIYUgf2a2Pn1jWr1iLf2wr/MS8+l36o3m155lee2Y9LX/8MAv/p08+VVdXjzMbr33n4w7Lp67qnUvP8mXyk+5sDjsXuTVECXf8pOA4Z9Zjy64vWstaJh3u2FzFWhrgDwR1Ykl7/ePIR51OazfsiavnD3pIAAAghCHtIQEA+E/ApxRCQQ8JAACEgIAEAABCQEACAAAhICABAIAQhiipYZ2ePV8Zt5ZZ87zHJw/Nq1fot5ZS7miqoo2R6vwlajcTyPb+E3Hned4te77ASft27apzksL5+eL8y2wr5beWM6Jrz6OpdpNoPj7+fL3nOWFbdz60/qsLdP42+atXTp6fNviZ9RJYxj0m6z9cISCpgQ89JAAAEAICEgAACAEBCQAAhICABAAAQkBAAgAAIQxVlt2F58vLlsFVg69bcnVeJ8vojT+wDq7KTQF6oecsu17T4b7h72xuT182g4jWAAAbzklEQVQyfPgXpIvxUulijy3o6gvhXKGuZhM/i0PHRc6z7Ma7qdG0as64o7Nl2XEmSudnrtpKLWvWEuu4x9wr2Sc6suy40EMCAAAhICABAIAQEJAAAEAICEgAACCEIUpqAADor2MK9LsJQhuigLT+X6bMuklLkW2S8h90M3id1UVOKT9DyXk6XM9TmNvaxk+C4pf2nEzIOc+pWXaWLV1NTM5tkAVnjnDbnoeOq+Brbap5y2HnorbHGvij59lWbclp5tNMWr+U8W6+d2sbrHfURUsuK3dKdd4qP63OvNrFPxToGzyyAwAAISAgAQCAEBCQAABACEP0DgkAoL8wHxIfekgAACCEYeohjZsWg4fmkklrss4Z5oTljbfGT2NzXEV3OzuvhZ9W11Wlzg/sbhLYw87ndSvLjjMPLH+KWBdHujvssNxVkt0xb9XihZaJVjnz5I53Pg+dMi+t9fu9aBkw0laLbc7cw86rPX/vICT0kAAAQAgISAAAIAQEJAAAEMIwvUMCAOgrZNnxoYcEAABCGKYe0vjzxUlrXtGLvafV2Tkfy86VPU+dR9WlWmy6mBi3i2y9UzOmOGdyK70q2HnVo7S61nlxDzsWupZEZqqCwt1k2dny387yn/vmASMPrUW21Z4/dhdTK4MQhO4hqaqqaZpto6Zpus4ZAhMAAM4lQXtIuq4vLi7KsqxpmizLKysrxvbZ2dnf//73/W0eAAC4TtCAVC6Xo9Eoi0PT09NGQMrn85Ik6boeCGBaEQA4Z5DUwCdoQIrH42zB/HSuWCyOjY1pmoZoBAAweAR9hxQKhUKhkKIoi4uL6XSaiFRVVVU1lUr1u2kAAOCJkZOTk9P36odcLtdoNDKZTCgUIqLl5eXR0VEiqlQqsVgsnU6z7Z1EIpGDgwPzlnA47GmDAUAo5r8A4XC4Xq/3sTFEFIlE/k/9jltnexq51vdP5DpBH9mVSqVGo7G1tWVsSaVSzWaTiKrVaiKRkCTp1JOsW1d3678xli3TJxP9mP7TvGpLCnee4hyJhOv1g9P3a6e3sVZPrdGLZ9bfj1z87/pXxmp3A606T2S2FkWS4fodx9fWjfzoyGq4/s+mGp2P+HlqAzqXRi6G66Y/o12lfXO+aU5RJGyp0b637cjxXltAltZHnobrf3PQtqjNqvMabatXzDU6PudA03W9VqsRkSzLAr77EDQgsYTvhYUFtrq7uzsxMcGWJUmKxWJ9axkAQK/6mNSQy+UqlUqj0TC2yLKcSqUSiYRtT+MPb1uZTMb4a+w6QQPS5uZmp6Ld3V0fGwIAcL6xX9GwjpFZrVa7fv26oijmv7e6rlerVc7Z2JMqjwgakAAAwBWFQoFFo3g8nkqlWP9mb2+vWCzWarVyuTwxMZFMJtnORtxaWlpqezb+y/szQkACABhYmqZtb28TUTQaNb+VTyQSsVhsdna20WgUCgUjIBmD4xi//vTTEAUk8zywQepi8KEzjPnWBd5pzzBsWpA3V2jPLvLGoLPpudTFOVid83/K0Zax7Dgjydk4fx3R3bU0rdtvn64axB+C8bPOpXz4aWk3FEVhC5lMxlYUCATm5+fz+Xyj0VBVlfWcjo6OiCgajfrcTkbQ3yEBAMDZsQBDRG0zEWRZZgvGmyGW9dCvBLwh6iEBAPSXi1l2LzjbbWxsLBqNdgowrQNVsy3sR5/+Q0ACABhYyWTSeD/Uij3QM/+WhqXYBYPBbDZbqVRYz2l0dNTJcARnh0d2AADDqFQqlctlIpqbm7MV5fP5crlsfo5XLpdnZ2dLpZKnTUIPCQBg6ORyOZZ9J8syGy+UTBkQRLS0tDQ1NSVJUrPZ3N/f397ebjaba2trsiwP3Q9jvWBODAucIZvKq6S7nk9jO7CriU175jwVjHOg86IeKjrjUa34o9TwL3vn1bO0jv/Nu3IeTgJea43Hzj/MRd4Ht38Qtz6YAA5HLhnL4ydf9HwsnWFwTkVRstksS16QZXlnZ8d4w2S8UtrY2DA/6IvFYhMTE9evXyeifD7v3egEQxSQAAD6K3jyzFjuNlyajyUiiny/29o1Tctms8ZADEtLS+l02pzvkEgkWkcSMooURSmXy9Vq1bsZ6fAOCQBg8BWLxdnZWRaN4vH4w4cPV1ZWuoorRuJD6yhEbkEPCQBgwGWzWZa/IMtyJpPpbXxqH36chIAEADDIcrmckU3HGbdaBAhIAAADyxjLzkk0Wl5e1nV9amqq7dzcxjB33v0aaYgCkjkdzpYa59HwdKc4S75Zb6ftKlvJozHozrJzb/xJO+wm6Y6fouhWumXPzOdxPOXe6aW/sq5uUUdvc2u92fnyvT3JbZAAXByp4aKz3SqVChFJktQ6ll1b1Wq1VqvNzc21PqDb29sjIlmWEZAAAKBr+/v7RCRJUqFQ4Ow2Pz8fCoWSySQbnWFtbW1jY8OcDp7P51lCRNvOk1sQkAAABhbLiGs0GuzBXSdTU1OhUCgWi83NzZXL5UqloihKPB5nP4w1xhCam5vrlBfuCgQkAICB1e0Er5ubm5IksXGDWCoEI0lSOp32tHtECEgAAAOsXq93e8jKyko6na5UKkdHR7VaTZblsbGxeDyOtG8AAPBbIBDgjBHuHQQkAAC/iDe8nlCGKCC5ldvdxXn8GXXU+azSLiZkn99/V12l3XJ29iybvOczuZVQzDkPJ1fbRTd7vQQ3OdOi9zgSKfgKY9kBAIAQEJAAAEAICEgAACCEIXqHBADQZ+f35asv0EMCAAAhDFEPKUi6abmLwVW7S89zvq9baXX+TGHOqdSf0Uvd4laWHd8ZMjGdXz/Ohe+q4b86fZf2+KOgctwlenu8Y+nNwx5PC+cdekgAACAEBCQAABDCED2yAwDoM8EfaPcbekgAACAEBCQAABDCED2yCziewvxMaXVeTEB9ljHovM8TE+IphPNB57q6IOYpsW0f85B7oK3UWumx6VQuTmHufNjC/+y1CntqXM8DHl5EKh20gR4SAAAIAQEJAACEMESP7AAA+kyE59sCQw8JAACEgIAEAABCGKJHdr3PGHuue9k9D3DmFucj3Z3avPGejuVm2R1zaz0ef15qv3/4yXtfWlc7z2TaVZadR6MG/qBzLfarc9h7g2yFP+58Vn7iqk4wsITuIamqqmkafwsAAAwGQXtIuq4vLi7KsqxpmizLKysrrVv63UYAgC6d68ct3hO0h1Qul6PR6Obm5u7ubqVSabsFAAAGiaA9pHg8zhZ0Xe+0BQAABomgASkUChGRoij5fD6dTrfdAgAAg2Tk5OSk321oL5fLNRqNTCbDQlHbLRyRSOTg4MC8JRwOe9JQABCS+S9AOByu1+t9bAyxP0q3XWtDeDXS90/kOkF7SKVSqdFobG1tcbacat26+mI9ZiyPWxNNJ61pufYEX8cjSEZ+FK7/5qBtURucfGPHA7ZGroXrdw46FjuvkVtLd5W6NUe4tSjy/XD9v51WyknmPiXP21SajHxTqv+NqQbugLyH1sv3kfW81knCD023my0hfDUc/vmB0y/UlYFz32mpkZP2bVu1DdL6dudabKs/CofvmCrtatza3rzvwTnBdYIGJJbevbCwwFZ3d3dbt/SpaQAAvUKWHZegAWlzc/PULQAAMEgETfsGAIBhg4AEAABCEPSRnRfM76ttb7Ztq6eMetfzLJk962oIM+e5Em7hv1vvubSl6Djo9C1+z0kNvZ2TiILBLu4Zt74Ht2Yn5qx2Nd3uTeuqbYg8TuYC/4PYfngY4O4M59oQBSQAgD5DUgMXHtkBAIAQEJAAAEAICEgAACAEBCQAABDCECU16CPvPV85ed1cdKakOy/w0+q6ynzyAb955Lj0DB/Eee7cf450nr3VJhzuYmeun1pXjzsst25x67t1Kx/PxpbwZmvtoeMzezSAv4jzAojYJoGghwQAAEJAQAIAACEM0SM7AIAhp+t6rVYjIlmWAwGnPzJmR4VCISdT/5wFAhIAwODL5XKVSqXRaBhbZFlOpVKJROLUY5eXl6vV6tLS0srKipdtREACABhouq4vLi6yjpFZrVa7fv26oij8uRSKxWK1WvWygc8hIJ2ZDyl4XY1l5xaPajnDabdGOAcLPSTLoXWVn2XXac9WznPw/MnW6+o76DndzHzg+RvXrh/3aaFQYNEoHo+nUqmJiQki2tvbKxaLtVqtXC5PTEwkk8m2x6qqms/nfWsqkhoAAAaWpmnb29tEFI1Gt7a2WDQiokQisbOzMzo6SkSFQqHT4Wtra/60k0FAAgAYWIqisIVMJmMrCgQC8/PzRNRoNFRVbT02l8vVarVoNBqNRr1uJ4OABAAwsI6OjtiC0Tcyk2WZLTSbTVuRoijb29uSJPk5WzfeIQEADKyxsbFoNNopw1vX27/L03V9eXmZiDKZjNep3mYISAAAfvE9qSGZTHZKWKBvH+hJkhSLxczb19bWms1mPB7nHOuFIQ1Ifxz5yLw6fjJpXrWPZWebDLSrodv6zqPENC8GR7NeyZt/Zy0N91pjP3zXusqfobXv/Jlk2ItR3DAy3FmUSqVyuUxEc3Nztu2VSkWSpI2NDZ+bNKQBCQBgmOVyOZZ9J8tyOp02tmuaxvK8NzY2nA/l4BYEJAAAv+yOPF9eOOn9WCIK9/jQQFGUbDbLhmyQZXlnZ8cceJaXl5vN5tLSkpMRHFyHgAQA4JdugxDnWCXS7Qk0Tctms8awC0tLS+l02hyNWJ736Oiouc/kJwQkAAC/9O/9YbFYLBQKLL07Ho+3ps+pqsoe4m1tbfn/sI5BQAIAGHDZbJblL8iynMlkbDl1jPFTpH/8x39se5Lt7W0WsbwbZRUBCQBgkOVyOSObzs9fufYAAYlIhDnLXeRWnrfzcTT5Y786Xr35M2cNE9J3uaV9H6hW8N8mgHeMsexOjUaxWKxer7ctWlhY8Gf6CQwdBAAwsCqVChFJktQ6lp2A0EMCABhY+/v7RCRJEmdIbyKan5/3c4igThCQAAD84vvbADYTUqPRYA/uOpmamhIhIOGRHQDAwGodxltk6CEBAAysTnkKXdnd3T37SZxAQCIi+u3IJ+bV8ZNx3t78pLKgsyIXOU+H82g01a7mtf6ym2PPLRen+nZFX9qDkU+hW3hkBwAAQkAPCQDALwP6SMAt6CEBAIAQEJAAAEAICEgAACAEvENq41cjh+bVt3uewaQvI4g5ny7bnynMbblWtgMPz9CGvuIPXmeDFwcATggdkFRVDQQC5t8Pq6oaCoX6NVcHAMCZuPjfJi+4dyphCBqQdF1fXFyUZVnTNFmW2RCzCwsLoVCoVqt1ms8DAADOL0HfIZXL5Wg0urm5ubu7y0ar3dvbC4VCm5ubW1tb/FECAQDgPBK0hxSPx9mCrv/lFYSqqhMTE0QUCoWMOeEBAGBgCBqQ2HsjRVHy+Xw6nTZvJKJoNOrkJOvW1XA43Ftj7ka62Dnyox5r6Vnkmt81ElHkfQ8qPe2Ud3v9BnvGqfGuZ5W+7/vH9L9G8uXbPDg4MJb78Amhe4IGJCLK5XKNRmNra8uIQ5qmsQWHPaR127rp7jyLt591LIp8P1z/Ta+1cAbBszGVRpLh+p2DTqX21a5eqHYeiSxyN1x/0/HHtFX6WcfSm7Yiq7vh8JsufYMO8Ws8yxSxnNL3w+F/9vdjnqXGnker68O36Wdl0CtBA1KpVGLRyNgyMTGhqioRsTSH/jUNAKBXLmbZXXTvVMIQNCCpqqpp2sLCAlvd3d1NJBLFYjGXy1Wr1VQq1dfWAQCA+wQNSJubm60b79+/ryiKIFPtAgCAuwQNSJ3g50cAAINK0N8hAQDAsDlnPaRB5nzQubOscvBTpno+7SGvkJ9ZJ5QfWFfP8nKaP8ah83RLH2DWV5chqYELPSQAABACAhIAAAgBAQkAAISAgAQAAEJAUgMAgG/6nqciNASkrt28YFl9+wtrsfP7zf/5ZLvKo+OXupUTeH781rra1YyxgkMqHQgCj+wAAEAICEgAACAEBCQAABACAhIAAAgBSQ0AAL45t1k9vkBAOqubl0wrYe79Zkurc35nnmXwOr1zUc/TmtpKuznPORq8rivn688M0upATHhkBwAAQkBAAgAAISAgAQCAEPAOCQDAN+frbaPf0EMCAAAhoIfko57/26jntDripsO5lWV3yDvw5qD8F6GLg9ed5cI7hDw6OI/QQwIAACEgIAEAgBDwyA4AwDeD8gjbG+ghAQCAEBCQAABACHhk57Kbf+90z7f/w/FJ/enl27LjPum8Z9hS+rYnremPFzsXneVL8OcLRGYdnErX9VqtRkSyLAcCAf7Oqqo2m01JkiYmJnxpHQISAMAQyOVylUql0WgYW2RZTqVSiUTCtqemaYVCoVwumzfOzc2l0+lQKORpI/HIDgBgkOm6Pjs7u729bY5GRFSr1a5fv57NZs0bVVWdnZ21RSMiKpfLs7Ozqqp62lQEJAAA3xy79z+nCoUCe0wXj8c//PDDer1er9ffe+89WZaJqFwul0olY+e1tbVms0lEmUzmd7/7Xb1e/93vfpfJZCRJajaby8vLbl8QCwQkAICBpWna9vY2EUWj0a2tLeNtUCKR2NnZGR0dJaJCocA2qqrKQlcmk0mlUuwlUyAQSKVS6XSaiBqNxt7ennetRUACABhYiqKwhUwmYysKBALz8/NE1Gg02LM4I9ikUinbznNzc2zB06d2CEgAAAPr6OiILbTNlGNP7YiIPaZjL5mi0Wjrnqem5LkCWXZ9c/MnHYve3rGu8x8Xc7N9b77nvEU9uul5DX3DyQIHOBfGxsai0WincKLrlj8fqVQqmUxKksTZMxgMut5IAwISAIBv/B46KJlMJpPJTqXsgZ4kSbFYjDr0oph8Ps8W4vG42218DgEJAGAYlUollt5tvB9q3YE98TN+wLSxseHpT5EQkAAAhk4ul2PZd7Isswy6Vnt7e9Vq1Vh97733Wn9F6y4EJAAA3yyalnc67nX6sUQU7q0FiqJks1nW45FleWdnp9MbJiPlgYWl69evq6q6srLSW71OICABAPim2yDEOfadbo/XNC2bzRqdnqWlpXQ6zUmfM2KPruv5fL5cLrNOlXcxCQFJRDet/yX09i+tpT8zrYTp5rr3DYKzsaUlYUqcIda3L79YLBYKBZbeHY/HM5mM87dBgUBgc3NT07RqtVoul70LSEL/DknXdU3TzFtUVbXlKQrl4OBg4GscnkrxMQes0mGWzWbz+Xyz2WTP6La2tnrITWBP8JrNpnd/hIUOSOVy+d69e2yZjQ9YKpXY//e3YQAA50UulzOy6e7fv88yvNuKRCKRSKTTH1gjKZwNL+QFcQPS7OyskflORJVKJRqNbm5u3r9/39PBlAAABoYxlt3c3Nzm5qaTQ4zBHWyMQYM4Ie2MxA1I9+/fX1paMlZlWa5Wq4qilMtlI/cDAAA4KpUKEUmS1DqWXSv2o1d2SKdTsfFYPXJukhpCoZAkSXt7e7VajfPDY7O74R7TInsWJrrrwWnvvt9Sjcc18g1JpeEe0pjcqPT90/dyucZh+DYFEQ77/PXS/v4+EUmSZAzp3db8/HwoFIrFYuw3sNls1tadMjLFO/1oyRUjJycn3p39jHK5HH2bYpjL5cbGxlgomp6e/vjjj/vcOAAA4b366qsss45vZ2eHPYhbXl42ekKxWIxNg6QoCotG8Xh8a2vLu9aemx4SAAB0y0k0MtvY2BgdHWXTy9rmjV1aWvL0V7F0jgLS/Pz8wsLC0dFRrVbzdHQ/AICBUa/Xu9o/EAisrKzMz89XKpXj4+NarSbL8tjYWCwW83QUO0boR3atFEUJhUI+XBcAAPDZOQtIAAAwqMRN+wYAgKGCgAQAAEJAQAIAACEgIAEAgBAQkAAAQAgISAAAIAQEJAAAEAICEgAACAEBCQAAhDCYAck207kPE5+3VqFpmm3+dR8qZfV6+mFbr63/H9P/L3R4biGv75/WSgf1FoIe/NX6+nq/2+AmXdf/6Z/+Sdf1f//3f5ck6eLFi2+++eZXX32Vz+f/+q//2piC17sav/Od7xBRNpt9/PjxRx999D//8z++Vcq2/8M//MO//Mu/uF5ja6XGtf3oo4/+/Oc/T01NeV0j+5gLCwuaphWLRY9GNdR13XzPhEIhH26h1ip8uIVaK/X0/mmt1Li2nt5CrR/T61sIencyWN5999179+6dnJwcHx//5Cc/uXfv3rvvvstWf/rTn/pQ48nJyf7+/i9+8Qu25V//9V/9qZT5xS9+8cMf/vD4+NiHSguFAru2JycnP/zhD32o8eTk5Ne//jW7tkdHRx59obZ7xodbqLUKH26htp/L0/untVIfbqHWj+nDLQQ9OzfTTzgnSRIRBQIBNnB6qVRSFIUt+1AjEe3v74+NjZVKJSLybjIrW6VEVCwWx8bGNE0LBAI+VGp8NE8ffdg+pqqq7D9yQ6FQtVr1okbbPePDLdRahQ+3UGulPtw/tkqNeWS8u4VaP6YPtxD0bNDeIU1NTRWLxb29veXlZUmSjInP9/b2xsbGfKiRbbx37x4RHR0dLSws+FOpqqqqqqZSKS+qa1spe9yhKMri4qJHsxq3vbbGM5ZoNOpFpbZ7xodbqG0VXt9Ctkp9uH9aK/XhFmp7bb2+haBnA9JDUhSF/UdlMpnc2NjY29tLJpOKohQKhUQiYUx8zhY8rZGVxuNxVtfs7KxbNfIrLRaLo6OjuVxO07RsNptOp115OG6uMRaL2SrN5XKNRmNra8vdB/H8a2u8AHfxP2/NNdrumaOjI69vobZVeH0L2SqdmJjw4v4h7rVNJpNe30Jtr60XtxC4YkACUiwWYxPC7+3tNZvNlZUVTdPYFv9rZP8M2G7uPovgVJpKpdhcxdVqNZFIGP0Jt2psrbRUKrE/Ja5U1LbS1o85MTGhqioRaZrm4gM088f0jVFpLpezFflwC9kq9ej+Ie619eEWar22Ht1C4IoBCUiGWCw2Ozt7dHRUrVY3NjYCgYDXE5/baiSiZDI5Ozuby+VqtZpHDyJaKzUSsSRJ8uhvq63SUqmkaZrxQGl3d9frGokokUgUi8VcLletVj16vjQ/P2++Z2yrPtRIvtxCtkp9uH9aK2UJ357eQq3X1odbCHo2gDPG6rrOXmAa72a9nvi8tcZ+Veo1/yvty7VtrcL/GvtVqQ/8r7Qv1xZ6MIABCQAAzqNBy7IDAIBzCgEJAACEgIAEAABCQEACAAAhICABAIAQEJAAAEAICEgAACAEBCQAABACAhIAAAgBAQkAAISAgAQAAEJAQAIAACEgIAEAgBAQkAAAQAgISAAAIAQEJAAAEAICEgAACAEBCQAAhICABAAAQkBAAkEpihKJRNw61cLCAhGpqnrqzrlcLpfLuVIvAHQFAQkGnyzLmUyGiPL5fL/bAgAdISDBeZLL5aanp6enp1knRlXVbDabzWYjkcjy8rKu68Zur7766sLCQrFYVFVV07S9vb1cLler1XK5nKqqRh+oVCqxbhM7ZHZ2ttFoGEXT09OvvvoqOkwA/vj/zKkO7nRE79QAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<IPython.core.display.Image object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "% рисуем\n",
    "figure\n",
    "axes('fontsize',18)\n",
    "% contourf(x,y,aam,30) % 1 контур\n",
    "imagesc(x,y,aam) % 2 рисунок\n",
    "% imagescn(x,y,aam) % 3 без нанов\n",
    "axis xy\n",
    "c=colorbar('fontsize',28)\n",
    "xlabel('Longitude')\n",
    "%ylabel\n",
    "title('Average SST in 2013')\n",
    "grid on\n",
    "colormap(jet)\n",
    "c1=prctile(aam(:),3); % автоподбор по шкале\n",
    "c2=prctile(aam(:),97); % автоподбор по шкале\n",
    "caxis([c1 c2])% ограничение шкалы\n",
    "savname='SST_Reinolds_7D_year2013.jpeg'; % имя рисунка\n",
    "print(gcf,'-r150','-djpeg',[picdir savname]) % сохраняем рисунок"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "'shading' is a function from the file /usr/share/octave/4.2.2/m/plot/appearance/shading.m\r\n",
       "\r\n",
       " -- shading (TYPE)\r\n",
       " -- shading (HAX, TYPE)\r\n",
       "     Set the shading of patch or surface graphic objects.\r\n",
       "\r\n",
       "     Valid arguments for TYPE are\r\n",
       "\r\n",
       "     \"flat\"\r\n",
       "          Single colored patches with invisible edges.\r\n",
       "\r\n",
       "     \"faceted\"\r\n",
       "          Single colored patches with black edges.\r\n",
       "\r\n",
       "     \"interp\"\r\n",
       "          Colors between patch vertices are interpolated and the patch\r\n",
       "          edges are invisible.\r\n",
       "\r\n",
       "     If the first argument HAX is an axes handle, then plot into this\r\n",
       "     axes, rather than the current axes returned by 'gca'.\r\n",
       "\r\n",
       "     See also: fill, mesh, patch, pcolor, surf, surface, hidden,\r\n",
       "     lighting.\r\n",
       "\r\n",
       "Additional help for built-in functions and operators is\r\n",
       "available in the online version of the manual.  Use the command\r\n",
       "'doc <topic>' to search the manual index.\r\n",
       "\r\n",
       "Help and information about Octave is also available on the WWW\r\n",
       "at http://www.octave.org and via the help@octave.org\r\n",
       "mailing list.\r\n"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "shading?"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "c = -7.1088\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjAAAAGkCAIAAACgjIjwAAAJMmlDQ1BkZWZhdWx0X3JnYi5pY2MAAEiJlZVnUJNZF8fv8zzphUASQodQQ5EqJYCUEFoo0quoQOidUEVsiLgCK4qINEWQRQEXXJUia0UUC4uCAhZ0gywCyrpxFVFBWXDfGZ33HT+8/5l7z2/+c+bec8/5cAEgiINlwct7YlK6wNvJjhkYFMwE3yiMn5bC8fR0A9/VuxEArcR7ut/P+a4IEZFp/OW4uLxy+SmCdACg7GXWzEpPWeGjy0wPj//CZ1dYsFzgMt9Y4eh/eexLzr8s+pLj681dfhUKABwp+hsO/4b/c++KVDiC9NioyGymT3JUelaYIJKZttIJHpfL9BQkR8UmRH5T8P+V/B2lR2anr0RucsomQWx0TDrzfw41MjA0BF9n8cbrS48hRv9/z2dFX73kegDYcwAg+7564ZUAdO4CQPrRV09tua+UfAA67vAzBJn/eqiVDQ0IgALoQAYoAlWgCXSBETADlsAWOAAX4AF8QRDYAPggBiQCAcgCuWAHKABFYB84CKpALWgATaAVnAad4Dy4Aq6D2+AuGAaPgRBMgpdABN6BBQiCsBAZokEykBKkDulARhAbsoYcIDfIGwqCQqFoKAnKgHKhnVARVApVQXVQE/QLdA66At2EBqGH0Dg0A/0NfYQRmATTYQVYA9aH2TAHdoV94fVwNJwK58D58F64Aq6HT8Id8BX4NjwMC+GX8BwCECLCQJQRXYSNcBEPJBiJQgTIVqQQKUfqkVakG+lD7iFCZBb5gMKgaCgmShdliXJG+aH4qFTUVlQxqgp1AtWB6kXdQ42jRKjPaDJaHq2DtkDz0IHoaHQWugBdjm5Et6OvoYfRk+h3GAyGgWFhzDDOmCBMHGYzphhzGNOGuYwZxExg5rBYrAxWB2uF9cCGYdOxBdhK7EnsJewQdhL7HkfEKeGMcI64YFwSLg9XjmvGXcQN4aZwC3hxvDreAu+Bj8BvwpfgG/Dd+Dv4SfwCQYLAIlgRfAlxhB2ECkIr4RphjPCGSCSqEM2JXsRY4nZiBfEU8QZxnPiBRCVpk7ikEFIGaS/pOOky6SHpDZlM1iDbkoPJ6eS95CbyVfJT8nsxmpieGE8sQmybWLVYh9iQ2CsKnqJO4VA2UHIo5ZQzlDuUWXG8uIY4VzxMfKt4tfg58VHxOQmahKGEh0SiRLFEs8RNiWkqlqpBdaBGUPOpx6hXqRM0hKZK49L4tJ20Bto12iQdQ2fRefQ4ehH9Z/oAXSRJlTSW9JfMlqyWvCApZCAMDQaPkcAoYZxmjDA+SilIcaQipfZItUoNSc1Ly0nbSkdKF0q3SQ9Lf5RhyjjIxMvsl+mUeSKLktWW9ZLNkj0ie012Vo4uZynHlyuUOy33SB6W15b3lt8sf0y+X35OQVHBSSFFoVLhqsKsIkPRVjFOsUzxouKMEk3JWilWqUzpktILpiSTw0xgVjB7mSJleWVn5QzlOuUB5QUVloqfSp5Km8oTVYIqWzVKtUy1R1WkpqTmrpar1qL2SB2vzlaPUT+k3qc+r8HSCNDYrdGpMc2SZvFYOawW1pgmWdNGM1WzXvO+FkaLrRWvdVjrrjasbaIdo12tfUcH1jHVidU5rDO4Cr3KfFXSqvpVo7okXY5upm6L7rgeQ89NL0+vU++Vvpp+sP5+/T79zwYmBgkGDQaPDamGLoZ5ht2GfxtpG/GNqo3uryavdly9bXXX6tfGOsaRxkeMH5jQTNxNdpv0mHwyNTMVmLaazpipmYWa1ZiNsulsT3Yx+4Y52tzOfJv5efMPFqYW6RanLf6y1LWMt2y2nF7DWhO5pmHNhJWKVZhVnZXQmmkdan3UWmijbBNmU2/zzFbVNsK20XaKo8WJ45zkvLIzsBPYtdvNcy24W7iX7RF7J/tC+wEHqoOfQ5XDU0cVx2jHFkeRk4nTZqfLzmhnV+f9zqM8BR6f18QTuZi5bHHpdSW5+rhWuT5z03YTuHW7w+4u7gfcx9aqr01a2+kBPHgeBzyeeLI8Uz1/9cJ4eXpVez33NvTO9e7zofls9Gn2eedr51vi+9hP0y/Dr8ef4h/i3+Q/H2AfUBogDNQP3BJ4O0g2KDaoKxgb7B/cGDy3zmHdwXWTISYhBSEj61nrs9ff3CC7IWHDhY2UjWEbz4SiQwNCm0MXwzzC6sPmwnnhNeEiPpd/iP8ywjaiLGIm0iqyNHIqyiqqNGo62ir6QPRMjE1MecxsLDe2KvZ1nHNcbdx8vEf88filhICEtkRcYmjiuSRqUnxSb7JicnbyYIpOSkGKMNUi9WCqSOAqaEyD0tandaXTlz/F/gzNjF0Z45nWmdWZ77P8s85kS2QnZfdv0t60Z9NUjmPOT5tRm/mbe3KVc3fkjm/hbKnbCm0N39qzTXVb/rbJ7U7bT+wg7Ijf8VueQV5p3tudATu78xXyt+dP7HLa1VIgViAoGN1tubv2B9QPsT8M7Fm9p3LP58KIwltFBkXlRYvF/OJbPxr+WPHj0t6ovQMlpiVH9mH2Je0b2W+z/0SpRGlO6cQB9wMdZcyywrK3BzcevFluXF57iHAo45Cwwq2iq1Ktcl/lYlVM1XC1XXVbjXzNnpr5wxGHh47YHmmtVagtqv14NPbogzqnuo56jfryY5hjmceeN/g39P3E/qmpUbaxqPHT8aTjwhPeJ3qbzJqamuWbS1rgloyWmZMhJ+/+bP9zV6tua10bo63oFDiVcerFL6G/jJx2Pd1zhn2m9az62Zp2WnthB9SxqUPUGdMp7ArqGjzncq6n27K7/Ve9X4+fVz5ffUHyQslFwsX8i0uXci7NXU65PHsl+spEz8aex1cDr97v9eoduOZ67cZ1x+tX+zh9l25Y3Th/0+LmuVvsW523TW939Jv0t/9m8lv7gOlAxx2zO113ze92D64ZvDhkM3Tlnv296/d5928Prx0eHPEbeTAaMip8EPFg+mHCw9ePMh8tPN4+hh4rfCL+pPyp/NP637V+bxOaCi+M24/3P/N59niCP/Hyj7Q/Fifzn5Ofl08pTTVNG02fn3Gcufti3YvJlykvF2YL/pT4s+aV5quzf9n+1S8KFE2+Frxe+rv4jcyb42+N3/bMec49fZf4bmG+8L3M+xMf2B/6PgZ8nFrIWsQuVnzS+tT92fXz2FLi0tI/QiyQvpTNDAsAAAAJcEhZcwAACxMAAAsTAQCanBgAAAAddEVYdFNvZnR3YXJlAEdQTCBHaG9zdHNjcmlwdCA5LjI2WJButwAAIABJREFUeJzt3U9sG2d+N/Af3xaLnob0JYdC1NIJljyMgHXjXRSkDnUhBUueNlFjUc1lJUBk+p7kBUgrL1AqsoSidihgI6DAwqQAqSjQhAzippeQhR1s3oOHeHftrAsMC1CLOFoNsYc9xBoe9sr38MSjh0POaPhnho/I7wdBMH+feTii9PPzzG+ex9dutwkAAGDc/te4KwAAAECEgAQAAIJAQAIAACEgIAEAgBAQkAAAQAgISAAAIAQEJAAAEAICEgAACAEBCQAAhICABBPo0aNHDjcCgDgQkGAEXrx44fP5fD7f8+fPx1uTzc1Nn8/3xhtv+Hy+xcXFFy9e9Ny4uLjo67S5uWkU8ujRI5/P5ySAbW5uXrly5cqVK8bp3VuIiF3x6dOnNscAALUBhnb//n32dbp79+7Yq3H37t2HDx/evXuXiG7fvt1z45MnTx4+fJhKpYxdX3/9tVHOt99++/Dhw2+//dbh5VixrBx+Cyv29u3bV69eZavtdpvtun///u3bt4moXC67e1MALg+0kGAEyuXy66+//vrrrxuR6ebNm6+++ipbfvXVV2/evEkvWypGc4Q1RFhzgV42I9gBRuOmUCj4fL5XX301nU4vLi6yAk3lGL7++msiWlxcXFxc3NzcvHv37muvvdZz4/Xr1xcXF9l12bJRWyJ6+vTpG2+88fTpU1bDdDrNLlcoFLovt7m5ubm5GQgEHj58yAIS20JEbJVd3TjrypUrqVQqnU6/++677FOP5EcAMAnGHRHh0mN/l41mwZMnT9ovWw9Pnjxhe+/fv8/2lsvlcrlMXHvi6tWrt2/fbr9sOhgH3L9//8mTJ0SUSqXYloWFhe7D+DbZkydPAoEAEb3++uupVOr+/ftWGxnWRmENFx6rmBFg3n777YcPH77++utWvy+snnfv3l1YWLh69Srb+Prrr7MK8wXyl2Ah/MJ2GMD0QAsJhvXJJ58Q0c2bN1kziMWJdDodCATK5TJ7EnPz5s2HDx9evXr1yZMnLEIYrYd79+7du3ePiDY3N9leFsy+/vprVtS9e/du3rz59ttvs+OtyiGi69evP3369O7du9evX3/06NG7776bTqd7buzrA7777rusgWX18RcXF1OpFGsVsVYXv2Dl1Vdf/eqrr5BqAWBAQIJhsbDx2muvvfbaa/QyPhHRzZs3Hz169PDhw7ffftv01zmdTi8vL7NlY9e9e/d+9KMfERHryzKwA/gutZ7lsEvfv39/c3OzUCg8f/58YWGhWCz23Diij06bm5vLy8vvvfce6827fv36V199xXZ98cUXPev89OnTp0+fLi4ufvLJJ1evXmV3DwAIAQmG9PTp06+++ur27dusg+v27dvffPMNSyd74403WAvgjTfeIKLr16+/ePHivffee++9954/f9797IQ1fdhetoVFuHv37j1//tx4fmNTzvPnzz/44IN79+69ePHi6dOnz58/v3r1as+NI/nshULhgw8+YN1uLOAZIZm1e1h8NSmXyz/60Y/YA6pvvvmmZ9ACmFLj7jOEy409hjFS1NgTI/ZMqN1us4c37DHJt99+yx7DEBH7I256smK0FRYWFoxC+FPYI5nucozKfPvtt0bPHhGxnr2eG/nKX/gMiR3ADuYPY0l6BlY9Vnlj1VQgu0VGRMQzJACer40pzMFDjx49unLlyvXr13vuZS0eY+/z58+NLXfv3n3+/LnxxMWmHHaWaW/PjS5hDUT7C9nfB4DphIAE4nr+/Pn169evXLmyuLhYLBbv37/fbz4CAFwiCEggtBcvXnzyyScvXrzgX2wCgImEgAQAAEJAlh0AAAgBAQkAAISAgAQAAEJAQAIAACGIFZAURdF1nd+iqqqmaeOqDwAAeObPtre3x10HIiJd19955512u/1P//RPr7zyyg9+8AO25Y9//ONnn332u9/9bn5+ftx1BAAAF4mS9s3Gu0ylUpqm5XK5o6OjYrF4dnaWzWaJaGFh4Ysvvhh3HQEAwEV/Pu4KfMcYFkxV1WAwSETxeJxtMXXiAQDARBIlIDH5fL5arbLhYVhYUhRlb29vgAFjfD5fz+3hcHjISgKA4I6Pj7s3jr03KBKJjLbARqMx2gLHTpQuO4Ou64uLi7/5zW+IKJ/PN5vNTCbDglNffD7hPpr3cBMIN4GIcBPEuAORSOSdXpFyMP8eDk9eQBIlyy6Xy1UqFSLy+/1sS6lUajab+/v7A0QjAAC4dETpsksmkxsbG6qq1ut11kHHEr5XV1fZAUdHR2Os3oW2LXoIx07YinkJN4Em6CbcoV/0e0q7fcuNmsDIiRKQ5ubmHjx4UK/XV1ZWWJNod3d33JUCAADviBKQiMjv98disXHXAgAAxkOggAQAMNn8466Aruv1ep2IZFk2Hthb0TSNDZTjWVNhYgMS0ruJ3YTRZfVcUrgJhJuAPwgvX6ppNpvGFlmWU6lUIpFwcvDy8nI6nXY7xUyULDsAAHCDrutLS0sHBwd8gCGier1+69atXC7n5OByuby0tKSqqqtVRUACAJhkhUKBddPF4/FPP/200Wg0Go0PP/xQlmUiKpfLpVKp++Dl5eVHjx6xg3d2diRJarVaa2trrg6dg4AEADCxNE07ODggomg0ur+/Pzc3x7YnEonDw8OZmRkiKhQKxvHs4Hg8vru7a3TQJZPJ/f19Imq1WtVq1b3aIiABAHgkMLr/HFIUhS1kMhnTLr/fv7KyQkTNZpP1xRkHG0OJGoy8htPT0wE+uEMISAAAE8uIH0bbiMd67Yio1Wrx220S8ExHjhYCEgDAxJqdnY1Go90tHsb0QMjoo3v8+LHpSGOi1NnZ2VHX8RwCEgDAxEomk0dHR+wJUDfWRydJEuuRCwaDLHSVy2Wj+46IdF3f2NhgR1rFtpGY2PeQAADARqlUKpfLRLS8vGxs3NnZIaJqtbq2tibLcjQardfr9Xq91WpJknR4eOjqq0gISAAgugEGVCV6c/T1mCD5fJ4l1MmyzE845/f7M5mMoiitVouFImNXJpPp+SBqhBCQAAA8ssYt/0ef577VuTrwyBOKouRyOfbeqyzLh4eHfApDqVTa29tj7aF4PC5JEhGxURu2trZUVXV12GsEJAAAj/QbhGzO3ey/BE3TcrlcrVZjq+vr6+l0mo9GmqaxaGQKVNlstlgs7u3tlctlSZKy2eyAn+EiSGoAAJh8xWJxaWmJRaN4PP7o0aNsNmtK7y4UCsazItOuVCrFHjUdHBy4N1gDAhIAwITL5XJ808dqJm42CkM8Hu/5HpLxbiz/YGm00GUHADDJ8vm8kU1n/wSIvfTKnht1u3C6iuEhIAHAJFoNjbsGPTgf8mdUjLHsLoxGRDQzM9NsNk3jfBtcHVaVQZcdAMDEYr1wkiR1j2XXjXXKKYpijMvQsyj35utDCwkAYGKxQYAkSeKH9O62srISDAaTyWS5XG61WhsbG6bnTGzKPup8i3bkEJAAACYWS0BoNpus487K/Px8MBicm5vLZDJ7e3v1en1xcTEej7P5KYzZY01v0Y4cAhIAwMTqd3DuVColSRJLyTNNfbS8vJzJZFxNbUBAAgDwiPdJDY1Go99TkslkPB6vVqunp6fNZnNmZiYQCMTjcVdHsWMQkAAAoIPf708mk95fFwEJACbGrfPFjfHVAgaFtG8AABACAhIAAAgBAQkAAISAZ0gAAB7xPsvuckELCQAAhIAWEkyyO/SLMP2STYD9Pv3c4SkjrIDNRe9Q21kZZz2X36erg1ZqknQ2ObhZywPXziwPA1GhhQQAAEJAQAIAACGgyw4AwCPoOrSHFhIAAAgBAQkAAISALjuYKDY5cvwuU/LbQJl1IW75rHPXecfMHfqGOyNkcZRtV84Zt+8ZX/JvO487MZbep7esi3PE5l6NBV+Hzh9WqOO4AL/jhNtxzY1awcihhQQAAEJACwkAwCNIarCHFhIAAAhBrICkKIqu6/wWVVVNWwAAYCKJEpB0XV9aWnr8+PHS0lKlUmEbV1dXS6XS2tqaoijjrR4AALhNlGdI5XI5kUikUqmVlZVcLpdIJCqVSjAY3N3d1TQtl8vFYrFx1xEuge6UsI8ofGGemHUSl8kNbtlZ4lbI+owfcss2zxZ+b3HYl5YVuEP/YSw7zLhzmJ1oIkACHh7KTBRRAlIqlWILqqoGg0G2MDc3R0TBYLBWq42zcgAA4D5fu+1wyGEv5PP5arWaTqeTyWQ+n5+fn2cNo9XV1aOjo76K8vl8PbeHw+Hh6wkAIjs+Pu7eOPa/dZFIpPHHHhUbsLRXwo1GY1SlCUKUFhKTzWbT6fTi4mIymSQiTdPY9gFaSOGwpz+tbYv4N14fhcN/3+s3c6r0exOE7rJ7xi1/aTqO33diLLEuuwtvwmCTboyly66zqjc69q2+vMUxunZ4fkN+S9cikYjrNYOhiZLUkMvlWC6D3+9nW+bm5k5PT4lI0zRZlsdZOQAAcJ8oLaRkMrmxsaGqar1eT6fTRJRIJIrFYj6fr9VqxhMmAACYVKIEpLm5uQcPHtTr9ZWVFZbUQEQPHjxQFIXfAq4a7WSpJo4z2Xqf4irr+tzoWAtxnW58/5tNh9v3+dM7d/GrdmPZWZRmOoVPuuNOuUO/IqIwvcsWrIs+sa6BJdGGvIPLTpSARER+v787txvZ3gAwOZCmbkuUZ0gAADDlEJAAAEAICEgAACAEBCQAABCCQEkN4JlRZ9PxD2rPrA4a4KIeJnHxH4FPpet8l/WnvY+ye1JtMy2sw+fb/B21KS3ELfMZdyfG8TeIiM744vhl/nzTrmfkgOnn697PqzNd84XVYWdi5g+MsFJ/Gl1RwkALCQAAhICABAAAQkBAAgAAISAgAQCAEBCQAABACMiym1h36Bdh+uVwCXV8SlDIYpkGShc7sVi21DUUG5dIFuAqEDKfGP5T5M61iybC4avGp8/9tPOwN/lrWuYTjpjF3T0Lde7gq80Pecdmr/icaJWIiJ5xZ51xyyem0q0+nfk4i8M8y5DsrMDJ+Sc6OQudbxcn4Q5ZdrbQQgIAACEgIAEAgBAQkAAAQAh4hgQAMC10Xa/X60Qky7IxPbc4EJAAALwywqSGP/R3eD6fr1arzWbT2CLLciqVSiQS/DEsXNlLJBLJZLK/yzuDgDRRHOfUOUyfCzk6yunvGJ/fxRVhStSyGmvNxKo+oY6jyE/0UVeynH1pfPretbPOo07Od9lV7pzDEdX4w2xK5nedBTpKPrvGJZiFQuc72EB0n7+8Cd0JeN+d03klPhnvxGbYPtNpHTWy3jUy79NVfvXOMy6d8ktux5s0zXRdX1tb64409Xr91q1biqLs7u4aW2q12oUFyrI8+loSEQISAMBkKxQKLBrF4/FUKjU3N0dElUqlWCzW6/VyuTw3N8daPIlEwj7YHBwcEBErwQ0ISAAAE0vTNBZFotHo/v6+sT2RSMRisaWlpWazWSgUWECy74jL5XJEtLy8zPfyjRay7AAAJpaiKGwhk8mYdvn9/pWVFSJqNpuqqtqXU6lUyuXyzMxMdzkjhBYSAIBXPB8z4vT0lC307GczOuharZZNIbqub21tEdHu7q6ruXkISAAAE2t2djYajVpFEV3XnRRSKBRarVY8Ho/FYiOtnRkCkqBs5t+8Q4fcHptZSENEv/wuwSgQ6jiKX7NKuDPt8ltsJ9tpTJ2wybKzYXXRUNdhH70cxs26nnz2Gp9KF+rMIhsgs44vYbA5TB1e9IT/5HyK3I2X2298V6Fz/Eywdv/qtcq469xlN67d+TiE79Pf2lxpWFbDJU6xZDJp82SIdehJkmQTaRRFOTg4kCTJ1c46Bs+QAACmUalUKpfLRLS8vGxzWKFQYMcEg0G3q4QWEgDA1Mnn8yz7TpbldDptdVipVKrVapIk2RwzQghIAAAe8R2dL7dXBz+XiMLhAeugKEoul2NDNsiyfHh4aJOnYDSPvBlnCAEJAMAj7VsjOzfyed8laJqWy+WMsRjW19fT6bRNpCmVSs1m07PmESEgAQBMg2KxyJLliCgej2cymQufCVUqFXawZ8OwIiABAEy4XC7H8hdkWc5kMk6yt1VVZW0pl8ZR7QkBSSA2Q6Peod9yayFusTOTmNtD14g+J7oVIuocUpOs075NrPKkHWYve/UOoFXSNtv1B6IboS9Nh5lyqa1Sq/3k6C0Nk+93jFp68fbB6J2Z2lYDv7LlPxFdo2fUOe5qx/CsDlP5n3UednJemk12+vvks9w3UvyF7gQumrd+auTzeSObzhhH9UKseTQzM+PeyHXdkPYNADCxjLHs+opGRFStVokoHo+7VbNe0EICAPCK50MHsbjS72utqqqyNDwvm0eEgAQAMMEeP35MRJIksQRuKysrK3yOA+uvIyL3BvbuCQEJAGBisZmQms0m67izMj8/zwckdlY0GnW7eiZ4hgQAMLHsh/G2wvLr3JsZ1gpaSEIJccudo1he43bd4Lb/sLMA/qQQ0edE7xOR3ViiDjmcV3sADk8f7DCWZccSzGxS5mwS8KwOczhS6vAznTu8EP/pui/6Py/T8PhdJ9yQuyf8F6yrBr2XbQ4zObHe5R4vplC/BBqNhmdnDQ8BCQDAKx69YHpZocsOAACEgIAEAABCECsgqaqqaZppi8M5DQEA4FIT5RmSrutra2uyLGuaJstyNps1tiiKkk6nvRxPCQAAvCdKQCqXy9FoNJvNEtHCwkI2m61Wq2yLrusbGxuTGpA6B6njkuRWO4/7G275xvliINSRS2Qa0OwPRDcCX5JtrpdNEteQ6XMOuZe/Ry9zz9j/HY5lN3zW3wCl8T8F00/EapepZKtd7LP/D9EP6b+J6PfmYQ1fnh7ovCg/RuI188HnHM1gjiw7cEqUgGSMmGR00MmyXCqVFEWp1+vep8MDAIye50MHXS6+dlugMXEVRdnb20smk8lkkjWMgsFgvV5nW/oqyufrPcBweOB5FgHgkjg+Pu7eOPa/dZFIpPF/elRswNL+OTyut4XcI0oLiYjy+Xyz2dzf32cjWBQKhUQiweLQwsJCvwEpHPb0p7VtEf8u5G6XXeTdv2zcp+nusvso8u9/33iHRt1lZ8ObLjubErq3GzeB77I74V7EftbZMXd2wpXGTzlhmn6Cn0/DZvqJzzyafoJ35/3z8NPepkgk4n0doF+iBCQ2V+7+/v64KwIAAOMhSkBiCd+rq6ts9ejoaGVlZXV19fT0tF6vezwnBwAAeE+UgNQ9c1QwGPziiy8URTGNi37Z3aFfdW6w6Kb7WedRN/gTzvtNTD1CpvlS/9BrBDPBDTZJq6HnJ2UbbXrSQhZ5YIPdt44LnTnrsgsM22V34SnsM/LzzNrckI4sO4f3wLpqfO/Z+3e8mj2Wv9C2ME/KkdRgS5SAZMXJ3O8AADABxBqpAQAAphYCEgAACAEBCQAAhCD6MyQAgMmBpAZbCEjeC3Ws3eCWf3q+GLjRkdvEZ9ZZpYRRr6wwtmWw1LUhc/M8e+X2wldru++YaYtVCQOOa8eXbTPdKnd7AgEu4a3zZechs+yMt19ZJb/PvcvKH2bzpbKbTLbjOEdHAdhAlx0AAAgBAQkAAISAgAQAAEJAQAIAACEgqQEAwCvIsrOFgOSFO/QLbq3zKxnilq/xm086jzpfdZj6xU+WajJkUplDrg6g5/AjWGXZDT/nRcdhpuL5VZvC+F3Wf6c6EvACjurWMxnPJo+OHM/H8eya9fSxpmtemtETQSDosgMAACEgIAEAgBAQkAAAQAh4hgQA4BUkNdhCCwkAAISAFpL3rEc0c8bhwHTd06SadjnfLiCrNLBeH+F7F36uQT64zSB1Nrs6r9r7sBPrw5x9WfhkvFCAFfdK99zBNveN/47Z3JyOpLuTzn0WJ935bcfkre//lUcTyMKlgBYSAAAIAQEJAACEgIAEAABCwDMkAACvIMvOFlpIAAAgBLSQvPA+/dxYvkNvduw76b18FnL0T6kLx2S7MMvOvcw6h9ObOuRqBiBfVafD3DlMpbO9Kl+0o9Ic5ux1bGbnvBI4OyPTyHjWXwOrXaafacfUtKGQZVUvTfImjBlaSAAAIAQEJAAAEAK67AAAvIKkBltoIQEAgBAQkAAAQAjosvPes8610PnyyfmiTYqawxy5sYxlZ5OuNmTSnel0q6o6PMwh8+kOM8cGuObQE8tefPqJuQQ+4+5aoONraXXffk/f7yzbskJngUDPZZN/bX/Tc/vPfFetToEJhoAEADAtdF2v1+tEJMuy3++/8HhVVVutVjAYDAaD7tcOAQkAwCujfTmvL/l8vlqtNptNY4ssy6lUKpFIdB+s63qhUCiXy61Wyzg4k8nEYjFXK4lnSAAAk0zX9aWlpYODAz4aEVG9Xr9161YulzMdr6oqO96IRuzgtbW1UqnkalURkAAAJlmhUGDddPF4/NNPP200Go1G48MPP5RlmYjK5TIfZnRd39raajabkiTt7Oywgw8PD2dmZohob29P0zT3qoqABAAwsTRNOzg4IKJoNLq/vz83N8e2JxIJI8wUCgXjeCN6HR4eJpNJtjEWix0dHRFRq9X6+OOP3astAhIAwMRSFIUtZDIZ0y6/37+yskJEzWZTVVW2sVqtEtHy8rIRuphgMBiPx4moVqu5V1skNXhvkETkfmbsvuD04Xc5TEm3ycAe/tFuPyX86cKDLW+Cw7nJhx88dIASbD5Tz9LObC/UWVro2gl3HjeCKoU6i7RO++Z2OfzCgBtOT0/ZginAMKzXjojY4yJVVdlzpp6ZDvv7+27V8iUEJAAAj3ifZTc7OxuNRq0yvHVd51crlQpbcDubzgoCEgDAxEomk8ajoG6sQ0+SJD4CGc0mVVVrtdrZ2dns7GwsFvPgVSQEJACAaVQqlcrlMhEtLy+zLSydQZIkVVW3trbYqmF9fT2dTjt5nXZgYgUkVVX9fj8fh1mKoTcvCQMATIl8Ps+y72RZTqfT/C5N09bW1lqtliRJrLVUr9dbrdbBwUGtVnvw4IF7tRIlIOm6vra2JsuypmmyLGezWSJib2xpmpZIJGxanQAAl8JV34mx/E07NPC5RBQOhwerg6IouVyOJS/Isnx4eGhq9LBdmUwmlUqxLbqu7+3tlcvler2ez+fZ32c3iBKQyuVyNBpln3NhYSGbzbLOzd3dXfai1gQFpGsda6Hey6GOuc3JT+fPHvvKsutrcFWHbE53OLiqtxl332XZDTBj94iz7EY6N7lTZ10LTpycL4ZC5yumr+XwWXY8/s7r7V/wu/y+nzspQXy/bZ//+vf7S8ifS0TJyJ/6vbqmablczsjbtumCW19fN6IREfn9/t3dXU3TarVauVye/IDEMtyJy/p4/Pjx7Owse4XYg3RDAIAJViwWC4UCS++Ox+OZTMbmUYipE49JJBK1Wq3VaimK4lIanigBid0aRVH29vaMe/Hxxx+n0+nT09PV1VX2nrBzx8fHPp+ve/vA7dwLOC42TJ1NPf5fOe/23kxEH1ksX+jdyB/6Odzke/2fYvOvNqtdw9TQkV43wfTRXrFYnhyRzZF88//HdnUo/A/JXG6fv7bHx8f86p1efwqmSi6XY/kL9mOkstbSzMxMz2bTdGXZ5fP5ZrO5v79vfOx4PM566paWlvotLRwONxqNEVfR2rbjb/wd+m3H+jWuGc71Uly70TE/zU/pP893cTMq2Xe+vRv5w/3GX3ZvH/51WisO+3AcnjK8AJ0ZN2GQLruTzuL4k04sttsYvsvOarl7tbPkyGa4ce/Y4ohep4e4AkLn+76kG/xRz7j+Z5upkgbosjN3WQ/XZbfdbkcikWFKuNTy+byRTbe7u2tzJBtJiB9T1WOiDB1UKpVM0Wh2dtbYa3p7CwAAnDDGsrswGtHL0RxarVbPEVSNjZIkjbqa3xGlhaSqqqZpq6urbPXo6CiZTC4tLeXz+Xq93rNDEwDgchlpf4CjpAY2Np0kSd1j2XVLJBJbW1tsBNXuzAU2jsPMzEzPUYhGQpSA1DN0P3jwQFGUlZWViXoPKdCZZfdTbvnG+eK1zpnO+R6MAXrSpmfEsBF8UocdazZZdgNUge8CMPXeOxt9zvIwU5adw9L4PSHLnrQBBqlzeIq/445QqL1qLJ/4jqxKgG6PHz8mIkmS+CG9uxl/ZpeXlw8ODg4ODmZnZ/ncZiM9j43H6hJRApKVcQ2pBAAwAdhoC81mk3XcWZmfn2cBKZvN1mq1er2+tbVVKpWi0SgRGVPNLi8v8+ngIyd6QAIAgIENkKFweHi4sbHBwhI/etD6+rp7byAxCEgAABNrgGRjv99/dHSkqmqlUmEDCAUCgXg8Pl1p3wAAIIi5uTn3khesICABAHhkpFl2rr9R7j0EJC/cCbXPV2507nuT3/OlsfxD+m/+KKucpcHGshvJwQbv5xwzuXDSW/7/PU+xHL9usLHsdOtdTphOsbq7J45LuPB4+xK4E/mMO7J9lbWzsPPPYEqfs2LzA7rWPv+deeb7zElpcFmI8mIsAABMOQQkAAAQAgISAAAIAc+QAAA8MvYHroJDCwkAAISAFpJb7tCL85UQt+OnHYeFrp0Yy39D/5c744Q/rK9ZYt3m8F95A0xF4ZzDtMMAnRF978Isu2HvqClxzGEyntU9cDhInU22Ws+L/nfXYfyFbErjDuNnj6WB5pgwHXbS8bvBX9PR/CA32h2lBZ5xZ3UMBgmXA1pIAAAgBAQkAAAQArrsAAA8gqQGe2ghAQCAEBCQAABACOiyc8v7dMVYvrNxPpZd4M2O9KE36XwwLn4sO7vB1vrwXYLZYBlQDjPZvOmFuHDAup6HBeiM6JXB8xIdnmczx6vD+WcHYHP6Sddh3+uVdcb/3ELWpfGHnXTs4ZPuhv+O8btsskw7dnXWp2O1O6sQhIcWEgAACAEBCQAAhIAuOwAAj+jIsrOFFhIAAAgBAQkAAISALju33PkVN0vstfNFPpWOrMevG+3krQ7z4ga7qFXv+mmFAAAaCUlEQVSilM1Fhx+OzyqzzmmWncOpYG1r4OgU6xS1QfSVWccvh3td3eYjOLsJgWvcnQ85PMd0HYtvhamwEwfL1JlZZ9oFlwFaSAAAIAS0kAAAPHJmfmcNOqCFBAAAQkBAAgAAISAgAQCAEPAMyTUhbjFwYiz/sHOMLavMOps8NJtBw6wS3mwOs2F12GAVGO1Etw6z7L5b4//f05BZdqb7wc++emZ9mNWFHFbAOg+t56fh9r+szlnv5T4uyq+GrM+yfhnUMrPOJsvu9xbb7UuAywABCQDAI5gPyR667AAAQAgISAAAIAQEJAAAEAICEgAACAFJDa4JnS9e42brtJkKc8jkN+oxMljvyVKHHFaunwqMjMMZY3ungfWVZTcAm/Q5h0PeDZDk1iudsK/C7Jw4K5pf/X3nrgFGJOCzE086d1nVByYLAhIAgEeQZWcPXXYAACAEBCQAABCCWAFJVVVN00wbNU3Tdb3n8QAAMDFEeYak6/ra2posy5qmybKczWaN7UtLS7/5zW/GWz0AAHCbKAGpXC5Ho1EWhxYWFoyAtLe3J0mSrut+/yWbR4TPprNaJtdT0b7LsrMZVm6ACgw/Sp4NpyU4zCpzmGVnVbJDAetV3fow9+pjbZBH6nwFTqx3OZw216Zwh+l8NmwGGBQDkhrsiRKQ4vE4W+B754rF4uzsrKZply4aAQBAv0R5hhQMBoPBoKIoa2tr6XSaiFRVVVU1lUqNu2oAAOAFX7vdHncdvpPP55vNZiaTCQaDRLSxsTEzM0NE1Wo1Foul02m23SGfz9dzezgcHkltAUBYx8fH3RvH/rcuEon8ZeP+qEr7Q+TdRqMxqtIEIUqXXalUajab+/v7xpZUKtVqtYioVqslEglJkvoqMBwOe/nT2u6Kf//a/sZYfpM+M5Z/Sv/JH8YP4jDa50lEFImEG41jGvVTH1f7wUf+DCmSDDdKxxccZlXyYE64ZZvJe/jVAQZx6Gekhp+Ew/91fEzWD1ac/kRtnpaFHJ/Fc+kZUqjzor9tRyIRZ0VMOF3X6/U6EcmyLOCjEFECEkv4Xl1dZatHR0dzc3NsWZKkWCw2tpoBAIzIGJMa8vl8tVptNpvGFlmWU6lUIpEwHWn8He4pk8kYf5xHTpSAtLu7a7Xr6OjIw4oAAEwU9lINaxjx6vX6rVu3FEXh//zqul6r1WxKYx1XLhElIE0ePr3bT45e7B3LyKQDXCfgdL7rgQwwsbfDtG+HpwzAs+E+LxpFtluga+Gigvs4LmDTydazEvaHDZ9wD70UCgUWjeLxeCqVYu2bSqVSLBbr9Xq5XJ6bm0smk+xgI26tr6/3LK2vZ/n9QkACAJhYmqYdHBwQUTQa5R/SJxKJWCy2tLTUbDYLhYIRkIyxcoyXQb0kSto3AACMnKIobCGTyZh2+f3+lZUVImo2m6qqso2np6dEFI1GPazjOQQkAICJxQIMEfXMRJBlmS0YT4ZY1sO4EvDQZQcA4JERZtl9z9lhs7Oz0WjUKsB0j1vNtrB3QL2HgAQAMLGSyaTxfKgb69DjX61hKXaBQCCXy1WrVdZympmZGWB0ggEgILlltAOYDpt0N9jZA0zFPbwhX18dbRLX8J+U/4epzVulzt4JHeYTnPWqQl8ldJ9+ZrnS+eGGvo0ByxUYXKlUKpfLRLS8vGzatbe3x682m81yuVytVjOZjE14Gx4CEgDA1Mnn8yz7TpZlNnwocRkQRLS+vj4/Py9JUqvVevz48cHBQavV2trakmV58l+MBQCYeCe+q8ZyiBtdrN9zaYhhORVFyeVyLHlBluXDw0PjCZPxSGlnZ4dvCcVisbm5uVu3bhHR3t6ee4MVICABAHgk0H5hLPfbkcmfS0QU+et+r65pWi6XMwZiWF9fT6fTfL5DIpHoHknI2KUoSrlcrtVq7k1Qh7RvAIDJVywWl5aWWDSKx+OPHj3KZrN9xRUj8aF7FKJRQQsJAGDC5XI5lr8gy3ImkxlsuGoPXk5CQHKLVZbdyOeYsGQzb/do880cZj25l+Tm3h0dbNLxAQq3TmJ0OFycQ0OWYDrd4bQSNtuPnF132+JK29YZidvXnBU9BfL5vJFNZzOMtQgQkAAAJpYxlp2TaLSxsaHr+vz8fM+puo1h7tx7GwkBCQDAIyMcqeEVZ4dVq1UikiSpeyy7nmq1Wr1eX15e7u6gq1QqRCTLMgISAAD07fHjx0QkSVKhULA5bGVlJRgMJpNJNjrD1tbWzs4Onw6+t7fHEiJ6Np5GBQEJAGBisYy4ZrPJOu6szM/PB4PBWCy2vLzMBmVQFCUej7MXY40xhJaXl63ywkcCAQkAYGL1O8Hr7u6uJEnlcrnVarFUCEaSpHQ67WrziBCQ3DNkNp3T04fPQxt+ms6x58J5Zsj+f9PAyt4OTzjwWHY8m9M/HK5ke+ZsOienPOOWR1eTS6fRaPR7SjabTafT1Wr19PS0Xq/Lsjw7OxuPx5H2DQAAXvP7/a4OomoFAQkAwCsT0E/gJgwdBAAAQkBAAgAAISAgAQCAEPAMyS1+LqfK4Vh2rmTW9TWWncMB61wdUs2zqWkH4N5cpf3/5B3eG5uMSIef5sjZYSZWo8/ZsT5l+2SgSsBlg4AEAOAV0f6NJRh02QEAgBAQkAAAQAgISAAAIAQEJAAAEAKSGtzicMbYQTLrRjv16gAj0Q2WbCbykHdWnyhgu+oE/3FOrA/jd3Ve5YwrYZhvQfdYdjYlfNZn4US0Hep1Pat69CxhGh74T8NnHAJaSAAAIAQEJAAAEAICEgAACAEBCQAAhICABAAAQkCWnVsGmTH2smTgmOrp3ghvNlexulc9KxOw3mVzVvd2bvXM2cfu+BrY5Oz9nlt+RlYcZtn1TKUb+Mt1w6Jk8+qJg0p04fe8aV2Y1Qe/LL8x5y5fjT2FFhIAAAgBAQkAAISAgAQAAEIQKyCpqqppmv0WAACYSKIkNei6vra2JsuypmmyLGez2e4t464jAMBwkNRgS5SAVC6Xo9EoizoLCwvZbLZ7y7jrOFZDzvHqMEVttIa/isO8uAsPC/VxIZv0OatdNqd0DGMYsh7G8Pu9KzMYm3EHHf5MQtzyiXVpNpl0n3Hr29aHOUyZ6yunEi4vUQJSPB5nC7quW20BAIAJ5mu32+OuwzlFUfb29pLJZDKZtNrikM/n67k9HA4PW0sAENvx8XH3xrH/rYtEIsf3GqMqLbwZaTRGVpogRGkhEVE+n282m/v7+8Fg0GqLc+Fw2Muf1nZX/Au1V8+XuW6Pa53vPXb02wzQZ2HbZRf5SbjxXz1+M0f8bqFVv4mrXXYODwtQ5K/Djf/X6yZYlOBil51p5pETbpWf7+FfO0o44b4v/FfH4d1lNdsMh+/1+httxWEXmU2/Kf+Btq1Ltpqdw3TYCQ3lVrsdiUSGKwO8IEpAKpVKLPbYbAEAgAkmSkBi6d2rq6ts9ejoqHvLmKoGADAiyLKzJUpA2t3dvXALAABMMFEC0uThnytYLZPNGKzDTzrukMNJrW2e07hXH4cVuHBXgMjxcKgDPEMarLRAwNGPfsi7O9jgqg4fWTrMvd/mlm9YFzhACjjSvieMWCM1AADA1EJAAgAAIaDLDgDAK0hqsIUWEgAACAEBCQAAhIAuO7ec+T48X2mfz85sSrviVweZ9XwwVpl1IiQtOazPAMl4dtd0dNpnPuvZxUdntXPVYR7akD+6wb55Awzte+KsNPRsTSe0kAAAQAhoIQEAeAVTF9hCCwkAAISAgAQAAEJAlx0AwLTQdb1erxORLMt+v7+vs4LB4AAzAfUFAemScC/ryOFYdsMb7TTU/Q/N9uEVz8YHHMpJ56o3k3wPP1Th8Dd3gJ/BaFMNJ1s+n69Wq81m09giy3IqlUokEheeu7GxUavV1tfXs9msm3VEQAIAmGi6rq+trbGGEa9er9+6dUtRFPupFYrFYq1Wc7OC5xCQAAC8Mo6meKFQYNEoHo+nUqm5uTkiqlQqxWKxXq+Xy+W5ublkMtnzXFVV9/b2PKsqkhoAACaWpmkHBwdEFI1G9/f3WTQiokQicXh4ODMzQ0SFQsHq9K2tLW/qySAgAQBMLEVR2EImkzHt8vv9KysrRNRsNlVV7T43n8/X6/VoNBqNRt2uJ4OABAAwsU5PT9mC0TbiybLMFlqtlmmXoigHBweSJHk5eTeeIXnhme8zYznUvsbv6hjLjp9C1CZ/aOwZYaPNVhvtMGqdiVbbf0UUpu0rA13CQ/x3wvTJxv7TdnWu4NF+urHfKwHNzs5Go1GrDG9d7z10hK7rGxsbRJTJZNxO9eYhIAEAeMXzmJlMJq0SFuhlh54kSbFYjN++tbXVarXi8bjNuW5Alx0AwDQqlUrlcpmIlpeXTdur1aokSTs7Ox5XCS0kAICpk8/nWfadLMvpdNrYrmkay/Pe2dlxPpTDqCAgAQB45ch3vrzaHvxcIgqHB6uCoii5XI4N2SDL8uHhIR94NjY2Wq3W+vq6kxEcRg4BCQDAK/0GIZtzlUi/BWialsvljGEX1tfX0+k0H41YnvfMzAzfZvISApLXTJOTejdL7JA8G6rMKp/QJvmMW97+ubOrCOCaxXaXEs/si8UocB4Z3697sVgsFAosvTsej3enz6mqyjrx9vf3ve+sYxCQAAAmXC6XY/kLsixnMhlTTh1jvIr0d3/3dz0LOTg4YBHLvVFWEZAAACZZPp83sum8fMt1AAhIAAATyxjL7sJoFIvFGo1Gz12rq6veTD+B95AAACZWtVolIkmSuseyExBaSAAAE+vx48dEJEmSzZDeRLSysuLlEEFWEJC89qXvS351tR3qfZwp7enMYperU2ZaZQQNP9qa1cdxftHfD1cBwYx9LlvPbuHl/1kNx/PPz2ZCajabrOPOyvz8vAgBCV12AAATq3sYb5GhhQQAMLGs8hT6cnR0NHwhTqCFBAAAQkBAAgAAIaDLDgDAK9Oe1HEBBKQxO/KdGMvbA4y76NkYZDbDylkdNkDJzg/jV0/6v+g4WA1eZ4I/WTC10GUHAABCQEACAAAhICABAIAQxHqGpKqq3+/nXxhWVTUYDI5rcg4AgFEa4RPC742uKGGIEpB0XV9bW5NlWdM0WZbZmLKrq6vBYLBer1tN4AEAABNDlIBULpej0SiLQwsLC9lstlKpBIPB3d1dNu3uNASkbR+3/KJzn834dU44HC9uMKOdMdaK3rl6MrqSBePqz8ozl7TaMF6iBKR4PM4WdP27Pzyqqs7NzRFRMBg0JoEHAIBJJUpAYs+NFEXZ29tLp9P8RiKKRqP9Fnh8fOzz+bq3h8PhIappbdTFfvTXoykn8hN3Pq+wen3cj1z6offpo7Fe/Zdi3ARvHB8f86t3ev0pAAGJEpCIKJ/PN5vN/f19Iw5pmsYWBmghhcPhkYwq6ND2qL/x5i47nuPekMhPwo3/OiYadIoHq5dhB5t+YrRdds96l7z9jEw+Cof/vvPP07jYvBjr9vQTvwyH/7e3N0GoLrvtdjsSiYy7FnAxUQJSqVRi0cjYMjc3p6oqEbE0h/FVDQBgREYYqF8ZXVHCECUgqaqqadrq6ipbPTo6SiQSxWIxn8/XarVUKjXW2gEAgOtECUi7u7vdGx88eKAoiiBz615uDruEBuuyG+CiuvUuKyeWe7q76cRxg1se7B/HDkcrHEsXmVD9cjABRAlIVqYh2xsAAAhDBwEAgCBEbyEBAEwOJDXYQgsJAACEgIAEAABCQEACAAAh4BmSoLavdK5+w6047IYey+zmPJvc7gESxwcbHmLcvuSWHU5hLppLcqdhEiAgAQB4BvHdDrrsAABACAhIAAAgBAQkAAAQAgISAAAIAUkNl8P2VW75V9bHdc90fuEz1OEHV9Wtd1mV5rA+1ttFHlDVIZGfbotcN5hgCEgAAJ5BrLeDLjsAABACAhIAAAgBAQkAAISAgAQAAEJAUsNk6X5i6kaWne7ssH5rwpz0PmX7cj4MHn78Ovc+9+W8o5cd7rodtJAAAEAICEgAACAEBCQAABACAhIAAAgBSQ0AAJ5BUoMdBKTLZ/tvnR32H0TUZ5bdaJlGn/vS0Unbo6+H60Kdqyfc8gB319W/WPhzCCJDlx0AAAgBLSQAgGmh63q9XiciWZb9fr/9waqqtlotSZLm5uY8qR0CEgDAFMjn89VqtdlsGltkWU6lUolEwnSkpmmFQqFcLvMbl5eX0+l0MBh0tZLosgMAmGS6ri8tLR0cHPDRiIjq9fqtW7dyuRy/UVXVpaUlUzQionK5vLS0pKqqq1VFQAIA8MzZ6P5zqlAosG66eDz+6aefNhqNRqPx4YcfyrJMROVyuVQqGQdvbW21Wi0iymQyv/71rxuNxq9//etMJiNJUqvV2tjYGPUN6YAuu4m1/RZRmLbf6rXrkFvpf/S57Q+HqdcFtl0s2yOhcVcAwKBp2sHBARFFo9H9/X1jeyKRiMViS0tLzWazUCgkk0kiUlWVha5MJpNKpdiRfr+fLe/t7TWbzUql0t3LNypoIQEATCxFUdhCJpMx7fL7/SsrK0TUbDZZX1ylUmG7jGhkWF5eZguu9tohIAEATKzT01O20DNTjvXaERHrpmMPmaLRaPeRF6bkjQS67AAAJtbs7Gw0GrUKJ7reMZdMKpVKJpOSJNkcGQgERl5JAwISAIBnvB4rI5lMsudDPbEOPUmSYrEYWbSimL29PbYQj8dHXcdzCEgAANOoVCqx9G7j+VD3AazHz3iBaWdnx9VXkRCQptH2Grf8C275597XBcxMHSIYfQ7ckM/nWfadLMvpdLrnMZVKpVarGasffvihe/l1DAISAIBnuH8M0qHlURefS0ThwWqgKEoul2MtHlmWDw8PrZ4wGSkPLCzdunVLVdVsNjvYdZ2Y2IB0fHw87iqMH24C4SYQEW6CQHeg3yBkc+4/93u+pmm5XM5o9Kyvr6fTaZv0OSP26Lq+t7dXLpdZo8q9mDSxAQkAQDxj64ItFouFQoGld8fj8Uwm4/xpkN/v393d1TStVquVy2X3ApJY7yHpuq5pGr9FVVVTYiIAAPQll8vt7e21Wi3WR7e/vz9AbgLrwWu1Wu79TRYrIJXL5Y8//pgtswEBS6US+/94KwYAcEnl83kjm+7Bgwcsw7unSCQSiUSs/t4aSeFseCE3CNRlt7S0VK/X19fX2Wq1Wo1Go9lsVtf1jY0Nm1R6EWy32+OuQg93fL6+KrZ9y726jE2/N2Ei/dznuzX1N2E6GWPZLS8v7+7uOjnFGNzBxBg0yCakDUmgFtKDBw+MaEREsizXajVFUcrlspHsAQAAzlWrVSKSJKl7LLtu7KVXdopVUTMzMyOtYAeBWkgmwWBQkqRKpVKv1wdrHvl8vu6N4fCAuZKXUTgcjkQi467FmOEm0PTdBGFy6szC4V96fMXHjx8TkSRJhULB5rCVlZVgMBiLxdg7sLlcztScMjLFrV5aGglfW6SGfD6fp5c5hfl8fnZ2loWihYWFL774YsyVAwC4bH784x+zzDp7h4eHrCNuY2PDaAnFYjE2DZKiKCwaxeNxfg6LkRO3hQQAAENyEo14Ozs7MzMzbHpZ07yx6+vrrr4VSyIHpJWVldXV1dPT03q97upwfgAAk6rRaPR1vN/vz2azKysr1Wr17OysXq/Lsjw7OxuLxVwdxY4Rq8uum6IowWDQgxsBAADjJXpAAgCAKSFQ2jcAAEwzBCQAABACAhIAAAgBAQkAAISAgAQAAEJAQAIAACEgIAEAgBAQkAAAQAgISAAAIAQEJAAAEAICEgAACAEBCQAAhICABAAAQkBAAgAAISAgAQCAEBCQAABACAhIAAAghIkKSLqua5rGb1FVVdd1q9XJNlUf1so0fwGYKfzI3bpvAm6LmP5se3t73HUYmX/7t3978uTJ/Pw8W11dXdU0rVgsBoPBYDBoWh1vVV2l6/o//MM/6Lr+L//yL5Ik/eAHP/jxj3+sKMpnn332u9/9zrg/k637JkzPF4DRdf2dd9754x//uLe39xd/8Rdzc3NENG3fhJ43Ydq+CZdJe1K89dZb4XD4gw8+YKuff/75P/7jP7bb7dPT05/97Gem1THW0wMffPDBxx9/3G63z87O3nrrrWn4yN1MN2GqvgDMxx9/zH4dzs7O2Eeens9u6L4JU/hNuET+fNwBcWQePHiQz+eNVVVV2b+GgsFgrVaTZZlfHVstvSJJEhH5/f56va5pmt/vz+VykiSl02m/3z/u2nmEvwmm78O4q+YFWZZLpZKiKPV6XZZlIprCb0L3TZjCb8IlMlHPkEyMxng0Gu1enWDz8/PFYrFSqWxsbEiSpOv6zMxMIpEIBAIbGxvjrp1HTDeBpukLwASDQUmSKpVKpVKZnZ0loin8JnTfBJq+b8IlcrlbSIqiPH78eHZ2NplMdu81EhxYC4lf9a6KXuFvRSwW29nZqVQqyWRSUZREIpFIJIgoFotVKpVx19RFNjeBOr8PY62mu4ybcHp6mkgk2K/GwsJCMpmcwm9C902gqfkmXEaXOyDFYrFYLNZz19zcnKqqRKRpGuuv41c9raUn+FtRqVRarVY2m9U0LRaLFYtFWZatbtQksbkJE/8FMBg3ge/BZqbwm9B9E6bnm3AZXe6AZCORSBSLxXw+X6vVUqmUaXXctXNXLBZbWlo6PT2t1Wo7OztEtLGxEY/H6/V6z6bkRDLdhLm5uen5AjArKyurq6unp6f1ej0ejxNRNBqdtm9C902Yqj8Fl46v3W6Puw4uUhSFz+w0rU4wXdfZU1zjwfX0fHYDbgL1+si4CT23gAgmPCABAMBlMclZdgAAcIkgIAEAgBAQkAAAQAgISAAAIAQEJAAAEAICEgAACAEBCQAAhICABAAAQkBAAgAAISAgAQCAEBCQAABACAhIAAAgBAQkAAAQAgISAAAI4f8DoHAaZXVQjysAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<IPython.core.display.Image object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "% рисуем\n",
    "figure\n",
    "axes('fontsize',18)\n",
    "% contourf(x,y,aam,30) % 1 контур\n",
    "% imagesc(x,y,aam) % 2 рисунок\n",
    "pcolor(x,y,aam) % 2 рисунок\n",
    "shading('flat')\n",
    "% imagescn(x,y,aam) % 3 без нанов\n",
    "axis xy\n",
    "c=colorbar('fontsize',28)\n",
    "xlabel('Longitude')\n",
    "%ylabel\n",
    "title('Average SST in 2013')\n",
    "grid on\n",
    "colormap(jet)\n",
    "c1=prctile(aam(:),3); % автоподбор по шкале\n",
    "c2=prctile(aam(:),97); % автоподбор по шкале\n",
    "caxis([c1 c2])% ограничение шкалы\n",
    "savname='SST_Reinolds_7D_year2013_pcolor.jpeg'; % имя рисунка\n",
    "% print(gcf,'-r150','-djpeg',[picdir savname]) % сохраняем рисунок"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Octave",
   "language": "octave",
   "name": "octave"
  },
  "language_info": {
   "file_extension": ".m",
   "help_links": [
    {
     "text": "GNU Octave",
     "url": "https://www.gnu.org/software/octave/support.html"
    },
    {
     "text": "Octave Kernel",
     "url": "https://github.com/Calysto/octave_kernel"
    },
    {
     "text": "MetaKernel Magics",
     "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md"
    }
   ],
   "mimetype": "text/x-octave",
   "name": "octave",
   "version": "4.2.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}