now this is readable by paraview

Author: kbarnhart

lithology/introduction_to_lithology.ipynb

@@ -80,7 +80,7 @@
     "ids = np.tile([0,1], 20) \n",
     "\n",
     "# Anticline\n",
-    "lith = LithoLayers(mg, z0s, ids, x0=6000, y0=6000, function = lambda x, y : ((0.002*x)**2+(0.001*y)**2), attrs=attrs)\n",
+    "lith = LithoLayers(mg, z0s, ids, x0=6000, y0=10000, function = lambda x, y : ((0.002*x)**2+(0.001*y)**2), attrs=attrs)\n",
     "\n",
     "# Shallow dips\n",
     "#lith = LithoLayers(mg, z0s, ids, function = lambda x, y : ((0.001*x)+(0.003*y)), attrs=attrs)\n",
@@ -138,11 +138,11 @@
     "out_fields = ['topographic__elevation',\n",
     "              'rock_type__id']\n",
     "\n",
-    "ds = xr.Dataset(data_vars={'topographic__elevation' : (('t', 'y', 'x'), \n",
+    "ds = xr.Dataset(data_vars={'topographic__elevation' : (('time', 'y', 'x'), \n",
     "                                                       np.empty((nts, mg.shape[0], mg.shape[1])),\n",
     "                                                      {'units' : 'meters',\n",
     "                                                       'long_name': 'Topographic Elevation'}),\n",
-    "                           'rock_type__id': (('t', 'y', 'x'), \n",
+    "                           'rock_type__id': (('time', 'y', 'x'), \n",
     "                                             np.empty((nts, mg.shape[0], mg.shape[1])),\n",
     "                                            {'units' : '-',\n",
     "                                             'long_name' : 'Rock Type ID Code'})\n",
@@ -153,9 +153,10 @@
     "                        'y': (('y'), \n",
     "                              mg.y_of_node.reshape(mg.shape)[:, 1],\n",
     "                              {'units' : 'meters'}),\n",
-    "                        'time': (('t'),  \n",
+    "                        'time': (('time'),  \n",
     "                                 dt*np.arange(nts)/1e6,\n",
-    "                                 {'units': 'millions of years since model start'})\n",
+    "                                 {'units': 'millions of years since model start',\n",
+    "                                  'standard_name' : 'time'})\n",
     "                       }\n",
     "                )\n"
    ]
@@ -174,7 +175,9 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "metadata": {},
+   "metadata": {
+    "collapsed": true
+   },
    "outputs": [],
    "source": [
     "for i in range(nts):\n",
@@ -272,7 +275,9 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "metadata": {},
+   "metadata": {
+    "collapsed": true
+   },
    "outputs": [],
    "source": [
     "mg2 = RasterModelGrid((100, 60), 200)\n",
@@ -310,11 +315,11 @@
     "dt = 1000\n",
     "\n",
     "\n",
-    "ds2 = xr.Dataset(data_vars={'topographic__elevation' : (('t', 'y', 'x'), \n",
+    "ds2 = xr.Dataset(data_vars={'topographic__elevation' : (('time', 'y', 'x'), \n",
     "                                                        np.empty((nts, mg2.shape[0], mg2.shape[1])),\n",
     "                                                       {'units' : 'meters',\n",
     "                                                        'long_name': 'Topographic Elevation'}),\n",
-    "                            'rock_type__id': (('t', 'y', 'x'), \n",
+    "                            'rock_type__id': (('time', 'y', 'x'), \n",
     "                                              np.empty((nts, mg2.shape[0], mg2.shape[1])),\n",
     "                                             {'units' : '-',\n",
     "                                              'long_name' : 'Rock Type ID Code'})\n",
@@ -325,9 +330,10 @@
     "                         'y': (('y'), \n",
     "                               mg2.y_of_node.reshape(mg2.shape)[:, 1],\n",
     "                               {'units' : 'meters'}),\n",
-    "                         'time': (('t'),  \n",
+    "                         'time': (('time'),  \n",
     "                                  dt*np.arange(nts)/1e6,\n",
-    "                                  {'units': 'millions of years since model start'})\n",
+    "                                  {'units': 'millions of years since model start',\n",
+    "                                   'standard_name' : 'time'})\n",
     "                        }\n",
     "                 )\n",
     "\n",
@@ -429,7 +435,9 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "metadata": {},
+   "metadata": {
+    "scrolled": true
+   },
    "outputs": [],
    "source": [
     "imshow_grid(mg2, 'topographic__elevation', cmap='viridis')"