# Millenium eruption

## Simulations with gravity current ON

| ID | Meteo data | Column height |
| -------- | ----------------------- | ----- |
| ->**plinian-w1**<- | 1996-11-29 / 1996-12-04 | 30 km |
| plinian-w2 | 1980-01-30 / 1980-02-04 | 44 km |
| coign-w1 | 1976-06-08 / 1976-06-13 | 50 km |
| coign-w2 | 1962-10-23 / 1962-10-28 | 20 km |

## Datasets
<details>
<summary>plinian-ds</summary>

label,latitude,longitude,thickness
P1,35.583,135.883,0.0009666666667
P2,42.283,126.600,0.000004141176471
P3,42.306,126.283,0.0001
P4,40.530,132.200,0.033
P5,41.350,134.530,0.1088
P6,41.940,138.720,0.0174
P7,41.650,140.010,0.0284
P8,45.007,145.001,0.000009381818182
P9,43.650,142.850,0.007333333333
P10,43.067,141.350,0.0008
P11,42.637,141.603,0.008571428571
P12,42.633,140.307,0.02702702703
P13,44.351,146.256,0.01
P14,43.078,145.101,0.01
P15,40.783,141.317,0.007692307692
P16,39.955,139.739,0.0009714285714

</details>

<details>
<summary>coign-ds</summary>

label,latitude,longitude,thickness
P1,35.583,135.883,0.00003333333333
P2,42.283,126.600,0.00002145882353
P4,40.530,132.200,0.017
P5,41.350,134.530,0.0512
P6,41.940,138.720,0.0126
P7,41.650,140.010,0.0116
P8,45.007,145.001,0.000005018181818
P9,43.650,142.850,0.002666666667
P10,43.067,141.350,0.0002
P11,42.637,141.603,0.001428571429
P12,42.633,140.307,0.02297297297
P15,40.783,141.317,0.01230769231
P16,39.955,139.739,0.001028571429

</details>

<details>
<summary>all-ds</summary>

label,latitude,longitude,thickness,load
P01,42.000,130.283,0.09,90
P02,41.333,130.417,0.1,100
P03,40.283,132.217,0.05,50
P04,40.867,132.650,0.04,40
P05,40.767,133.200,0.09,90
P06,41.350,134.435,0.16,160
P07,40.017,136.398,0.1,100
P08,40.745,136.938,0.001,1
P09,41.005,137.018,0.08,80
P10,41.873,137.657,0.11,110
P11,41.147,138.002,0.09,90
P12,41.423,138.260,0.074,74
P13,40.750,138.450,0.02,20
P14,40.952,138.667,0.11,110
P15,41.998,138.668,0.03,30
P16,40.898,139.015,0.06,60
P17,41.833,139.163,0.032,32
P18,40.895,142.278,0.002,2
P19,40.678,142.783,0.002,2
P20,42.633,141.567,0.02,20
P21,42.767,141.600,0.01,10
P22,42.466,140.550,0.01,10
P23,42.183,140.667,0.04,40
P24,41.883,141.083,0.05,50
P25,41.516,140.000,0.04,40
P26,40.300,141.333,0.05,50
P27,40.483,140.483,0.04,40
P28,46.710,151.660,0.015,15
P29,42.996,131.847,0.03,30
P30,43.393,132.795,0.02,20
P31,43.504,133.382,0.015,15
P32,43.522,134.144,0.01,10
P33,35.583,153.883,0.001,1
P34,43.078,145.101,0.01,10
P35,45.007,145.001,0.0000144,0.0144
P36,42.283,126.600,0.0000256,0.0256
P37,39.955,139.739,0.002,2
P38,42.633,140.307,0.01,10
P39,42.082,140.688,0.05,50
P40,41.998,128.539,0.72,720
P41,42.011,128.423,0.9,900
P42,42.033,128.421,0.5,500
P43,42.074,128.416,0.5,500
P44,42.113,128.376,0.2,200
P45,41.630,140.427,0.03,30
P46,41.781,141.100,0.04,40
P47,41.838,141.139,0.04,40
P48,42.301,140.271,0.02,20
P49,42.319,140.274,0.02,20.02
P50,42.484,140.353,0.01,10.02
P51,42.537,140.399,0.01,10.03
P52,42.555,140.42,0.01,10.04
P53,42.503,140.798,0.02,20.1
P54,42.395,140.911,0.01,10.06
P55,42.386,141.076,0.01,10.07
P56,42.440,141.151,0.01,10.08
P57,42.560,141.376,0.01,10.09
P58,42.617,141.782,0.01,10.1
P59,42.581,141.905,0.03,30.33
P60,42.555,141.939,0.01,10.12
P61,42.484,142.046,0.03,30.39
P62,42.473,142.118,0.02,20.28
P63,42.448,142.196,0.02,20.3
P64,42.309,142.432,0.03,30.48
P65,42.304,142.454,0.02,20.34
P66,42.266,142.513,0.02,20.36
P67,42.228,142.614,0.02,20.38
P68,42.217,142.642,0.01,10.2
P69,42.192,142.723,0.03,30.63
P70,42.137,142.86,0.04,40.88
P71,42.137,142.91,0.01,10.23
P72,41.990,143.245,0.04,40.96
P73,41.998,143.253,0.05,51.25

</details>

## Results: Deposit for each phase (2-step optimization)

| plinian-dataset | coign-dataset |
| -------- | -------- |
|  |  |

## Correction factor
Each deposit field was corrected using the correction factor:
$$K=\exp\left[\frac{1}{n}\sum_i \ln(Obs_i/Mod_i)\right]$$

| plinian-dataset | coign-dataset |
| -------- | -------- |
| $$k_1=3.74$$ | $$k_2=4.34$$ |

| plinian-dataset | coign-dataset |
| -------- | -------- |
|  |  |

## Total deposit

$$ deposit(total) = k_1 *deposit(plinian) + k_2*deposit(coign) $$