workout-depths.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "92778a84",
   "metadata": {},
   "outputs": [],
   "source": [
    "import src.depthdistr as dd"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "8c73150c",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX4AAAEGCAYAAABiq/5QAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Z1A+gAAAACXBIWXMAAAsTAAALEwEAmpwYAAAT/ElEQVR4nO3dcbBc5X3e8e+DVIQTx1CL29RFkCsHOYxIbMWWcdKQ1DVjVxjHshNRQ9wx6TCjemqatCmTym1DiCaZgUxr6qlJXdVQE9IEUbl2bo0SNRO5E6fjyLpgKBZUzbUsB9EkCKHIxg5gwa9/7CHe7Fzlru5dsdK+38+M5p7znvec/b2wes7RObvvTVUhSWrHWeMuQJL00jL4JakxBr8kNcbgl6TGGPyS1Jjl4y5g0Pnnn1/T09PjLkOSzij333//k1U1NUzf0y74p6enmZ2dHXcZknRGSfKVYft6q0eSGmPwS1JjDH5JaozBL0mNMfglqTEGvyQ1xuCXpMYY/JLUGINfkhpz2n1zVydnest9Y3vtg7dcNbbXlrR4XvFLUmMMfklqjMEvSY0x+CWpMQa/JDXG4Jekxhj8ktQYg1+SGmPwS1JjDH5JaozBL0mNMfglqTFDBX+SDUn2J5lLsmWe7SuSbO+270ky3bfttUk+l2RfkoeTnDPC+iVJJ2nB4E+yDLgduBJYC1ybZO1At+uBo1V1MXAbcGu373Lg14D3V9WlwJuBb46seknSSRvmiv8yYK6qDlTVc8A9wMaBPhuBu7rlHcAVSQK8DfjfVfUQQFUdqarnR1O6JGkxhgn+C4DH+tYPdW3z9qmq48AxYCXwGqCS7EryQJKfne8FkmxOMptk9vDhwyc7BknSSTjVD3eXA5cD7+1+vjvJFYOdqmpbVa2vqvVTU1OnuCRJatswwf84cGHf+qqubd4+3X39c4Ej9P518HtV9WRVfQPYCbx+qUVLkhZvmODfC6xJsjrJ2cA1wMxAnxngum55E7C7qgrYBXxfkm/rTgh/B3hkNKVLkhZjwd+5W1XHk9xAL8SXAXdW1b4kW4HZqpoB7gDuTjIHPEXv5EBVHU3yIXonjwJ2VtX4fkmsJGm4X7ZeVTvp3abpb7upb/kZ4OoT7Ptr9D7SKUk6DfjNXUlqjMEvSY0x+CWpMQa/JDXG4Jekxhj8ktQYg1+SGmPwS1JjDH5JaozBL0mNMfglqTEGvyQ1xuCXpMYY/JLUGINfkhpj8EtSYwx+SWqMwS9JjTH4JakxBr8kNcbgl6TGGPyS1Jihgj/JhiT7k8wl2TLP9hVJtnfb9ySZ7tqnk/x5kge7Px8dcf2SpJO0fKEOSZYBtwNvBQ4Be5PMVNUjfd2uB45W1cVJrgFuBd7TbftSVa0bbdmSpMUa5or/MmCuqg5U1XPAPcDGgT4bgbu65R3AFUkyujIlSaOy4BU/cAHwWN/6IeBNJ+pTVceTHANWdttWJ/kC8FXgX1fVZ5dW8ulpest94y5BkoYyTPAvxR8DF1XVkSRvAD6V5NKq+mp/pySbgc0AF1100SkuSZLaNsytnseBC/vWV3Vt8/ZJshw4FzhSVc9W1RGAqrof+BLwmsEXqKptVbW+qtZPTU2d/CgkSUMbJvj3AmuSrE5yNnANMDPQZwa4rlveBOyuqkoy1T0cJsmrgTXAgdGULklajAVv9XT37G8AdgHLgDural+SrcBsVc0AdwB3J5kDnqJ3cgD4EWBrkm8CLwDvr6qnTsVAJEnDGeoef1XtBHYOtN3Ut/wMcPU8+30C+MQSa5QkjZDf3JWkxhj8ktQYg1+SGmPwS1JjDH5JaozBL0mNMfglqTEGvyQ1xuCXpMYY/JLUGINfkhpj8EtSYwx+SWqMwS9JjTH4JakxBr8kNcbgl6TGGPyS1BiDX5IaY/BLUmMMfklqzPJxF6Az1/SW+8byugdvuWosrytNCq/4JakxQwV/kg1J9ieZS7Jlnu0rkmzvtu9JMj2w/aIkTye5cUR1S5IWacHgT7IMuB24ElgLXJtk7UC364GjVXUxcBtw68D2DwG/tfRyJUlLNcwV/2XAXFUdqKrngHuAjQN9NgJ3dcs7gCuSBCDJu4AvA/tGUrEkaUmGCf4LgMf61g91bfP2qarjwDFgZZKXA/8C+IW/6gWSbE4ym2T28OHDw9YuSVqEU/1w92bgtqp6+q/qVFXbqmp9Va2fmpo6xSVJUtuG+Tjn48CFfeururb5+hxKshw4FzgCvAnYlOSXgfOAF5I8U1UfWWrhkqTFGSb49wJrkqymF/DXAD8x0GcGuA74HLAJ2F1VBfzwix2S3Aw8behL0ngtGPxVdTzJDcAuYBlwZ1XtS7IVmK2qGeAO4O4kc8BT9E4OkqTT0FDf3K2qncDOgbab+pafAa5e4Bg3L6I+SdKI+c1dSWqMwS9JjTH4JakxBr8kNcbgl6TGGPyS1BiDX5IaY/BLUmMMfklqjMEvSY0x+CWpMQa/JDXG4Jekxhj8ktQYg1+SGmPwS1JjDH5JaozBL0mNMfglqTEGvyQ1xuCXpMYY/JLUGINfkhozVPAn2ZBkf5K5JFvm2b4iyfZu+54k0137ZUke7P48lOTdI65fknSSFgz+JMuA24ErgbXAtUnWDnS7HjhaVRcDtwG3du1fBNZX1TpgA/AfkywfUe2SpEUY5or/MmCuqg5U1XPAPcDGgT4bgbu65R3AFUlSVd+oquNd+zlAjaJoSdLiDRP8FwCP9a0f6trm7dMF/TFgJUCSNyXZBzwMvL/vRPAXkmxOMptk9vDhwyc/CknS0E75w92q2lNVlwJvBD6Y5Jx5+myrqvVVtX5qaupUlyRJTRsm+B8HLuxbX9W1zdunu4d/LnCkv0NVPQo8DXzvYouVJC3dMMG/F1iTZHWSs4FrgJmBPjPAdd3yJmB3VVW3z3KAJN8FXAIcHEnlkqRFWfATNlV1PMkNwC5gGXBnVe1LshWYraoZ4A7g7iRzwFP0Tg4AlwNbknwTeAH4x1X15KkYiCRpOEN9tLKqdgI7B9pu6lt+Brh6nv3uBu5eYo2SpBHym7uS1BiDX5IaY/BLUmMMfklqjMEvSY0x+CWpMQa/JDXG4Jekxhj8ktQYg1+SGmPwS1JjDH5JaozBL0mNMfglqTEGvyQ1xuCXpMYY/JLUGINfkhpj8EtSYwx+SWqMwS9JjTH4JakxBr8kNWao4E+yIcn+JHNJtsyzfUWS7d32PUmmu/a3Jrk/ycPdz7eMuH5J0klaMPiTLANuB64E1gLXJlk70O164GhVXQzcBtzatT8J/GhVfR9wHXD3qAqXJC3OMFf8lwFzVXWgqp4D7gE2DvTZCNzVLe8ArkiSqvpCVf2/rn0f8LIkK0ZRuCRpcYYJ/guAx/rWD3Vt8/apquPAMWDlQJ8fBx6oqmcHXyDJ5iSzSWYPHz48bO2SpEV4SR7uJrmU3u2ffzTf9qraVlXrq2r91NTUS1GSJDVrmOB/HLiwb31V1zZvnyTLgXOBI936KuCTwPuq6ktLLViStDTDBP9eYE2S1UnOBq4BZgb6zNB7eAuwCdhdVZXkPOA+YEtV/a8R1SxJWoIFg7+7Z38DsAt4FLi3qvYl2ZrknV23O4CVSeaAnwFe/MjnDcDFwE1JHuz+/I2Rj0KSNLTlw3Sqqp3AzoG2m/qWnwGunme/XwR+cYk1SpJGyG/uSlJjhrril04n01vuG9trH7zlqrG9tjQqXvFLUmMMfklqjMEvSY0x+CWpMQa/JDXG4Jekxhj8ktQYg1+SGmPwS1JjDH5JaozBL0mNMfglqTEGvyQ1ZuJm5xznzI2SdCbwil+SGmPwS1JjDH5JaozBL0mNMfglqTEGvyQ1xuCXpMYMFfxJNiTZn2QuyZZ5tq9Isr3bvifJdNe+Mslnkjyd5CMjrl2StAgLBn+SZcDtwJXAWuDaJGsHul0PHK2qi4HbgFu79meAnwNuHFnFkqQlGeaK/zJgrqoOVNVzwD3AxoE+G4G7uuUdwBVJUlVfr6rfp3cCkCSdBoYJ/guAx/rWD3Vt8/apquPAMWDlsEUk2ZxkNsns4cOHh91NkrQIp8XD3araVlXrq2r91NTUuMuRpIk2TPA/DlzYt76qa5u3T5LlwLnAkVEUKEkarWGCfy+wJsnqJGcD1wAzA31mgOu65U3A7qqq0ZUpSRqVBadlrqrjSW4AdgHLgDural+SrcBsVc0AdwB3J5kDnqJ3cgAgyUHgFcDZSd4FvK2qHhn5SCRJQxlqPv6q2gnsHGi7qW/5GeDqE+w7vYT6JEkjdlo83JUkvXQMfklqjMEvSY0x+CWpMQa/JDVmqE/1SOqZ3nLfWF734C1XjeV1NZm84pekxhj8ktQYg1+SGmPwS1JjDH5JaozBL0mNMfglqTEGvyQ1xuCXpMYY/JLUGINfkhpj8EtSY5ykTToDjGtyOHCCuEnkFb8kNcbgl6TGGPyS1BiDX5IaM9TD3SQbgA8Dy4CPVdUtA9tXAL8KvAE4Arynqg522z4IXA88D/xUVe0aWfWSTjl/69jkWfCKP8ky4HbgSmAtcG2StQPdrgeOVtXFwG3Ard2+a4FrgEuBDcCvdMeTJI3JMFf8lwFzVXUAIMk9wEbgkb4+G4Gbu+UdwEeSpGu/p6qeBb6cZK473udGU76kSTXOj7COy0v1r5xhgv8C4LG+9UPAm07Up6qOJzkGrOza/2Bg3wsGXyDJZmBzt/p0kv1DVT9a5wNPjuF1x6m1Mbc2XmhvzGf0eHPronZ7cczfNewOp8UXuKpqG7BtnDUkma2q9eOs4aXW2phbGy+0N+bWxguLG/Mwn+p5HLiwb31V1zZvnyTLgXPpPeQdZl9J0ktomODfC6xJsjrJ2fQe1s4M9JkBruuWNwG7q6q69muSrEiyGlgDfH40pUuSFmPBWz3dPfsbgF30Ps55Z1XtS7IVmK2qGeAO4O7u4e1T9E4OdP3upfcg+Djwgap6/hSNZanGeqtpTFobc2vjhfbG3Np4YRFjTu/CXJLUCr+5K0mNMfglqTFNBn+SO5M8keSLfW2vTPI7Sf6w+/nXx1njKCW5MMlnkjySZF+Sn+7aJ3nM5yT5fJKHujH/Qte+OsmeJHNJtncfWJgYSZYl+UKST3frkz7eg0keTvJgktmubZLf1+cl2ZHk/yR5NMkPLma8TQY/8HF6U0j02wL8blWtAX63W58Ux4F/XlVrgR8APtBNpzHJY34WeEtVvQ5YB2xI8gP0phO5rZte5Ci96UYmyU8Dj/atT/p4Af5uVa3r+yz7JL+vPwz8dlVdAryO3v/rkx9vVTX5B5gGvti3vh94Vbf8KmD/uGs8hWP/TeCtrYwZ+DbgAXrfOH8SWN61/yCwa9z1jXCcq7q/+G8BPg1kksfbjekgcP5A20S+r+l9P+rLdB/KWcp4W73in893VtUfd8t/AnznOIs5VZJMA98P7GHCx9zd9ngQeAL4HeBLwJ9V1fGuy7xTiJzB/h3ws8AL3fpKJnu8AAX8jyT3d1O/wOS+r1cDh4H/3N3O+1iSb2cR4zX451G9U+fEfc41ycuBTwD/tKq+2r9tEsdcVc9X1Tp6V8KXAZeMt6JTJ8k7gCeq6v5x1/ISu7yqXk9v9uAPJPmR/o0T9r5eDrwe+A9V9f3A1xm4rTPseA3+b/nTJK8C6H4+MeZ6RirJX6MX+v+lqv5b1zzRY35RVf0Z8Bl6tzrO66YVgcmaQuSHgHcmOQjcQ+92z4eZ3PECUFWPdz+fAD5J7wQ/qe/rQ8ChqtrTre+gdyI46fEa/N/SP+3EdfTug0+EborsO4BHq+pDfZsmecxTSc7rll9G75nGo/ROAJu6bhMz5qr6YFWtqqppet+c311V72VCxwuQ5NuTfMeLy8DbgC8yoe/rqvoT4LEk39M1XUFvVoSTHm+T39xN8hvAm+lNZ/qnwM8DnwLuBS4CvgL8/ap6akwljlSSy4HPAg/zrfu//5Leff5JHfNrgbvoTTNyFnBvVW1N8mp6V8SvBL4A/IPq/b6IiZHkzcCNVfWOSR5vN7ZPdqvLgV+vql9KspLJfV+vAz4GnA0cAP4h3fubkxhvk8EvSS3zVo8kNcbgl6TGGPyS1BiDX5IaY/BLUmMMfk2UJDcnuXGR+65L8vaTPVZ6did5RZLp/llfF1HDv0nylsXuLw3D4Je+ZR3w9oU6zePtwEOD02As0r9nsmaT1GnI4NcZL8m/SvJ/k/w+8D197d+d5Le7Cbw+m+SSrv3jST6aZLbb7x3dPPVbgfd0c7u/pzvM2iT/M8mBJD91ghLeyzzflkzy6m4yrTcm+ckkn+rmSz+Y5IYkP9Nt/4MkrwSoqq8AK5P8zRH+J5L+EoNfZ7Qkb6A3RcE6elfeb+zbvA34J1X1BuBG4Ff6tk3Tm9flKuCj9P4u3ARsr97c7tu7fpcAf6/r+/PdnEeDfgj4S5OjdV+r/wTwk1W1t2v+XuDHuhp/CfhGN9nW54D39e3+QHdM6ZRYvnAX6bT2w8Anq+obAElmup8vB/428F97UxUBsKJvv3ur6gXgD5Mc4MQzd97XTXHwbJIn6E15e2igzyur6mt961P0/gXwY1X1SF/7Z7p+X0tyDPjvXfvDwGv7+j0B/K0Fxi0tmsGvSXUWvbno151g++BcJSeau6R/Xpvnmf/vzPEkZ3UnEoBjwB8Bl9ObRGu+Y73Qt/7CwHHPAf78BPVIS+atHp3pfg94V5KXdTM1/ihA96D1y0muhr/45M3r+va7OslZSb4beDW932L0NeA7FlHD/u4YL3oOeDfwviQ/sYjjvYbeLJPSKWHw64xWVQ8A24GHgN8C9vZtfi9wfZKHgH3Axr5tfwR8vtvn/VX1DL0pjNcOPNwdxn30Znvtr+vrwDuAf5bkncMeqHuGcDEwexKvL50UZ+dUc5J8HPh0Ve0Y0fFeBfxqVb11BMd6N/D6qvq5pVcmzc8rfmmJut93+p+SvGIEh1sO/NsRHEc6Ia/4JakxXvFLUmMMfklqjMEvSY0x+CWpMQa/JDXm/wODAMr4le4R1QAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "d = dd.sample_depths(1000, mean_depth = 20, isaftershock=False, doplot=True)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "39134779",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYAAAAEGCAYAAABsLkJ6AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Z1A+gAAAACXBIWXMAAAsTAAALEwEAmpwYAAAYcUlEQVR4nO3df5RV5X3v8fdHCDRrpRqRSS8F2sFImlLTi3Egrrba1kSDMQXaguL1Rrz1hqa59La1tiFNg13UrKVJe13NXTSRVOKPapCQGudexkvMUtPee4uZ0aCIljgikUFbJ2KMrRE64ds/9kOyPT0zZ5+Zw5wZns9rrbNm72c/z3Oeh+Gcz+y9z9lbEYGZmeXnpHYPwMzM2sMBYGaWKQeAmVmmHABmZplyAJiZZWpquwfQjJkzZ0ZnZ2e7h2FmNqk8/PDD346IjtrySRUAnZ2d9PX1tXsYZmaTiqRv1Sv3ISAzs0w5AMzMMuUAMDPLlAPAzCxTDgAzs0w5AMzMMuUAMDPLlAPAzCxTDgAzs0xV+iawpCXAXwBTgL+KiOtrtl8N/FdgCBgEfiMivpW2rQb+OFW9LiJuTeVnA7cAbwR6gN8J353Gxqhz3fa2PO/+6y9uy/OajUXDPQBJU4CNwEXAAuAySQtqqn0D6IqInwW2AZ9MbWcA1wLvAhYD10o6NbX5DPBBYH56LBnzbMzMrLIqh4AWA/0RsS8ijgBbgGXlChHxQES8mlZ3AnPS8nuB+yLiUES8BNwHLJE0Czg5Inamv/pvA5aPfTpmZlZVlQCYDRworQ+ksuFcBdzboO3stNywT0lrJPVJ6hscHKwwXDMzq6KlJ4El/WegC/hUq/qMiE0R0RURXR0d/+5qpmZmNkpVAuAgMLe0PieVvY6k9wAfA5ZGxOEGbQ/yw8NEw/ZpZmbHT5UA6AXmS5onaRqwCuguV5B0FnATxZv/C6VNO4ALJZ2aTv5eCOyIiOeB70o6R5KAK4B7WjAfMzOrqOHHQCNiSNJaijfzKcDmiNgjaQPQFxHdFId83gR8sXg/59mIWBoRhyT9KUWIAGyIiENp+cP88GOg9/LD8wZmZjYOKn0PICJ6KD6rXy5bX1p+zwhtNwOb65T3AWdWHqmZmbWUvwlsZpYpB4CZWaYcAGZmmXIAmJllygFgZpYpB4CZWaYcAGZmmXIAmJllygFgZpYpB4CZWaYcAGZmmXIAmJllqtLF4Mya0a4bs5tZc7wHYGaWKQeAmVmmHABmZpmqFACSlkjaK6lf0ro628+T9IikIUkrSuW/LGlX6fGapOVp2y2SniltW9iqSZmZWWMNTwJLmgJsBC4ABoBeSd0R8USp2rPAlcA15bYR8QCwMPUzA+gHvlKq8gcRsW0M4zczs1Gq8imgxUB/ROwDkLQFWAb8IAAiYn/adnSEflYA90bEq6MerZmZtUyVQ0CzgQOl9YFU1qxVwBdqyj4h6TFJN0qaXq+RpDWS+iT1DQ4OjuJpzcysnnE5CSxpFvAOYEep+KPA24FFwAzgI/XaRsSmiOiKiK6Ojo7jPlYzs1xUCYCDwNzS+pxU1oxLgLsj4l+PFUTE81E4DHye4lCTmZmNkyoB0AvMlzRP0jSKQzndTT7PZdQc/kl7BUgSsBx4vMk+zcxsDBoGQEQMAWspDt88CWyNiD2SNkhaCiBpkaQBYCVwk6Q9x9pL6qTYg/haTdd3SNoN7AZmAte1YD5mZlZRpWsBRUQP0FNTtr603EtxaKhe2/3UOWkcEec3M1AzM2stfxPYzCxTDgAzs0w5AMzMMuUAMDPLlAPAzCxTDgAzs0w5AMzMMuUAMDPLlAPAzCxTDgAzs0w5AMzMMuUAMDPLlAPAzCxTDgAzs0w5AMzMMuUAMDPLlAPAzCxTlQJA0hJJeyX1S1pXZ/t5kh6RNCRpRc2270valR7dpfJ5kh5Kfd6V7jdsZmbjpGEASJoCbAQuAhYAl0laUFPtWeBK4M46XXwvIhamx9JS+Q3AjRFxBvAScNUoxm9mZqNUZQ9gMdAfEfsi4giwBVhWrhAR+yPiMeBolSeVJOB8YFsquhVYXnXQZmY2dlUCYDZwoLQ+QJ2bvI/gRyT1SdopaXkqOw34TkQMNepT0prUvm9wcLCJpzUzs5FMHYfn+MmIOCjpdOB+SbuBl6s2johNwCaArq6uOE5jNDPLTpU9gIPA3NL6nFRWSUQcTD/3AQ8CZwEvAm+WdCyAmurTzMzGrkoA9ALz06d2pgGrgO4GbQCQdKqk6Wl5JvDzwBMREcADwLFPDK0G7ml28GZmNnoNAyAdp18L7ACeBLZGxB5JGyQtBZC0SNIAsBK4SdKe1PyngT5Jj1K84V8fEU+kbR8BrpbUT3FO4OZWTszMzEZW6RxARPQAPTVl60vLvRSHcWrb/X/gHcP0uY/iE0ZmZtYG/iawmVmmHABmZplyAJiZZcoBYGaWKQeAmVmmHABmZplyAJiZZcoBYGaWqfG4GJzZCa9z3fa2Pff+6y9u23Pb5OY9ADOzTDkAzMwy5QAwM8uUA8DMLFMOADOzTDkAzMwy5QAwM8tUpQCQtETSXkn9ktbV2X6epEckDUlaUSpfKOnvJe2R9JikS0vbbpH0jKRd6bGwJTMyM7NKGn4RTNIUYCNwATAA9ErqLt3aEeBZ4ErgmprmrwJXRMRTkn4ceFjSjoj4Ttr+BxGxbYxzMDOzUajyTeDFQH+6hSOStgDLgB8EQETsT9uOlhtGxDdLy89JegHoAL4z1oGbmdnYVDkENBs4UFofSGVNkbQYmAY8XSr+RDo0dKOk6cO0WyOpT1Lf4OBgs09rZmbDGJeTwJJmAbcD/yUiju0lfBR4O7AImAF8pF7biNgUEV0R0dXR0TEewzUzy0KVADgIzC2tz0lllUg6GdgOfCwidh4rj4jno3AY+DzFoSYzMxsnVQKgF5gvaZ6kacAqoLtK56n+3cBttSd7014BkgQsBx5vYtxmZjZGDQMgIoaAtcAO4Elga0TskbRB0lIASYskDQArgZsk7UnNLwHOA66s83HPOyTtBnYDM4HrWjkxMzMbWaX7AURED9BTU7a+tNxLcWiott1fA389TJ/nNzVSMzNrKX8T2MwsUw4AM7NMOQDMzDLlADAzy5QDwMwsU5U+BWSTU+e67e0egplNYN4DMDPLlAPAzCxTDgAzs0w5AMzMMuUAMDPLlAPAzCxTDgAzs0w5AMzMMuUAMDPLlAPAzCxTDgAzs0xVCgBJSyTtldQvaV2d7edJekTSkKQVNdtWS3oqPVaXys+WtDv1+el0b2AzMxsnDQNA0hRgI3ARsAC4TNKCmmrPAlcCd9a0nQFcC7wLWAxcK+nUtPkzwAeB+emxZNSzMDOzplXZA1gM9EfEvog4AmwBlpUrRMT+iHgMOFrT9r3AfRFxKCJeAu4DlkiaBZwcETsjIoDbgOVjnIuZmTWhSgDMBg6U1gdSWRXDtZ2dlhv2KWmNpD5JfYODgxWf1szMGpnwJ4EjYlNEdEVEV0dHR7uHY2Z2wqgSAAeBuaX1OamsiuHaHkzLo+nTzMxaoEoA9ALzJc2TNA1YBXRX7H8HcKGkU9PJ3wuBHRHxPPBdSeekT/9cAdwzivGbmdkoNQyAiBgC1lK8mT8JbI2IPZI2SFoKIGmRpAFgJXCTpD2p7SHgTylCpBfYkMoAPgz8FdAPPA3c29KZmZnZiCrdEzgieoCemrL1peVeXn9Ip1xvM7C5TnkfcGYzgzUzs9aZ8CeBzczs+HAAmJllygFgZpYpB4CZWaYcAGZmmXIAmJllygFgZpYpB4CZWaYcAGZmmXIAmJllygFgZpYpB4CZWaYcAGZmmXIAmJllygFgZpYpB4CZWaYqBYCkJZL2SuqXtK7O9umS7krbH5LUmcovl7Sr9DgqaWHa9mDq89i2t7RyYmZmNrKGASBpCrARuAhYAFwmaUFNtauAlyLiDOBG4AaAiLgjIhZGxELgA8AzEbGr1O7yY9sj4oUxz8bMzCqrsgewGOiPiH0RcQTYAiyrqbMMuDUtbwPenW72XnZZamtmZhNAlQCYDRworQ+ksrp10k3kXwZOq6lzKfCFmrLPp8M/H68TGABIWiOpT1Lf4OBgheGamVkV43ISWNK7gFcj4vFS8eUR8Q7g3PT4QL22EbEpIroioqujo2McRmtmlocqAXAQmFtan5PK6taRNBU4BXixtH0VNX/9R8TB9PMV4E6KQ01mZjZOqgRALzBf0jxJ0yjezLtr6nQDq9PyCuD+iAgASScBl1A6/i9pqqSZafkNwPuBxzEzs3EztVGFiBiStBbYAUwBNkfEHkkbgL6I6AZuBm6X1A8cogiJY84DDkTEvlLZdGBHevOfAnwV+FxLZmRmZpU0DACAiOgBemrK1peWXwNWDtP2QeCcmrJ/Ac5ucqxmZtZC/iawmVmmHABmZplyAJiZZcoBYGaWKQeAmVmmKn0KyMwmrs5129vyvPuvv7gtz2ut4z0AM7NMOQDMzDLlADAzy5QDwMwsUw4AM7NMOQDMzDLlADAzy5QDwMwsUw4AM7NMOQDMzDLlADAzy1SlAJC0RNJeSf2S1tXZPl3SXWn7Q5I6U3mnpO9J2pUeny21OVvS7tTm05LUslmZmVlDDQNA0hRgI3ARsAC4TNKCmmpXAS9FxBnAjcANpW1PR8TC9PhQqfwzwAeB+emxZPTTMDOzZlXZA1gM9EfEvog4AmwBltXUWQbcmpa3Ae8e6S96SbOAkyNiZ0QEcBuwvNnBm5nZ6FUJgNnAgdL6QCqrWycihoCXgdPStnmSviHpa5LOLdUfaNAnAJLWSOqT1Dc4OFhhuGZmVsXxPgn8PPATEXEWcDVwp6STm+kgIjZFRFdEdHV0dByXQZqZ5ahKABwE5pbW56SyunUkTQVOAV6MiMMR8SJARDwMPA28LdWf06BPMzM7jqoEQC8wX9I8SdOAVUB3TZ1uYHVaXgHcHxEhqSOdREbS6RQne/dFxPPAdyWdk84VXAHc04L5mJlZRQ1vCRkRQ5LWAjuAKcDmiNgjaQPQFxHdwM3A7ZL6gUMUIQFwHrBB0r8CR4EPRcShtO3DwC3AG4F708PMzMZJpXsCR0QP0FNTtr60/Bqwsk67LwFfGqbPPuDMZgZrZmat428Cm5llqtIegI1e57rt7R6CmVld3gMwM8uUA8DMLFMOADOzTDkAzMwy5QAwM8uUA8DMLFP+GKiZjUo7P+K8//qL2/bcJxLvAZiZZcoBYGaWKQeAmVmmHABmZplyAJiZZcoBYGaWKQeAmVmmKgWApCWS9krql7Suzvbpku5K2x+S1JnKL5D0sKTd6ef5pTYPpj53pcdbWjYrMzNrqOEXwdI9fTcCFwADQK+k7oh4olTtKuCliDhD0irgBuBS4NvAr0TEc5LOpLit5OxSu8vTncHMzGycVdkDWAz0R8S+iDgCbAGW1dRZBtyalrcB75akiPhGRDyXyvcAb5Q0vRUDNzOzsakSALOBA6X1AV7/V/zr6kTEEPAycFpNnV8HHomIw6Wyz6fDPx+XpHpPLmmNpD5JfYODgxWGa2ZmVYzLSWBJP0NxWOg3S8WXR8Q7gHPT4wP12kbEpojoioiujo6O4z9YM7NMVAmAg8Dc0vqcVFa3jqSpwCnAi2l9DnA3cEVEPH2sQUQcTD9fAe6kONRkZmbjpEoA9ALzJc2TNA1YBXTX1OkGVqflFcD9ERGS3gxsB9ZFxP87VlnSVEkz0/IbgPcDj49pJmZm1pSGAZCO6a+l+ATPk8DWiNgjaYOkpanazcBpkvqBq4FjHxVdC5wBrK/5uOd0YIekx4BdFHsQn2vhvMzMrIFK9wOIiB6gp6ZsfWn5NWBlnXbXAdcN0+3Z1YdpZmat5m8Cm5llygFgZpYpB4CZWaYcAGZmmXIAmJllygFgZpYpB4CZWaYqfQ/gRNC5bnu7h2BmLdKu1/P+6y9uy/MeL94DMDPLlAPAzCxTDgAzs0w5AMzMMuUAMDPLlAPAzCxTDgAzs0w5AMzMMpXNF8HMzMbqRPsCWqU9AElLJO2V1C9pXZ3t0yXdlbY/JKmztO2jqXyvpPdW7dPMzI6vhgEgaQqwEbgIWABcJmlBTbWrgJci4gzgRuCG1HYBxU3kfwZYAvylpCkV+zQzs+Ooyh7AYqA/IvZFxBFgC7Csps4y4Na0vA14tySl8i0RcTgingH6U39V+jQzs+OoyjmA2cCB0voA8K7h6kTEkKSXgdNS+c6atrPTcqM+AZC0BliTVv9Z0t4KY55oZgLfbvcgxshzmDhOhHl4Dk3QDWPu4ifrFU74k8ARsQnY1O5xjIWkvojoavc4xsJzmDhOhHl4DhNDlUNAB4G5pfU5qaxuHUlTgVOAF0doW6VPMzM7jqoEQC8wX9I8SdMoTup219TpBlan5RXA/RERqXxV+pTQPGA+8PWKfZqZ2XHU8BBQOqa/FtgBTAE2R8QeSRuAvojoBm4GbpfUDxyieEMn1dsKPAEMAf8tIr4PUK/P1k9vwpjUh7ASz2HiOBHm4TlMACr+UDczs9z4UhBmZplyAJiZZcoB0GKSNkt6QdLjNeW/LekfJO2R9Ml2ja+KenOQtFDSTkm7JPVJWtzOMTYiaa6kByQ9kf7NfyeVz5B0n6Sn0s9T2z3W4Ywwh0+l/0uPSbpb0pvbPNRhDTeH0vbflxSSZrZrjFWMNI/J9Nr+dyLCjxY+gPOAdwKPl8p+GfgqMD2tv6Xd4xzFHL4CXJSW3wc82O5xNpjDLOCdaflHgW9SXHbkk8C6VL4OuKHdYx3FHC4EpqbyGybjHNL6XIoPgnwLmNnusY7ydzGpXtu1D+8BtFhE/C3FJ6HKfgu4PiIOpzovjPvAmjDMHAI4OS2fAjw3roNqUkQ8HxGPpOVXgCcpvoVevmzJrcDytgywguHmEBFfiYihVG0nxfdoJqQRfg9QXDfsDyn+b01oI8xjUr22azkAxsfbgHPTlVK/JmlRuwc0Cr8LfErSAeDPgI+2dzjVpavTngU8BPxYRDyfNv0j8GPtGlczauZQ9hvAveM+oFEoz0HSMuBgRDza3lE1r+Z3Malf2xP+UhAniKnADOAcYBGwVdLpkfYZJ4nfAn4vIr4k6RKK7368p81jakjSm4AvAb8bEd8trlFYiIiQNOF/B7VzKJV/jOL7NXe0a2xVledAMeY/ojiUNanU+f80qV/b3gMYHwPA30Th68BRigtJTSargb9Jy1+kuKLrhCbpDRQv1jsi4tjY/0nSrLR9FjChd9mHmQOSrgTeD1w+0d9s6szhrcA84FFJ+ykOYT0i6T+0b5SNDfO7mNSvbQfA+PgyxckiJL0NmMbkuxLic8AvpuXzgafaOJaG0uXIbwaejIj/UdpUvmzJauCe8R5bVcPNQdISimPnSyPi1XaNr4p6c4iI3RHxlojojIhOijfRd0bEP7ZxqCMa4f/Tl5nEr21/E7jFJH0B+CWKvwL+CbgWuB3YDCwEjgDXRMT9bRpiQ8PMYS/wFxSHs14DPhwRD7drjI1I+gXg74DdFH+VQXHY4SFgK/ATFJ8+uSQiak94TwgjzOHTwHSKCy4C7IyID43/CBsbbg4R0VOqsx/oiogJ+8Y5wu/iq0yi13YtB4CZWaZ8CMjMLFMOADOzTDkAzMwy5QAwM8uUA8DMLFMOAMuGpD+RdM0o2y6U9L5m+1LhfkknS+qsvUpsk2P4M0nnj7a9WS0HgFk1Cymugtqs9wGPli/hMAb/k+IKpmYt4QCwE5qkj0n6pqT/C/xUqfytkv6PpIcl/Z2kt6fyWyR9Nt3z4JuS3i9pGrABuDTdD+HS1M0CSQ9K2ifpvw8zhMup821jSadL+oakRZKulPTldH+C/ZLWSro6bd8paQZARHwLOG2iXzLBJg8HgJ2wJJ0NrOKHf72Xr9S4CfjtiDgbuAb4y9K2ToprHV0MfJbidbIeuCsiFkbEXane24H3prrXpmvF1Pp54HXfmJb0UxTXlLkyInpT8ZnAr6UxfgJ4NSLOAv4euKLU/JHUp9mY+WqgdiI7F7j72PVyJHWnn28Cfg74YunqoNNL7bZGxFHgKUn7KN7o69mergN/WNILFJeWHqipMyNdP/6YDoo9gl+LiCdK5Q+keq9Iehn4X6l8N/CzpXovAD/eYN5mlTgALEcnAd+JiIXDbK+9Pspw10s5XFr+PvVfT0OSTkqBAvAy8CzwC0A5AMp9HS2tH63p90eA7w0zHrOm+BCQncj+Flgu6Y2SfhT4FYB0QvYZSSvhB5/U+Y+ldislnSTprcDpFBfCe4XiVoDN2pv6OOYI8KvAFZL+0yj6exsw6k8SmZU5AOyElW7hdxfwKMVds3pLmy8HrpL0KLCH4laRxzwLfD21+VBEvAY8QHHSt3wSuIrtFFdWLY/rXyiu5f97kpZW7SidYzgD6Gvi+c2G5auBmpVIugX43xGxrUX9zQJui4gLWtDXr1JcN//jYx+ZmfcAzI6rdP/hz0k6uQXdTQX+vAX9mAHeAzAzy5b3AMzMMuUAMDPLlAPAzCxTDgAzs0w5AMzMMvVvzjqsXfaMWsoAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "d = dd.sample_depths(1000, mean_depth = 20, isaftershock=True,  doplot=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "77153340",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "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.8.12"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}