changing the model grid size for memory reasons

Author: kbarnhart

rockblock/introduction_to_rockblock.ipynb

@@ -23,9 +23,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "metadata": {
-    "collapsed": true
-   },
+   "metadata": {},
    "outputs": [],
    "source": [
     "import os\n",
@@ -73,7 +71,7 @@
     "attrs = {'K_sp': {0: 0.0003,\n",
     "                  1: 0.0001}}\n",
     "\n",
-    "mg = RasterModelGrid((100, 100), 200)\n",
+    "mg = RasterModelGrid((70, 70), 200)\n",
     "z = mg.add_zeros('node', 'topographic__elevation') \n",
     "random_field = 0.01*np.random.randn(mg.size('node'))\n",
     "z += random_field - random_field.min()\n",
@@ -83,7 +81,7 @@
     "ids = np.tile([0,1], 20) \n",
     "\n",
     "# Anticline\n",
-    "rb = LayeredRockBlock(mg, z0s, ids, x0=10000, y0=10000, function = lambda x, y : ((0.0015*x)**2+(0.0007*y)**2), attrs=attrs)\n",
+    "rb = LayeredRockBlock(mg, z0s, ids, x0=7000, y0=7000, function = lambda x, y : ((0.002*x)**2+(0.001*y)**2), attrs=attrs)\n",
     "\n",
     "# Shallow dips\n",
     "#rb = LayeredRockBlock(mg, z0s, ids, x0=5000, y0=5000, function = lambda x, y : ((0.001*x)+(0.003*y)), attrs=attrs)\n",
@@ -104,9 +102,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "metadata": {
-    "collapsed": true
-   },
+   "metadata": {},
    "outputs": [],
    "source": [
     "imshow_grid(mg, 'K_sp')"
@@ -191,9 +187,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "metadata": {
-    "collapsed": true
-   },
+   "metadata": {},
    "outputs": [],
    "source": [
     "imshow_grid(mg, 'topographic__elevation', cmap='viridis')"
@@ -213,9 +207,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "metadata": {
-    "collapsed": true
-   },
+   "metadata": {},
    "outputs": [],
    "source": [
     "hvds = hv.Dataset(ds)\n",
@@ -232,9 +224,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "metadata": {
-    "collapsed": true
-   },
+   "metadata": {},
    "outputs": [],
    "source": [
     "%opts Image (cmap='viridis')\n",
@@ -272,7 +262,7 @@
    },
    "outputs": [],
    "source": [
-    "mg2 = RasterModelGrid((100, 100), 200)\n",
+    "mg2 = RasterModelGrid((70, 70), 200)\n",
     "mg2.set_closed_boundaries_at_grid_edges(True, False, True, False)\n",
     "z2 = mg2.add_zeros('node', 'topographic__elevation') \n",
     "random_field = 0.01*np.random.randn(mg2.size('node'))\n",
@@ -331,9 +321,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "metadata": {
-    "collapsed": true
-   },
+   "metadata": {},
    "outputs": [],
    "source": [
     "imshow_grid(mg2, 'topographic__elevation')"
@@ -351,9 +339,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "metadata": {
-    "collapsed": true
-   },
+   "metadata": {},
    "outputs": [],
    "source": [
     "volcanic_deposits = np.zeros(mg2.size('node'))\n",
@@ -408,9 +394,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "metadata": {
-    "collapsed": true
-   },
+   "metadata": {},
    "outputs": [],
    "source": [
     "imshow_grid(mg2, 'topographic__elevation')"
@@ -426,9 +410,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "metadata": {
-    "collapsed": true
-   },
+   "metadata": {},
    "outputs": [],
    "source": [
     "%opts Image (cmap='viridis')\n",