lhs_dict.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "traditional-carnival",
   "metadata": {},
   "outputs": [],
   "source": [
    "#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++\n",
    "# Script to set ensemsble dictionary with Latin Hypercube Sampling\n",
    "# Mikayla Pascual 07/22/22\n",
    "#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++\n",
    "\n",
    "#import \n",
    "from MatlabFuncs import *\n",
    "from model import *\n",
    "from loadmodel import *\n",
    "from os.path import exists\n",
    "import lhsmdu\n",
    "import numpy as np\n",
    "import pickle"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "detected-while",
   "metadata": {},
   "outputs": [],
   "source": [
    "### frontal forcings ###\n",
    "#from frontalforcingsrignot import frontalforcingsrignot\n",
    "#md.frontalforcings = frontalforcingsrignot()\n",
    "#md.frontalforcings.subglacial_discharge = 2.6 * np.ones(md.mesh.numberofvertices)\n",
    "#md.frontalforcings.thermalforcing = 15.7 * np.ones(md.mesh.numberofvertices)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "6e9df5af-c066-4ad6-969c-aecd788f7727",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Opening ./Models/SAtoES_friction_coefficient.nc for reading \n",
      "group toolkits is empty\n",
      "empty init\n",
      "Opening ./Models/SAtoES_inversion.nc for reading \n",
      "empty init\n"
     ]
    }
   ],
   "source": [
    "# Load a separate model to extract the friction coefficient and rheology\n",
    "md_friction = loadmodel(\"./Models/SAtoES_friction_coefficient.nc\")\n",
    "friction_coefficient = md_friction.results.StressbalanceSolution.FrictionCoefficient\n",
    "rheology_B = md_friction.materials.rheology_B\n",
    "md_friction = None\n",
    "\n",
    "# Load the original model\n",
    "md = loadmodel(\"./Models/SAtoES_inversion.nc\")\n",
    "md.friction.coefficient = friction_coefficient\n",
    "md.materials.rheology_B = rheology_B"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "instrumental-franklin",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Basal shear stress parameters: Sigma_b = coefficient^2 * Neff ^r * |u_b|^(s - 1) * u_b,\n",
       "(effective stress Neff = rho_ice * g * thickness + rho_water * g * base, r = q / p and s = 1 / p)\n",
       "         coefficient            : (12112,)        -- friction coefficient [SI]\n",
       "         p                      : (24063,)        -- p exponent\n",
       "         q                      : (24063,)        -- q exponent\n",
       "         coupling               : 0               -- Coupling flag 0: uniform sheet (negative pressure ok, default), 1: ice pressure only, 2: water pressure assuming uniform sheet (no negative pressure), 3: use provided effective_pressure, 4: used coupled model (not implemented yet)\n",
       "         effective_pressure     : N/A             -- Effective Pressure for the forcing if not coupled [Pa]\n",
       "         effective_pressure_l...: 0               -- Neff do not allow to fall below a certain limit: effective_pressure_limit * rho_ice * g * thickness (default 0)"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "### basal sliding ###\n",
    "md.friction"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "id": "egyptian-booth",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "   Materials:\n",
       "         rho_ice                : 917.0           -- ice density [kg/m^3]\n",
       "         rho_water              : 1023.0          -- water density [kg/m^3]\n",
       "         rho_freshwater         : 1000.0          -- fresh water density [kg/m^3]\n",
       "         mu_water               : 0.001787        -- water viscosity [Ns/m^2]\n",
       "         heatcapacity           : 2093.0          -- heat capacity [J/kg/K]\n",
       "         thermalconductivity    : 2.4             -- ice thermal conductivity [W/m/K]\n",
       "         temperateiceconducti...: 0.24            -- temperate ice thermal conductivity [W/m/K]\n",
       "         effectiveconductivit...: 1               -- computation of effectiveconductivity: (0) arithmetic mean, (1) harmonic mean, (2) geometric mean (default)\n",
       "         meltingpoint           : 273.15          -- melting point of ice at 1atm in K\n",
       "         latentheat             : 334000.0        -- latent heat of fusion [J/m^3]\n",
       "         beta                   : 9.8e-08         -- rate of change of melting point with pressure [K/Pa]\n",
       "         mixed_layer_capacity   : 3974.0          -- mixed layer capacity [W/kg/K]\n",
       "         thermal_exchange_vel...: 0.0001          -- thermal exchange velocity [m/s]\n",
       "         rheology_B             : 210000000.0     -- flow law parameter [Pa s^(1/n)]\n",
       "         rheology_n             : 3               -- Glen's flow law exponent\n",
       "         rheology_law           : 'Paterson'      -- law for the temperature dependance of the rheology: 'None', 'BuddJacka', 'Cuffey', 'CuffeyTemperate', 'Paterson', 'Arrhenius', 'LliboutryDuval', 'NyeCO2', or 'NyeH2O'\n",
       "         earth_density          : 5512            -- Mantle density [kg/m^-3]"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "### ice temperature ###\n",
    "md.materials.rheology_B"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "elect-astrology",
   "metadata": {},
   "outputs": [],
   "source": [
    "### calving forcings ###\n",
    "# We set the calving model (no moving front ... spclevelset is actually ignored)\n",
    "from calvingvonmises import calvingvonmises\n",
    "md.calving = calvingvonmises()\n",
    "md.calving.stress_threshold_groundedice = md.calving.stress_threshold_groundedice * 0.5"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "id": "olympic-novelty",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([[0.53037254, 0.32026347, 0.10739125, 0.6831801 , 0.407963  ,\n",
       "        0.92431733, 0.59452998, 0.30512883, 0.24713031, 0.45273503,\n",
       "        0.31121587, 0.06556057, 0.89326978, 0.11442491, 0.81646947,\n",
       "        0.82494126, 0.13620279, 0.145442  , 0.41739796, 0.55705937,\n",
       "        0.9453419 , 0.88678903, 0.2987356 , 0.42432271, 0.77412763,\n",
       "        0.99764246, 0.52582554, 0.56071007, 0.65892112, 0.19030783,\n",
       "        0.80748396, 0.83758794, 0.50953111, 0.12054134, 0.61596857,\n",
       "        0.05427841, 0.60330633, 0.01090319, 0.93562654, 0.90479041,\n",
       "        0.33104554, 0.3861249 , 0.57737953, 0.04767278, 0.49905373,\n",
       "        0.09409471, 0.97953073, 0.62865429, 0.39296502, 0.35159283,\n",
       "        0.96435389, 0.74401899, 0.58575719, 0.91583386, 0.43749987,\n",
       "        0.1550989 , 0.73104823, 0.47858291, 0.46309765, 0.16971297,\n",
       "        0.2111786 , 0.02958251, 0.84091012, 0.87231036, 0.6304658 ,\n",
       "        0.95705113, 0.66850299, 0.22434071, 0.67928522, 0.75590299,\n",
       "        0.20630007, 0.54845419, 0.032968  , 0.51354266, 0.18905952,\n",
       "        0.69900963, 0.44976483, 0.23901838, 0.28127372, 0.85823724,\n",
       "        0.70409838, 0.79990801, 0.26472281, 0.64553836, 0.25529775,\n",
       "        0.00321737, 0.27075022, 0.34895651, 0.76802676, 0.86139527,\n",
       "        0.07776768, 0.78625847, 0.48409144, 0.72254021, 0.08859729,\n",
       "        0.17700188, 0.3654774 , 0.71609569, 0.98091066, 0.3747392 ],\n",
       "       [0.561937  , 0.11208384, 0.05661304, 0.46911256, 0.43403757,\n",
       "        0.89207749, 0.02648747, 0.40871187, 0.03748319, 0.27868757,\n",
       "        0.77536555, 0.74751876, 0.01140065, 0.37348275, 0.44389549,\n",
       "        0.50511296, 0.19929076, 0.24874297, 0.95134525, 0.6417279 ,\n",
       "        0.42357641, 0.06290744, 0.04614771, 0.35902776, 0.79962292,\n",
       "        0.07653603, 0.99838159, 0.20525202, 0.15208327, 0.59524688,\n",
       "        0.30555083, 0.48742458, 0.58219883, 0.8836831 , 0.84364027,\n",
       "        0.9295944 , 0.65443945, 0.82755094, 0.25750811, 0.98265544,\n",
       "        0.1757187 , 0.32894303, 0.47475198, 0.61945118, 0.78214458,\n",
       "        0.91074159, 0.93891454, 0.51847681, 0.6999259 , 0.70679393,\n",
       "        0.26936688, 0.60826949, 0.62844426, 0.29569142, 0.16382365,\n",
       "        0.87934558, 0.71432117, 0.68538819, 0.13446114, 0.67564818,\n",
       "        0.31470505, 0.83047968, 0.00882048, 0.49444394, 0.18588944,\n",
       "        0.81139188, 0.57003388, 0.6624529 , 0.63704586, 0.34099195,\n",
       "        0.52908887, 0.39215023, 0.10204594, 0.23181508, 0.72013982,\n",
       "        0.36771805, 0.38875035, 0.54614163, 0.731882  , 0.90285278,\n",
       "        0.94475927, 0.09954346, 0.96247351, 0.22731646, 0.86439448,\n",
       "        0.53297847, 0.85455916, 0.28071207, 0.41798687, 0.12243302,\n",
       "        0.33725412, 0.45068348, 0.55955099, 0.97193015, 0.80882742,\n",
       "        0.76809473, 0.2196786 , 0.08477891, 0.75727397, 0.14370711],\n",
       "       [0.02184368, 0.01443076, 0.16165749, 0.92982045, 0.12096246,\n",
       "        0.66589212, 0.00626464, 0.9953808 , 0.17662186, 0.86116181,\n",
       "        0.38555019, 0.03287213, 0.13940176, 0.27762072, 0.22972514,\n",
       "        0.77037556, 0.94285133, 0.43195423, 0.70011155, 0.96242326,\n",
       "        0.48678776, 0.34231089, 0.06509003, 0.93887818, 0.18693758,\n",
       "        0.44394629, 0.85372469, 0.91455736, 0.82177041, 0.58634001,\n",
       "        0.51364153, 0.98852365, 0.14756553, 0.90610513, 0.49364202,\n",
       "        0.11209932, 0.54776152, 0.81486409, 0.5022475 , 0.71195175,\n",
       "        0.46165401, 0.72777161, 0.79228987, 0.53526642, 0.65451546,\n",
       "        0.60579611, 0.80877155, 0.26867747, 0.40451076, 0.52015366,\n",
       "        0.30426444, 0.67606194, 0.15557765, 0.68651336, 0.37154076,\n",
       "        0.08646535, 0.47182652, 0.74750524, 0.97339088, 0.28709895,\n",
       "        0.42047866, 0.9540786 , 0.89705035, 0.39226871, 0.05795287,\n",
       "        0.5563078 , 0.83587562, 0.23353549, 0.09482451, 0.21597349,\n",
       "        0.78472462, 0.59698313, 0.04433697, 0.62202273, 0.84213934,\n",
       "        0.10005804, 0.32300958, 0.73748918, 0.76794163, 0.64202468,\n",
       "        0.45959044, 0.87502403, 0.56916979, 0.75706765, 0.33722433,\n",
       "        0.35085044, 0.36219991, 0.29008815, 0.24716128, 0.25344438,\n",
       "        0.88786653, 0.63338936, 0.20819334, 0.57130698, 0.41173373,\n",
       "        0.69340877, 0.07408487, 0.19348703, 0.61763647, 0.31096951]])"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Set Latin Hypercube sampling parameters\n",
    "k = lhsmdu.sample(3, 100) # Latin Hypercube Sampling with multi-dimensional uniformity, \n",
    "                        # 3 variables with 10 samples each\n",
    "k = np.array(k)\n",
    "k"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "id": "ee27d845-6df4-491c-92c8-a03eeb879cdd",
   "metadata": {},
   "outputs": [],
   "source": [
    "k[0] = np.interp(k[0], (k.min(), k.max()), (+1, +20))\n",
    "k[1] = np.interp(k[1], (k.min(), k.max()), (-15, -1))\n",
    "k[2] = np.interp(k[2], (k.min(), k.max()), (-0.4, +0.2))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "id": "focal-guatemala",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "{'A000': {'basal': 11.06461846505946,\n",
       "  'ice_temp': -14.610632648514423,\n",
       "  'calve': -0.14239633532102136},\n",
       " 'A001': {'basal': 7.053147637341737,\n",
       "  'ice_temp': -14.925851016542094,\n",
       "  'calve': -0.14252347176810665},\n",
       " 'A002': {'basal': 2.988921635135776,\n",
       "  'ice_temp': -14.964720175353104,\n",
       "  'calve': -0.13999843766960296},\n",
       " 'A003': {'basal': 13.982070221507596,\n",
       "  'ice_temp': -14.67567602519086,\n",
       "  'calve': -0.1268239438366477},\n",
       " 'A004': {'basal': 8.727535648744237,\n",
       "  'ice_temp': -14.700253562140894,\n",
       "  'calve': -0.1406963837952911},\n",
       " 'A005': {'basal': 18.585940912033244,\n",
       "  'ice_temp': -14.379298617215296,\n",
       "  'calve': -0.13135048617401351},\n",
       " 'A006': {'basal': 12.289533283379336,\n",
       "  'ice_temp': -14.985829581277255,\n",
       "  'calve': -0.14266352604560384},\n",
       " 'A007': {'basal': 6.764192013370105,\n",
       "  'ice_temp': -14.717999629191503,\n",
       "  'calve': -0.1256995412481035},\n",
       " 'A008': {'basal': 5.656865510146335,\n",
       "  'ice_temp': -14.978124728348485,\n",
       "  'calve': -0.13974178902280743},\n",
       " 'A009': {'basal': 9.582337791557848,\n",
       "  'ice_temp': -14.809109462698375,\n",
       "  'calve': -0.12800148409949402},\n",
       " 'A010': {'basal': 6.880407909612024,\n",
       "  'ice_temp': -14.461080300491247,\n",
       "  'calve': -0.13615853246803172},\n",
       " 'A011': {'basal': 2.1902767725827914,\n",
       "  'ice_temp': -14.48059292939334,\n",
       "  'calve': -0.1422071903518407},\n",
       " 'A012': {'basal': 17.99317091966237,\n",
       "  'ice_temp': -14.996401128458075,\n",
       "  'calve': -0.14038013789219533},\n",
       " 'A013': {'basal': 3.1232106900724776,\n",
       "  'ice_temp': -14.742685166402106,\n",
       "  'calve': -0.1380095929866884},\n",
       " 'A014': {'basal': 16.526874360772112,\n",
       "  'ice_temp': -14.693345976655733,\n",
       "  'calve': -0.1388310331631077},\n",
       " 'A015': {'basal': 16.68862054569352,\n",
       "  'ice_temp': -14.650450047222797,\n",
       "  'calve': -0.1295585272274602},\n",
       " 'A016': {'basal': 3.5390009854058198,\n",
       "  'ice_temp': -14.864743910446679,\n",
       "  'calve': -0.1266004556827141},\n",
       " 'A017': {'basal': 3.715399032403361,\n",
       "  'ice_temp': -14.830092064948872,\n",
       "  'calve': -0.13536267301894556},\n",
       " 'A018': {'basal': 8.907671018788983,\n",
       "  'ice_temp': -14.337768872219936,\n",
       "  'calve': -0.13076360073778975},\n",
       " 'A019': {'basal': 11.57413207998271,\n",
       "  'ice_temp': -14.554722060005835,\n",
       "  'calve': -0.12626478444209088},\n",
       " 'A020': {'basal': 18.987348874856615,\n",
       "  'ice_temp': -14.70758384282586,\n",
       "  'calve': -0.1344222422816922},\n",
       " 'A021': {'basal': 17.869438402188496,\n",
       "  'ice_temp': -14.960309602434634,\n",
       "  'calve': -0.13690011449078066},\n",
       " 'A022': {'basal': 6.6421303803306495,\n",
       "  'ice_temp': -14.972053377652367,\n",
       "  'calve': -0.1416546322823628},\n",
       " 'A023': {'basal': 9.039880588969368,\n",
       "  'ice_temp': -14.752813975839567,\n",
       "  'calve': -0.12666859779721162},\n",
       " 'A024': {'basal': 15.718470206757152,\n",
       "  'ice_temp': -14.444082824977194,\n",
       "  'calve': -0.13956486771592697},\n",
       " 'A025': {'basal': 19.98588817805362,\n",
       "  'ice_temp': -14.950759859329965,\n",
       "  'calve': -0.13515700131487124},\n",
       " 'A026': {'basal': 10.977805649251234,\n",
       "  'ice_temp': -14.304809856852417,\n",
       "  'calve': -0.12812903551186972},\n",
       " 'A027': {'basal': 11.643832466475223,\n",
       "  'ice_temp': -14.860566773803871,\n",
       "  'calve': -0.12708571571498983},\n",
       " 'A028': {'basal': 13.518909920046104,\n",
       "  'ice_temp': -14.897822857933798,\n",
       "  'calve': -0.12867707209097357},\n",
       " 'A029': {'basal': 4.571992115657638,\n",
       "  'ice_temp': -14.587291951639875,\n",
       "  'calve': -0.13271485654427734},\n",
       " 'A030': {'basal': 16.355320195963067,\n",
       "  'ice_temp': -14.790286007332332,\n",
       "  'calve': -0.13396168270537023},\n",
       " 'A031': {'basal': 16.93007516298551,\n",
       "  'ice_temp': -14.662844542293987,\n",
       "  'calve': -0.12581714568793595},\n",
       " 'A032': {'basal': 10.66670712291882,\n",
       "  'ice_temp': -14.596434904647618,\n",
       "  'calve': -0.1402401240012065},\n",
       " 'A033': {'basal': 3.23998758524581,\n",
       "  'ice_temp': -14.38518068223505,\n",
       "  'calve': -0.12723067689153117},\n",
       " 'A034': {'basal': 12.698845792313001,\n",
       "  'ice_temp': -14.41323924933937,\n",
       "  'calve': -0.13430468722792466},\n",
       " 'A035': {'basal': 1.9748739470092094,\n",
       "  'ice_temp': -14.353010001509219,\n",
       "  'calve': -0.1408483925184777},\n",
       " 'A036': {'basal': 12.457094193053301,\n",
       "  'ice_temp': -14.545814898025013,\n",
       "  'calve': -0.13337650269869356},\n",
       " 'A037': {'basal': 1.1467400933406036,\n",
       "  'ice_temp': -14.424513270084702,\n",
       "  'calve': -0.1287955200252235},\n",
       " 'A038': {'basal': 18.801860095539475,\n",
       "  'ice_temp': -14.823950203731181,\n",
       "  'calve': -0.13415709768590117},\n",
       " 'A039': {'basal': 18.213126565861888,\n",
       "  'ice_temp': -14.315829388921598,\n",
       "  'calve': -0.13056053360951891},\n",
       " 'A040': {'basal': 7.259002459552609,\n",
       "  'ice_temp': -14.881261185322828,\n",
       "  'calve': -0.13485330243913163},\n",
       " 'A041': {'basal': 8.310595389105302,\n",
       "  'ice_temp': -14.773894765514592,\n",
       "  'calve': -0.13028921263142973},\n",
       " 'A042': {'basal': 11.962091300238171,\n",
       "  'ice_temp': -14.671724408219145,\n",
       "  'calve': -0.1291826825516788},\n",
       " 'A043': {'basal': 1.8487572433734991,\n",
       "  'ice_temp': -14.570331668446103,\n",
       "  'calve': -0.1335908017142462},\n",
       " 'A044': {'basal': 10.466669585304622,\n",
       "  'ice_temp': -14.456330136006104,\n",
       "  'calve': -0.1315456031759913},\n",
       " 'A045': {'basal': 2.735059769642265,\n",
       "  'ice_temp': -14.366220424021613,\n",
       "  'calve': -0.1323811717645022},\n",
       " 'A046': {'basal': 19.640093148032786,\n",
       "  'ice_temp': -14.346479247885231,\n",
       "  'calve': -0.12890001107909366},\n",
       " 'A047': {'basal': 12.941045801273514,\n",
       "  'ice_temp': -14.641085808558822,\n",
       "  'calve': -0.13816297540898703},\n",
       " 'A048': {'basal': 8.4411892039855,\n",
       "  'ice_temp': -14.513941910682384,\n",
       "  'calve': -0.13583334627842802},\n",
       " 'A049': {'basal': 7.6512979435831205,\n",
       "  'ice_temp': -14.509129382282392,\n",
       "  'calve': -0.13384999537039644},\n",
       " 'A050': {'basal': 19.350332087012305,\n",
       "  'ice_temp': -14.815640601900459,\n",
       "  'calve': -0.13755263585003813},\n",
       " 'A051': {'basal': 15.14362615908376,\n",
       "  'ice_temp': -14.578166830465381,\n",
       "  'calve': -0.13117606718478142},\n",
       " 'A052': {'basal': 12.12204035356174,\n",
       "  'ice_temp': -14.564030086468849,\n",
       "  'calve': -0.14010271092011178},\n",
       " 'A053': {'basal': 18.423971812993464,\n",
       "  'ice_temp': -14.797194633216337,\n",
       "  'calve': -0.13099681844318622},\n",
       " 'A054': {'basal': 9.291463168898611,\n",
       "  'ice_temp': -14.889596212881003,\n",
       "  'calve': -0.1363988032598482},\n",
       " 'A055': {'basal': 3.899771707539529,\n",
       "  'ice_temp': -14.388220044974595,\n",
       "  'calve': -0.14128803171458676},\n",
       " 'A056': {'basal': 14.89598416457146,\n",
       "  'ice_temp': -14.503854942749427,\n",
       "  'calve': -0.13467883733175584},\n",
       " 'A057': {'basal': 10.07583403924337,\n",
       "  'ice_temp': -14.524128682761116,\n",
       "  'calve': -0.12995076809138922},\n",
       " 'A058': {'basal': 9.780184259508786,\n",
       "  'ice_temp': -14.910170926908927,\n",
       "  'calve': -0.12607668265262334},\n",
       " 'A059': {'basal': 4.178788397689978,\n",
       "  'ice_temp': -14.53095364356135,\n",
       "  'calve': -0.13784703504966694},\n",
       " 'A060': {'basal': 4.970463666956434,\n",
       "  'ice_temp': -14.78387151475672,\n",
       "  'calve': -0.13555948649892785},\n",
       " 'A061': {'basal': 1.503371881866497,\n",
       "  'ice_temp': -14.422461064111035,\n",
       "  'calve': -0.1264079007326236},\n",
       " 'A062': {'basal': 16.993503243312844,\n",
       "  'ice_temp': -14.998209085059347,\n",
       "  'calve': -0.12738597236221638},\n",
       " 'A063': {'basal': 17.593006933139897,\n",
       "  'ice_temp': -14.65792597736387,\n",
       "  'calve': -0.13604330548582777},\n",
       " 'A064': {'basal': 12.97563165281408,\n",
       "  'ice_temp': -14.874134406682707,\n",
       "  'calve': -0.1417770392826609},\n",
       " 'A065': {'basal': 19.210905362884564,\n",
       "  'ice_temp': -14.435836149246482,\n",
       "  'calve': -0.13322992846596632},\n",
       " 'A066': {'basal': 13.701850054945162,\n",
       "  'ice_temp': -14.604959055240583,\n",
       "  'calve': -0.1284351586889244},\n",
       " 'A067': {'basal': 5.221758960174272,\n",
       "  'ice_temp': -14.54019976126117,\n",
       "  'calve': -0.13876568329242},\n",
       " 'A068': {'basal': 13.90770797250201,\n",
       "  'ice_temp': -14.558002827865915,\n",
       "  'calve': -0.14114466662524688},\n",
       " 'A069': {'basal': 15.370519433776364,\n",
       "  'ice_temp': -14.765451922293964,\n",
       "  'calve': -0.1390668829272172},\n",
       " 'A070': {'basal': 4.877321185653662,\n",
       "  'ice_temp': -14.633649788873784,\n",
       "  'calve': -0.12931243151716365},\n",
       " 'A071': {'basal': 11.409839287422072,\n",
       "  'ice_temp': -14.729604603945946,\n",
       "  'calve': -0.13253232002515208},\n",
       " 'A072': {'basal': 1.5680087677513148,\n",
       "  'ice_temp': -14.932884713087104,\n",
       "  'calve': -0.14201056083414465},\n",
       " 'A073': {'basal': 10.7432968748266,\n",
       "  'ice_temp': -14.841953670015432,\n",
       "  'calve': -0.13210287463932963},\n",
       " 'A074': {'basal': 4.548158938920953,\n",
       "  'ice_temp': -14.499777737036881,\n",
       "  'calve': -0.12832773181906287},\n",
       " 'A075': {'basal': 14.284292923821198,\n",
       "  'ice_temp': -14.74672456586348,\n",
       "  'calve': -0.1410549083054265},\n",
       " 'A076': {'basal': 9.525629746997925,\n",
       "  'ice_temp': -14.731986947399086,\n",
       "  'calve': -0.1372311445342989},\n",
       " 'A077': {'basal': 5.50198986929816,\n",
       "  'ice_temp': -14.621700688019098,\n",
       "  'calve': -0.130122550086966},\n",
       " 'A078': {'basal': 6.3087425531300765,\n",
       "  'ice_temp': -14.491549826371088,\n",
       "  'calve': -0.12960027077037434},\n",
       " 'A079': {'basal': 17.32431836294291,\n",
       "  'ice_temp': -14.371748223632967,\n",
       "  'calve': -0.1317598283096827},\n",
       " 'A080': {'basal': 14.381448995636955,\n",
       "  'ice_temp': -14.342383765543982,\n",
       "  'calve': -0.13488869398965436},\n",
       " 'A081': {'basal': 16.210677686808147,\n",
       "  'ice_temp': -14.934638233252501,\n",
       "  'calve': -0.12776373801594482},\n",
       " 'A082': {'basal': 5.992747165874494,\n",
       "  'ice_temp': -14.329971146560792,\n",
       "  'calve': -0.1330093369772839},\n",
       " 'A083': {'basal': 13.2634019126666,\n",
       "  'ice_temp': -14.845105916000133,\n",
       "  'calve': -0.1297867664947827},\n",
       " 'A084': {'basal': 5.8128008359964,\n",
       "  'ice_temp': -14.398696487560274,\n",
       "  'calve': -0.13698735221401587},\n",
       " 'A085': {'basal': 1.0,\n",
       "  'ice_temp': -14.630924292719733,\n",
       "  'calve': -0.13675365571105297},\n",
       " 'A086': {'basal': 6.107824467859663,\n",
       "  'ice_temp': -14.405588235035179,\n",
       "  'calve': -0.1365590048391917},\n",
       " 'A087': {'basal': 7.600964498884227,\n",
       "  'ice_temp': -14.807690870995453,\n",
       "  'calve': -0.13779576843026586},\n",
       " 'A088': {'basal': 15.601990356584823,\n",
       "  'ice_temp': -14.711500505462546,\n",
       "  'calve': -0.13853199217682122},\n",
       " 'A089': {'basal': 17.38461240550493,\n",
       "  'ice_temp': -14.9185992004821,\n",
       "  'calve': -0.13842423281806582},\n",
       " 'A090': {'basal': 2.423338887165583,\n",
       "  'ice_temp': -14.768071067587934,\n",
       "  'calve': -0.1275434808437324},\n",
       " 'A091': {'basal': 15.950076084157008,\n",
       "  'ice_temp': -14.688589536235948,\n",
       "  'calve': -0.13190792954247954},\n",
       " 'A092': {'basal': 10.181004586641913,\n",
       "  'ice_temp': -14.612304564056814,\n",
       "  'calve': -0.1392003176491986},\n",
       " 'A093': {'basal': 14.733546367951384,\n",
       "  'ice_temp': -14.32334475014043,\n",
       "  'calve': -0.13297268269351792},\n",
       " 'A094': {'basal': 2.630101241377056,\n",
       "  'ice_temp': -14.43763310010725,\n",
       "  'calve': -0.13570946763502606},\n",
       " 'A095': {'basal': 4.317950607087835,\n",
       "  'ice_temp': -14.466175061380715,\n",
       "  'calve': -0.1308785577783022},\n",
       " 'A096': {'basal': 7.9163866595419785,\n",
       "  'ice_temp': -14.850457870553528,\n",
       "  'calve': -0.1415003649414105},\n",
       " 'A097': {'basal': 14.610505372088026,\n",
       "  'ice_temp': -14.944983960752632,\n",
       "  'calve': -0.13945254041151178},\n",
       " 'A098': {'basal': 19.666439300743352,\n",
       "  'ice_temp': -14.47375731776849,\n",
       "  'calve': -0.13217810187402196},\n",
       " 'A099': {'basal': 8.093215938503246,\n",
       "  'ice_temp': -14.903692151734218,\n",
       "  'calve': -0.137437639442064}}"
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Create 10 ensembles with parameters from [0], [1], [2], ..., [9]\n",
    "lhs_params_for_ens = {\n",
    "    'basal' : k[0],\n",
    "    'ice_temp' : k[1],\n",
    "    'calve' : k[2]\n",
    "} \n",
    "lhs_params_for_ens\n",
    "\n",
    "# ens needs to match LHS sample number\n",
    "ens = list(range(0, 100))\n",
    "\n",
    "ens_dict = {}\n",
    "for i in range(len(ens)):\n",
    "    trial_name = 'A{:03d}'.format(i)\n",
    "    res = dict()\n",
    "    for sub in lhs_params_for_ens:\n",
    "        res[sub] = lhs_params_for_ens[sub][i]\n",
    "    ens_dict[trial_name] = res\n",
    "ens_dict"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "id": "aging-archives",
   "metadata": {},
   "outputs": [],
   "source": [
    "# save output\n",
    "with open('ens_dict.pickle', 'wb') as f:\n",
    "    # Pickle the 'data' dictionary \n",
    "    pickle.dump(ens_dict, f)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "controlling-hungary",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.9.13"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}