new_test_dataset.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "288d0819-0f40-4b95-bcbb-65af51a05adc",
   "metadata": {},
   "outputs": [],
   "source": [
    "import pandas as pd\n",
    "import numpy as np"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "eb7b2411",
   "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>Unnamed: 0</th>\n",
       "      <th>Unnamed: 1</th>\n",
       "      <th>Unnamed: 2</th>\n",
       "      <th>Unnamed: 3</th>\n",
       "      <th>Unnamed: 4</th>\n",
       "      <th>Unnamed: 5</th>\n",
       "      <th>Unnamed: 6</th>\n",
       "      <th>Unnamed: 7</th>\n",
       "      <th>Unnamed: 8</th>\n",
       "      <th>Unnamed: 9</th>\n",
       "      <th>Unnamed: 10</th>\n",
       "      <th>Unnamed: 11</th>\n",
       "      <th>Unnamed: 12</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>LNG FLOW</td>\n",
       "      <td>LNG RUNDOWN TEMPERATURE</td>\n",
       "      <td>HMR FLOW</td>\n",
       "      <td>LMR/HMR RATIO</td>\n",
       "      <td>PERCENTAGE NITROGEN IN MIXED REFERIGERANT</td>\n",
       "      <td>PERCENTAGE METHANE IN MIXED REFERIGERANT</td>\n",
       "      <td>PERCENTAGE ETHANE IN MIXED REFERIGERANT</td>\n",
       "      <td>PERCENTAGE PROPANE IN MIXED REFERIGERANT</td>\n",
       "      <td>PERCENTAGE BUTHANE IN MIXED REFERIGERANT</td>\n",
       "      <td>AMBIENT TEMPERATURE</td>\n",
       "      <td>LOW PRESSSURE PROPANE MIXED REFRIGERANT KETTLE...</td>\n",
       "      <td>PR PEAK MARGIN</td>\n",
       "      <td>MR PEAK MARGIN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>9469.942383</td>\n",
       "      <td>-160.132996</td>\n",
       "      <td>16287.15625</td>\n",
       "      <td>0.336019</td>\n",
       "      <td>9.436035</td>\n",
       "      <td>38.877541</td>\n",
       "      <td>45.310646</td>\n",
       "      <td>6.139355</td>\n",
       "      <td>6.139355</td>\n",
       "      <td>26.593405</td>\n",
       "      <td>57.273438</td>\n",
       "      <td>28.285618</td>\n",
       "      <td>7.644921</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>9507.515625</td>\n",
       "      <td>-160.059555</td>\n",
       "      <td>16293.646484</td>\n",
       "      <td>0.333414</td>\n",
       "      <td>9.390137</td>\n",
       "      <td>38.877541</td>\n",
       "      <td>45.310646</td>\n",
       "      <td>6.160547</td>\n",
       "      <td>6.160547</td>\n",
       "      <td>26.280834</td>\n",
       "      <td>57.09375</td>\n",
       "      <td>28.944149</td>\n",
       "      <td>7.525976</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>9519.392578</td>\n",
       "      <td>-159.999573</td>\n",
       "      <td>16333.569336</td>\n",
       "      <td>0.33673</td>\n",
       "      <td>9.42041</td>\n",
       "      <td>38.982128</td>\n",
       "      <td>45.207718</td>\n",
       "      <td>6.139355</td>\n",
       "      <td>6.139355</td>\n",
       "      <td>25.690422</td>\n",
       "      <td>56.572266</td>\n",
       "      <td>29.707829</td>\n",
       "      <td>6.763085</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "    Unnamed: 0               Unnamed: 1    Unnamed: 2     Unnamed: 3  \\\n",
       "0          NaN                      NaN           NaN            NaN   \n",
       "1     LNG FLOW  LNG RUNDOWN TEMPERATURE      HMR FLOW  LMR/HMR RATIO   \n",
       "2  9469.942383              -160.132996   16287.15625       0.336019   \n",
       "3  9507.515625              -160.059555  16293.646484       0.333414   \n",
       "4  9519.392578              -159.999573  16333.569336        0.33673   \n",
       "\n",
       "                                   Unnamed: 4  \\\n",
       "0                                         NaN   \n",
       "1  PERCENTAGE NITROGEN IN MIXED REFERIGERANT    \n",
       "2                                    9.436035   \n",
       "3                                    9.390137   \n",
       "4                                     9.42041   \n",
       "\n",
       "                                  Unnamed: 5  \\\n",
       "0                                        NaN   \n",
       "1  PERCENTAGE METHANE IN MIXED REFERIGERANT    \n",
       "2                                  38.877541   \n",
       "3                                  38.877541   \n",
       "4                                  38.982128   \n",
       "\n",
       "                                 Unnamed: 6  \\\n",
       "0                                       NaN   \n",
       "1  PERCENTAGE ETHANE IN MIXED REFERIGERANT    \n",
       "2                                 45.310646   \n",
       "3                                 45.310646   \n",
       "4                                 45.207718   \n",
       "\n",
       "                                  Unnamed: 7  \\\n",
       "0                                        NaN   \n",
       "1  PERCENTAGE PROPANE IN MIXED REFERIGERANT    \n",
       "2                                   6.139355   \n",
       "3                                   6.160547   \n",
       "4                                   6.139355   \n",
       "\n",
       "                                  Unnamed: 8            Unnamed: 9  \\\n",
       "0                                        NaN                   NaN   \n",
       "1  PERCENTAGE BUTHANE IN MIXED REFERIGERANT   AMBIENT TEMPERATURE    \n",
       "2                                   6.139355             26.593405   \n",
       "3                                   6.160547             26.280834   \n",
       "4                                   6.139355             25.690422   \n",
       "\n",
       "                                         Unnamed: 10     Unnamed: 11  \\\n",
       "0                                                NaN             NaN   \n",
       "1  LOW PRESSSURE PROPANE MIXED REFRIGERANT KETTLE...  PR PEAK MARGIN   \n",
       "2                                          57.273438       28.285618   \n",
       "3                                           57.09375       28.944149   \n",
       "4                                          56.572266       29.707829   \n",
       "\n",
       "      Unnamed: 12  \n",
       "0             NaN  \n",
       "1  MR PEAK MARGIN  \n",
       "2        7.644921  \n",
       "3        7.525976  \n",
       "4        6.763085  "
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "data = pd.read_excel(\"./T3 LNG Validation DATA.xlsx\")\n",
    "data.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "b25f8ae9",
   "metadata": {},
   "outputs": [],
   "source": [
    "data.drop([0], inplace=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "5d73e2c8",
   "metadata": {},
   "outputs": [],
   "source": [
    "headers = data.iloc[0, :]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "8c79927b",
   "metadata": {},
   "outputs": [],
   "source": [
    "data.drop([1], inplace=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "9fbbbd7f",
   "metadata": {},
   "outputs": [],
   "source": [
    "data.columns = headers\n",
    "data = data.reset_index(drop=True)\n",
    "data = data.drop(columns=['PR PEAK MARGIN', 'MR PEAK MARGIN'])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "71b6da0c",
   "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>1</th>\n",
       "      <th>LNG FLOW</th>\n",
       "      <th>LNG RUNDOWN TEMPERATURE</th>\n",
       "      <th>HMR FLOW</th>\n",
       "      <th>LMR/HMR RATIO</th>\n",
       "      <th>PERCENTAGE NITROGEN IN MIXED REFERIGERANT</th>\n",
       "      <th>PERCENTAGE METHANE IN MIXED REFERIGERANT</th>\n",
       "      <th>PERCENTAGE ETHANE IN MIXED REFERIGERANT</th>\n",
       "      <th>PERCENTAGE PROPANE IN MIXED REFERIGERANT</th>\n",
       "      <th>PERCENTAGE BUTHANE IN MIXED REFERIGERANT</th>\n",
       "      <th>AMBIENT TEMPERATURE</th>\n",
       "      <th>LOW PRESSSURE PROPANE MIXED REFRIGERANT KETTLE LEVEL</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>9469.942383</td>\n",
       "      <td>-160.132996</td>\n",
       "      <td>16287.15625</td>\n",
       "      <td>0.336019</td>\n",
       "      <td>9.436035</td>\n",
       "      <td>38.877541</td>\n",
       "      <td>45.310646</td>\n",
       "      <td>6.139355</td>\n",
       "      <td>6.139355</td>\n",
       "      <td>26.593405</td>\n",
       "      <td>57.273438</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>9507.515625</td>\n",
       "      <td>-160.059555</td>\n",
       "      <td>16293.646484</td>\n",
       "      <td>0.333414</td>\n",
       "      <td>9.390137</td>\n",
       "      <td>38.877541</td>\n",
       "      <td>45.310646</td>\n",
       "      <td>6.160547</td>\n",
       "      <td>6.160547</td>\n",
       "      <td>26.280834</td>\n",
       "      <td>57.09375</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>9519.392578</td>\n",
       "      <td>-159.999573</td>\n",
       "      <td>16333.569336</td>\n",
       "      <td>0.33673</td>\n",
       "      <td>9.42041</td>\n",
       "      <td>38.982128</td>\n",
       "      <td>45.207718</td>\n",
       "      <td>6.139355</td>\n",
       "      <td>6.139355</td>\n",
       "      <td>25.690422</td>\n",
       "      <td>56.572266</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>9497.821289</td>\n",
       "      <td>-160.071472</td>\n",
       "      <td>16283.046875</td>\n",
       "      <td>0.338297</td>\n",
       "      <td>9.451172</td>\n",
       "      <td>39.033596</td>\n",
       "      <td>45.129688</td>\n",
       "      <td>6.139355</td>\n",
       "      <td>6.139355</td>\n",
       "      <td>25.238932</td>\n",
       "      <td>57.736328</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>9516.681641</td>\n",
       "      <td>-160.009735</td>\n",
       "      <td>16368.020508</td>\n",
       "      <td>0.337543</td>\n",
       "      <td>9.466309</td>\n",
       "      <td>39.060158</td>\n",
       "      <td>45.129688</td>\n",
       "      <td>6.139355</td>\n",
       "      <td>6.139355</td>\n",
       "      <td>25.690422</td>\n",
       "      <td>56.693359</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "1     LNG FLOW LNG RUNDOWN TEMPERATURE      HMR FLOW LMR/HMR RATIO  \\\n",
       "0  9469.942383             -160.132996   16287.15625      0.336019   \n",
       "1  9507.515625             -160.059555  16293.646484      0.333414   \n",
       "2  9519.392578             -159.999573  16333.569336       0.33673   \n",
       "3  9497.821289             -160.071472  16283.046875      0.338297   \n",
       "4  9516.681641             -160.009735  16368.020508      0.337543   \n",
       "\n",
       "1 PERCENTAGE NITROGEN IN MIXED REFERIGERANT   \\\n",
       "0                                   9.436035   \n",
       "1                                   9.390137   \n",
       "2                                    9.42041   \n",
       "3                                   9.451172   \n",
       "4                                   9.466309   \n",
       "\n",
       "1 PERCENTAGE METHANE IN MIXED REFERIGERANT   \\\n",
       "0                                 38.877541   \n",
       "1                                 38.877541   \n",
       "2                                 38.982128   \n",
       "3                                 39.033596   \n",
       "4                                 39.060158   \n",
       "\n",
       "1 PERCENTAGE ETHANE IN MIXED REFERIGERANT   \\\n",
       "0                                45.310646   \n",
       "1                                45.310646   \n",
       "2                                45.207718   \n",
       "3                                45.129688   \n",
       "4                                45.129688   \n",
       "\n",
       "1 PERCENTAGE PROPANE IN MIXED REFERIGERANT   \\\n",
       "0                                  6.139355   \n",
       "1                                  6.160547   \n",
       "2                                  6.139355   \n",
       "3                                  6.139355   \n",
       "4                                  6.139355   \n",
       "\n",
       "1 PERCENTAGE BUTHANE IN MIXED REFERIGERANT  AMBIENT TEMPERATURE   \\\n",
       "0                                  6.139355            26.593405   \n",
       "1                                  6.160547            26.280834   \n",
       "2                                  6.139355            25.690422   \n",
       "3                                  6.139355            25.238932   \n",
       "4                                  6.139355            25.690422   \n",
       "\n",
       "1 LOW PRESSSURE PROPANE MIXED REFRIGERANT KETTLE LEVEL  \n",
       "0                                          57.273438    \n",
       "1                                           57.09375    \n",
       "2                                          56.572266    \n",
       "3                                          57.736328    \n",
       "4                                          56.693359    "
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "data.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "9ce0aa0c",
   "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>Unnamed: 0</th>\n",
       "      <th>LNG Flow</th>\n",
       "      <th>LNG Rundown Temperature</th>\n",
       "      <th>MR IGV Position</th>\n",
       "      <th>HMR Flow</th>\n",
       "      <th>HMR/MR Ratio</th>\n",
       "      <th>Percentage Nitrogen in MR</th>\n",
       "      <th>Percentage Methane in MR</th>\n",
       "      <th>Percentage Ethane in MR</th>\n",
       "      <th>Percentage Propane in MR</th>\n",
       "      <th>Percentage Butane in MR</th>\n",
       "      <th>Ambient Condition</th>\n",
       "      <th>Low Pressure Propane/Mixed Refrigerant kettle level</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>18778</td>\n",
       "      <td>13575.379883</td>\n",
       "      <td>-145.473633</td>\n",
       "      <td>67.037796</td>\n",
       "      <td>21107.593750</td>\n",
       "      <td>5.086884</td>\n",
       "      <td>3.413086</td>\n",
       "      <td>42.969826</td>\n",
       "      <td>42.611233</td>\n",
       "      <td>10.676367</td>\n",
       "      <td>0.249805</td>\n",
       "      <td>25.657589</td>\n",
       "      <td>59.044922</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>10519</td>\n",
       "      <td>11672.666016</td>\n",
       "      <td>-145.039062</td>\n",
       "      <td>52.585159</td>\n",
       "      <td>17607.126953</td>\n",
       "      <td>6.215993</td>\n",
       "      <td>2.857910</td>\n",
       "      <td>42.204494</td>\n",
       "      <td>43.970901</td>\n",
       "      <td>10.588183</td>\n",
       "      <td>0.281172</td>\n",
       "      <td>27.051105</td>\n",
       "      <td>60.626953</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>31757</td>\n",
       "      <td>13388.155273</td>\n",
       "      <td>-145.786133</td>\n",
       "      <td>60.087402</td>\n",
       "      <td>20145.201172</td>\n",
       "      <td>4.819857</td>\n",
       "      <td>3.526367</td>\n",
       "      <td>43.268654</td>\n",
       "      <td>42.958206</td>\n",
       "      <td>9.767187</td>\n",
       "      <td>0.307461</td>\n",
       "      <td>29.943972</td>\n",
       "      <td>60.833984</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>6400</td>\n",
       "      <td>12687.996094</td>\n",
       "      <td>-145.610352</td>\n",
       "      <td>77.094765</td>\n",
       "      <td>19162.550781</td>\n",
       "      <td>5.156615</td>\n",
       "      <td>2.703125</td>\n",
       "      <td>44.138577</td>\n",
       "      <td>42.788868</td>\n",
       "      <td>9.953125</td>\n",
       "      <td>0.336055</td>\n",
       "      <td>32.572247</td>\n",
       "      <td>60.548828</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>17354</td>\n",
       "      <td>11994.082031</td>\n",
       "      <td>-145.502930</td>\n",
       "      <td>67.006332</td>\n",
       "      <td>17399.679688</td>\n",
       "      <td>6.411177</td>\n",
       "      <td>3.293945</td>\n",
       "      <td>43.354980</td>\n",
       "      <td>42.325684</td>\n",
       "      <td>10.676367</td>\n",
       "      <td>0.288477</td>\n",
       "      <td>28.709211</td>\n",
       "      <td>57.835938</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   Unnamed: 0      LNG Flow  LNG Rundown Temperature  MR IGV Position  \\\n",
       "0       18778  13575.379883              -145.473633        67.037796   \n",
       "1       10519  11672.666016              -145.039062        52.585159   \n",
       "2       31757  13388.155273              -145.786133        60.087402   \n",
       "3        6400  12687.996094              -145.610352        77.094765   \n",
       "4       17354  11994.082031              -145.502930        67.006332   \n",
       "\n",
       "      HMR Flow   HMR/MR Ratio  Percentage Nitrogen in MR  \\\n",
       "0  21107.593750      5.086884                   3.413086   \n",
       "1  17607.126953      6.215993                   2.857910   \n",
       "2  20145.201172      4.819857                   3.526367   \n",
       "3  19162.550781      5.156615                   2.703125   \n",
       "4  17399.679688      6.411177                   3.293945   \n",
       "\n",
       "   Percentage Methane in MR  Percentage Ethane in MR  \\\n",
       "0                 42.969826                42.611233   \n",
       "1                 42.204494                43.970901   \n",
       "2                 43.268654                42.958206   \n",
       "3                 44.138577                42.788868   \n",
       "4                 43.354980                42.325684   \n",
       "\n",
       "   Percentage Propane in MR  Percentage Butane in MR  Ambient Condition  \\\n",
       "0                 10.676367                 0.249805          25.657589   \n",
       "1                 10.588183                 0.281172          27.051105   \n",
       "2                  9.767187                 0.307461          29.943972   \n",
       "3                  9.953125                 0.336055          32.572247   \n",
       "4                 10.676367                 0.288477          28.709211   \n",
       "\n",
       "   Low Pressure Propane/Mixed Refrigerant kettle level  \n",
       "0                                          59.044922    \n",
       "1                                          60.626953    \n",
       "2                                          60.833984    \n",
       "3                                          60.548828    \n",
       "4                                          57.835938    "
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df = pd.read_csv(\"inputs.csv\")\n",
    "df.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "7b04585e",
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "c:\\Users\\owner\\Downloads\\lng\\conda-env\\lib\\site-packages\\seaborn\\_decorators.py:36: FutureWarning: Pass the following variable as a keyword arg: x. From version 0.12, the only valid positional argument will be `data`, and passing other arguments without an explicit keyword will result in an error or misinterpretation.\n",
      "  warnings.warn(\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAj8AAAJNCAYAAADeRhT9AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAAsTAAALEwEAmpwYAAAaq0lEQVR4nO3df7DldX3f8deHXSzLwpQIzmoA3dg1kU0IEretbdP0Qm1nF4JmWieQKqyZamYggws2zVi8M64zV5JUh0Bo0xlL1N2aMbZWB+ws2I64NZPM2O6Kilmw3BgQt0LMVRRWglz49I9z7u3l7t794b17z919Px7/cM+53+/n+zl8v9+7z/2e797Teu8BAKjilFFPAABgOYkfAKAU8QMAlCJ+AIBSxA8AUIr4AQBKWX0sC59zzjl9/fr1x2kqAwcOHMjatWuP6zZYOezvOuzrWuzvOlbyvt67d+9f9d5fMv/5Y4qf9evXZ8+ePUs3q0PYvXt3xsbGjus2WDns7zrs61rs7zpW8r5urT1yqOe97QUAlCJ+AIBSxA8AUIr4AQBKET8AQCniBwAoRfwAAKWIHwCgFPEDAJQifgCAUsQPAFCK+AEAShE/AEAp4gcAKEX8AACliB8AoBTxAwCUIn4AgFLEDwBQivgBAEoRPwBAKeIHAChF/AAApYgfAKAU8QMAlCJ+AIBSxA8AUIr4AQBKET8AQCniBwAoRfwAAKWIHwCgFPEDAJQifgCAUsQPAFCK+AEAShE/AEAp4gcAKEX8AACliB8AoBTxAwCUIn4AgFLEDwBQivgBAEoRPwBAKeIHAChF/AAApYgfAKAU8QMAlCJ+AIBSxA8AUIr4AQBKWT3qCQCw9G6//fZMTk6ObPubNm3Ktm3bDnp+//79SZJzzz13uae0oA0bNuT6668f9TRYRuIH4CQ0OTmZL331gTx3+otHsv2NFz6bvV9//KDnV/3ge0mSx55ZGX/8rPrBd0Y9BUZgZRx9ACy5505/cZ5+9WUj2fbzp6095LbXPLgrSUY2r/lm5kMt7vkBAEoRPwBAKeIHAChF/AAApYgfAKAU8QMAlCJ+AIBSxA8AUIr4AQBKET8AQCniBwAoRfwAAKWIHwCgFPEDAJQifgCAUsQPAFCK+AEAShE/AEAp4gcAKEX8AACliB8AoBTxAwCUIn4AgFLEDwBQivgBAEoRPwBAKeIHAChF/AAApYgfAKAU8QMAlCJ+AIBSxA8AUIr4AQBKET8AQCniBwAoRfwAAKWIHwCgFPEDAJQifgCAUsQPAFCK+AEAShE/AEAp4gcAKEX8AACliB8AoBTxAwCUIn4AgFLEDwBQivgBAEoRPwBAKeIHAChF/AAApYgfAKAU8QMAlCJ+AIBSxA8AUIr4AQBKET8AQCniBwAoRfwAAKWIHwCgFPEDAJQifgCAUsQPAFCK+AEAShE/AEAp4gcAKEX8AACliB8AoBTxAwCUIn4AgFLEDwBQivgBAEoRPwBAKeIHAChF/AAApYgfAKAU8QMAlCJ+AIBSxA8AUIr4AQBKET8AQCniBwAoRfwAAKWIHwCgFPEDAJQifgCAUsQPAFCK+AEAShE/AEAp4gcAKEX8AACliB8AoBTxAwCUIn4AgFLEDwBQivgBAEoRPwBAKeIHAChF/AAApYgfAKAU8QMAlCJ+AIBSxA8AUIr4AQBKET/Agm6//fbcfvvto54GcBJZCT9XVo9068CKNjk5OeopACeZlfBzxZUfAKAU8QMAlCJ+AIBSxA8AUIr4AQBKET8AQCniBwAoRfwAAKWIHwCgFPEDAJQifgCAUsQPAFCK+AEAShE/AEAp4gcAKEX8AACliB8AoBTxAwCUIn4AgFLEDwBQivgBAEoRPwBAKeIHAChF/AAApYgfAKAU8QMAlCJ+AIBSxA8AUIr4AQBKET8AQCniBwAoRfwAAKWIHwCgFPEDAJQifgCAUsQPAFCK+AEAShE/AEAp4gcAKEX8AACliB8AoBTxAwCUIn4AgFLEDwBQivgBAEoRPwBAKeIHAChF/AAApYgfAKAU8QMAlCJ+AIBSxA8AUIr4AQBKET8AQCniBwAoRfwAAKWIHwCgFPEDAJQifgCAUsQPAFCK+AEAShE/AEAp4gcAKEX8AACliB8AoBTxAwCUIn4AgFLEDwBQivgBAEoRPwBAKeIHAChF/AAApYgfAKAU8QMAlCJ+AIBSxA8AUIr4AQBKET8AQCniBwAoRfwAAKWIHwCgFPEDAJQifgCAUsQPAFCK+AEAShE/AEAp4gcAKEX8AACliB8AoBTxAwCUIn4AgFLEDwBQivgBAEoRPwBAKeIHAChF/AAApYgfAKAU8QMAlCJ+AIBSxA8AUIr4AQBKET8AQCniBwAoRfwAAKWsqPgZGxvL9u3bMzY2lsnJybz97W/Pli1bMjk5mcnJyVx++eW59957c/nll+fOO+/MJZdckquuuipjY2O55JJLcvXVV+e6667L1NRUpqamct111+Xaa6/N5ORk3vGOd2RqaipJMjU1lSuvvDJjY2MZGxvLFVdcMfv161//+kxOTi7r656amnrB/EY15mLmcSzrziw7OTmZD3/4w0dcZzHzuvPOOzM2NpZPf/rTSzbm/HVnjs09e/YsOObRbG/Pnj259NJL87nPfW522ampqVx77bW57rrrDjqOF9ru4daZ2cYdd9yRsbGxXH311Zmampoda7mPfYBRWFHxM9fExEQeeuihPP3005mYmMjExEQOHDiQm2++OQcOHMitt96a3nsee+yxJEnvPY8++mj27duXnTt3ZseOHdm3b18eeOCBTExM5P7778/OnTuTJDt27Mjjjz8+u60nn3xy9uvp6elMTEws62vdsWPHC+Y3qjEXM49jWXdm2YmJiXzjG9844jqLmdett96aJLnllluWbMz5684cm9u3b19wzKPZ3vbt2/P888/nfe973+yyO3bsyAMPPJB9+/YddBwvtN3DrTOzjY9+9KNJkkcffTQ7d+6cHWu5j32AUVgx8TM2NvaCxw8//PALvp55PD09nWQQOwvZtWtX7r777hes33vPPffck8nJyYOuAsz38MMPL9vfgKempnLPPffMzm8prv78KGMuZh7Hsu7cZeful4XWWcy87rzzztnjpPc+u9+X8rXu3bt39th86qmnDjnm0Wxvz549eeqpp5IMjvHee+6+++7s2rVrdpm5/78W2u7k5OSCx/699947u4257rrrrtmxlvPYBxiV1aOewPHw7LPPHvL55557LhMTE4cNpxkTExP5yEc+ssQzO9iOHTvy/PPPJxnMb+fOnbnxxhuXfczFzONY1p277IzDrbOYec1c9Zlxyy235IorrljS1/qe97znoGXmj3k029u+fftB4zz77LOHPFYPt92JiYnZvyDM/97NN998yNc0fxtzj/39+/fn6aefzrZt2w657rF44okn8qlPfWrR43B0Jicnc8oPj/yzrrpT/vr7mZx8ckmO8aqO9dyenJzMmjVrjuOMjuyI8dNa+7Ukv5Yk69aty+7du4/3nI6b6enpF1xROpyHH354WV7rZz7zmdk/rKanp3PPPffk4osvXvYxFzOPY1l37rIzDrfOYuY1/w/13nt27969pK/1UFdS5o95NNs71DgLRfrhtrvQ8X2oIFrI3GP/mWeeyfT0dJ544omjXn8hzz333JKMw9E588wzc+5ZP5aLL1o7ku2fc/qqvO0Q27730VVJkktHNK/57nvyFfn+E991bC7CsZ7b09PTeeaZZ0baE0eMn977B5N8MEk2bdrU5789dSJZvXp1zjvvvKMKoPXr1x/0Vtzx8MUvfjG7du3K9PR0Vq9enc2bNy96uz/KmIuZx7GsO3fZGYdbZzHzaq29ICBaaxkbG1vS13raaacdFCLzxzya7Z1xxhkHjTN//nPHX2i75513Xh555JGD1lu9enCqH00EzT32Z/42d9tttx1xvSPZvXv3spxTDGzbti17v/549n75wEi2/7aL1uaOQ2x7zYHnkuSQ3xuFNfu/m9e+ct2SHONVHeu5PXOVbZQ/D1bMPT9L6dRTT82pp5560POrVq3K+Ph4WmtHHGN8fPx4TO0gW7duzSmnDHbDqlWrcs0114xkzMXM41jWnbvsjMOts5h53XDDDS94/M53vnPRY85f973vfe9By8wf82i2d6i3vU499dTZaJk//kLbHR8fX3Cdm2666ZCvaf75sFzHPsCorJj4mX/5a/369S/4eubxzA/2wwXMZZddli1btrxg/dZaNm/enA0bNuSKK6447FzWr1+fDRs2HNsL+BGdffbZ2bx58+z8zj777JGMuZh5HMu6c5edu18WWmcx83rjG984e5y01mb3+1K+1te+9rWzx+YZZ5xxyDGPZnubNm3KGWeckWRwjLfWsmXLllx22WWzy8z9/7XQdjds2LDgsX/ppZfObmOuN7zhDbNjLeexDzAqKyZ+5hsfH8+rXvWqrFmzJuPj4xkfH8/atWtz0003Ze3atbnhhhvSWstLX/rSJIM/3M4///xs3Lgx11xzTbZu3ZqNGzfmggsuyPj4eC688MLZv3Fv3bo169atm93WmWeeOfv16tWrl/1vvlu3bn3B/EY15mLmcSzrziw7Pj6el7/85UdcZzHzmrn6M3PVZynGnL/uzLG5ffv2Bcc8mu1t3749p5xySt797nfPLrt169ZccMEF2bhx40HH8ULbPdw6M9t4y1vekiQ5//zzc80118yO5aoPUEE7mn/5NGPTpk19z549x3E67guoxv5e2Wbem3fPz4ln5p6fp1992ZEXPg4WvOfnwcGvbxjVvOZb8+Au9/ws0o96z89y/D9vre3tvW+a//yKvfIDAHA8iB8AoBTxAwCUIn4AgFLEDwBQivgBAEoRPwBAKeIHAChF/AAApYgfAKAU8QMAlCJ+AIBSxA8AUIr4AQBKET8AQCniBwAoRfwAAKWIHwCgFPEDAJQifgCAUsQPAFCK+AEAShE/AEAp4gcAKEX8AACliB8AoBTxAwCUIn4AgFLEDwBQivgBAEoRPwBAKeIHAChF/AAApYgfAKAU8QMAlCJ+AIBSxA8AUIr4AQBKET8AQCniBwAoRfwAAKWIHwCgFPEDAJQifgCAUsQPAFCK+AEAShE/AEAp4gcAKEX8AACliB8AoBTxAwCUIn4AgFLEDwBQivgBAEoRPwBAKeIHAChF/AAApYgfAKAU8QMAlCJ+AIBSxA8AUIr4AQBKET8AQCniBwAoRfwAAKWIHwCgFPEDAJQifgCAUsQPAFCK+AEAShE/AEAp4gcAKEX8AACliB8AoBTxAwCUIn4AgFLEDwBQivgBAEoRPwBAKeIHAChF/AAApYgfAKAU8QMAlCJ+AIBSxA8AUIr4AQBKET8AQCniBwAoRfwAAKWIHwCgFPEDAJQifgCAUsQPAFCK+AEAShE/AEAp4gcAKEX8AACliB8AoBTxAwCUIn4AgFLEDwBQivgBAEoRPwBAKatHPQFg5dqwYcOopwCcZFbCzxXxAyzo+uuvH/UUgJPMSvi54m0vAKAU8QMAlCJ+AIBSxA8AUIr4AQBKET8AQCniBwAoRfwAAKWIHwCgFPEDAJQifgCAUsQPAFCK+AEAShE/AEAp4gcAKEX8AACliB8AoBTxAwCUIn4AgFLEDwBQivgBAEoRPwBAKeIHAChF/AAApYgfAKAU8QMAlCJ+AIBSxA8AUIr4AQBKET8AQCniBwAoRfwAAKWIHwCgFPEDAJQifgCAUsQPAFCK+AEAShE/AEAp4gcAKEX8AACliB8AoBTxAwCUIn4AgFLEDwBQivgBAEoRPwBAKeIHAChF/AAApYgfAKAU8QMAlCJ+AIBSxA8AUIr4AQBKET8AQCniBwAoRfwAAKWIHwCgFPEDAJQifgCAUsQPAFCK+AEAShE/AEAp4gcAKEX8AACliB8AoBTxAwCUIn4AgFLEDwBQivgBAEoRPwBAKeIHAChF/AAApYgfAKAU8QMAlCJ+AIBSxA8AUIr4AQBKET8AQCniBwAoRfwAAKWIHwCgFPEDAJQifgCAUsQPAFCK+AEAShE/AEAp4gcAKEX8AACliB8AoBTxAwCUIn4AgFLEDwBQivgBAEoRPwBAKeIHAChF/AAApYgfAKAU8QMAlCJ+AIBSxA8AUIr4AQBKET8AQCniBwAoRfwAAKWIHwCglNWjngAAx8eqH3wnax7cNZJtn/JT/zRrHvzvBz2/6gdTSTKyec236gffSbJu1NNgmYkfgJPQhg0bRrr90//GqXntKw+Oiv37p5Mk5567UoJj3cj/X7H8xA/ASej6668f6fZ3796dq6++eqRzgIW45wcAKEX8AACliB8AoBTxAwCUIn4AgFLEDwBQivgBAEoRPwBAKeIHAChF/AAApYgfAKAU8QMAlCJ+AIBSxA8AUIr4AQBKET8AQCniBwAoRfwAAKWIHwCgFPEDAJQifgCAUsQPAFCK+AEAShE/AEAp4gcAKEX8AACliB8AoBTxAwCUIn4AgFLEDwBQivgBAEoRPwBAKeIHAChF/AAApYgfAKAU8QMAlCJ+AIBSxA8AUIr4AQBKET8AQCniBwAoRfwAAKWIHwCgFPEDAJQifgCAUsQPAFCK+AEAShE/AEAp4gcAKEX8AACliB8AoBTxAwCUIn4AgFJa7/3oF27t20keOX7TSZKck+SvjvM2WDns7zrs61rs7zpW8r5+Re/9JfOfPKb4WQ6ttT29902jngfLw/6uw76uxf6u40Tc1972AgBKET8AQCkrMX4+OOoJsKzs7zrs61rs7zpOuH294u75AQA4nlbilR8AgONmRcVPa21za+1rrbXJ1tq7Rj0flk5r7fzW2udaa/taa3/WWts2fP7FrbX/0Vp7aPjfHxv1XFkarbVVrbX7Wmv/bfj4J1prXxie3x9vrb1o1HNkabTWzmqtfaK19mBr7YHW2t9zbp+8Wms3Dn+Of7W19rHW2mkn2vm9YuKntbYqyb9PsiXJxiS/0lrbONpZsYSmk/yr3vvGJK9L8uvD/fuuJJ/tvb8qyWeHjzk5bEvywJzHv5Pkd3vvG5J8N8m/HMmsOB5uS3JP7/3VSS7KYL87t09CrbVzk7wjyabe+88kWZXkqpxg5/eKiZ8kfyfJZO/96733Hyb5oyRvHPGcWCK992/13r84/PrJDH44npvBPt4xXGxHkl8ayQRZUq2185JcnuSO4eOW5NIknxguYl+fJFprfzPJLyT5gyTpvf+w9/5EnNsns9VJ1rTWVic5Pcm3coKd3yspfs5N8uicx98cPsdJprW2PsnFSb6QZF3v/VvDbz2WZN2o5sWSujXJbyZ5fvj47CRP9N6nh4+d3yePn0jy7SQfHr7NeUdrbW2c2yel3vv+JB9I8o0Moud7SfbmBDu/V1L8UEBr7Ywk/zXJDb3378/9Xh/800P//PAE11r7xSR/2XvfO+q5sCxWJ/m5JP+h935xkgOZ9xaXc/vkMbx3640ZRO+PJ1mbZPNIJ/UjWEnxsz/J+XMenzd8jpNEa+3UDMLnD3vvnxw+/Xhr7WXD778syV+Oan4smX+Q5A2ttYczePv60gzuCTlreJk8cX6fTL6Z5Ju99y8MH38igxhybp+cXp/kL3rv3+69P5vkkxmc8yfU+b2S4ud/J3nV8I7xF2VwA9VdI54TS2R4z8cfJHmg937LnG/dlWTr8OutSe5c7rmxtHrv/6b3fl7vfX0G5/G9vfc3J/lckjcNF7OvTxK998eSPNpa+6nhU/84yb44t09W30jyutba6cOf6zP7+4Q6v1fULzlsrV2Wwb0Cq5J8qPf+vtHOiKXSWvv5JH+c5P78//tAbsrgvp//nOTlSR5J8su99++MZJIsudbaWJLf6L3/YmvtlRlcCXpxkvuSvKX3/swIp8cSaa29JoOb21+U5OtJfjWDv1w7t09CrbX3Jrkyg3/Fe1+St2Vwj88Jc36vqPgBADjeVtLbXgAAx534AQBKET8AQCniBwAoRfwAAKWIHyiktdZbax+d83h1a+3bcz55/a3Dx18afkL3jQuM89bW2r+b8/gtrbWvDD/p+cvDjzg4q7X2ntbab81b9zWttQcOMebu1trXhuv/yZzfG3Msr2/XcLtntdaum/P8j7fWPnG4dYE6xA/UciDJz7TW1gwf/5Mc/JtYP957f00Gv7X13a2183MYrbXNSW5MsqX3/tMZ/HbfP83gs5w+lsHvA5nrquHzh/Lm3vtFGXww4vuP6hXN0Xu/bPihmmcluW7O8/+39/6mhdYDahE/UM+uDD5xPUl+JQuESO99KslkkpcdYbx3Z/CLDPcP13uu9/6h3vvXeu//J8l3W2t/d87yv7zQNuf4fJINbeD9rbWvttbub61dmQw+LqG19vnhFaqvttb+4fD5h1tr5yT57SR/a/j997fW1rfWvjpc5rTW2oeH493XWrtk+PxbW2ufbK3d01p7qLX2b48wR+AEJX6gnj9KclVr7bQkP5vBb9k+SGvt5UlOS/KVI4z300m+eJjvfyyDqz1prb0uyXd67w8dYcwrMvht4P8syWuSXJTBZwq9f/g5Uf8iyWeGV6guSvKleeu/K8mf995f03v/1/O+9+sZfNbmhRnE347h/4sMt3VlkguTXHmkq17AiUn8QDG9968kWZ/BH/y7DrHIla21r2Rw1ef3e+9/fbRjt9YuHF5t+fOZqzRJPp7kTa21U3L4t7yS5A9ba1/K4C2330jy80k+Nrya9HiS/5nkb2fwWYC/2lrbnuTC3vuTRzvH4ZgfTZLe+4MZfPTCTw6/99ne+/eGr3lfklccw7jACUL8QE13JflADh0iH++9/2ySv5/kt1trLz3CWH+WwX0+6b3fP7wac3eSNcPnHk3yF0n+UZJ/nkEMLeTNw6s1vzRc75B6759P8gsZ3K/0kdbaNUeY49Ga+1lEzyVZvdCCwIlL/EBNH0ry3t77/Qst0Hvfk+Q/Jdl2hLF+K8kHWmvnzXluzbxlPpbkd5N8vff+zWOY5x9ncCVqVWvtJRkEz/9qrb0iyeO99/+YwQdq/ty89Z5McuZhxnxzkrTWfjKDD9782jHMCTjBiR8oqPf+zd777x3For+TwdtLC4VEeu+7kvxekrtba/taa3+awVWTz8xZ7L9kcG/QkW50nu9TGdxz9OUk9yb5zd77Y0nGkny5tXZfBvfo3DZvTlNJ/mR4M/T8fzX2+0lOaa3dn8FVqLeu5E+fBpaeT3UHAEpx5QcAKEX8AACliB8AoBTxAwCUIn4AgFLEDwBQivgBAEoRPwBAKf8Pj+ku58oieNIAAAAASUVORK5CYII=",
      "text/plain": [
       "<Figure size 720x720 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "import seaborn as sns\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "fig, ax = plt.subplots(nrows=1, figsize=(10, 10))\n",
    "sns.boxplot(df['MR IGV Position'], ax=ax)\n",
    "ax.grid()\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "id": "be906e49",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(59.49700975316839, 13.074677920582042)"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "mean = np.mean(df['MR IGV Position'])\n",
    "std = np.std(df['MR IGV Position'])\n",
    "mean, std"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "id": "b8e33be3",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([70.94910264, 41.63668832, 65.22293378, ..., 51.66631638,\n",
       "       63.03461718, 67.29762789])"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "size = len(data)\n",
    "random_MR_IGV_position = np.random.normal(mean, std, size)\n",
    "random_MR_IGV_position"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "id": "7b7dfd60",
   "metadata": {},
   "outputs": [],
   "source": [
    "data['MR IGV Position'] = random_MR_IGV_position"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "id": "7b1e50e3",
   "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>1</th>\n",
       "      <th>LNG FLOW</th>\n",
       "      <th>LNG RUNDOWN TEMPERATURE</th>\n",
       "      <th>HMR FLOW</th>\n",
       "      <th>LMR/HMR RATIO</th>\n",
       "      <th>PERCENTAGE NITROGEN IN MIXED REFERIGERANT</th>\n",
       "      <th>PERCENTAGE METHANE IN MIXED REFERIGERANT</th>\n",
       "      <th>PERCENTAGE ETHANE IN MIXED REFERIGERANT</th>\n",
       "      <th>PERCENTAGE PROPANE IN MIXED REFERIGERANT</th>\n",
       "      <th>PERCENTAGE BUTHANE IN MIXED REFERIGERANT</th>\n",
       "      <th>AMBIENT TEMPERATURE</th>\n",
       "      <th>LOW PRESSSURE PROPANE MIXED REFRIGERANT KETTLE LEVEL</th>\n",
       "      <th>MR IGV Position</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>9469.942383</td>\n",
       "      <td>-160.132996</td>\n",
       "      <td>16287.15625</td>\n",
       "      <td>0.336019</td>\n",
       "      <td>9.436035</td>\n",
       "      <td>38.877541</td>\n",
       "      <td>45.310646</td>\n",
       "      <td>6.139355</td>\n",
       "      <td>6.139355</td>\n",
       "      <td>26.593405</td>\n",
       "      <td>57.273438</td>\n",
       "      <td>70.949103</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>9507.515625</td>\n",
       "      <td>-160.059555</td>\n",
       "      <td>16293.646484</td>\n",
       "      <td>0.333414</td>\n",
       "      <td>9.390137</td>\n",
       "      <td>38.877541</td>\n",
       "      <td>45.310646</td>\n",
       "      <td>6.160547</td>\n",
       "      <td>6.160547</td>\n",
       "      <td>26.280834</td>\n",
       "      <td>57.09375</td>\n",
       "      <td>41.636688</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>9519.392578</td>\n",
       "      <td>-159.999573</td>\n",
       "      <td>16333.569336</td>\n",
       "      <td>0.33673</td>\n",
       "      <td>9.42041</td>\n",
       "      <td>38.982128</td>\n",
       "      <td>45.207718</td>\n",
       "      <td>6.139355</td>\n",
       "      <td>6.139355</td>\n",
       "      <td>25.690422</td>\n",
       "      <td>56.572266</td>\n",
       "      <td>65.222934</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>9497.821289</td>\n",
       "      <td>-160.071472</td>\n",
       "      <td>16283.046875</td>\n",
       "      <td>0.338297</td>\n",
       "      <td>9.451172</td>\n",
       "      <td>39.033596</td>\n",
       "      <td>45.129688</td>\n",
       "      <td>6.139355</td>\n",
       "      <td>6.139355</td>\n",
       "      <td>25.238932</td>\n",
       "      <td>57.736328</td>\n",
       "      <td>63.372071</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>9516.681641</td>\n",
       "      <td>-160.009735</td>\n",
       "      <td>16368.020508</td>\n",
       "      <td>0.337543</td>\n",
       "      <td>9.466309</td>\n",
       "      <td>39.060158</td>\n",
       "      <td>45.129688</td>\n",
       "      <td>6.139355</td>\n",
       "      <td>6.139355</td>\n",
       "      <td>25.690422</td>\n",
       "      <td>56.693359</td>\n",
       "      <td>54.960265</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "1     LNG FLOW LNG RUNDOWN TEMPERATURE      HMR FLOW LMR/HMR RATIO  \\\n",
       "0  9469.942383             -160.132996   16287.15625      0.336019   \n",
       "1  9507.515625             -160.059555  16293.646484      0.333414   \n",
       "2  9519.392578             -159.999573  16333.569336       0.33673   \n",
       "3  9497.821289             -160.071472  16283.046875      0.338297   \n",
       "4  9516.681641             -160.009735  16368.020508      0.337543   \n",
       "\n",
       "1 PERCENTAGE NITROGEN IN MIXED REFERIGERANT   \\\n",
       "0                                   9.436035   \n",
       "1                                   9.390137   \n",
       "2                                    9.42041   \n",
       "3                                   9.451172   \n",
       "4                                   9.466309   \n",
       "\n",
       "1 PERCENTAGE METHANE IN MIXED REFERIGERANT   \\\n",
       "0                                 38.877541   \n",
       "1                                 38.877541   \n",
       "2                                 38.982128   \n",
       "3                                 39.033596   \n",
       "4                                 39.060158   \n",
       "\n",
       "1 PERCENTAGE ETHANE IN MIXED REFERIGERANT   \\\n",
       "0                                45.310646   \n",
       "1                                45.310646   \n",
       "2                                45.207718   \n",
       "3                                45.129688   \n",
       "4                                45.129688   \n",
       "\n",
       "1 PERCENTAGE PROPANE IN MIXED REFERIGERANT   \\\n",
       "0                                  6.139355   \n",
       "1                                  6.160547   \n",
       "2                                  6.139355   \n",
       "3                                  6.139355   \n",
       "4                                  6.139355   \n",
       "\n",
       "1 PERCENTAGE BUTHANE IN MIXED REFERIGERANT  AMBIENT TEMPERATURE   \\\n",
       "0                                  6.139355            26.593405   \n",
       "1                                  6.160547            26.280834   \n",
       "2                                  6.139355            25.690422   \n",
       "3                                  6.139355            25.238932   \n",
       "4                                  6.139355            25.690422   \n",
       "\n",
       "1 LOW PRESSSURE PROPANE MIXED REFRIGERANT KETTLE LEVEL  MR IGV Position  \n",
       "0                                          57.273438          70.949103  \n",
       "1                                           57.09375          41.636688  \n",
       "2                                          56.572266          65.222934  \n",
       "3                                          57.736328          63.372071  \n",
       "4                                          56.693359          54.960265  "
      ]
     },
     "execution_count": 13,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "data.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "id": "8af2c91d",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "['LNG FLOW',\n",
       " 'LNG RUNDOWN TEMPERATURE',\n",
       " 'HMR FLOW',\n",
       " 'LMR/HMR RATIO',\n",
       " 'PERCENTAGE NITROGEN IN MIXED REFERIGERANT ',\n",
       " 'PERCENTAGE METHANE IN MIXED REFERIGERANT ',\n",
       " 'PERCENTAGE ETHANE IN MIXED REFERIGERANT ',\n",
       " 'PERCENTAGE PROPANE IN MIXED REFERIGERANT ',\n",
       " 'PERCENTAGE BUTHANE IN MIXED REFERIGERANT ',\n",
       " 'AMBIENT TEMPERATURE ',\n",
       " 'LOW PRESSSURE PROPANE MIXED REFRIGERANT KETTLE LEVEL',\n",
       " 'PR PEAK MARGIN',\n",
       " 'MR PEAK MARGIN']"
      ]
     },
     "execution_count": 14,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "list(headers)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "id": "afe2e0fd",
   "metadata": {},
   "outputs": [],
   "source": [
    "headers = [header.strip() for header in data.columns]\n",
    "data.columns = headers"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "id": "0321e19b",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Index(['LNG FLOW', 'LNG RUNDOWN TEMPERATURE', 'HMR FLOW', 'LMR/HMR RATIO',\n",
       "       'PERCENTAGE NITROGEN IN MIXED REFERIGERANT',\n",
       "       'PERCENTAGE METHANE IN MIXED REFERIGERANT',\n",
       "       'PERCENTAGE ETHANE IN MIXED REFERIGERANT',\n",
       "       'PERCENTAGE PROPANE IN MIXED REFERIGERANT',\n",
       "       'PERCENTAGE BUTHANE IN MIXED REFERIGERANT', 'AMBIENT TEMPERATURE',\n",
       "       'LOW PRESSSURE PROPANE MIXED REFRIGERANT KETTLE LEVEL',\n",
       "       'MR IGV Position'],\n",
       "      dtype='object')"
      ]
     },
     "execution_count": 16,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "data.columns"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "id": "12548dab",
   "metadata": {},
   "outputs": [],
   "source": [
    "new_header_order = ['LNG FLOW', 'LNG RUNDOWN TEMPERATURE','MR IGV Position', 'HMR FLOW', 'LMR/HMR RATIO',\n",
    "       'PERCENTAGE NITROGEN IN MIXED REFERIGERANT',\n",
    "       'PERCENTAGE METHANE IN MIXED REFERIGERANT',\n",
    "       'PERCENTAGE ETHANE IN MIXED REFERIGERANT',\n",
    "       'PERCENTAGE PROPANE IN MIXED REFERIGERANT',\n",
    "       'PERCENTAGE BUTHANE IN MIXED REFERIGERANT', 'AMBIENT TEMPERATURE',\n",
    "       'LOW PRESSSURE PROPANE MIXED REFRIGERANT KETTLE LEVEL']"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "id": "72a42fec",
   "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>LNG FLOW</th>\n",
       "      <th>LNG RUNDOWN TEMPERATURE</th>\n",
       "      <th>MR IGV Position</th>\n",
       "      <th>HMR FLOW</th>\n",
       "      <th>LMR/HMR RATIO</th>\n",
       "      <th>PERCENTAGE NITROGEN IN MIXED REFERIGERANT</th>\n",
       "      <th>PERCENTAGE METHANE IN MIXED REFERIGERANT</th>\n",
       "      <th>PERCENTAGE ETHANE IN MIXED REFERIGERANT</th>\n",
       "      <th>PERCENTAGE PROPANE IN MIXED REFERIGERANT</th>\n",
       "      <th>PERCENTAGE BUTHANE IN MIXED REFERIGERANT</th>\n",
       "      <th>AMBIENT TEMPERATURE</th>\n",
       "      <th>LOW PRESSSURE PROPANE MIXED REFRIGERANT KETTLE LEVEL</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>9469.942383</td>\n",
       "      <td>-160.132996</td>\n",
       "      <td>70.949103</td>\n",
       "      <td>16287.15625</td>\n",
       "      <td>0.336019</td>\n",
       "      <td>9.436035</td>\n",
       "      <td>38.877541</td>\n",
       "      <td>45.310646</td>\n",
       "      <td>6.139355</td>\n",
       "      <td>6.139355</td>\n",
       "      <td>26.593405</td>\n",
       "      <td>57.273438</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>9507.515625</td>\n",
       "      <td>-160.059555</td>\n",
       "      <td>41.636688</td>\n",
       "      <td>16293.646484</td>\n",
       "      <td>0.333414</td>\n",
       "      <td>9.390137</td>\n",
       "      <td>38.877541</td>\n",
       "      <td>45.310646</td>\n",
       "      <td>6.160547</td>\n",
       "      <td>6.160547</td>\n",
       "      <td>26.280834</td>\n",
       "      <td>57.09375</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>9519.392578</td>\n",
       "      <td>-159.999573</td>\n",
       "      <td>65.222934</td>\n",
       "      <td>16333.569336</td>\n",
       "      <td>0.33673</td>\n",
       "      <td>9.42041</td>\n",
       "      <td>38.982128</td>\n",
       "      <td>45.207718</td>\n",
       "      <td>6.139355</td>\n",
       "      <td>6.139355</td>\n",
       "      <td>25.690422</td>\n",
       "      <td>56.572266</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>9497.821289</td>\n",
       "      <td>-160.071472</td>\n",
       "      <td>63.372071</td>\n",
       "      <td>16283.046875</td>\n",
       "      <td>0.338297</td>\n",
       "      <td>9.451172</td>\n",
       "      <td>39.033596</td>\n",
       "      <td>45.129688</td>\n",
       "      <td>6.139355</td>\n",
       "      <td>6.139355</td>\n",
       "      <td>25.238932</td>\n",
       "      <td>57.736328</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>9516.681641</td>\n",
       "      <td>-160.009735</td>\n",
       "      <td>54.960265</td>\n",
       "      <td>16368.020508</td>\n",
       "      <td>0.337543</td>\n",
       "      <td>9.466309</td>\n",
       "      <td>39.060158</td>\n",
       "      <td>45.129688</td>\n",
       "      <td>6.139355</td>\n",
       "      <td>6.139355</td>\n",
       "      <td>25.690422</td>\n",
       "      <td>56.693359</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "      LNG FLOW LNG RUNDOWN TEMPERATURE  MR IGV Position      HMR FLOW  \\\n",
       "0  9469.942383             -160.132996        70.949103   16287.15625   \n",
       "1  9507.515625             -160.059555        41.636688  16293.646484   \n",
       "2  9519.392578             -159.999573        65.222934  16333.569336   \n",
       "3  9497.821289             -160.071472        63.372071  16283.046875   \n",
       "4  9516.681641             -160.009735        54.960265  16368.020508   \n",
       "\n",
       "  LMR/HMR RATIO PERCENTAGE NITROGEN IN MIXED REFERIGERANT  \\\n",
       "0      0.336019                                  9.436035   \n",
       "1      0.333414                                  9.390137   \n",
       "2       0.33673                                   9.42041   \n",
       "3      0.338297                                  9.451172   \n",
       "4      0.337543                                  9.466309   \n",
       "\n",
       "  PERCENTAGE METHANE IN MIXED REFERIGERANT  \\\n",
       "0                                38.877541   \n",
       "1                                38.877541   \n",
       "2                                38.982128   \n",
       "3                                39.033596   \n",
       "4                                39.060158   \n",
       "\n",
       "  PERCENTAGE ETHANE IN MIXED REFERIGERANT  \\\n",
       "0                               45.310646   \n",
       "1                               45.310646   \n",
       "2                               45.207718   \n",
       "3                               45.129688   \n",
       "4                               45.129688   \n",
       "\n",
       "  PERCENTAGE PROPANE IN MIXED REFERIGERANT  \\\n",
       "0                                 6.139355   \n",
       "1                                 6.160547   \n",
       "2                                 6.139355   \n",
       "3                                 6.139355   \n",
       "4                                 6.139355   \n",
       "\n",
       "  PERCENTAGE BUTHANE IN MIXED REFERIGERANT AMBIENT TEMPERATURE  \\\n",
       "0                                 6.139355           26.593405   \n",
       "1                                 6.160547           26.280834   \n",
       "2                                 6.139355           25.690422   \n",
       "3                                 6.139355           25.238932   \n",
       "4                                 6.139355           25.690422   \n",
       "\n",
       "  LOW PRESSSURE PROPANE MIXED REFRIGERANT KETTLE LEVEL  \n",
       "0                                          57.273438    \n",
       "1                                           57.09375    \n",
       "2                                          56.572266    \n",
       "3                                          57.736328    \n",
       "4                                          56.693359    "
      ]
     },
     "execution_count": 18,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "data = data.reindex(columns=new_header_order)\n",
    "data.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "id": "4031162c",
   "metadata": {},
   "outputs": [],
   "source": [
    "data.to_csv('new_inputs.csv')"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "id": "6d6e6472",
   "metadata": {},
   "source": [
    "# Generate optimization results"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "id": "387d376f",
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "c:\\Users\\owner\\Downloads\\lng\\conda-env\\lib\\site-packages\\sklearn\\base.py:329: UserWarning: Trying to unpickle estimator StandardScaler from version 1.2.0 when using version 1.1.0. This might lead to breaking code or invalid results. Use at your own risk. For more info please refer to:\n",
      "https://scikit-learn.org/stable/modules/model_persistence.html#security-maintainability-limitations\n",
      "  warnings.warn(\n",
      "c:\\Users\\owner\\Downloads\\lng\\conda-env\\lib\\site-packages\\sklearn\\base.py:329: UserWarning: Trying to unpickle estimator KNeighborsRegressor from version 1.2.0 when using version 1.1.0. This might lead to breaking code or invalid results. Use at your own risk. For more info please refer to:\n",
      "https://scikit-learn.org/stable/modules/model_persistence.html#security-maintainability-limitations\n",
      "  warnings.warn(\n",
      "c:\\Users\\owner\\Downloads\\lng\\conda-env\\lib\\site-packages\\sklearn\\base.py:329: UserWarning: Trying to unpickle estimator RegressorChain from version 1.2.0 when using version 1.1.0. This might lead to breaking code or invalid results. Use at your own risk. For more info please refer to:\n",
      "https://scikit-learn.org/stable/modules/model_persistence.html#security-maintainability-limitations\n",
      "  warnings.warn(\n",
      "c:\\Users\\owner\\Downloads\\lng\\conda-env\\lib\\site-packages\\sklearn\\base.py:329: UserWarning: Trying to unpickle estimator Pipeline from version 1.2.0 when using version 1.1.0. This might lead to breaking code or invalid results. Use at your own risk. For more info please refer to:\n",
      "https://scikit-learn.org/stable/modules/model_persistence.html#security-maintainability-limitations\n",
      "  warnings.warn(\n"
     ]
    }
   ],
   "source": [
    "from input_search_algorithm import Optimize\n",
    "import pickle\n",
    "import csv\n",
    "\n",
    "with open('model.pkl', 'rb') as obj:\n",
    "    model = pickle.load(obj)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "id": "20bba186",
   "metadata": {},
   "outputs": [],
   "source": [
    "spaces = [(0.0, 82.99989318847656), \n",
    "        (0.0, 22607.90234375), \n",
    "        (0.0, 755.8787231445312), \n",
    "        (0.0, 24.98193359375), \n",
    "        (-0.00830078125, 84.9219741821289), \n",
    "        (-0.0166015625, 84.89873504638672), \n",
    "        (0.0, 20.802440643310547), \n",
    "        (-0.0007421874906867, 0.9607030749320984), \n",
    "        (35, 37.24669647216797), \n",
    "        (50, 90)]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "id": "ad2d9b98",
   "metadata": {},
   "outputs": [],
   "source": [
    "headers = ['MR IGV Position', 'HMR FLOW', 'LMR/HMR RATIO',\n",
    "       'PERCENTAGE NITROGEN IN MIXED REFERIGERANT',\n",
    "       'PERCENTAGE METHANE IN MIXED REFERIGERANT',\n",
    "       'PERCENTAGE ETHANE IN MIXED REFERIGERANT',\n",
    "       'PERCENTAGE PROPANE IN MIXED REFERIGERANT',\n",
    "       'PERCENTAGE BUTHANE IN MIXED REFERIGERANT', 'AMBIENT TEMPERATURE',\n",
    "       'LOW PRESSSURE PROPANE MIXED REFRIGERANT KETTLE LEVEL', 'LNG FLOW', 'LNG RUNDOWN TEMPERATURE',]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "id": "fd1b8681",
   "metadata": {},
   "outputs": [],
   "source": [
    "from tqdm import tqdm"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "52ddd757",
   "metadata": {},
   "outputs": [],
   "source": [
    "with open(\"new_test_dataset_analysis.csv\", 'w', newline='') as f_obj:\n",
    "    writer = csv.writer(f_obj,)\n",
    "    writer.writerow(headers)\n",
    "    for i, row in tqdm(list(data.iterrows())[:50]):\n",
    "        row = list(row)\n",
    "        actual_output = np.array(row[:2])\n",
    "        actual_input = np.array(row[2:])\n",
    "        actual_prediction = model.predict(actual_input.reshape(1, -1))\n",
    "        p = Optimize(model.predict, actual_input.copy(), spaces, actual_output.reshape(1, -1))\n",
    "        result = p.maximize()\n",
    "        optimized_input = result[0]\n",
    "        predicted_output = result[1]\n",
    "        \n",
    "        writer.writerow(list(actual_input) + list(actual_output) + [actual_prediction])\n",
    "        writer.writerow(list(optimized_input) + list(predicted_output[0]))\n",
    "\n",
    "        writer.writerow(list())"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3.10.0 ('conda-env')",
   "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.10.0"
  },
  "vscode": {
   "interpreter": {
    "hash": "72c4305ed16563862b7b8f497eeac66b171a2e02864ff61cfb80a89c78fe124b"
   }
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}