examples/mkGridsExample12.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "9b16579e-353e-43bc-9bc5-256c613511ff",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Example 12:\n",
    "#\n",
    "# This example creates two simple grids in cartesian and lat/lon space and\n",
    "# demonstrates the use of the extendGrid() routine on two simple 5x5 grids.\n",
    "\n",
    "# NOTE: The i grid lines are shown in yellow.  This helps see that for the\n",
    "# stereographic grid, that the grid is actually rotated 90 degrees from\n",
    "# the nominal view of the grid where (J,I) (0,0) is at the bottom left\n",
    "# and J grows upwards and I grows rightwards."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "c3caefa2-f550-4f8d-a909-a49cdd62f601",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Initialize library\n",
    "# conda: gridTools\n",
    "import os\n",
    "from gridtools.gridutils import GridUtils\n",
    "%matplotlib inline\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "c3c6e2d5-d15f-4060-86c1-fd4c2a339960",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Create a grid in spherical space using the LCC projection\n",
    "# NOTE: Specification of the grid is in degrees, but\n",
    "# the computation is spherical coordinates using meters.\n",
    "grd = GridUtils()\n",
    "grd.setGridParameters({\n",
    "    'projection': {\n",
    "        'name': 'LambertConformalConic',\n",
    "        'lon_0': 230.0,\n",
    "        'lat_0': 40.0,\n",
    "        'ellps': 'WGS84'\n",
    "    },\n",
    "    'centerX': 230.0,\n",
    "    'centerY': 40.0,\n",
    "    'centerUnits': 'degrees',\n",
    "    'dx': 20.0,\n",
    "    'dxUnits': 'degrees',\n",
    "    'dy': 20.0,\n",
    "    'dyUnits': 'degrees',\n",
    "    'tilt': 0.0,\n",
    "    'gridResolution': 4.0,\n",
    "    'gridResolutionUnits': 'degrees',\n",
    "    'gridMode': 2,\n",
    "    'gridType': 'MOM6',\n",
    "    'ensureEvenI': True,\n",
    "    'ensureEvenJ': True,\n",
    "    'tileName': 'tile1',\n",
    "})\n",
    "grd.makeGrid()\n",
    "\n",
    "# Define plot parameters so we can see what the grid looks like\n",
    "grd.setPlotParameters(\n",
    "    {\n",
    "        'figsize': (8,8),\n",
    "        'projection': {\n",
    "            'name': 'NearsidePerspective',\n",
    "            'lat_0': 40.0,\n",
    "            'lon_0': 230.0\n",
    "        },\n",
    "        'extent': [-160.0 ,-100.0, 20.0, 60.0],\n",
    "        'iLinewidth': 1.0,\n",
    "        'jLinewidth': 1.0,\n",
    "        'showGridCells': True,\n",
    "        'title': \"Lambert Conformal Conic: 5x5\",\n",
    "        'iColor': 'y',\n",
    "        'jColor': 'k'\n",
    "    }\n",
    ")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "d362a8d3-1d3a-4d83-b4d1-62a22200610b",
   "metadata": {},
   "outputs": [],
   "source": [
    "(figure, axes) = grd.plotGrid()\n",
    "figure"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "56ec0b9d-896d-465a-aae0-4a65fc835bc6",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Extend this grid one regular cell (two cells in the supergrid) using meters\n",
    "extGrid = grd.extendGrid(2,2,2,2)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "2f8c73e2-5a73-488e-a980-78079d18a421",
   "metadata": {},
   "outputs": [],
   "source": [
    "eGridLCC = GridUtils()\n",
    "eGridLCC.readGrid(local=extGrid)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "2dd52472-7420-42cf-a36e-fc9560c13e10",
   "metadata": {},
   "outputs": [],
   "source": [
    "eGridLCC.setPlotParameters(\n",
    "    {\n",
    "        'figsize': (8,8),\n",
    "        'projection': {\n",
    "            'name': 'NearsidePerspective',\n",
    "            'lat_0': 40.0,\n",
    "            'lon_0': 230.0\n",
    "        },\n",
    "        'extent': [-160.0 ,-100.0, 20.0, 60.0],\n",
    "        'iLinewidth': 1.0,\n",
    "        'jLinewidth': 1.0,\n",
    "        'showGridCells': True,\n",
    "        'title': \"Lambert Conformal Conic: 5x5 extended by 1 regular grid point to 7x7\",\n",
    "        'iColor': 'y',\n",
    "        'jColor': 'k'\n",
    "    }\n",
    ")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "7cd5b8aa-ff3c-403e-b04e-163f558e4a3e",
   "metadata": {},
   "outputs": [],
   "source": [
    "(figure, axes) = eGridLCC.plotGrid()\n",
    "figure"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "0331708d-f6ee-4667-8cb8-54405c30cafd",
   "metadata": {},
   "outputs": [],
   "source": [
    "eGridLCC.grid"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "5e7b2710-3ee6-42a8-b03b-cfa00ff79a2d",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Create a mercator grid in lat/lon space\n",
    "grd = GridUtils()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "e9b9b370-c9e7-49f6-857c-187d65e65950",
   "metadata": {},
   "outputs": [],
   "source": [
    "grd.setGridParameters({\n",
    "    'projection': {\n",
    "        'name': 'Mercator',\n",
    "        'lon_0': 230.0,\n",
    "        'lat_0': 40.0,\n",
    "        'ellps': 'WGS84'\n",
    "    },\n",
    "    'centerX': 230.0,\n",
    "    'centerY': 40.0,\n",
    "    'centerUnits': 'degrees',\n",
    "    'dx': 20.0,\n",
    "    'dxUnits': 'degrees',\n",
    "    'dy': 20.0,\n",
    "    'dyUnits': 'degrees',\n",
    "    'tilt': 0.0,\n",
    "    'gridResolution': 4.0,\n",
    "    'gridResolutionUnits': 'degrees',\n",
    "    'gridMode': 2,\n",
    "    'gridType': 'MOM6',\n",
    "    'ensureEvenI': True,\n",
    "    'ensureEvenJ': True,\n",
    "    'tileName': 'tile1',\n",
    "})\n",
    "grd.makeGrid()\n",
    "\n",
    "# Define plot parameters so we can see what the grid looks like\n",
    "grd.setPlotParameters(\n",
    "    {\n",
    "        'figsize': (8,8),\n",
    "        'projection': {\n",
    "            'name': 'NearsidePerspective',\n",
    "            'lat_0': 40.0,\n",
    "            'lon_0': 230.0\n",
    "        },\n",
    "        'extent': [-160.0 ,-100.0, 20.0, 60.0],\n",
    "        'iLinewidth': 1.0,\n",
    "        'jLinewidth': 1.0,\n",
    "        'showGridCells': True,\n",
    "        'title': \"Mercator: 5x5\",\n",
    "        'iColor': 'y',\n",
    "        'jColor': 'k'\n",
    "    }\n",
    ")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "3cac5524-2a33-44f9-b600-fd7b53e23188",
   "metadata": {},
   "outputs": [],
   "source": [
    "(figure, axes) = grd.plotGrid()\n",
    "figure"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "77db29f7-71e4-44d2-bb4e-19e8ac19da96",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Extend this grid one regular cell (two cells in the supergrid) using degrees\n",
    "extGrid = grd.extendGrid(2,2,2,2)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "1198cab1-7c98-432b-9168-e23899ed7355",
   "metadata": {},
   "outputs": [],
   "source": [
    "eGridMercator = GridUtils()\n",
    "eGridMercator.readGrid(local=extGrid)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "a6a9b540-093b-424f-b46c-837566927fa2",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Define plot parameters so we can see what the grid looks like\n",
    "eGridMercator.setPlotParameters(\n",
    "    {\n",
    "        'figsize': (8,8),\n",
    "        'projection': {\n",
    "            'name': 'NearsidePerspective',\n",
    "            'lat_0': 40.0,\n",
    "            'lon_0': 230.0\n",
    "        },\n",
    "        'extent': [-160.0 ,-100.0, 20.0, 60.0],\n",
    "        'iLinewidth': 1.0,\n",
    "        'jLinewidth': 1.0,\n",
    "        'showGridCells': True,\n",
    "        'title': \"Mercator: 5x5 extended by 1 regular grid point to 7x7\",\n",
    "        'iColor': 'y',\n",
    "        'jColor': 'k'\n",
    "    }\n",
    ")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "c6759833-fcbc-487c-8034-12a57d803261",
   "metadata": {},
   "outputs": [],
   "source": [
    "(figure, axes) = eGridMercator.plotGrid()\n",
    "figure"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "4427c871-7ed1-4bd1-bb31-eef72cf145d5",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Create a 5x5 Stereographic grid for testing\n",
    "grd = GridUtils()\n",
    "grd.setGridParameters({\n",
    "    'projection': {\n",
    "        'name': \"Stereographic\",\n",
    "        'ellps': 'WGS84',\n",
    "        'lon_0': 0.0,\n",
    "        'lat_0': 90.0,\n",
    "        'lat_ts': 75.0,\n",
    "    },\n",
    "    'centerX': 0.0,\n",
    "    'centerY': 90.0,\n",
    "    'cneterUnits': 'degrees',\n",
    "    'dx': 2500000.0,\n",
    "    'dy': 2500000.0,\n",
    "    'dxUnits': 'meters',\n",
    "    'dyUnits': 'meters',\n",
    "    'gridResolution': 500000.0,\n",
    "    'gridResolutionUnits': 'meters',\n",
    "    'tilt': 0.0,\n",
    "    'gridMode': 2,\n",
    "    'gridType': 'MOM6',\n",
    "    'ensureEvenI': True,\n",
    "    'ensureEvenJ': True\n",
    "})\n",
    "\n",
    "grd.makeGrid()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "becd19ed-ffc7-4c3b-8327-ba54eb291360",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Define plot parameters so we can see what the grid looks like\n",
    "grd.setPlotParameters(\n",
    "    {\n",
    "        'figsize': (8,8),\n",
    "        'projection': {\n",
    "            'name': 'Stereographic',\n",
    "            'lon_0': 0.0,\n",
    "            'lat_0': 90.0\n",
    "        },\n",
    "        'extent': [-180, 180, 60, 90],\n",
    "        'iLinewidth': 1.0,\n",
    "        'jLinewidth': 1.0,\n",
    "        'showGrid': False,\n",
    "        'showGridCells': True,\n",
    "        'title': 'Stereographic: 5x5',\n",
    "        'iColor': 'y',\n",
    "        'jColor': 'k'\n",
    "    }\n",
    ")\n",
    "\n",
    "(figure, axes) = grd.plotGrid()\n",
    "figure"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "36218729-8e74-4b2d-925f-fde34cc5dcfb",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Extend this grid one regular cell (two cells in the supergrid) using meters\n",
    "extGrid = grd.extendGrid(2,2,2,2)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "0a4dd38a-9251-4cd5-abde-bdd09531e1c1",
   "metadata": {},
   "outputs": [],
   "source": [
    "eGridStereographic = GridUtils()\n",
    "eGridStereographic.readGrid(local=extGrid)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "3bbf6fbf-b1b6-45e6-94a5-a5509aeeb658",
   "metadata": {},
   "outputs": [],
   "source": [
    "eGridStereographic.setPlotParameters(\n",
    "    {\n",
    "        'figsize': (8,8),\n",
    "        'projection': {\n",
    "            'name': 'Stereographic',\n",
    "            'lon_0': 0.0,\n",
    "            'lat_0': 90.0\n",
    "        },\n",
    "        'extent': [-180, 180, 60, 90],\n",
    "        'iLinewidth': 1.0,\n",
    "        'jLinewidth': 1.0,\n",
    "        'showGrid': False,\n",
    "        'showGridCells': True,\n",
    "        'title': 'Stereographic: 5x5 extended by 1 regular grid point to 7x7',\n",
    "        'iColor': 'y',\n",
    "        'jColor': 'k'\n",
    "    }\n",
    ")\n",
    "\n",
    "(figure, axes) = eGridStereographic.plotGrid()\n",
    "figure"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "9744662d-a47b-430c-ac8d-fabb623a5301",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Grid clipping tests"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "41d54e1f-665d-48cc-bedc-58ba6070b597",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Repeat last example, but only keep sides of the grid (jstart)\n",
    "# Row along the bottom\n",
    "extGrid = grd.extendGrid(2,0,0,0)\n",
    "eGridClip = GridUtils()\n",
    "eGridClip.readGrid(local=extGrid)\n",
    "eGridClip.setPlotParameters(\n",
    "    {\n",
    "        'figsize': (8,8),\n",
    "        'projection': {\n",
    "            'name': 'Stereographic',\n",
    "            'lon_0': 0.0,\n",
    "            'lat_0': 90.0\n",
    "        },\n",
    "        'extent': [-180, 180, 60, 90],\n",
    "        'iLinewidth': 1.0,\n",
    "        'jLinewidth': 1.0,\n",
    "        'showGrid': False,\n",
    "        'showGridCells': True,\n",
    "        'title': 'Stereographic: 5x5 extended by 1 regular grid point (jstart)',\n",
    "        'iColor': 'y',\n",
    "        'jColor': 'k'\n",
    "    }\n",
    ")\n",
    "\n",
    "(figure, axes) = eGridClip.plotGrid()\n",
    "figure"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "0d9de7be-1ccf-4da0-a264-713281aaba54",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Repeat last example, but only keep sides of the grid (jend)\n",
    "# Row across the top\n",
    "extGrid = grd.extendGrid(0,2,0,0)\n",
    "eGridClip = GridUtils()\n",
    "eGridClip.readGrid(local=extGrid)\n",
    "eGridClip.setPlotParameters(\n",
    "    {\n",
    "        'figsize': (8,8),\n",
    "        'projection': {\n",
    "            'name': 'Stereographic',\n",
    "            'lon_0': 0.0,\n",
    "            'lat_0': 90.0\n",
    "        },\n",
    "        'extent': [-180, 180, 60, 90],\n",
    "        'iLinewidth': 1.0,\n",
    "        'jLinewidth': 1.0,\n",
    "        'showGrid': False,\n",
    "        'showGridCells': True,\n",
    "        'title': 'Stereographic: 5x5 extended by 1 regular grid point (jend)',\n",
    "        'iColor': 'y',\n",
    "        'jColor': 'k'\n",
    "    }\n",
    ")\n",
    "\n",
    "(figure, axes) = eGridClip.plotGrid()\n",
    "figure"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "44c5e41c-0b7b-43b7-807d-16d89c7abc07",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Repeat last example, but only keep sides of the grid (istart)\n",
    "# Row on the left\n",
    "extGrid = grd.extendGrid(0,0,2,0)\n",
    "eGridClip = GridUtils()\n",
    "eGridClip.readGrid(local=extGrid)\n",
    "eGridClip.setPlotParameters(\n",
    "    {\n",
    "        'figsize': (8,8),\n",
    "        'projection': {\n",
    "            'name': 'Stereographic',\n",
    "            'lon_0': 0.0,\n",
    "            'lat_0': 90.0\n",
    "        },\n",
    "        'extent': [-180, 180, 60, 90],\n",
    "        'iLinewidth': 1.0,\n",
    "        'jLinewidth': 1.0,\n",
    "        'showGrid': False,\n",
    "        'showGridCells': True,\n",
    "        'title': 'Stereographic: 5x5 extended by 1 regular grid point (istart)',\n",
    "        'iColor': 'y',\n",
    "        'jColor': 'k'\n",
    "    }\n",
    ")\n",
    "\n",
    "(figure, axes) = eGridClip.plotGrid()\n",
    "figure"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "7ac69401-617c-4a9d-b4b0-f5023d7687d3",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Repeat last example, but only keep sides of the grid (iend)\n",
    "# Row along the right\n",
    "extGrid = grd.extendGrid(0,0,0,2)\n",
    "eGridClip = GridUtils()\n",
    "eGridClip.readGrid(local=extGrid)\n",
    "eGridClip.setPlotParameters(\n",
    "    {\n",
    "        'figsize': (8,8),\n",
    "        'projection': {\n",
    "            'name': 'Stereographic',\n",
    "            'lon_0': 0.0,\n",
    "            'lat_0': 90.0\n",
    "        },\n",
    "        'extent': [-180, 180, 60, 90],\n",
    "        'iLinewidth': 1.0,\n",
    "        'jLinewidth': 1.0,\n",
    "        'showGrid': False,\n",
    "        'showGridCells': True,\n",
    "        'title': 'Stereographic: 5x5 extended by 1 regular grid point (iend)',\n",
    "        'iColor': 'y',\n",
    "        'jColor': 'k'\n",
    "    }\n",
    ")\n",
    "\n",
    "(figure, axes) = eGridClip.plotGrid()\n",
    "figure"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "81d368de-613c-44ba-bfb0-42d70acfab2f",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Repeat last example, but only keep two sides of the grid (i/j start)\n",
    "# Row on the left and bottom\n",
    "extGrid = grd.extendGrid(2,0,2,0)\n",
    "eGridClip = GridUtils()\n",
    "eGridClip.readGrid(local=extGrid)\n",
    "eGridClip.setPlotParameters(\n",
    "    {\n",
    "        'figsize': (8,8),\n",
    "        'projection': {\n",
    "            'name': 'Stereographic',\n",
    "            'lon_0': 0.0,\n",
    "            'lat_0': 90.0\n",
    "        },\n",
    "        'extent': [-180, 180, 60, 90],\n",
    "        'iLinewidth': 1.0,\n",
    "        'jLinewidth': 1.0,\n",
    "        'showGrid': False,\n",
    "        'showGridCells': True,\n",
    "        'title': 'Stereographic: 5x5 extended by 1 regular grid point (i/j start)',\n",
    "        'iColor': 'y',\n",
    "        'jColor': 'k'\n",
    "    }\n",
    ")\n",
    "\n",
    "(figure, axes) = eGridClip.plotGrid()\n",
    "figure"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "3b4f6931-5a74-48ad-a241-fe070f7ab33c",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Repeat last example, but only keep two sides of the grid (i/j end)\n",
    "# Row on the top and right\n",
    "extGrid = grd.extendGrid(0,2,0,2)\n",
    "eGridClip = GridUtils()\n",
    "eGridClip.readGrid(local=extGrid)\n",
    "eGridClip.setPlotParameters(\n",
    "    {\n",
    "        'figsize': (8,8),\n",
    "        'projection': {\n",
    "            'name': 'Stereographic',\n",
    "            'lon_0': 0.0,\n",
    "            'lat_0': 90.0\n",
    "        },\n",
    "        'extent': [-180, 180, 60, 90],\n",
    "        'iLinewidth': 1.0,\n",
    "        'jLinewidth': 1.0,\n",
    "        'showGrid': False,\n",
    "        'showGridCells': True,\n",
    "        'title': 'Stereographic: 5x5 extended by 1 regular grid point (i/j end)',\n",
    "        'iColor': 'y',\n",
    "        'jColor': 'k'\n",
    "    }\n",
    ")\n",
    "\n",
    "(figure, axes) = eGridClip.plotGrid()\n",
    "figure"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "065adba5-881b-4900-bb18-259a5db7506c",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Repeat last example, but only keep two sides of the grid (jstart and iend)\n",
    "# Row along bottom and right\n",
    "extGrid = grd.extendGrid(2,0,0,2)\n",
    "eGridClip = GridUtils()\n",
    "eGridClip.readGrid(local=extGrid)\n",
    "eGridClip.setPlotParameters(\n",
    "    {\n",
    "        'figsize': (8,8),\n",
    "        'projection': {\n",
    "            'name': 'Stereographic',\n",
    "            'lon_0': 0.0,\n",
    "            'lat_0': 90.0\n",
    "        },\n",
    "        'extent': [-180, 180, 60, 90],\n",
    "        'iLinewidth': 1.0,\n",
    "        'jLinewidth': 1.0,\n",
    "        'showGrid': False,\n",
    "        'showGridCells': True,\n",
    "        'title': 'Stereographic: 5x5 extended by 1 regular grid point (istart and jend)',\n",
    "        'iColor': 'y',\n",
    "        'jColor': 'k'\n",
    "    }\n",
    ")\n",
    "\n",
    "(figure, axes) = eGridClip.plotGrid()\n",
    "figure"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "d612e99a-d29e-439e-98ba-1c8608f30b13",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Create a 1x3 mercator grid for an asymmetric test\n",
    "grd = GridUtils()\n",
    "grd.setGridParameters({\n",
    "    'projection': {\n",
    "        'name': 'Mercator',\n",
    "        'lon_0': 230.0,\n",
    "        'lat_0': 40.0,\n",
    "        'ellps': 'WGS84'\n",
    "    },\n",
    "    'centerX': 230.0,\n",
    "    'centerY': 40.0,\n",
    "    'centerUnits': 'degrees',\n",
    "    'dx': 12.0,\n",
    "    'dxUnits': 'degrees',\n",
    "    'dy': 4.0,\n",
    "    'dyUnits': 'degrees',\n",
    "    'tilt': 0.0,\n",
    "    'gridResolution': 4.0,\n",
    "    'gridResolutionUnits': 'degrees',\n",
    "    'gridMode': 2,\n",
    "    'gridType': 'MOM6',\n",
    "    'ensureEvenI': True,\n",
    "    'ensureEvenJ': True,\n",
    "    'tileName': 'tile1',\n",
    "})\n",
    "grd.makeGrid()\n",
    "\n",
    "# Define plot parameters so we can see what the grid looks like\n",
    "grd.setPlotParameters(\n",
    "    {\n",
    "        'figsize': (8,8),\n",
    "        'projection': {\n",
    "            'name': 'Mercator',\n",
    "            'lat_0': 40.0,\n",
    "            'lon_0': 230.0\n",
    "        },\n",
    "        'extent': [-160.0 ,-100.0, 20.0, 60.0],\n",
    "        'iLinewidth': 1.0,\n",
    "        'jLinewidth': 1.0,\n",
    "        'showGridCells': True,\n",
    "        'title': \"Mercator: 1x3\",\n",
    "        'iColor': 'y',\n",
    "        'jColor': 'k'\n",
    "    }\n",
    ")\n",
    "\n",
    "(figure, axes) = grd.plotGrid()\n",
    "figure"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "3d01b956-4454-4ed7-959e-a7e2ef5ff73c",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Create a 3x1 mercator grid for an asymmetric test\n",
    "grd = GridUtils()\n",
    "grd.setGridParameters({\n",
    "    'projection': {\n",
    "        'name': 'Mercator',\n",
    "        'lon_0': 230.0,\n",
    "        'lat_0': 40.0,\n",
    "        'ellps': 'WGS84'\n",
    "    },\n",
    "    'centerX': 230.0,\n",
    "    'centerY': 40.0,\n",
    "    'centerUnits': 'degrees',\n",
    "    'dx': 4.0,\n",
    "    'dxUnits': 'degrees',\n",
    "    'dy': 12.0,\n",
    "    'dyUnits': 'degrees',\n",
    "    'tilt': 0.0,\n",
    "    'gridResolution': 4.0,\n",
    "    'gridResolutionUnits': 'degrees',\n",
    "    'gridMode': 2,\n",
    "    'gridType': 'MOM6',\n",
    "    'ensureEvenI': True,\n",
    "    'ensureEvenJ': True,\n",
    "    'tileName': 'tile1',\n",
    "})\n",
    "grd.makeGrid()\n",
    "\n",
    "# Define plot parameters so we can see what the grid looks like\n",
    "grd.setPlotParameters(\n",
    "    {\n",
    "        'figsize': (8,8),\n",
    "        'projection': {\n",
    "            'name': 'Mercator',\n",
    "            'lat_0': 40.0,\n",
    "            'lon_0': 230.0\n",
    "        },\n",
    "        'extent': [-160.0 ,-100.0, 20.0, 60.0],\n",
    "        'iLinewidth': 1.0,\n",
    "        'jLinewidth': 1.0,\n",
    "        'showGridCells': True,\n",
    "        'title': \"Mercator: 3x1\",\n",
    "        'iColor': 'y',\n",
    "        'jColor': 'k'\n",
    "    }\n",
    ")\n",
    "\n",
    "(figure, axes) = grd.plotGrid()\n",
    "figure"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "ba51367e-03d2-4148-b89c-654c59c93c83",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Repeat last example, but only keep sides of the grid (jstart)\n",
    "# Row along the bottom\n",
    "extGrid = grd.extendGrid(2,0,0,0)\n",
    "eGridClip = GridUtils()\n",
    "eGridClip.readGrid(local=extGrid)\n",
    "eGridClip.setPlotParameters(\n",
    "    {\n",
    "        'figsize': (8,8),\n",
    "        'projection': {\n",
    "            'name': 'Mercator',\n",
    "            'lat_0': 40.0,\n",
    "            'lon_0': 230.0\n",
    "        },\n",
    "        'extent': [-160.0 ,-100.0, 20.0, 60.0],\n",
    "        'iLinewidth': 1.0,\n",
    "        'jLinewidth': 1.0,\n",
    "        'showGrid': False,\n",
    "        'showGridCells': True,\n",
    "        'title': 'Mercator: 3x1 extended by 1 regular grid point (jstart)',\n",
    "        'iColor': 'y',\n",
    "        'jColor': 'k'\n",
    "    }\n",
    ")\n",
    "\n",
    "(figure, axes) = eGridClip.plotGrid()\n",
    "figure"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "1f23a3e6-c9c9-4776-9a5b-23ffebb242a4",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.10"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}