{
  "nbformat": 4,
  "nbformat_minor": 0,
  "metadata": {
    "colab": {
      "provenance": []
    },
    "kernelspec": {
      "name": "python3",
      "display_name": "Python 3"
    },
    "language_info": {
      "name": "python"
    }
  },
  "cells": [
    {
      "cell_type": "code",
      "execution_count": 9,
      "metadata": {
        "id": "EHpYKhbDBXUA"
      },
      "outputs": [],
      "source": [
        "import numpy as np\n",
        "import matplotlib.pyplot as plt"
      ]
    },
    {
      "cell_type": "code",
      "source": [
        "dx = 2000.0\n",
        "a = 125000.0  #espessura da litosfera\n",
        "T1 = 1300.0   #temperatura basal\n",
        "\n",
        "x = np.arange(0.0,1500.0E3,dx)\n",
        "z = np.arange(0.0, a, dx)\n",
        "\n",
        "X, Z = np.meshgrid(x,z)\n",
        "\n",
        "beta = np.copy(x)*0 +1.0\n",
        "betamax = 3.0 \n",
        "beta[x>700.0E3] = betamax\n",
        "\n",
        "T = np.copy(Z)*T1*beta/a\n",
        "T[T>T1] = T1\n",
        "\n",
        "plt.figure(figsize=(10,3))\n",
        "plt.contourf(X/1000.0,Z/1000.0,T,100)\n",
        "plt.ylim(a/1000.0,0)"
      ],
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 232
        },
        "id": "UyPs2vXbBdE6",
        "outputId": "c32a6b16-96ce-40c7-d050-3adf5ecdd45c"
      },
      "execution_count": 16,
      "outputs": [
        {
          "output_type": "execute_result",
          "data": {
            "text/plain": [
              "(125.0, 0.0)"
            ]
          },
          "metadata": {},
          "execution_count": 16
        },
        {
          "output_type": "display_data",
          "data": {
            "text/plain": [
              "<Figure size 720x216 with 1 Axes>"
            ],
            "image/png": "iVBORw0KGgoAAAANSUhEUgAAAlYAAADGCAYAAAAQXM51AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAUaUlEQVR4nO3df4xd5X3n8fdn7njG/DY/stSx0Zps2VZpqjQsoqyy6qLQTWiK4v0DIZooCwmV1Z/bX6sAjbTpSq0UtlVSVrtKZIW0pKIxhKYFddOmlA1KV1poDPlFIE5cfgRbJk5TIFFpA/Z8949zDDdm7BnfueNzzvj9kkZzznPOPfeZR8+c+z3f8zznpqqQJEnSys10XQFJkqS1wsBKkiRpSgysJEmSpsTASpIkaUoMrCRJkqbEwEqSJGlKVi2wSnJ5kl1Jdie5YbXeR5IkqS+yGs+xSjICvgb8B2AP8DngZ6rqkam/mSRJUk+sVsbqYmB3VT1WVS8AO4Ctq/RekiRJvbBagdUm4Kmx9T1tmSRJ0po129UbJ9kGbAMYMfo3J3N6V1WR1DOZcV6NpP76zsK3/76qXrXYttUKrPYC542tb27LXlJV24HtAGfMnF2XzF2+SlWRNDQ56aSuqyBJR/TpZ2958kjbViuw+hxwQZLzaQKqq4G3r9J7SVpjsn6+6ypI0kRWJbCqqgNJfgn4NDACPlpVXzniCxKY7eyupKS+8XwgaaBW7exVVZ8CPrVax5e0NmXdOjBjJWmg+nFZmDQnU0kCas7zgaRh6kdgJUmtzM5Ss84KlDRM/Qis0pxMJQlgYd6MlaRhMpqR1C+zIzNWkgarJ4FVYHbUdSUk9cHsLAvzng8kDVNPAitJelnNpusqSNJE+hFY+RwrSYfMjjg4561AScPUj2gmeCtQUmN2loPzBlaShqkfgZUkjamRtwIlDVM/AitvBUo6ZGaGg/MGVpKGyWhGUu8sODJA0kD1JLAKzDimQhLU7IgFnw8qaaB6ElhJUmt2hoNz3gqUNEz9CKzSXKVKEkB5OpA0UP0IrCSpVbMzHPRWoKSB6kdglYDfDSapZcZK0lD1I7CSpFaNRizMdV0LSZrMxIFVkvOAjwHnAgVsr6qbk5wF3A5sAZ4ArqqqZ45+MPw2e0lAcy4wsJI0VCvJWB0AfqOqHkpyGvBgknuAa4F7q+r9SW4AbgCuX3lVJZ0oyussSQM1cWBVVfuAfe3yd5M8CmwCtgKXtrvdCtzHEoFVEWrkoApJULPhoBkrSQM1lTFWSbYAbwAeAM5tgy6Ap2luFUrSstWouq6CJE1kxYFVklOBPwF+taq+k7z8YL+qqiSLniGTbAO2AayfO8MxVpIAqNEMNWdgJWmYVhRYJVlHE1TdVlWfbIu/mWRjVe1LshHYv9hrq2o7sB3g9NM2Vc36pGVJza3AhVkDK0nDtJJZgQFuAR6tqg+MbbobuAZ4f/v7rhXVUNKJxyGXkgZqJRmrNwLvBL6c5Att2W/SBFR3JLkOeBK4aulDhRp5K1ASLIwCcwtdV0OSJrKSWYH/FzjS/bvLJj2uJOHgdUkD1Y8nr6cZVyFJNQoz6w50XQ1Jmkg/AitJatVsmJ072HU1JGkivQisKu24CkkCRjOOsZI0TL0IrCTpkIVRmJ/zVqCkYepNYFVmrCTRnAvWzXgrUNIw9SawkqRDTpp7sesqSNJE+hFYOStQUmthFk6dNbCSNEz9CKwkacyctwIlDVQvAqtKnBUoCWjGWJ267ntdV0OSJtKLwEqSxpmxkjRU/Qis4qxASY2FEZw++89dV0OSJtKPwEqSWjUD8yOfYyVpmHoTWC30piaSunb66J+6roIkTaQf4Yy3AiW1agRnzD7fdTUkaSL9CKwkqbUwCuvjc6wkDVMvAqvmS5i7roWkvtgwMmMlaZh6EVhJ0iE1gvUzZqwkDdOKA6skI2AnsLeqrkhyPrADOBt4EHhnVb2w1HFqZqU1kbRWnDbj4HVJwzSNjNWvAI8Cp7frNwEfrKodST4MXAd8aArvI+kEUCPYYGAlaaBWFFgl2Qz8NPA7wK8nCfAm4O3tLrcCv8VSgVWak6kk1QjWxyevSxqmlWasfh94D3Bau3428GxVHXq63x5g0wrfQ9IJZsNMdV0FSZrIxIFVkiuA/VX1YJJLJ3j9NmAbwLpTz/RLmCUBsDAD83HQpaRhWknG6o3A25K8FVhPM8bqZmBDktk2a7UZ2LvYi6tqO7Ad4ORXneflqaSXrI9jAyQN08SBVVXdCNwI0Gas/ktVvSPJJ4AraWYGXgPcteTBHGMlqVUjmM+6rqshSRNZjXz79TQD2XfTjLm6ZRXeQ5IkqXem8oDQqroPuK9dfgy4+JiPYcZKEp4LJA2bI0QlSZKmpB9faeMYK0ktzwWShqwXgVXRTLGWJEkaMsMZSb1ixkrSkPUiY+WtQEmStBaYsZLUK15kSRoyM1aSJElTYsZKUq+UZyVJA9aPjBVmrCRJ0vB5bSipV7zIkjRk/chYOcZKkiStAWasJPXKghdZkgasFxmrwoyVJEkaPjNWknqlRtV1FSRpYr3IWBGnWEuSpOEznJHUKw4LkDRk/chY4clUkiQN34oCqyQbgI8Ar6MZg/5uYBdwO7AFeAK4qqqeOfqBnAkkqeUYK0kDttJbgTcDf1lVPwy8HngUuAG4t6ouAO5t1yVJkta8iTNWSc4AfgK4FqCqXgBeSLIVuLTd7VbgPuD6ox2r4kwgSY1K1zWQpMmtJGN1PvAt4A+SfD7JR5KcApxbVfvafZ4Gzl1pJSVJkoZgJWOsZoELgV+uqgeS3Mxht/2qqpIsmopKsg3YBjA680wHr0sCzF5LGraVZKz2AHuq6oF2/U6aQOubSTYCtL/3L/biqtpeVRdV1UWjU05ZQTUkSZL6YeKMVVU9neSpJD9UVbuAy4BH2p9rgPe3v+9a8mDBmUCSGmavJQ3YSp9j9cvAbUnmgMeAd9Fkwe5Ich3wJHDVCt9DkiRpEFYUWFXVF4CLFtl02TEfy5lAknCMlaRh8yttJEmSpqQfX2njc6wkteK5QNKAmbGSJEmakn5krMCZQJIaZqwkDZgZK0mSpCnpR8Yq5RgrSQDMzCx0XQVJmpgZK0mSpCnpR8YKZwJJangukDRk/Qis/EobSZK0BngrUFKvzIwcYyVpuHqTsYoDViVJ0sCZsZLUKzMOC5A0YP3IWFEOWJUkSYNnxkpSr4wcFiBpwHqRsUocsCpJkobPjJWkXpn1IkvSgPUiYwUOWJUkScNnxkpSr4zMWEkasBVlrJL8GvCzQAFfBt4FbAR2AGcDDwLvrKoXjn6ccsCqJEkavIkzVkk2Af8ZuKiqXgeMgKuBm4APVtUPAs8A102jopJODOtmDnZdBUma2ErHWM0CJyV5ETgZ2Ae8CXh7u/1W4LeADy15INP/kiRp4CbOWFXVXuD3gG/QBFTP0dz6e7aqDrS77QE2Lfb6JNuS7Eyy8+B3np+0GpLWGC+yJA3ZxBmrJGcCW4HzgWeBTwCXL/f1VbUd2A5w8gWvLgesSpKkoVvJrcCfBB6vqm8BJPkk8EZgQ5LZNmu1Gdi71IFCOa5CEgDrRp4LJA3XSh638A3gkiQnJwlwGfAI8Bngynafa4C7VlZFSZKkYZg4Y1VVDyS5E3gIOAB8nubW3v8GdiT57bbslqWOlTiuQlLD7LWkIVvRrMCqeh/wvsOKHwMuXslxJUmShqgXX2mTlOMqJAEw68OCJQ2YX2kjSZI0Jf3IWDkrUFJrznOBpAEzYyVJkjQlPclYOa5CUmPO8ZaSBsyMlSRJ0pT0I2OVclyFJADmZg4svZMk9ZQZK0mSpCnpR8YKx1VIapi9ljRkZqwkSZKmpB8Zq5TjKiQBMO+5QNKAmbGSJEmakl5krGZwVqAkSRq+XgRWiel/SQ3PBZKGzFuBknplfmRgJWm4+pGxorxKlSRJg2fGSlKvzOfFrqsgSRNbMmOV5KPAFcD+qnpdW3YWcDuwBXgCuKqqnkkS4GbgrcDzwLVV9dBS7zGTMv0vSZIGbzkZqz8ELj+s7Abg3qq6ALi3XQf4KeCC9mcb8KHpVFOSJKn/lsxYVdVnk2w5rHgrcGm7fCtwH3B9W/6xqirg/iQbkmysqn1He49Qpv8lAc4KlDRsk46xOncsWHoaOLdd3gQ8NbbfnrbsFZJsS7Izyc7nn/nehNWQJEnqjxXPCqyqSlITvG47sB1g049sKK9SJQGsnzF7LWm4Js1YfTPJRoD29/62fC9w3th+m9sySZKkNW/SjNXdwDXA+9vfd42V/1KSHcCPA88tNb4Kmi9h9ipVEsC85wJJA7acxy18nGag+jlJ9gDvowmo7khyHfAkcFW7+6doHrWwm+ZxC+9ahTpLkiT10nJmBf7METZdtsi+BfzisVZihvIqVRIA650hLGnAfPK6JEnSlPTmuwK9SpUEzgqUNGy9CKxmHLwuSZLWAG8FSuoVs9eShqwfGStvBUqSpDXAjJWkXnFYgKQh60XGKjjGSpIkDZ8ZK0m9sj4vdF0FSZpYLzJWzRgrT6aSJGnYzFhJ6pV1Odh1FSRpYr3IWCXlyVSSJA2eGStJveKjVyQNWT8yVj7HSpIkrQFmrCT1yvoc6LoKkjSxXmSsmlmBnkwlSdKwmbGS1CvrncgiacB6krHyZCpJkoZvyYxVko8m2Z/k4bGy303y1SRfSvKnSTaMbbsxye4ku5K8ZbUqLmltWp+FrqsgSRNbTsbqD4H/CXxsrOwe4MaqOpDkJuBG4PokrwWuBn4EeDXw10n+dVUdNR3VzAr0ZCpJkoZtycCqqj6bZMthZX81tno/cGW7vBXYUVXfAx5Pshu4GPh/R3uPmYT1yTFUW5IkqX+mMXj93cBftMubgKfGtu1py14hybYkO5Ps/Pa3zVZJkqThW9Hg9STvBQ4Atx3ra6tqO7Ad4MdeP1fzcYKiJEkatokDqyTXAlcAl1VVtcV7gfPGdtvclkmSJK15EwVWSS4H3gP8+6p6fmzT3cAfJ/kAzeD1C4C/Xep4M4T1GU1SFUmSpN5YMrBK8nHgUuCcJHuA99HMApwH7kkz6Pz+qvq5qvpKkjuAR2huEf7iUjMCJUmS1oq8fBevOxe+fr7+5i9/oOtqSOqB+azrugqSdFTrNv7dg1V10WLb+vHk9XU/yikbd3ZdDUmSpGU48iOiepGxSvJdYFfX9RiIc4C/77oSA2A7LZ9ttTy20/LZVstjOy1f39rqX1bVqxbb0IuMFbDrSCk1fb8kO22rpdlOy2dbLY/ttHy21fLYTss3pLby4VGSJElTYmAlSZI0JX0JrLZ3XYEBsa2Wx3ZaPttqeWyn5bOtlsd2Wr7BtFUvBq9LkiStBX3JWEmSJA1e54FVksuT7EqyO8kNXdenS0nOS/KZJI8k+UqSX2nLz0pyT5Kvt7/PbMuT5H+0bfelJBd2+xccX0lGST6f5M/b9fOTPNC2x+1J5try+XZ9d7t9S5f1Pt6SbEhyZ5KvJnk0yb+1T71Skl9r/+8eTvLxJOvtU40kH02yP8nDY2XH3IeSXNPu//Uk13Txt6y2I7TV77b/f19K8qdJNoxtu7Ftq11J3jJWvqY/Gxdrp7Ftv5GkkpzTrg+rT1VVZz/ACPg74DXAHPBF4LVd1qnj9tgIXNgunwZ8DXgt8N+BG9ryG4Cb2uW3An9B86SyS4AHuv4bjnN7/Trwx8Cft+t3AFe3yx8Gfr5d/gXgw+3y1cDtXdf9OLfTrcDPtstzwAb71CvaaBPwOHDSWF+61j71Uvv8BHAh8PBY2TH1IeAs4LH295nt8pld/23Hqa3eDMy2yzeNtdVr28+9eeD89vNwdCJ8Ni7WTm35ecCngSeBc4bYp7rOWF0M7K6qx6rqBWAHsLXjOnWmqvZV1UPt8neBR2lO+FtpPhxpf//Hdnkr8LFq3A9sSLLxOFe7E0k2Az8NfKRdD/Am4M52l8Pb6VD73Qlc1u6/5iU5g+YEdgtAVb1QVc9in1rMLHBSklngZGAf9ikAquqzwD8cVnysfegtwD1V9Q9V9QxwD3D56tf++Fqsrarqr6rqQLt6P7C5Xd4K7Kiq71XV48Bums/FNf/ZeIQ+BfBB4D3A+ADwQfWprgOrTcBTY+t72rITXntr4Q3AA8C5VbWv3fQ0cG67fCK33+/T/PMttOtnA8+OnbzG2+Kldmq3P9fufyI4H/gW8AftbdOPJDkF+9T3qaq9wO8B36AJqJ4DHsQ+dTTH2odOyL61iHfTZF/Atvo+SbYCe6vqi4dtGlQ7dR1YaRFJTgX+BPjVqvrO+LZq8p8n9FTOJFcA+6vqwa7rMgCzNOn2D1XVG4B/pLlt8xL7FLTjg7bSBKKvBk6hB1e+Q2EfWp4k7wUOALd1XZe+SXIy8JvAf+26LivVdWC1l+Z+6iGb27ITVpJ1NEHVbVX1ybb4m4dux7S/97flJ2r7vRF4W5InaFLkbwJupkkPH/qapvG2eKmd2u1nAN8+nhXu0B5gT1U90K7fSRNo2ae+308Cj1fVt6rqReCTNP3MPnVkx9qHTtS+BUCSa4ErgHe0gSjYVuP+Fc2FzRfbc/tm4KEkP8DA2qnrwOpzwAXtzJs5mkGgd3dcp860YzRuAR6tqg+MbbobODTb4RrgrrHy/9TOmLgEeG4sNb9mVdWNVbW5qrbQ9Jn/U1XvAD4DXNnudng7HWq/K9v9T4ir66p6GngqyQ+1RZcBj2CfOtw3gEuSnNz+Hx5qJ/vUkR1rH/o08OYkZ7YZwje3ZWtekstphi68raqeH9t0N3B1O8v0fOAC4G85AT8bq+rLVfUvqmpLe27fQzOZ62mG1qe6Hj1PM9r/azQzIN7bdX06bot/R5NO/xLwhfbnrTRjN+4Fvg78NXBWu3+A/9W23ZeBi7r+Gzpos0t5eVbga2hOSruBTwDzbfn6dn13u/01Xdf7OLfRjwE72371ZzSzZ+xTr2yn/wZ8FXgY+COamVr2qebv/TjN2LMXaT7wrpukD9GML9rd/ryr67/rOLbVbpqxQIfO6x8e2/+9bVvtAn5qrHxNfzYu1k6HbX+Cl2cFDqpP+eR1SZKkKen6VqAkSdKaYWAlSZI0JQZWkiRJU2JgJUmSNCUGVpIkSVNiYCVJkjQlBlaSJElTYmAlSZI0Jf8fE3sUwEAxsAEAAAAASUVORK5CYII=\n"
          },
          "metadata": {
            "needs_background": "light"
          }
        }
      ]
    },
    {
      "cell_type": "code",
      "source": [],
      "metadata": {
        "id": "beLuDODECOE7"
      },
      "execution_count": 4,
      "outputs": []
    }
  ]
}