DFN_mesh_node.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import pandas as pd"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [],
   "source": [
    "data=pd.read_csv('./data/n40e_node.csv', float_precision='round_trip')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Id</th>\n",
       "      <th>X</th>\n",
       "      <th>Y</th>\n",
       "      <th>Z</th>\n",
       "      <th>pressure</th>\n",
       "      <th>temperature</th>\n",
       "      <th>depth</th>\n",
       "      <th>disp_k_bc</th>\n",
       "      <th>disp_j_bc</th>\n",
       "      <th>disp_i_bc</th>\n",
       "      <th>vertical</th>\n",
       "      <th>h_max</th>\n",
       "      <th>h_min</th>\n",
       "      <th>depth_below_wt</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>count</th>\n",
       "      <td>127449.000000</td>\n",
       "      <td>127449.000000</td>\n",
       "      <td>127449.000000</td>\n",
       "      <td>127449.000000</td>\n",
       "      <td>1.274490e+05</td>\n",
       "      <td>127449.000000</td>\n",
       "      <td>127449.000000</td>\n",
       "      <td>1.274490e+05</td>\n",
       "      <td>1.274490e+05</td>\n",
       "      <td>1.274490e+05</td>\n",
       "      <td>1.274490e+05</td>\n",
       "      <td>1.274490e+05</td>\n",
       "      <td>1.274490e+05</td>\n",
       "      <td>127449.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>mean</th>\n",
       "      <td>63724.000000</td>\n",
       "      <td>2000.000000</td>\n",
       "      <td>2000.000000</td>\n",
       "      <td>-376.000000</td>\n",
       "      <td>1.873044e+07</td>\n",
       "      <td>450.416959</td>\n",
       "      <td>2054.224202</td>\n",
       "      <td>-4.428907e+07</td>\n",
       "      <td>-3.492181e+07</td>\n",
       "      <td>-2.834829e+07</td>\n",
       "      <td>-2.555864e+07</td>\n",
       "      <td>-1.619137e+07</td>\n",
       "      <td>-9.617856e+06</td>\n",
       "      <td>1949.007647</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>std</th>\n",
       "      <td>36791.501566</td>\n",
       "      <td>1177.572735</td>\n",
       "      <td>1177.572735</td>\n",
       "      <td>1131.375288</td>\n",
       "      <td>1.087857e+07</td>\n",
       "      <td>62.995666</td>\n",
       "      <td>1132.635458</td>\n",
       "      <td>2.441962e+07</td>\n",
       "      <td>1.925480e+07</td>\n",
       "      <td>1.563037e+07</td>\n",
       "      <td>1.354379e+07</td>\n",
       "      <td>8.379729e+06</td>\n",
       "      <td>4.756801e+06</td>\n",
       "      <td>1132.375333</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>min</th>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>-2296.000000</td>\n",
       "      <td>1.013250e+05</td>\n",
       "      <td>323.150000</td>\n",
       "      <td>35.271242</td>\n",
       "      <td>-8.867295e+07</td>\n",
       "      <td>-6.991838e+07</td>\n",
       "      <td>-5.675727e+07</td>\n",
       "      <td>-4.989358e+07</td>\n",
       "      <td>-3.127636e+07</td>\n",
       "      <td>-1.836012e+07</td>\n",
       "      <td>-30.064704</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25%</th>\n",
       "      <td>31862.000000</td>\n",
       "      <td>960.000000</td>\n",
       "      <td>960.000000</td>\n",
       "      <td>-1336.000000</td>\n",
       "      <td>9.273462e+06</td>\n",
       "      <td>397.980000</td>\n",
       "      <td>1074.295128</td>\n",
       "      <td>-6.541298e+07</td>\n",
       "      <td>-5.157795e+07</td>\n",
       "      <td>-4.186916e+07</td>\n",
       "      <td>-3.726533e+07</td>\n",
       "      <td>-2.343566e+07</td>\n",
       "      <td>-1.371555e+07</td>\n",
       "      <td>965.041706</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>50%</th>\n",
       "      <td>63724.000000</td>\n",
       "      <td>2000.000000</td>\n",
       "      <td>2000.000000</td>\n",
       "      <td>-376.000000</td>\n",
       "      <td>1.869320e+07</td>\n",
       "      <td>457.780000</td>\n",
       "      <td>2056.037810</td>\n",
       "      <td>-4.432818e+07</td>\n",
       "      <td>-3.495264e+07</td>\n",
       "      <td>-2.837332e+07</td>\n",
       "      <td>-2.555744e+07</td>\n",
       "      <td>-1.618654e+07</td>\n",
       "      <td>-9.612952e+06</td>\n",
       "      <td>1945.305898</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>75%</th>\n",
       "      <td>95586.000000</td>\n",
       "      <td>3040.000000</td>\n",
       "      <td>3040.000000</td>\n",
       "      <td>584.000000</td>\n",
       "      <td>2.820853e+07</td>\n",
       "      <td>507.140000</td>\n",
       "      <td>3033.997263</td>\n",
       "      <td>-2.316180e+07</td>\n",
       "      <td>-1.826302e+07</td>\n",
       "      <td>-1.482527e+07</td>\n",
       "      <td>-1.385670e+07</td>\n",
       "      <td>-8.953497e+06</td>\n",
       "      <td>-5.519811e+06</td>\n",
       "      <td>2935.517430</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>max</th>\n",
       "      <td>127448.000000</td>\n",
       "      <td>4000.000000</td>\n",
       "      <td>4000.000000</td>\n",
       "      <td>1544.000000</td>\n",
       "      <td>3.878688e+07</td>\n",
       "      <td>535.880000</td>\n",
       "      <td>4112.845681</td>\n",
       "      <td>-7.604480e+05</td>\n",
       "      <td>-5.996111e+05</td>\n",
       "      <td>-4.867431e+05</td>\n",
       "      <td>-6.591230e+05</td>\n",
       "      <td>-4.982861e+05</td>\n",
       "      <td>-3.854181e+05</td>\n",
       "      <td>4036.351535</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                  Id              X              Y              Z  \\\n",
       "count  127449.000000  127449.000000  127449.000000  127449.000000   \n",
       "mean    63724.000000    2000.000000    2000.000000    -376.000000   \n",
       "std     36791.501566    1177.572735    1177.572735    1131.375288   \n",
       "min         0.000000       0.000000       0.000000   -2296.000000   \n",
       "25%     31862.000000     960.000000     960.000000   -1336.000000   \n",
       "50%     63724.000000    2000.000000    2000.000000    -376.000000   \n",
       "75%     95586.000000    3040.000000    3040.000000     584.000000   \n",
       "max    127448.000000    4000.000000    4000.000000    1544.000000   \n",
       "\n",
       "           pressure    temperature          depth     disp_k_bc     disp_j_bc  \\\n",
       "count  1.274490e+05  127449.000000  127449.000000  1.274490e+05  1.274490e+05   \n",
       "mean   1.873044e+07     450.416959    2054.224202 -4.428907e+07 -3.492181e+07   \n",
       "std    1.087857e+07      62.995666    1132.635458  2.441962e+07  1.925480e+07   \n",
       "min    1.013250e+05     323.150000      35.271242 -8.867295e+07 -6.991838e+07   \n",
       "25%    9.273462e+06     397.980000    1074.295128 -6.541298e+07 -5.157795e+07   \n",
       "50%    1.869320e+07     457.780000    2056.037810 -4.432818e+07 -3.495264e+07   \n",
       "75%    2.820853e+07     507.140000    3033.997263 -2.316180e+07 -1.826302e+07   \n",
       "max    3.878688e+07     535.880000    4112.845681 -7.604480e+05 -5.996111e+05   \n",
       "\n",
       "          disp_i_bc      vertical         h_max         h_min  depth_below_wt  \n",
       "count  1.274490e+05  1.274490e+05  1.274490e+05  1.274490e+05   127449.000000  \n",
       "mean  -2.834829e+07 -2.555864e+07 -1.619137e+07 -9.617856e+06     1949.007647  \n",
       "std    1.563037e+07  1.354379e+07  8.379729e+06  4.756801e+06     1132.375333  \n",
       "min   -5.675727e+07 -4.989358e+07 -3.127636e+07 -1.836012e+07      -30.064704  \n",
       "25%   -4.186916e+07 -3.726533e+07 -2.343566e+07 -1.371555e+07      965.041706  \n",
       "50%   -2.837332e+07 -2.555744e+07 -1.618654e+07 -9.612952e+06     1945.305898  \n",
       "75%   -1.482527e+07 -1.385670e+07 -8.953497e+06 -5.519811e+06     2935.517430  \n",
       "max   -4.867431e+05 -6.591230e+05 -4.982861e+05 -3.854181e+05     4036.351535  "
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "data.describe()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "All is good.  BUT...because the file is so big, it can't be edited in Excel, so need to add some additional columns to the file now, and then later add the 10 row header that is usually written by Leapfrog"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Add the disp_k_bc\tdisp_j_bc\tdisp_i_bc\tvertical\th_max\th_min columns"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "data.insert(loc=8, column='disp_k_bc', value=1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "data.insert(loc=9, column='disp_j_bc', value=1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "data.insert(loc=10, column='disp_i_bc', value=1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [],
   "source": [
    "data.insert(loc=11, column='vertical', value=1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [],
   "source": [
    "data.insert(loc=12, column='h_max', value=1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [],
   "source": [
    "data.insert(loc=13, column='h_min', value=1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [],
   "source": [
    "data.rename(columns={'temperature_C': 'temperature'}, inplace=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [],
   "source": [
    "data.rename(columns={'pressure_MPa': 'pressure'}, inplace=True)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "now check the file"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Id</th>\n",
       "      <th>X</th>\n",
       "      <th>Y</th>\n",
       "      <th>Z</th>\n",
       "      <th>dX</th>\n",
       "      <th>dY</th>\n",
       "      <th>dZ</th>\n",
       "      <th>pressure</th>\n",
       "      <th>disp_k_bc</th>\n",
       "      <th>disp_j_bc</th>\n",
       "      <th>disp_i_bc</th>\n",
       "      <th>vertical</th>\n",
       "      <th>h_max</th>\n",
       "      <th>h_min</th>\n",
       "      <th>sigma_h_max_MPa</th>\n",
       "      <th>sigma_h_min_MPa</th>\n",
       "      <th>Sigma_V_MPa</th>\n",
       "      <th>temperature</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>0</td>\n",
       "      <td>50</td>\n",
       "      <td>420</td>\n",
       "      <td>-1475</td>\n",
       "      <td>10</td>\n",
       "      <td>10</td>\n",
       "      <td>10</td>\n",
       "      <td>27.667349</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>54.383797</td>\n",
       "      <td>47.417642</td>\n",
       "      <td>78.659853</td>\n",
       "      <td>226.074346</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>1</td>\n",
       "      <td>60</td>\n",
       "      <td>420</td>\n",
       "      <td>-1475</td>\n",
       "      <td>10</td>\n",
       "      <td>10</td>\n",
       "      <td>10</td>\n",
       "      <td>27.666658</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>54.381446</td>\n",
       "      <td>47.412370</td>\n",
       "      <td>78.656734</td>\n",
       "      <td>226.210827</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>2</td>\n",
       "      <td>70</td>\n",
       "      <td>420</td>\n",
       "      <td>-1475</td>\n",
       "      <td>10</td>\n",
       "      <td>10</td>\n",
       "      <td>10</td>\n",
       "      <td>27.666128</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>54.381283</td>\n",
       "      <td>47.408446</td>\n",
       "      <td>78.656918</td>\n",
       "      <td>226.350794</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>3</td>\n",
       "      <td>80</td>\n",
       "      <td>420</td>\n",
       "      <td>-1475</td>\n",
       "      <td>10</td>\n",
       "      <td>10</td>\n",
       "      <td>10</td>\n",
       "      <td>27.665727</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>54.383207</td>\n",
       "      <td>47.405856</td>\n",
       "      <td>78.660337</td>\n",
       "      <td>226.494008</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>4</td>\n",
       "      <td>90</td>\n",
       "      <td>420</td>\n",
       "      <td>-1475</td>\n",
       "      <td>10</td>\n",
       "      <td>10</td>\n",
       "      <td>10</td>\n",
       "      <td>27.665387</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>54.386726</td>\n",
       "      <td>47.404308</td>\n",
       "      <td>78.666232</td>\n",
       "      <td>226.640055</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   Id   X    Y     Z  dX  dY  dZ   pressure  disp_k_bc  disp_j_bc  disp_i_bc  \\\n",
       "0   0  50  420 -1475  10  10  10  27.667349          1          1          1   \n",
       "1   1  60  420 -1475  10  10  10  27.666658          1          1          1   \n",
       "2   2  70  420 -1475  10  10  10  27.666128          1          1          1   \n",
       "3   3  80  420 -1475  10  10  10  27.665727          1          1          1   \n",
       "4   4  90  420 -1475  10  10  10  27.665387          1          1          1   \n",
       "\n",
       "   vertical  h_max  h_min  sigma_h_max_MPa  sigma_h_min_MPa  Sigma_V_MPa  \\\n",
       "0         1      1      1        54.383797        47.417642    78.659853   \n",
       "1         1      1      1        54.381446        47.412370    78.656734   \n",
       "2         1      1      1        54.381283        47.408446    78.656918   \n",
       "3         1      1      1        54.383207        47.405856    78.660337   \n",
       "4         1      1      1        54.386726        47.404308    78.666232   \n",
       "\n",
       "   temperature  \n",
       "0   226.074346  \n",
       "1   226.210827  \n",
       "2   226.350794  \n",
       "3   226.494008  \n",
       "4   226.640055  "
      ]
     },
     "execution_count": 13,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "data.head()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "now sort so values are in correct order for em2ex"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [],
   "source": [
    "a = data.sort_values([\"Z\",\"Y\",\"X\"], ascending = [True, True, True])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Id</th>\n",
       "      <th>X</th>\n",
       "      <th>Y</th>\n",
       "      <th>Z</th>\n",
       "      <th>pressure</th>\n",
       "      <th>temperature_K</th>\n",
       "      <th>depth</th>\n",
       "      <th>disp_k_bc</th>\n",
       "      <th>disp_j_bc</th>\n",
       "      <th>disp_i_bc</th>\n",
       "      <th>vertical</th>\n",
       "      <th>h_max</th>\n",
       "      <th>h_min</th>\n",
       "      <th>depth_below_wt</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>-2296</td>\n",
       "      <td>36616764.63</td>\n",
       "      <td>508.64</td>\n",
       "      <td>3888.058474</td>\n",
       "      <td>-8.382654e+07</td>\n",
       "      <td>-6.765222e+07</td>\n",
       "      <td>-5.676565e+07</td>\n",
       "      <td>-4.720978e+07</td>\n",
       "      <td>-3.103545e+07</td>\n",
       "      <td>-2.014889e+07</td>\n",
       "      <td>3810.519152</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>2499</td>\n",
       "      <td>80</td>\n",
       "      <td>0</td>\n",
       "      <td>-2296</td>\n",
       "      <td>36622110.35</td>\n",
       "      <td>509.96</td>\n",
       "      <td>3890.393152</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>3811.075453</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>4998</td>\n",
       "      <td>160</td>\n",
       "      <td>0</td>\n",
       "      <td>-2296</td>\n",
       "      <td>36627815.37</td>\n",
       "      <td>511.28</td>\n",
       "      <td>3892.429949</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>3811.669145</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>7497</td>\n",
       "      <td>240</td>\n",
       "      <td>0</td>\n",
       "      <td>-2296</td>\n",
       "      <td>36633797.36</td>\n",
       "      <td>512.59</td>\n",
       "      <td>3894.237728</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>3812.291661</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19</th>\n",
       "      <td>9996</td>\n",
       "      <td>320</td>\n",
       "      <td>0</td>\n",
       "      <td>-2296</td>\n",
       "      <td>36640112.56</td>\n",
       "      <td>513.88</td>\n",
       "      <td>3896.918477</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>3812.948851</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "      Id    X  Y     Z     pressure  temperature_K        depth     disp_k_bc  \\\n",
       "0      0    0  0 -2296  36616764.63         508.64  3888.058474 -8.382654e+07   \n",
       "2   2499   80  0 -2296  36622110.35         509.96  3890.393152  0.000000e+00   \n",
       "6   4998  160  0 -2296  36627815.37         511.28  3892.429949  0.000000e+00   \n",
       "12  7497  240  0 -2296  36633797.36         512.59  3894.237728  0.000000e+00   \n",
       "19  9996  320  0 -2296  36640112.56         513.88  3896.918477  0.000000e+00   \n",
       "\n",
       "       disp_j_bc     disp_i_bc      vertical         h_max         h_min  \\\n",
       "0  -6.765222e+07 -5.676565e+07 -4.720978e+07 -3.103545e+07 -2.014889e+07   \n",
       "2   0.000000e+00  0.000000e+00  0.000000e+00  0.000000e+00  0.000000e+00   \n",
       "6   0.000000e+00  0.000000e+00  0.000000e+00  0.000000e+00  0.000000e+00   \n",
       "12  0.000000e+00  0.000000e+00  0.000000e+00  0.000000e+00  0.000000e+00   \n",
       "19  0.000000e+00  0.000000e+00  0.000000e+00  0.000000e+00  0.000000e+00   \n",
       "\n",
       "    depth_below_wt  \n",
       "0      3810.519152  \n",
       "2      3811.075453  \n",
       "6      3811.669145  \n",
       "12     3812.291661  \n",
       "19     3812.948851  "
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "a.head()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "now we have all the needed columns in the dataframe.  Now can populate the new columns with data, and also fix the pressure and temperature"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [],
   "source": [
    "a['temperature'] =+ a['temperature'] + 273.15"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [],
   "source": [
    "a['pressure'] = a['pressure']* 1e6"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [],
   "source": [
    "a['disp_k_bc'] = (a['Sigma_V_MPa'] * -1e6)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [],
   "source": [
    "a['disp_j_bc'] = (a['sigma_h_max_MPa'] * -1e6)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [],
   "source": [
    "a['disp_i_bc'] = (a['sigma_h_min_MPa'] * -1e6)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "to get IC for effective stress, must convert total to negtive and Pa, the add pressure.  Note that alpha (0.6) was already factored into the calculation of total stress, so we don't need to worry about it converting back into MOOSE conventions"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {},
   "outputs": [],
   "source": [
    "a['vertical'] = (a['Sigma_V_MPa'] * -1e6) + (a['pressure'])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {},
   "outputs": [],
   "source": [
    "a['h_max'] = (a['sigma_h_max_MPa'] * -1e6) + (a['pressure'])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {},
   "outputs": [],
   "source": [
    "a['h_min'] = (a['sigma_h_min_MPa'] * -1e6) + (a['pressure'])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Id</th>\n",
       "      <th>X</th>\n",
       "      <th>Y</th>\n",
       "      <th>Z</th>\n",
       "      <th>dX</th>\n",
       "      <th>dY</th>\n",
       "      <th>dZ</th>\n",
       "      <th>pressure</th>\n",
       "      <th>disp_k_bc</th>\n",
       "      <th>disp_j_bc</th>\n",
       "      <th>disp_i_bc</th>\n",
       "      <th>vertical</th>\n",
       "      <th>h_max</th>\n",
       "      <th>h_min</th>\n",
       "      <th>sigma_h_max_MPa</th>\n",
       "      <th>sigma_h_min_MPa</th>\n",
       "      <th>Sigma_V_MPa</th>\n",
       "      <th>temperature</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>0</td>\n",
       "      <td>50</td>\n",
       "      <td>420</td>\n",
       "      <td>-1475</td>\n",
       "      <td>10</td>\n",
       "      <td>10</td>\n",
       "      <td>10</td>\n",
       "      <td>27667348.9</td>\n",
       "      <td>-78659853.2</td>\n",
       "      <td>-54383797.3</td>\n",
       "      <td>-47417642.2</td>\n",
       "      <td>-50992504.3</td>\n",
       "      <td>-26716448.4</td>\n",
       "      <td>-19750293.3</td>\n",
       "      <td>54.383797</td>\n",
       "      <td>47.417642</td>\n",
       "      <td>78.659853</td>\n",
       "      <td>499.224346</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>1</td>\n",
       "      <td>60</td>\n",
       "      <td>420</td>\n",
       "      <td>-1475</td>\n",
       "      <td>10</td>\n",
       "      <td>10</td>\n",
       "      <td>10</td>\n",
       "      <td>27666657.9</td>\n",
       "      <td>-78656733.8</td>\n",
       "      <td>-54381446.3</td>\n",
       "      <td>-47412369.6</td>\n",
       "      <td>-50990075.9</td>\n",
       "      <td>-26714788.4</td>\n",
       "      <td>-19745711.7</td>\n",
       "      <td>54.381446</td>\n",
       "      <td>47.412370</td>\n",
       "      <td>78.656734</td>\n",
       "      <td>499.360827</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>2</td>\n",
       "      <td>70</td>\n",
       "      <td>420</td>\n",
       "      <td>-1475</td>\n",
       "      <td>10</td>\n",
       "      <td>10</td>\n",
       "      <td>10</td>\n",
       "      <td>27666128.4</td>\n",
       "      <td>-78656918.4</td>\n",
       "      <td>-54381283.4</td>\n",
       "      <td>-47408446.2</td>\n",
       "      <td>-50990790.0</td>\n",
       "      <td>-26715155.0</td>\n",
       "      <td>-19742317.8</td>\n",
       "      <td>54.381283</td>\n",
       "      <td>47.408446</td>\n",
       "      <td>78.656918</td>\n",
       "      <td>499.500794</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>3</td>\n",
       "      <td>80</td>\n",
       "      <td>420</td>\n",
       "      <td>-1475</td>\n",
       "      <td>10</td>\n",
       "      <td>10</td>\n",
       "      <td>10</td>\n",
       "      <td>27665726.5</td>\n",
       "      <td>-78660336.9</td>\n",
       "      <td>-54383207.2</td>\n",
       "      <td>-47405856.2</td>\n",
       "      <td>-50994610.4</td>\n",
       "      <td>-26717480.7</td>\n",
       "      <td>-19740129.7</td>\n",
       "      <td>54.383207</td>\n",
       "      <td>47.405856</td>\n",
       "      <td>78.660337</td>\n",
       "      <td>499.644008</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>4</td>\n",
       "      <td>90</td>\n",
       "      <td>420</td>\n",
       "      <td>-1475</td>\n",
       "      <td>10</td>\n",
       "      <td>10</td>\n",
       "      <td>10</td>\n",
       "      <td>27665386.7</td>\n",
       "      <td>-78666231.6</td>\n",
       "      <td>-54386726.3</td>\n",
       "      <td>-47404308.1</td>\n",
       "      <td>-51000844.9</td>\n",
       "      <td>-26721339.6</td>\n",
       "      <td>-19738921.4</td>\n",
       "      <td>54.386726</td>\n",
       "      <td>47.404308</td>\n",
       "      <td>78.666232</td>\n",
       "      <td>499.790055</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   Id   X    Y     Z  dX  dY  dZ    pressure   disp_k_bc   disp_j_bc  \\\n",
       "0   0  50  420 -1475  10  10  10  27667348.9 -78659853.2 -54383797.3   \n",
       "1   1  60  420 -1475  10  10  10  27666657.9 -78656733.8 -54381446.3   \n",
       "2   2  70  420 -1475  10  10  10  27666128.4 -78656918.4 -54381283.4   \n",
       "3   3  80  420 -1475  10  10  10  27665726.5 -78660336.9 -54383207.2   \n",
       "4   4  90  420 -1475  10  10  10  27665386.7 -78666231.6 -54386726.3   \n",
       "\n",
       "    disp_i_bc    vertical       h_max       h_min  sigma_h_max_MPa  \\\n",
       "0 -47417642.2 -50992504.3 -26716448.4 -19750293.3        54.383797   \n",
       "1 -47412369.6 -50990075.9 -26714788.4 -19745711.7        54.381446   \n",
       "2 -47408446.2 -50990790.0 -26715155.0 -19742317.8        54.381283   \n",
       "3 -47405856.2 -50994610.4 -26717480.7 -19740129.7        54.383207   \n",
       "4 -47404308.1 -51000844.9 -26721339.6 -19738921.4        54.386726   \n",
       "\n",
       "   sigma_h_min_MPa  Sigma_V_MPa  temperature  \n",
       "0        47.417642    78.659853   499.224346  \n",
       "1        47.412370    78.656734   499.360827  \n",
       "2        47.408446    78.656918   499.500794  \n",
       "3        47.405856    78.660337   499.644008  \n",
       "4        47.404308    78.666232   499.790055  "
      ]
     },
     "execution_count": 24,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "a.head()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "now we can delete some columns"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {},
   "outputs": [],
   "source": [
    "del a['sigma_h_max_MPa']"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {},
   "outputs": [],
   "source": [
    "del a['sigma_h_min_MPa']"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "metadata": {},
   "outputs": [],
   "source": [
    "del a['Sigma_V_MPa']"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "and now the dataframe should be ready for em2ex"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Id</th>\n",
       "      <th>X</th>\n",
       "      <th>Y</th>\n",
       "      <th>Z</th>\n",
       "      <th>dX</th>\n",
       "      <th>dY</th>\n",
       "      <th>dZ</th>\n",
       "      <th>pressure</th>\n",
       "      <th>disp_k_bc</th>\n",
       "      <th>disp_j_bc</th>\n",
       "      <th>disp_i_bc</th>\n",
       "      <th>vertical</th>\n",
       "      <th>h_max</th>\n",
       "      <th>h_min</th>\n",
       "      <th>temperature</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>0</td>\n",
       "      <td>50</td>\n",
       "      <td>420</td>\n",
       "      <td>-1475</td>\n",
       "      <td>10</td>\n",
       "      <td>10</td>\n",
       "      <td>10</td>\n",
       "      <td>27667348.9</td>\n",
       "      <td>-78659853.2</td>\n",
       "      <td>-54383797.3</td>\n",
       "      <td>-47417642.2</td>\n",
       "      <td>-50992504.3</td>\n",
       "      <td>-26716448.4</td>\n",
       "      <td>-19750293.3</td>\n",
       "      <td>499.224346</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>1</td>\n",
       "      <td>60</td>\n",
       "      <td>420</td>\n",
       "      <td>-1475</td>\n",
       "      <td>10</td>\n",
       "      <td>10</td>\n",
       "      <td>10</td>\n",
       "      <td>27666657.9</td>\n",
       "      <td>-78656733.8</td>\n",
       "      <td>-54381446.3</td>\n",
       "      <td>-47412369.6</td>\n",
       "      <td>-50990075.9</td>\n",
       "      <td>-26714788.4</td>\n",
       "      <td>-19745711.7</td>\n",
       "      <td>499.360827</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>2</td>\n",
       "      <td>70</td>\n",
       "      <td>420</td>\n",
       "      <td>-1475</td>\n",
       "      <td>10</td>\n",
       "      <td>10</td>\n",
       "      <td>10</td>\n",
       "      <td>27666128.4</td>\n",
       "      <td>-78656918.4</td>\n",
       "      <td>-54381283.4</td>\n",
       "      <td>-47408446.2</td>\n",
       "      <td>-50990790.0</td>\n",
       "      <td>-26715155.0</td>\n",
       "      <td>-19742317.8</td>\n",
       "      <td>499.500794</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>3</td>\n",
       "      <td>80</td>\n",
       "      <td>420</td>\n",
       "      <td>-1475</td>\n",
       "      <td>10</td>\n",
       "      <td>10</td>\n",
       "      <td>10</td>\n",
       "      <td>27665726.5</td>\n",
       "      <td>-78660336.9</td>\n",
       "      <td>-54383207.2</td>\n",
       "      <td>-47405856.2</td>\n",
       "      <td>-50994610.4</td>\n",
       "      <td>-26717480.7</td>\n",
       "      <td>-19740129.7</td>\n",
       "      <td>499.644008</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>4</td>\n",
       "      <td>90</td>\n",
       "      <td>420</td>\n",
       "      <td>-1475</td>\n",
       "      <td>10</td>\n",
       "      <td>10</td>\n",
       "      <td>10</td>\n",
       "      <td>27665386.7</td>\n",
       "      <td>-78666231.6</td>\n",
       "      <td>-54386726.3</td>\n",
       "      <td>-47404308.1</td>\n",
       "      <td>-51000844.9</td>\n",
       "      <td>-26721339.6</td>\n",
       "      <td>-19738921.4</td>\n",
       "      <td>499.790055</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   Id   X    Y     Z  dX  dY  dZ    pressure   disp_k_bc   disp_j_bc  \\\n",
       "0   0  50  420 -1475  10  10  10  27667348.9 -78659853.2 -54383797.3   \n",
       "1   1  60  420 -1475  10  10  10  27666657.9 -78656733.8 -54381446.3   \n",
       "2   2  70  420 -1475  10  10  10  27666128.4 -78656918.4 -54381283.4   \n",
       "3   3  80  420 -1475  10  10  10  27665726.5 -78660336.9 -54383207.2   \n",
       "4   4  90  420 -1475  10  10  10  27665386.7 -78666231.6 -54386726.3   \n",
       "\n",
       "    disp_i_bc    vertical       h_max       h_min  temperature  \n",
       "0 -47417642.2 -50992504.3 -26716448.4 -19750293.3   499.224346  \n",
       "1 -47412369.6 -50990075.9 -26714788.4 -19745711.7   499.360827  \n",
       "2 -47408446.2 -50990790.0 -26715155.0 -19742317.8   499.500794  \n",
       "3 -47405856.2 -50994610.4 -26717480.7 -19740129.7   499.644008  \n",
       "4 -47404308.1 -51000844.9 -26721339.6 -19738921.4   499.790055  "
      ]
     },
     "execution_count": 28,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "a.head()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "now can write the data to a CSV file."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [],
   "source": [
    "a.to_csv(\"processed_n40e_node.csv\", encoding='utf-8', index=False)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "metadata": {},
   "outputs": [],
   "source": [
    "a.describe().to_csv(\"processed_local_large_well_3_1b_node_Summary.csv\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "now need to read into leapfrog and verify all the mpaaings"
   ]
  }
 ],
 "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.8.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}