


"""Author: SKHEO, KHU"""


import matplotlib.pyplot as plt
import numpy as np
from matplotlib.font_manager import FontProperties
import matplotlib.gridspec as gridspec
import random


class influent(object):
    def buffer_tank(t1,t2):
        switch = np.random.choice(2,1)
        #print(switch)
        switch = 1

        if switch == 0:


            si_m =np.array([30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30,	30
            ])
            ss_m =np.array( [60.6354685714286,	58.76627,	58.6790042857143,	59.0787371428571,	61.3495721428571,	64.4120114285714,	68.4505514285714,	68.7160178571429,	66.2135607142857,	64.0918514285714,	59.6748057142857,	59.0760535714286,	50.9444857142857,	51.5911435714286,	47.4257928571428,	43.8942921428571,	42.0521914285714,	40.3570792857143,	37.9440521428572,	38.3863864285714,	38.73851,	38.5068207142857,	38.8947528571429,	38.5647657142857,	38.5663271428571,	38.6095528571429,	38.3977407142857,	38.2504228571429,	38.1998464285714,	38.517515,	38.7972178571429,	39.1536492857143,	39.4473192857143,	40.1436171428571,	40.7646014285714,	42.3019328571429,	47.1787142857143,	57.9955635714286,	69.5699585714286,	78.1772535714286,	92.812395,	103.705670714286,	113.199675,	118.184842142857,	117.334249285714,	110.729128571429,	101.336726428571,	93.8272685714286,	87.5917528571428,	80.7261778571429,	75.3976914285714,	70.1927371428572,	68.27595,	67.03307,	66.76507,	67.07967,	68.08952,	68.72261,	68.31479,	65.16874,	61.64516,	59.76219,	58.94655,	59.38933,	60.21274,	62.42662,	65.86786,	65.61928,	64.00353,	59.54469,	58.34064,	57.33857,	58.6902,	60.46908,	63.25004,	62.37225,	62.21689,	62.01492,	64.7641378571428,	65.1182907142857,	65.89489,	68.7446878571428,	72.4965928571428,	73.8964728571429,	76.2051228571429,	78.787575,	80.96422,	81.6603207142857,	79.9468214285714,	79.4172964285714,	76.744505,	74.60987,	72.4581742857143,	68.7893885714286,	66.4230035714286,	64.49190214
            ])
            xi_m = np.array([49.4308707142857,	48.6576828571429,	44.2617214285714,	47.5107407142857,	45.6272685714286,	46.2123042857143,	45.9966964285714,	37.4277971428571,	36.085415,	34.7306935714286,	32.8243407142857,	29.9764685714286,	26.6990357142857,	24.8555064285714,	22.4881757142857,	20.7251442857143,	19.8658014285714,	18.2974835714286,	18.4372371428571,	18.2232457142857,	17.4806371428571,	17.87371,	18.0521878571429,	17.8038764285714,	17.1174485714286,	17.2007414285714,	16.0687707142857,	16.0351614285714,	15.9654492857143,	16.0657157142857,	15.7817571428571,	16.1133721428571,	16.0612071428571,	15.8060992857143,	16.9256157142857,	17.1571664285714,	19.4147192857143,	21.4793057142857,	30.1330971428571,	36.7191007142857,	41.0205742857143,	49.4635685714286,	58.20672,	68.8762578571428,	75.8283385714286,	81.5284757142857,	86.0908792857143,	83.8933014285714,	83.3657107142857,	79.6946985714286,	76.6667735714286,	75.3878621428571,	69.5391428571428,	74.6537142857143,	70.989,	70.7332857142857,	67.836,	62.0564285714286,	57.0402857142857,	64.5568571428571,	60.4957142857143,	57.1474285714286,	58.1817142857143,	50.9847142857143,	48.5514285714286,	49.035,	47.1648571428571,	46.5202857142857,	44.6344285714286,	43.1411428571428,	41.9254285714286,	42.4255714285714,	45.1215714285714,	44.4241428571429,	45.2405714285714,	49.5551428571429,	48.1884285714286,	45.9525714285714,	50.4756271428571,	54.1749914285714,	64.5585035714286,	78.3004821428571,	73.5002542857143,	70.7127142857143,	61.7074421428571,	52.7604064285714,	63.6332557142857,	58.4051,	57.5799785714286,	58.665235,	50.398865,	65.6329171428571,	61.0644014285714,	61.726605,	58.548715,	53.98449214
            ])
            xs_m = np.array([216.725581428571,	212.693755,	201.949778571429,	199.152942857143,	198.874110714286,	198.924661428571,	200.957670714286,	200.463356428571,	196.688628571429,	188.630685714286,	182.280029285714,	176.809364285714,	171.548662142857,	164.180526428571,	156.899692857143,	149.356193571429,	142.01543,	135.406609285714,	129.827009285714,	125.17339,	121.874577857143,	119.186505,	117.272132857143,	116.175762857143,	115.350113571429,	114.426378571429,	113.10468,	111.670890714286,	110.121382857143,	109.188875714286,	108.466118571429,	107.993685714286,	107.632940714286,	107.460777142857,	107.729895,	108.390417857143,	109.337793571429,	112.799132142857,	119.018597857143,	129.590735714286,	144.404012142857,	158.307502142857,	172.399188571429,	186.794237857143,	203.212808571429,	218.839162142857,	232.302918571429,	246.42998,	255.36492,	260.716212857143,	262.784623571429,	261.943842142857,	262.4,	256.213857142857,	254.890428571429,	253.932857142857,	252.894714285714,	252.568428571429,	248.933857142857,	239.285285714286,	237.559571428571,	237.319,	233.425714285714,	232.687142857143,	228.195571428571,	221.154142857143,	216.265285714286,	212.509,	209.343428571429,	206.297285714286,	202.508714285714,	199.170142857143,	196.757714285714,	197.271,	198.217714285714,	197.498285714286,	200.136857142857,	202.351,	208.636487142857,	213.972012857143,	216.035177857143,	217.007077857143,	217.190538571429,	219.88186,	216.422905714286,	210.827879285714,	206.039118571429,	211.184556428571,	215.944957857143,	220.015212857143,	221.504147857143,	213.589959285714,	216.737107857143,	219.749815,	221.829720714286,	222.8158679
            ])
            xbh_m = np.array([29.5727378571429,	29.0391064285714,	27.3568278571429,	27.4070085714286,	27.1666078571429,	27.2373778571429,	27.4392157142857,	26.43243,	25.8637892857143,	24.8177992857143,	23.9004921428571,	22.9760521428571,	22.0274885714286,	21.0041392857143,	19.932075,	18.8978714285714,	17.9867792857143,	17.0781171428571,	16.4738314285714,	15.9328485714286,	15.4839492857143,	15.2291292857143,	15.0362,	14.8866371428571,	14.7184807142857,	14.6249764285714,	14.3525107142857,	14.1897885714286,	14.0098178571429,	13.9171285714286,	13.8051657142857,	13.7896192857143,	13.7435721428571,	13.6963578571429,	13.8503442857143,	13.9498392857143,	14.3057271428571,	14.9196342857143,	16.5725028571429,	18.4787778571429,	20.60287,	23.0857221428571,	25.6227992857143,	28.4079457142857,	31.0044892857143,	33.3741064285714,	35.376935,	36.702555,	37.6368121428571,	37.8234507142857,	37.7170285714286,	37.4812835714286,	36.8821428571429,	36.7632857142857,	36.2087142857143,	36.0738571428571,	35.6367142857143,	34.9584285714286,	33.9971428571429,	33.7601428571429,	33.1175714285714,	32.7185714285714,	32.4008571428571,	31.519,	30.7497142857143,	30.0212857142857,	29.27,	28.781,	28.2197142857143,	27.7152857142857,	27.1594285714286,	26.8441428571429,	26.8755714285714,	26.8548571428571,	27.0508571428571,	27.4502857142857,	27.5914285714286,	27.5892857142857,	28.7903142857143,	29.7940671428571,	31.1771042857143,	32.8120114285714,	32.2989928571429,	32.2880685714286,	30.9035092857143,	29.2877142857143,	29.9635492857143,	29.9543714285714,	30.3915278571429,	30.9642664285714,	30.2114871428571,	31.0244807142857,	30.8668485714286,	31.2752228571429,	31.1530642857143,	30.75556857
            ])
            xba_m = np.array([0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0
            ])
            xp_m = np.array([0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0
            ])
            so_m = np.array([0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0
            ])
            sno_m = np.array([0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0
            ])
            snh_m = np.array([28.8220571428571,	28.7887985714286,	29.5180928571429,	30.1037135714286,	31.5028664285714,	33.3389335714286,	35.9290307142857,	35.9040571428571,	34.9745692857143,	34.6332135714286,	32.656025,	32.8984735714286,	30.3977371428571,	28.3304792857143,	25.6398392857143,	23.3915521428571,	22.1718892857143,	21.0394971428571,	19.2716914285714,	19.523765,	19.680425,	19.7460357142857,	19.6819092857143,	19.6431385714286,	19.5588578571429,	19.4719685714286,	19.4346278571429,	19.2216121428571,	19.1036192857143,	19.0949664285714,	19.2159414285714,	19.5348928571429,	19.9514428571429,	20.2981207142857,	20.7682164285714,	21.3070535714286,	23.8743442857143,	28.7198535714286,	33.3418528571428,	37.0684792857143,	42.4418878571429,	45.7022585714286,	48.3885121428571,	49.2385892857143,	48.570855,	45.1777735714286,	40.4114407142857,	38.6294614285714,	36.5216007142857,	34.011615,	32.001965,	30.7008664285714,	29.79539,	29.23881,	28.52221,	27.95867,	28.12564,	27.55166,	27.4647,	26.45589,	26.1776,	26.28891,	26.78984,	27.29076,	27.12379,	27.42295,	28.33436,	28.2787,	27.7082,	26.68548,	26.57068,	26.62286,	27.4647,	28.33436,	29.26664,	30.63027,	30.58853,	28.16043,	29.0562521428571,	28.9043385714286,	29.1635307142857,	30.7279942857143,	32.366355,	32.7331092857143,	33.9136007142857,	34.3976907142857,	34.7011935714286,	34.374275,	34.2266692857143,	34.5052828571429,	33.8789728571428,	33.5312385714286,	32.7141457142857,	31.3006457142857,	30.090835,	29.31547357
            ])
            snd_m = np.array([6.06355142857143,	5.87662428571429,	5.86789785714286,	5.90787571428571,	6.13496071428572,	6.4412,	6.84505142857143,	6.87160285714286,	6.62135571428571,	6.40919,	5.96747785714286,	5.90760642857143,	5.09445142857143,	5.15911785714286,	4.74258,	4.38943214285714,	4.20522,	4.03570714285714,	3.79440214285714,	3.83863571428572,	3.87384857142857,	3.85068142857143,	3.88947642857143,	3.85647571428571,	3.85663714285714,	3.86095785714286,	3.83977857142857,	3.82504714285714,	3.81998214285714,	3.85175,	3.87971928571429,	3.91536928571429,	3.94472928571429,	4.01436642857143,	4.076465,	4.23019428571429,	4.71787071428572,	5.79955857142857,	6.95699214285714,	7.81772785714286,	9.28124214285714,	10.3705671428571,	11.3199671428571,	11.8184842857143,	11.733425,	11.0729128571429,	10.1336728571429,	9.38272428571429,	8.75917428571428,	8.07261785714286,	7.53976928571428,	7.01927214285714,	6.8276,	6.70331,	6.67651,	6.70797,	6.80895,	6.87226,	6.83148,	6.51687,	6.16452,	5.97622,	5.89466,	5.93893,	6.02127,	6.24266,	6.58679,	6.56193,	6.40035,	5.95447,	5.83406,	5.73386,	5.86902,	6.04691,	6.325,	6.23723,	6.22169,	6.20149,	6.47641214285714,	6.51182642857143,	6.58949357142857,	6.87447071428572,	7.24965714285714,	7.38964571428571,	7.62051714285715,	7.87875928571429,	8.09642142857143,	8.16602928571429,	7.99468214285714,	7.94172571428571,	7.67444928571429,	7.46099285714286,	7.24581714285714,	6.87894357142857,	6.64230142857143,	6.449195
            ])
            xnd_m = np.array([11.1176378571429,	10.9169378571429,	10.2845228571429,	10.3034121428571,	10.2130814285714,	10.2396371428571,	10.3155478571429,	9.93697571428571,	9.72322642857143,	9.33004357142857,	8.98514571428571,	8.63766571428571,	8.28102,	7.89624785714286,	7.49323214285714,	7.10448214285714,	6.76195428571429,	6.420385,	6.19316285714286,	5.98983071428572,	5.82102071428571,	5.72516142857143,	5.65265142857143,	5.59647571428571,	5.53331,	5.49819785714286,	5.39571357142857,	5.33442857142857,	5.26679142857143,	5.23202142857143,	5.18996428571428,	5.18408357142857,	5.16683142857143,	5.14899285714286,	5.20698642857143,	5.24426428571429,	5.37812785714286,	5.60895071428571,	6.23023142857143,	6.94693928571429,	7.74539285714286,	8.67882714285714,	9.63265928571428,	10.6796421428571,	11.65585,	12.5466828571429,	13.2996528571429,	13.7979642857143,	14.1491485714286,	14.2193378571429,	14.1792635714286,	14.090715,	13.86546,	13.8207085714286,	13.6123371428571,	13.5616542857143,	13.3972728571429,	13.1422285714286,	12.7808742857143,	12.6918157142857,	12.4501071428571,	12.3001871428571,	12.1807642857143,	11.8492814285714,	11.5600285714286,	11.2861071428571,	11.0037657142857,	10.8199342857143,	10.6089314285714,	10.4193128571429,	10.2102871428571,	10.0917228571429,	10.1035657142857,	10.0958657142857,	10.1695157142857,	10.31969,	10.37281,	10.3719114285714,	10.8233992857143,	11.2007928571429,	11.7207078571429,	12.3353228571429,	12.1424714285714,	12.1384457142857,	11.6178142857143,	11.0103757142857,	11.26451,	11.2610514285714,	11.4254307142857,	11.64078,	11.3576864285714,	11.6634342857143,	11.6040742857143,	11.75758,	11.7117114285714,	11.56225357
            ])
            salk_m = np.array([7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7,	7
            ])
            q_m = np.array([21622.2142857143,	21712.8571428571,	21056.5714285714,	22513.3571428571,	21852.2857142857,	22043,	21792.2857142857,	18642.7142857143,	18018.2142857143,	18009.6428571429,	17610.7142857143,	16626.7857142857,	15364.1428571429,	14909.3571428571,	14131.5,	13646.0714285714,	13702.7142857143,	13197.7857142857,	13751.2857142857,	14020.3571428571,	13753,	14274.1428571429,	14563.0714285714,	14430.1428571429,	13942.5,	14056.1428571429,	13239.3571428571,	13316.4285714286,	13381.1428571429,	13514.4285714286,	13304.3571428571,	13579.5714285714,	13515.5,	13255.3571428571,	14099.5714285714,	14131.9285714286,	15696.7857142857,	16761.2857142857,	21770.8571428571,	23846.2142857143,	23835.2857142857,	25733.4285714286,	27343.8571428571,	29255.5714285714,	29436.3571428571,	29300.1428571429,	29092.5,	27096.7857142857,	26090.9285714286,	24624.9285714286,	23587.5714285714,	23246.6428571429,	21537.5714285714,	23435.5714285714,	22513.7142857143,	22524,	21785.7142857143,	20117.5714285714,	18841.7142857143,	21912.1428571429,	20802.4285714286,	19755.2857142857,	20390.7142857143,	18081.4285714286,	17573.4285714286,	18296.7142857143,	18018.2857142857,	18088.4285714286,	17634.8571428571,	17306.2857142857,	17140.2857142857,	17623.7142857143,	18903.5714285714,	18585.7142857143,	18803.7142857143,	20603.1428571429,	19821.2857142857,	18723.7142857143,	19962.8571428571,	20031,	22057.2142857143,	24478.5714285714,	24203,	23325.2857142857,	22184.7857142857,	21406.1428571429,	25812.2857142857,	23460.7142857143,	22922.7142857143,	22712.8571428571,	20166,	26437.8571428571,	24812,	24957.5714285714,	23928.6428571429,	22512.57143
            ])

            si_s = np.array([0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0
            ])
            ss_s = np.array([11.2215351807508,	10.8764741913177,	11.3773538300105,	11.5297034283505,	12.0694187981056,	12.8589485959011,	13.8066729953971,	14.0486967582358,	13.2856498769012,	12.7489438687301,	11.8368240118816,	11.9690647401725,	10.4672277420269,	10.6925957842695,	9.88507395966846,	9.20551681769492,	8.62825658867592,	8.31482596722311,	7.69388723343822,	7.49208361489501,	7.62781758182509,	7.3321117259359,	7.26119930061931,	7.11530138865928,	7.25478503081339,	7.16571220410885,	7.15894782207598,	6.98353090594681,	6.91115388902453,	6.82704677619464,	6.68152302997339,	6.51048652559908,	6.42895581114131,	6.4395153028093,	6.20871383951853,	6.15336938154324,	5.97017249066192,	7.07516142273649,	7.07503046472795,	6.84130783606311,	7.81435791084271,	7.74683703068595,	7.52344404902104,	6.83285703559384,	6.12079290876657,	5.17011527254747,	4.18531906618582,	3.68276904473096,	2.94789149829404,	2.09990097644253,	1.40090859695719,	0.675679181353315,	1.47473008855391E-14,	1.47473008855391E-14,	1.47473008855391E-14,	1.47473008855391E-14,	0,	1.47473008855391E-14,	0,	1.47473008855391E-14,	1.47473008855391E-14,	1.47473008855391E-14,	0,	2.21209513283087E-14,	2.21209513283087E-14,	7.37365044276957E-15,	1.47473008855391E-14,	0,	1.47473008855391E-14,	2.21209513283087E-14,	0,	0,	0,	1.47473008855391E-14,	7.37365044276957E-15,	0,	1.47473008855391E-14,	1.47473008855391E-14,	6.16431372939959,	7.74575195004377,	9.63574060186352,	11.7614964308488,	11.7762198526655,	13.266970988704,	14.5109972366585,	14.2800541249098,	10.7426351065928,	9.57680682540928,	8.72127756431797,	9.64802660115403,	12.7108030168968,	11.7602910826476,	11.8473059977894,	11.5361337338307,	11.3510152346201,	11.25209116
            ])
            xi_s  = np.array([12.537949247238,	12.7622124355925,	11.6403925175808,	13.4345579644887,	12.4658087722062,	13.8442565821225,	13.9752812776766,	9.40145739589788,	9.29024560992272,	8.69554952786202,	8.0640599128178,	7.26815502078321,	6.37133414498782,	5.87700877198284,	5.242647255988,	4.76220091190022,	4.53287521503795,	4.12200005821598,	3.97586752057466,	4.16474704771022,	3.72454386141738,	3.9236717006298,	3.88616002820694,	3.8987244128849,	3.42029465736249,	3.34270499317741,	3.14512238270008,	3.0852753202513,	2.93504425439174,	2.87637283970554,	2.82580235487041,	2.81797795079593,	2.76135304640102,	2.7115732128741,	2.69752381720351,	2.86212376339833,	4.36392659814258,	3.77397577720348,	5.18029909945422,	8.06112437627031,	8.34018067183172,	10.376989237123,	12.6249200640541,	15.433534698335,	17.0910599125408,	18.4195057056107,	19.4851982338312,	19.1039186439617,	19.2856985490663,	17.1705411555248,	16.4831795067931,	16.0747146247641,	14.8755959718085,	16.0974184525087,	15.5620763495008,	14.936276991728,	13.961012207077,	12.739215781408,	11.8860988969894,	13.5081087149514,	11.7764363907736,	10.8757200837765,	11.3462653932307,	9.51205946662993,	9.18480748690078,	8.51587541911254,	7.8153377490591,	7.61629400519871,	7.13269481617283,	6.77437805101123,	6.43415606460842,	6.5344796710548,	8.12514845984406,	7.49263906323189,	8.300060293704,	7.93958734499414,	8.01959661093986,	8.84779676802948,	15.0699263933226,	33.4218283920202,	59.2959573123551,	92.0810493139303,	72.7007932764721,	67.3813043899916,	43.9908807084441,	19.8843353165704,	13.2926272132075,	10.4570328024879,	10.2203665519158,	10.920218208622,	10.9444543180446,	12.5964406844812,	11.4906026765968,	11.7721119607596,	11.8452420085774,	12.01217523
            ])
            xs_s  = np.array([44.1864737364508,	43.5072527617609,	43.9014451147492,	42.2212188872211,	42.0875963860641,	42.4484825060493,	42.9375499294695,	46.3033897362579,	46.503431431627,	42.8783209731347,	40.7033176621116,	39.9977827829263,	39.4458860028374,	37.9232850207006,	36.3802551940233,	34.5068277581962,	31.6693717457926,	29.8987419060294,	27.6705966775005,	25.216224412564,	24.6254933645434,	23.2348636988844,	22.2589955462998,	21.9179680996683,	22.4066557626343,	21.9792485594299,	21.6623080372635,	20.8466692166577,	20.2332026879611,	19.5738409757106,	18.7175280566285,	17.86333360046,	17.4451877580558,	17.2189791259956,	16.372258986249,	15.6686558600234,	13.5637132896097,	14.2894183920453,	15.823914292451,	18.2158797171636,	18.3940469824987,	18.5890639472014,	19.7248184180229,	20.0388873069665,	20.3077011195059,	20.3353975943519,	20.998464879947,	22.7934636417359,	23.5779172089296,	26.1143785619018,	28.2689269405191,	28.6544036263123,	28.7814299318692,	28.1551533439999,	29.3167346374417,	28.2664616241093,	27.9433579267016,	28.3271579060269,	28.0782798093999,	26.0540586397663,	25.1930911216497,	25.8559532915611,	25.1730539800267,	24.9953548203144,	24.4975675643814,	22.106098632442,	20.8681884637417,	20.4311817723641,	19.8879204687526,	19.4206172421459,	18.5326868882173,	17.0261986671995,	13.7582795404109,	13.1690252077774,	13.1362941468318,	15.0205572744544,	17.113258417681,	19.4613078460036,	34.8176919071434,	48.3915579565359,	57.3074253897461,	59.6798842965361,	42.9889284094536,	37.0305627097475,	22.6122230947671,	24.4669375616162,	18.5525435879968,	20.5372703817109,	21.0702380330887,	22.7843835543391,	28.2878719903873,	27.8144289968934,	30.9819829995996,	34.1927056989203,	37.129255959787,	39.85847733
            ])
            xbh_s = np.array([6.17883042534275,	6.15161756264479,	6.13306291400733,	6.04207206697892,	5.96348207178965,	6.08541724620448,	6.14135056424856,	6.14338757947738,	6.17753345152287,	5.70764899581183,	5.39255522293492,	5.24405869044518,	5.08099952858892,	4.86506191462614,	4.62068080667898,	4.36186602454086,	4.00629229267098,	3.7655263742525,	3.51373311244904,	3.22671650030433,	3.13967360341918,	2.99598087857322,	2.88556043989857,	2.84710955779428,	2.86850397349845,	2.81258654596578,	2.75538668054041,	2.65711415428571,	2.57383264709141,	2.49329656677461,	2.38923900787778,	2.28943410777044,	2.23738208244909,	2.20651626965426,	2.11308135510759,	2.03954209983313,	1.87418694392415,	1.96462104433442,	2.23505316717841,	2.78749187596242,	2.85274068270389,	3.08014009086045,	3.43198638974452,	3.75681510259729,	4.00073583107771,	4.20580970436696,	4.43443405034312,	4.50030902836811,	4.609219080719,	4.79580903716256,	4.95594832769804,	4.9508450159476,	4.79827075994199,	4.85767572513326,	4.88942352160522,	4.78003702684634,	4.63536730647457,	4.47840135821646,	4.32508238884881,	4.27252854079032,	4.08029674098659,	4.07095992675867,	4.04166977787729,	3.75814193042659,	3.65018805476717,	3.38870243895503,	3.17688613873102,	3.09720226109856,	2.98901433273888,	2.90226021811041,	2.76671401080459,	2.60877356558278,	2.32452919018969,	2.24334382484797,	2.33168413407539,	2.46357836766488,	2.68559667383349,	2.95995502638839,	5.47657413954515,	9.02882013188233,	12.8971918384744,	16.8111115703133,	12.7529449540377,	11.4317444980591,	6.62402903959218,	3.17757186999935,	3.05919028668715,	3.22605845610367,	3.34288150167621,	3.46769862255292,	3.4848283307248,	3.93992195566684,	4.1885465951199,	4.80042931861373,	5.22127163226136,	5.564300999
            ])
            xba_s = np.array([0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0
            ])
            xp_s  = np.array([0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0
            ])
            so_s  = np.array([0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0
            ])
            sno_s = np.array([0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0
            ])
            snh_s = np.array([5.33396779473915,	5.32823771371251,	5.72330328150751,	5.87498739935487,	6.19761982536678,	6.65564679653142,	7.24698902038684,	7.34043156839673,	7.01759335122857,	6.88912621746612,	6.47750114312232,	6.66369135991297,	6.24562261702653,	5.87167335995519,	5.34417450156776,	4.90567898196204,	4.5492219731738,	4.33479560617505,	3.90770679637414,	3.8105600075514,	3.87517843062094,	3.75985787734819,	3.67438505957217,	3.62421210662795,	3.67925455939982,	3.61388513224046,	3.62342903118915,	3.50936727836578,	3.45624643392323,	3.38449213569461,	3.3093035305074,	3.24827175330424,	3.25160182837848,	3.25606241850598,	3.16313568469304,	3.09938853327207,	3.02114780470225,	3.50367461391304,	3.39075406659261,	3.24386996064996,	3.57340307192798,	3.41396837806512,	3.21598392266881,	2.84673044561777,	2.53371942431478,	2.10942082187945,	1.66903843790583,	1.51622647761756,	1.22913177894282,	0.884733597389632,	0.59460548361918,	0.295529463469501,	0,	0,	1.10604756641544E-14,	0,	7.37365044276957E-15,	3.68682522138479E-15,	3.68682522138479E-15,	7.37365044276957E-15,	0,	1.10604756641544E-14,	7.37365044276957E-15,	3.68682522138479E-15,	1.10604756641544E-14,	1.10604756641544E-14,	0,	1.10604756641544E-14,	7.37365044276957E-15,	1.10604756641544E-14,	0,	0,	3.68682522138479E-15,	0,	0,	0,	3.68682522138479E-15,	0,	2.76560062302283,	3.43813948478368,	4.26455133396326,	5.25723709218605,	5.2575337521647,	5.87672459742735,	6.45783606094241,	6.23449653152569,	4.60428313423264,	4.03128037679714,	3.73373469337331,	4.19188032967477,	5.61120322619981,	5.28532090348625,	5.34893976951105,	5.24918718357969,	5.14221587150429,	5.114757452
            ])
            snd_s = np.array([1.12215512164253,	1.0876463500868,	1.13773377943361,	1.15297034283504,	1.20694375063925,	1.2858953941126,	1.38066623049475,	1.40486940856233,	1.32856632399634,	1.27489492139548,	1.18368133214319,	1.19690794099184,	1.04672330872516,	1.06926144925565,	0.988508465011808,	0.920551949030727,	0.862826193390069,	0.831481794938575,	0.769388990605059,	0.749208094228253,	0.76277962009258,	0.733209836287384,	0.726119395539446,	0.711529604343447,	0.725480641171264,	0.71657282397834,	0.7158966530363,	0.698353625117156,	0.691113518073767,	0.682702806790762,	0.668150432168645,	0.651051057911076,	0.642894244807904,	0.643952599325888,	0.620871918474321,	0.61533693815434,	0.59701618002123,	0.707515340489904,	0.707501977427827,	0.684132654435005,	0.78143712739048,	0.774683435807353,	0.752345741208313,	0.683285436298142,	0.612079023615415,	0.517011527254747,	0.418531104834857,	0.368275300905644,	0.294789149829404,	0.209989830383011,	0.140090325173235,	0.0675663145678801,	2.76511891603859E-15,	9.21706305346196E-16,	9.21706305346196E-16,	1.84341261069239E-15,	0,	1.84341261069239E-15,	0,	1.84341261069239E-15,	0,	9.21706305346196E-16,	9.21706305346196E-16,	0,	0,	0,	1.84341261069239E-15,	2.76511891603859E-15,	1.84341261069239E-15,	0,	9.21706305346196E-16,	9.21706305346196E-16,	9.21706305346196E-16,	1.84341261069239E-15,	1.84341261069239E-15,	1.84341261069239E-15,	0,	9.21706305346196E-16,	0.616430036633749,	0.774573858698165,	0.963575663753794,	1.17614991034608,	1.17762251978903,	1.32669549530296,	1.45110025818834,	1.42800621427469,	1.0742619070918,	0.957679880757189,	0.872127772873999,	0.964802392854191,	1.27108110347339,	1.17602964278723,	1.18473166882391,	1.15361444242805,	1.13510125620076,	1.125209918
            ])
            xnd_s = np.array([2.32286554016032,	2.31265687080459,	2.30568458190231,	2.27150169119354,	2.24191498542879,	2.28778440311356,	2.30877797724311,	2.30955159785056,	2.32240682731914,	2.1457070932539,	2.02730345998633,	1.97141966620978,	1.91018955137965,	1.82897035030508,	1.73712888494948,	1.63978618236636,	1.50610031779843,	1.41558895588074,	1.32093716547952,	1.21306262200879,	1.18034537845408,	1.12629072954515,	1.08478408454867,	1.07034162015049,	1.07837639351588,	1.05737323433103,	1.03588033228232,	0.998898034119121,	0.967600304598617,	0.93731853296729,	0.898229936466629,	0.860709074315678,	0.841133397056582,	0.829522739683931,	0.794406894477873,	0.766754732339384,	0.704603827455033,	0.738599491599991,	0.840281071966684,	1.04792553948973,	1.07243646112255,	1.1579856356584,	1.29018541235487,	1.41237539480985,	1.50401706849972,	1.58117912543106,	1.66705573576901,	1.69187408389775,	1.73270215879863,	1.80291990308107,	1.86314549222382,	1.86115656390051,	1.80384922422111,	1.82621503572853,	1.83813321149238,	1.7970185812477,	1.74252192616344,	1.68364744302755,	1.62593896862333,	1.60615487829554,	1.53392100498354,	1.53046880924884,	1.51935090049216,	1.41279595449184,	1.37223814706462,	1.27397271181637,	1.19431866864438,	1.16432329290263,	1.12371321174903,	1.09106167169291,	1.04006813016063,	0.980768023191155,	0.873938042951948,	0.843326277908491,	0.876523927688793,	0.92610769583074,	1.00958981579345,	1.11288986198343,	2.05887518022106,	3.39429728388304,	4.84858491052514,	6.31996174876134,	4.7943280411093,	4.29763396814271,	2.49023272137895,	1.19464746823045,	1.15005309182463,	1.2128286755658,	1.25669738488481,	1.30372148745869,	1.31012220179062,	1.48118289120405,	1.57464683386815,	1.804707393986,	1.96285593772817,	2.091872613
            ])
            salk_s = np.array([0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0,	0
            ])
            q_s = np.array([11247.1320253967,	11211.6866421889,	11307.9456320581,	11084.2307363427,	11182.0225872986,	11264.5789023005,	11351.9493441474,	11970.7153402895,	10655.2613235728,	10429.8896860244,	10226.6541519139,	9822.45376815144,	9637.0414775262,	9291.57989883763,	9099.41821259427,	8814.4009566887,	8741.63779258231,	8485.83927383253,	8134.16211399203,	8218.33637111497,	7763.66170546625,	7649.86932876979,	7401.55046609815,	7191.38852385956,	6680.09188962593,	6398.09749194896,	6202.91800572233,	5993.6117365368,	5659.45109946652,	5407.50907697642,	5322.95958173564,	5140.8199093782,	4895.27196934411,	4662.73318009847,	4374.42619910641,	4317.21246362093,	5041.69062271425,	4089.24733395301,	3984.16760647485,	4472.73975886698,	4022.68326882923,	4127.03142919732,	4225.75238375264,	4456.17496746474,	4382.71825744888,	4265.35380184868,	4116.69706570319,	3849.66939996294,	3712.21595575982,	2957.00197612602,	2685.30584747911,	2586.62526223896,	2437.052940124,	2768.87663181259,	2777.46440615762,	2421.77284077341,	2265.47904353655,	2246.55105625304,	2241.19622690038,	2731.7547715467,	2089.59759520586,	1795.98113012922,	1982.6317175594,	1844.71651832118,	1889.26480264295,	1518.19870132447,	1415.55139724482,	1460.47872627098,	1346.13765767346,	1258.78851578513,	1220.89963227446,	1349.90435151597,	2341.28320221349,	2017.30206077215,	2244.91985221235,	2147.4966770689,	2144.24594702162,	2553.7911314078,	3401.75228278113,	5707.18016051845,	8117.01608852161,	10462.8476246948,	10551.2847267342,	10745.2121962694,	11064.5378587392,	11229.6079242273,	7725.17494378528,	5821.48354642686,	4737.45550705425,	5994.01944065458,	9235.42029363039,	10021.1510843288,	10420.1320898168,	10345.6161857927,	10602.9340685792,	11009.51788
            ])



            rnd = np.random.randn(96)


            si = si_m + si_s*rnd
            ss =  ss_m + ss_s*rnd
            xi  = xi_m+ xi_s*rnd
            xs  = xs_m+ xs_s*rnd
            xbh = xbh_m+ xbh_s*rnd
            xba = xba_m+ xba_s*rnd
            xp  = xp_m+ xp_s*rnd
            so  = so_m+ so_s*rnd
            sno = sno_m+ sno_s*rnd
            snh = snh_m+ snh_s*rnd
            snd = snd_m+ snd_s*rnd
            xnd = xnd_m+ xnd_s*rnd
            salk = salk_m+ salk_s*rnd
            q = q_m+ q_s*rnd
            aa = q*(12.3*1)/sum(q)


            start_t = t1%24
            end_t = t2%24

            t_1 = start_t/24*96
            t_2 = end_t/24*96


            t_1 = int(t_1)
            t_2 = int(t_2)


            si_2 = si[t_1:t_2]
            ss_2 = ss[t_1:t_2]
            xi_2 = xi[t_1:t_2]
            xs_2 = xs[t_1:t_2]
            xbh_2 = xbh[t_1:t_2]
            xba_2 = xba[t_1:t_2]
            xp_2 = xp[t_1:t_2]
            so_2 = so[t_1:t_2]
            sno_2 = sno[t_1:t_2]
            snh_2 = snh[t_1:t_2]
            snd_2 = snd[t_1:t_2]
            xnd_2 = xnd[t_1:t_2]
            salk_2 = salk[t_1:t_2]
            q_2 = q[t_1:t_2]

            si_3 = sum(si_2*q_2)/sum(q_2)
            ss_3 = sum(ss_2*q_2)/sum(q_2)
            xi_3 = sum(xi_2*q_2)/sum(q_2)
            xs_3 = sum(xs_2*q_2)/sum(q_2)
            xbh_3 = sum(xbh_2*q_2)/sum(q_2)
            xba_3 = sum(xba_2*q_2)/sum(q_2)
            xp_3 = sum(xp_2*q_2)/sum(q_2)
            so_3 = sum(so_2*q_2)/sum(q_2)
            sno_3 = sum(sno_2*q_2)/sum(q_2)
            snh_3 = sum(snh_2*q_2)/sum(q_2)
            snd_3 = sum(snd_2*q_2)/sum(q_2)
            xnd_3 = sum(xnd_2*q_2)/sum(q_2)
            salk_3 = sum(salk_2*q_2)/sum(q_2)

            """
             List of variables :
                        0=V, 1=Si, 2=Ss, 3=Xi, 4=Xs, 5=Xbh, 6=Xba, 7=Xp, 8=So, 9=Sno, 10=Snh, 11=Snd, 12=Xnd, 13=Salk
                        (ref. BSM1 report Tbl. 1)
            """

            influent_mixed_instant=[0.66, si_3, ss_3, xi_3, xs_3, xbh_3, xba_3, xp_3, so_3, sno_3, snh_3, snd_3, xnd_3, salk_3]
            influent_var_instant=[si_2, ss_2, xi_2, xs_2, xbh_2, xba_2, xp_2, so_2, sno_2, snh_2, snd_2, xnd_2, salk_2]

        else:

            si_m = np.array(
                [30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30,
                 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30,
                 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30,
                 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30
                 ])
            ss_m = np.array([31.336885,	30.472625,	29.936155,	31.60591,	32.195675,	34.011635,	35.88137,	32.01173,	30.53369,	29.60085,	27.15013,	25.92058,	21.047395,	20.94314,	18.440965,	16.654585,	15.96361,	14.892095,	14.626475,	14.90654,	14.97495,	15.326255,	15.78325,	15.56881,	15.34235,	15.52757,	14.751975,	14.7944,	14.93187,	15.22528,	15.10848,	15.504365,	15.611045,	15.693865,	16.75311,	17.364545,	20.4694,	28.8483,	46.954685,	55.833565,	68.53352,	83.194635,	91.46498,	96.34346,	95.64263,	90.13225,	82.252335,	74.6787,	69.087805,	63.54216,	58.83903,	54.621745,	51.92047,	52.35794,	51.9488,	51.67747,	52.14088,	51.542385,	50.46358,	50.52671,	47.090415,	45.09842,	44.852595,	43.67163,	43.95967,	46.318605,	48.85126,	48.84758,	47.394335,	43.886985,	42.895165,	42.28631,	43.328285,	44.643205,	46.689405,	47.67464,	47.347405,	47.03185,	50.53279,	49.69281,	51.397005,	55.32019,	57.56891,	59.05613,	60.6225,	61.194075,	65.346985,	64.78882,	63.69962,	62.34435,	58.530195,	60.462215,	57.84723,	55.290275,	52.815575,	50.379985])
            xi_m = np.array(
                [28.973625,	28.52347,	25.373575,	28.606505,	26.556165,	27.439635,	26.819805,	19.34035,	19.213905,	18.25209,	16.81522,	14.81433,	12.38247,	11.350755,	9.79714,	8.76527,	8.315995,	7.35767,	7.6762,	7.51999,	7.108705,	7.44989,	7.56829,	7.395845,	7.003385,	7.108385,	6.32663,	6.38508,	6.36109,	6.42048,	6.1871,	6.419005,	6.357355,	6.190095,	6.940395,	6.949845,	8.39018,	10.259045,	18.433375,	22.51819,	27.29527,	36.77564,	44.193365,	53.719505,	60.691865,	67.36506,	73.0365,	73.206645,	74.30241,	70.90843,	68.94167,	67.26015,	62.14878,	65.74405,	61.219,	62.258385,	59.05968,	54.574425,	50.528135,	55.767135,	51.564335,	48.857805,	49.654615,	43.68107,	42.012365,	41.967925,	39.876975,	38.958755,	37.07527,	35.733645,	34.64554,	35.086375,	37.86842,	36.67548,	38.05183,	40.6201,	39.050535,	36.750755,	38.142655,	38.29649,	40.94412,	44.254515,	45.798975,	43.550095,	41.3662,	41.07767,	52.763825,	47.784705,	49.23367,	49.914685,	42.443925,	57.31898,	51.79366,	53.892525,	52.152845,	49.28295])
            xs_m = np.array(
                [122.918985,	120.277905,	111.938175,	113.367435,	110.17373,	111.22695,	110.088475,	101.701325,	102.813095,	97.254215,	91.82352,	86.139255,	79.15453,	74.655125,	68.529465,	63.501375,	59.826065,	55.01498,	54.227935,	51.72856,	49.73609,	49.601145,	48.96033,	48.071475,	47.061315,	47.029735,	44.634605,	44.536075,	43.79216,	43.450005,	42.547435,	42.955685,	42.53411,	42.08851,	43.798695,	43.5786,	46.345665,	52.883925,	68.790255,	76.171195,	91.27118,	111.03403,	122.526375,	135.734895,	151.713255,	169.239475,	184.2759,	196.273805,	207.125095,	219.123715,	224.056185,	223.11788,	218.53166,	219.04101,	216.662425,	211.79499,	208.87754,	206.7055,	203.698945,	204.618895,	196.67911,	195.425015,	193.403225,	186.37359,	183.798555,	181.013705,	176.475835,	173.640445,	169.54831,	166.062385,	162.32231,	158.244215,	151.474115,	151.166975,	153.352755,	158.29421,	160.77408,	165.586095,	164.908485,	155.45769,	154.883355,	158.636955,	159.61191,	168.39986,	175.193215,	170.795665,	166.384265,	166.638535,	176.184645,	178.14725,	178.611575,	179.84816,	174.99666,	184.830305,	187.34076,	188.083195])
            xbh_m = np.array(
                [16.876875,	16.533505,	15.25683,	15.77475,	15.1921,	15.407345,	15.21194,	13.449125,	13.55848,	12.834135,	12.071,	11.217155,	10.17062,	9.556255,	8.702875,	8.029645,	7.57146,	6.930355,	6.878315,	6.58312,	6.316005,	6.33913,	6.281065,	6.16312,	6.00723,	6.01528,	5.66234,	5.65807,	5.57254,	5.541145,	5.41478,	5.48601,	5.43226,	5.364185,	5.637615,	5.61436,	6.081725,	7.01578,	9.691805,	10.965345,	13.17393,	16.423495,	18.524515,	21.05068,	23.60041,	26.28937,	28.590255,	29.94201,	31.27007,	32.225905,	32.55552,	32.26436,	31.186765,	31.642865,	30.87546,	30.450185,	29.77103,	29.03124,	28.247335,	28.93185,	27.5828,	27.142595,	27.006395,	25.56152,	25.090215,	24.77582,	24.03937,	23.622215,	22.95829,	22.42161,	21.88526,	21.48101,	21.0379,	20.87135,	21.26737,	22.101585,	22.20245,	22.481765,	22.56132,	21.527965,	21.758555,	22.543525,	22.823655,	23.550045,	24.062045,	23.54115,	24.34941,	23.82489,	25.04662,	25.33996,	24.56156,	26.35159,	25.198815,	26.524895,	26.61043,	26.374045])
            xba_m = np.array(
                [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
                 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
                 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0
                 ])
            xp_m = np.array(
                [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
                 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
                 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0
                 ])
            so_m = np.array(
                [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
                 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
                 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0
                 ])
            sno_m = np.array(
                [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
                 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
                 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0
                 ])
            snh_m = np.array(
                [14.895465,	14.92812,	15.05919,	16.10487,	16.532405,	17.60404,	18.833785,	16.726105,	16.12816,	15.995365,	14.857445,	14.44261,	12.55864,	11.5006,	9.96975,	8.875335,	8.416765,	7.76375,	7.42875,	7.58164,	7.60776,	7.859195,	7.986795,	7.930045,	7.78085,	7.831025,	7.466565,	7.43449,	7.46738,	7.5479,	7.48311,	7.735575,	7.895665,	7.935405,	8.535155,	8.746345,	10.35835,	14.285905,	22.503335,	26.47401,	31.33948,	36.66321,	39.097765,	40.138955,	39.591545,	36.774195,	32.800895,	30.74584,	28.806335,	26.77163,	24.97377,	23.890435,	22.657915,	22.83774,	22.19266,	21.539065,	21.53776,	20.663915,	20.28795,	20.511815,	19.99693,	19.838435,	20.384465,	20.068115,	19.802335,	20.346975,	21.014335,	21.050915,	20.5178,	19.66834,	19.53619,	19.63395,	20.27593,	20.918735,	21.60381,	23.41245,	23.278045,	21.35675,	22.671395,	22.057365,	22.747105,	24.727415,	25.70184,	26.159445,	26.97886,	26.716585,	28.00766,	27.27235,	27.27095,	27.087415,	25.83824,	27.172985,	26.11745,	25.15826,	23.92642,	22.900755])
            snd_m = np.array([3.13369,	3.04726,	2.99362,	3.16059,	3.219565,	3.401165,	3.58814,	3.20117,	3.053365,	2.960085,	2.715015,	2.592055,	2.10474,	2.094315,	1.8441,	1.66546,	1.59636,	1.48921,	1.46265,	1.490655,	1.497495,	1.532625,	1.578325,	1.55688,	1.534235,	1.552755,	1.4752,	1.47944,	1.493185,	1.522525,	1.51085,	1.550435,	1.561105,	1.569385,	1.675315,	1.736455,	2.04694,	2.88483,	4.695465,	5.583355,	6.85335,	8.319465,	9.1465,	9.634345,	9.564265,	9.013225,	8.225235,	7.46787,	6.90878,	6.354215,	5.883905,	5.462175,	5.19205,	5.235795,	5.19488,	5.16775,	5.214085,	5.15424,	5.04636,	5.05267,	4.709045,	4.509845,	4.48526,	4.36716,	4.395965,	4.63186,	4.88513,	4.88476,	4.73943,	4.3887,	4.289515,	4.228635,	4.33283,	4.46432,	4.66894,	4.767465,	4.73474,	4.703185,	5.05328,	4.96928,	5.139705,	5.53202,	5.75689,	5.90561,	6.06225,	6.11941,	6.5347,	6.47888,	6.36996,	6.234435,	5.85302,	6.046225,	5.784725,	5.52903,	5.28156,	5.038])
            xnd_m = np.array(
                [6.34487,	6.215585,	5.73547,	5.93024,	5.71118,	5.79202,	5.718675,	5.05601,	5.097165,	4.824635,	4.53812,	4.216845,	3.82365,	3.592635,	3.271635,	3.01873,	2.846295,	2.605465,	2.585795,	2.4748,	2.374435,	2.38302,	2.36116,	2.31703,	2.258515,	2.261285,	2.12864,	2.12704,	2.0949,	2.083115,	2.035645,	2.06238,	2.042315,	2.01678,	2.119315,	2.11073,	2.28653,	2.637555,	3.643355,	4.122355,	4.95248,	6.174415,	6.964025,	7.913625,	8.872165,	9.8833,	10.748065,	11.256495,	11.755585,	12.11517,	12.23857,	12.12918,	11.72421,	11.89591,	11.60736,	11.44745,	11.19197,	10.91409,	10.619475,	10.87641,	10.369465,	10.20422,	10.152575,	9.60963,	9.43247,	9.31412,	9.0372,	8.880635,	8.630745,	8.42924,	8.22723,	8.075295,	7.90872,	7.8462,	7.99527,	8.309065,	8.346955,	8.451665,	8.481735,	8.09329,	8.17985,	8.47518,	8.580305,	8.85335,	9.045625,	8.850435,	9.15397,	8.956785,	9.416105,	9.526055,	9.23384,	9.906805,	9.47312,	9.97191,	10.00418,	9.914965])
            salk_m = np.array(
                [7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7,
                 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7,
                 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7
                 ])
            q_m = np.array([39466,	39585,	38839.5,	40780,	39937.5,	40241,	39948.5,	35852.5,	35574.5,	35630.5,	35230.5,	34275,	32884,	32502,	31666.5,	31217,	31265.5,	30777,	31532.5,	31756.5,	31623,	32226.5,	32630,	32541,	32195.5,	32423.5,	31545,	31688,	31891,	32124.5,	31868,	32224.5,	32228,	32004.5,	33056,	33087.5,	34655.5,	33638.75,	34908.5,	36432.333335,	35828.166665,	36578,	38483,	40694.5,	40963,	40864.5,	40725,	38816.5,	37801,	35960.5,	34865.5,	34513,	32862,	34668,	33565,	33839.5,	33019,	31298,	29924.5,	32868,	31785.5,	30689,	31336,	29008.5,	28470,	29121,	28763.5,	28783,	28290,	27946.5,	27773.5,	28350.5,	30006.5,	29589,	29876,	31497.5,	30576,	29216.5,	29851,	29881.5,	31370,	32966,	33037.5,	31513.5,	29980,	29579,	35512.5,	33393.5,	32755.5,	32638,	28989,	35506,	33685,	33654.5,	32524.5,	31002.5])

            si_s = np.array(
                [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
                 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
                 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0
                 ])
            ss_s = np.array(
                [1.79504006151673,	1.59286409063988,	0.32230634193264,	2.04819136024933,	1.63789265046583,	2.80046105155026,	3.21351747778038,	1.60469398708913,	0.844497628771093,	0.865682547935443,	1.05223145894808,	0.328182399284301,	0.75786997700793,	0.230347105039331,	0.554378787518064,	0.317540442227442,	1.04995457511266,	0.988528209030982,	0.389764328857836,	1.52573844390184,	0.807728076149392,	1.15233656596066,	1.06841006210163,	1.15020817454929,	0.206135768851503,	0.016899852070358,	0.254381664531861,	0.359379950470251,	0.0957846845795286,	0.19183806973591,	0.507589531806951,	0.66909979169777,	0.641423632282129,	0.68865836526539,	0.38658941941031,	0.982206674407172,	2.97947927534326,	1.59851387382156,	13.6929318988466,	15.6393847354188,	22.6582751379005,	31.9350140612533,	33.5810950082454,	33.4709560560078,	32.9900244498516,	31.0825049126514,	28.5713071031489,	28.4721212349519,	27.2827322736571,	25.0955590755655,	23.946977246613,	22.276090993737,	23.1301416351219,	20.7537678755883,	20.9533699777816,	21.7820001303829,	22.554782989406,	24.2965071996213,	25.2454232867702,	20.7069574066737,	20.583517775882,	20.7377024095197,	19.9318623084761,	22.2281845093116,	22.9853120241993,	22.7801732759093,	24.065106505478,	23.7187656040529,	23.4889488290994,	22.143338766637,	21.8432002222946,	21.287110236366,	21.7250285370227,	22.3811670614213,	23.4202746185105,	20.7855593964704,	21.0286273525033,	21.1892607999855,	22.4560567406168,	24.742543086126,	24.1450732113624,	23.4305347379056,	25.5619242820293,	26.0018578370123,	27.5217980053084,	30.2783194414758,	26.1464399605616,	27.4795978726072,	27.8280176679691,	27.7914036788393,	30.5631561950733,	24.4527777982635,	25.1408563359684,	23.4508782000003,	23.5340905260103,	24.21015532])
            xi_s = np.array(
                [5.15419325195029,	4.86591288832835,	0.678207327039452,	4.60458743520527,	3.24066330711014,	6.65143529257031,	7.42444442576345,	1.10519375685895,	0.136054415768105,	0.712014102247982,	1.25157900270019,	0.318311188618936,	1.07907323236192,	0.347309637715397,	0.69855078913419,	0.423415540574504,	1.31355691207119,	1.17662568389441,	0.7062441109135,	2.0441184251897,	1.18129965971806,	1.80272631222823,	1.76144541834256,	1.92544469410316,	0.945465406056721,	0.810492863663833,	0.896922525528264,	0.912153605595028,	0.448786531883478,	0.354854467070657,	0.898761003159351,	0.987863528656666,	0.990239407441448,	1.04966466133238,	0.901879344064382,	1.47125586641141,	3.48940225952813,	1.17014151471093,	3.30426049101006,	2.07797469787291,	6.08488012628019,	10.3723513833268,	12.767003853154,	15.221880607535,	17.8645366387278,	20.6743032155185,	23.9122300193855,	32.0819367082483,	35.0577743836685,	29.6356513012587,	30.8345649908962,	29.0420775848595,	32.8001662791029,	24.4332970064418,	22.265378326002,	27.9061317540868,	26.1096633368721,	29.3499872327272,	29.8976980032987,	19.5230273097297,	23.2421261360584,	22.7846704750376,	20.9408046902799,	25.2478557340975,	25.1215735340454,	22.1268914629066,	21.1708123840643,	19.5801332533783,	19.0999865357492,	18.6681210691395,	18.2128927944794,	18.6586034118647,	23.0991703581059,	21.7004282921052,	23.9174201831594,	18.0013830140909,	17.4251829113054,	13.4749733046137,	13.3549065731682,	21.4515125629919,	20.3192507585491,	17.6019878908733,	23.5481053824135,	17.9504783969133,	16.9250250717687,	22.7452634141221,	20.9093950070596,	20.7426875123271,	22.813145665116,	24.8029229341715,	26.8772348189178,	25.7938694483941,	27.4828929902076,	25.783199207066,	27.6000110863755,	29.45813921])
            xs_s = np.array(
                [7.3608896682704,	7.49227010821487,	0.552540309886976,	2.21434024062474,	1.59773605636225,	7.61732194246771,	10.0209405262805,	3.16607768383057,	4.21198760815485,	1.16200978672728,	0.759262977366862,	0.0753422275354256,	2.59888612143742,	1.06967578323995,	1.99247841695964,	1.41214173950422,	4.08876718046529,	3.97437851647273,	2.98036315881974,	6.74841498761005,	4.53228576682892,	6.68322681345246,	6.55147160691398,	7.37633288037149,	5.44216955965633,	5.43665412676308,	5.26380894516984,	4.84832714268024,	3.48817189372886,	3.18993546662781,	4.16012839682262,	3.93572098874526,	4.10458515916042,	4.51015332457778,	4.16043952380635,	5.51977452889156,	8.68879377693189,	1.01449317003615,	11.3778077286992,	13.2482607284295,	22.7773094934454,	33.0854838653117,	35.6532077909021,	38.5544386298394,	44.9673880183188,	53.2840668516514,	60.2145264800779,	67.2348169539864,	74.6420574717635,	89.6247733526978,	98.6524480545133,	98.5227717419114,	100.524195059653,	94.8017920115447,	98.9391939964202,	91.2295168372616,	88.0976701927083,	94.1802592930175,	98.768752977883,	92.9380211708968,	95.739274182042,	91.9776004564179,	87.7233490759516,	92.5438869221312,	92.3021059004401,	90.9880893689611,	93.0300299052646,	91.3137261838331,	90.2009310290764,	87.7075664525954,	85.3148231710773,	81.3480884789092,	74.9955825780856,	75.9079483157149,	80.1141305771963,	79.6368971395057,	84.9912935344133,	85.8101019738239,	80.2007794421629,	71.9867371581196,	64.1906515503002,	67.3378424120052,	70.5269576347669,	76.0537749471991,	85.5737586032801,	83.3429074193146,	64.2742174297008,	71.8158365201746,	77.4357656535535,	78.4602148870687,	90.9133198979185,	75.2727047967003,	80.9896631225936,	84.9612768515519,	90.1790814295378,	94.09300939])
            xbh_s = np.array(
                [1.39051841413553,	1.37293973955524,	0.136485750904627,	0.757848763804494,	0.537669854278628,	1.58549603797991,	1.93825039791045,	0.228911678273521,	0.452845324807489,	0.0499288098195823,	0.22360130634681,	0.0436779858738925,	0.408566298169587,	0.157324187746195,	0.299113239509721,	0.204148798796368,	0.600319515091755,	0.572311015488956,	0.40945018164607,	0.976769023259864,	0.634720259996481,	0.942898608441017,	0.92362994865368,	1.03347898711101,	0.709793786955056,	0.694096016412715,	0.684436937781706,	0.640058916194439,	0.43750110765574,	0.393950401002461,	0.562093322500811,	0.546918810976547,	0.566152115424821,	0.617777981519251,	0.56251051550171,	0.776728514862175,	1.35297104405453,	0.24282046865946,	1.6314296944858,	1.70272020016502,	3.20708280607159,	4.82883927979073,	5.380354284845,	5.97550484936628,	6.9813925298754,	8.21747175193199,	9.34759102282777,	11.0350942850616,	12.1890076991443,	13.2513154297243,	14.3874733925314,	14.1741534187831,	14.8142194060453,	13.2485435711421,	13.4669052162403,	13.2367773143031,	12.6896958687669,	13.7260174251383,	14.2958111404862,	12.496203169163,	13.2203512237762,	12.7511222348486,	12.0739967811843,	13.0878111287106,	13.0472302705383,	12.5683704872509,	12.6892150361557,	12.3217387130409,	12.1448559517271,	11.8191342840413,	11.5028454208165,	11.1114618174298,	10.8994853465657,	10.8451088350924,	11.5592583998196,	10.8490191355924,	11.3795401403132,	11.0311981266973,	10.3954313486647,	10.3817912588556,	9.38959316563023,	9.43771885315779,	10.4529403431786,	10.4453177320774,	11.3885981781802,	11.7872579103454,	9.46475154540257,	10.284316589069,	11.1388834184042,	11.4735712000754,	13.0877545601681,	11.229435512803,	12.0521896080526,	12.3052206986323,	13.0865241943688,	13.72770741])
            xba_s = np.array(
                [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
                 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
                 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0
                 ])
            xp_s = np.array(
                [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
                 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
                 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0
                 ])
            so_s = np.array(
                [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
                 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
                 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0
                 ])
            sno_s = np.array(
                [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
                 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
                 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0
                 ])
            snh_s = np.array(
                [0.853244539654371,	0.780320617310602,	0.162139584926076,	1.0436613247601,	0.841054018746716,	1.44948404861868,	1.6867395869102,	0.838451865791948,	0.446071241843722,	0.467786491093959,	0.575818265123641,	0.182857813614842,	0.45220892870442,	0.126487261018649,	0.29971428027373,	0.169217723805752,	0.553586827923136,	0.515353564264379,	0.197961614460986,	0.776007265945365,	0.410348207258178,	0.590907923834162,	0.540646773827422,	0.58585918141649,	0.104538666530619,	0.00852063671329751,	0.128757073786259,	0.180595071915044,	0.0478994133575765,	0.0951058620695907,	0.251404744983065,	0.333832182465981,	0.324413520140576,	0.348214734395315,	0.196950451763889,	0.494727259457168,	1.50773550525282,	0.791598970470527,	6.56242469095456,	7.41555609352395,	10.3613346596758,	14.073503839549,	14.3546282825453,	13.9447891921696,	13.6563179097167,	12.6817409174786,	11.3937741998185,	11.7222040864762,	11.375601555542,	10.5732828462687,	10.1641084362673,	9.74310412981663,	10.0939139460989,	9.05248000769955,	8.95133545371862,	9.07869245607814,	9.31666924328648,	9.74074239316747,	10.1494571837611,	8.40619148076286,	8.74078733855253,	9.1223492287486,	9.05856819708558,	10.2143625152062,	10.3541009573043,	10.006939812013,	10.3520786319101,	10.2216315729168,	10.1687611988875,	9.92373455707073,	9.94827116237792,	9.88381130820496,	10.1664560307808,	10.487277448473,	10.8368781121594,	10.2075389347678,	10.3385870345251,	9.62185653004657,	10.0748503472781,	10.9825774644229,	10.6860451646645,	10.4731494549849,	11.412208473604,	11.5177583028318,	12.2480369732378,	13.2191420711803,	11.2063414053294,	11.5673194171251,	11.9136744606859,	12.0748311671862,	13.4921489283287,	10.9895919636923,	11.3508457471327,	10.6706514500381,	10.6613742090689,	11.00497861])
            snd_s = np.array(
                [0.179506127472017,	0.159282873530082,	0.032229927086483,	0.204820550238495,	0.16378714372624,	0.280049640688932,	0.321351747778038,	0.160470812922475,	0.0844497628771093,	0.086571083220669,	0.10522456010837,	0.0328168257148679,	0.0757877048075742,	0.0230304678632457,	0.0554371716450253,	0.0317490944752759,	0.104991214870579,	0.0988535280098793,	0.0389757257790025,	0.152572430176621,	0.0807728076149394,	0.115237192129972,	0.106836763569475,	0.115017989027804,	0.0206121626715879,	0.00168998520703577,	0.0254417019870921,	0.0359351666199004,	0.00958129688507768,	0.0191838069735909,	0.0507561247535703,	0.066913514703683,	0.0641416561214318,	0.0688651294197578,	0.0386575277274686,	0.0982242029746232,	0.297946513320764,	0.159848558955031,	1.36929106856431,	1.56394342328935,	2.26582468536292,	3.19349928480499,	3.3581066723974,	3.34709701981434,	3.29899961655803,	3.10825049126514,	2.85712858899456,	2.8472078808545,	2.72827252025893,	2.50955732177012,	2.39469489623417,	2.22760556383979,	2.31301699193932,	2.07537961598595,	2.09534124041885,	2.17820001303829,	2.25547971315417,	2.42964718442823,	2.52454091446346,	2.07069149802669,	2.05835248469498,	2.07376741252484,	1.99319259480864,	2.22281845093116,	2.2985283739928,	2.27801520627058,	2.4065106505478,	2.37187656040528,	2.34889559001672,	2.21433316955692,	2.18431648669555,	2.12870960942304,	2.17250073238193,	2.23812024167604,	2.34202251210359,	2.07856159650129,	2.10286485657067,	2.11892325157142,	2.24560143142099,	2.47425148018547,	2.41450802824302,	2.34305488800412,	2.55619101398937,	2.6001871979148,	2.75218687159866,	3.0278312370408,	2.61464046052225,	2.74795837304716,	2.78280459522403,	2.77913471102968,	3.05631349818699,	2.44528131536026,	2.51408280516971,	2.34509135553393,	2.35340693128069,	2.421020482])
            xnd_s = np.array(
                [0.522665048381848,	0.516449579775219,	0.0514632315347568,	0.284921606411308,	0.202175970876858,	0.595949595184023,	0.728638182673678,	0.0862670273047587,	0.170235957570661,	0.0188443957186215,	0.0839901434693382,	0.0162846691707264,	0.153456313653105,	0.0590929137037597,	0.112408765005225,	0.0766079486737506,	0.225545849995073,	0.215151380311631,	0.153944217332123,	0.367271262148294,	0.238797031074509,	0.354628192900674,	0.347161145291346,	0.388498607719517,	0.266826743880744,	0.260887046918777,	0.257401010487527,	0.240670864044653,	0.164260905269635,	0.148004520370156,	0.211248150879481,	0.205782215460909,	0.212789643662468,	0.232298719755405,	0.211375430100095,	0.292063384901292,	0.508692618385605,	0.0912379879765001,	0.613167645305918,	0.640292261432231,	1.20532007707497,	1.81526331545868,	2.02229003885447,	2.24630146713437,	2.62454702652668,	3.0896323697165,	3.51422877861558,	4.14860255628929,	4.58263884071721,	4.98203396393481,	5.40856077434654,	5.3285021445806,	5.56887602377716,	4.98098744589865,	5.06237543641323,	4.97598112988785,	4.77018477229132,	5.16033800987881,	5.37475399913801,	4.69743762664285,	4.97030306243493,	4.79387284946107,	4.53881236241927,	4.92057224251408,	4.90524216749795,	4.72471780626103,	4.76985950317197,	4.63206560472215,	4.56544200379875,	4.4431196017213,	4.32433980461758,	4.17716967024922,	4.09737266999232,	4.07689485760916,	4.34549643951068,	4.07849999000245,	4.27805966578892,	4.14694792642131,	3.90784683943089,	3.90281931021665,	3.5300891837176,	3.54800726955287,	3.92966815469831,	3.92677608796325,	4.28135478338925,	4.43153012157765,	3.55820374933757,	3.86676393544395,	4.18733786576269,	4.31327358349201,	4.92027525766598,	4.22170325532835,	4.53097054821591,	4.626019841743,	4.91979442505477,	5.160514787])
            salk_s = np.array(
                [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
                 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
                 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0
                 ])
            q_s = np.array(
                [2193.44523524067,	2020.91118063115,	393.858477120907,	2733.67481606719,	2020.20407384997,	3330.47293938864,	3540.48365340104,	1421.99173696615,	765.796644025031,	813.879905145716,	1038.73986156304,	309.712770159708,	762.261110119098,	223.445742854949,	555.07882323144,	333.75440072005,	1156.11958724001,	1098.84393796389,	484.368145112785,	1901.41013461063,	989.949493661167,	1477.1460658987,	1391.58614537513,	1503.3090168026,	263.750829382582,	21.920310216783,	313.955410846827,	449.719912834644,	121.622366364086,	245.366053071732,	634.98188950552,	849.235244205044,	808.93015767741,	842.164176393178,	497.803173955329,	1222.58762467154,	3647.9638841414,	5083.0370965595,	9751.70961934368,	12707.6516646068,	12654.1472562111,	14515.4880041974,	14193.0473119764,	14159.8132932606,	14194.4615255388,	14346.4894834939,	14125.1650609825,	11704.7385489809,	11329.2648481709,	13926.4680554691,	13975.9655301521,	14374.0666479601,	12562.4590745602,	15570.4913217278,	16259.2133266035,	13735.5492245487,	13681.1020023973,	12307.900633333,	11800.9050712223,	16461.4458660228,	14946.1160339401,	14235.4737188476,	14832.271842169,	12411.8453301675,	12129.709724474,	13624.5334599024,	14063.6467710192,	14625.7966620625,	14488.6179465124,	14258.8082426267,	14140.0143033874,	13777.9756314199,	11594.4298911158,	12190.5209076561,	11762.014198257,	15470.7892655805,	15932.5299936953,	17215.9288015489,	16932.3789822931,	12051.2208717623,	12351.7412537666,	14741.7621741771,	12462.7570184129,	14999.8561493102,	15798.1797052698,	12732.164702045,	13905.2548520335,	14070.7178388311,	13599.7847225609,	12850.9586412843,	11446.6445738478,	13381.2887271742,	12522.8610948138,	13701.6080990517,	13029.8566569245,	11916.87058])

            rnd = np.random.randn(96)

            si = si_m + si_s * rnd
            ss = ss_m + ss_s * rnd
            xi = xi_m + xi_s * rnd
            xs = xs_m + xs_s * rnd
            xbh = xbh_m + xbh_s * rnd
            xba = xba_m + xba_s * rnd
            xp = xp_m + xp_s * rnd
            so = so_m + so_s * rnd
            sno = sno_m + sno_s * rnd
            snh = snh_m + snh_s * rnd
            snd = snd_m + snd_s * rnd
            xnd = xnd_m + xnd_s * rnd
            salk = salk_m + salk_s * rnd
            q = q_m + q_s *(-rnd)
            #aa = q * (12.3 * 1) / sum(q)

            t1 = 0
            t2 = 12

            start_t = t1 % 24
            end_t = t2 % 24

            t_1 = start_t / 24 * 96
            t_2 = end_t / 24 * 96

            t_1 = int(t_1)
            t_2 = int(t_2)

            si_2 = si[t_1:t_2]
            ss_2 = ss[t_1:t_2]
            xi_2 = xi[t_1:t_2]
            xs_2 = xs[t_1:t_2]
            xbh_2 = xbh[t_1:t_2]
            xba_2 = xba[t_1:t_2]
            xp_2 = xp[t_1:t_2]
            so_2 = so[t_1:t_2]
            sno_2 = sno[t_1:t_2]
            snh_2 = snh[t_1:t_2]
            snd_2 = snd[t_1:t_2]
            xnd_2 = xnd[t_1:t_2]
            salk_2 = salk[t_1:t_2]
            q_2 = q[t_1:t_2]

            si_3 = sum(si_2 * q_2) / sum(q_2)
            ss_3 = sum(ss_2 * q_2) / sum(q_2)
            xi_3 = sum(xi_2 * q_2) / sum(q_2)
            xs_3 = sum(xs_2 * q_2) / sum(q_2)
            xbh_3 = sum(xbh_2 * q_2) / sum(q_2)
            xba_3 = sum(xba_2 * q_2) / sum(q_2)
            xp_3 = sum(xp_2 * q_2) / sum(q_2)
            so_3 = sum(so_2 * q_2) / sum(q_2)
            sno_3 = sum(sno_2 * q_2) / sum(q_2)
            snh_3 = sum(snh_2 * q_2) / sum(q_2)
            snd_3 = sum(snd_2 * q_2) / sum(q_2)
            xnd_3 = sum(xnd_2 * q_2) / sum(q_2)
            salk_3 = sum(salk_2 * q_2) / sum(q_2)

            """
             List of variables :
                        0=V, 1=Si, 2=Ss, 3=Xi, 4=Xs, 5=Xbh, 6=Xba, 7=Xp, 8=So, 9=Sno, 10=Snh, 11=Snd, 12=Xnd, 13=Salk
                        (ref. BSM1 report Tbl. 1)
            """

            influent_mixed_instant = [0.66, si_3, ss_3, xi_3, xs_3, xbh_3, xba_3, xp_3, so_3, sno_3, snh_3, snd_3, xnd_3, salk_3]
            influent_var_instant = [si_2, ss_2, xi_2, xs_2, xbh_2, xba_2, xp_2, so_2, sno_2, snh_2, snd_2, xnd_2, salk_2]


        influent_mixed = influent_mixed_instant
        influent_var = influent_var_instant

        return switch, influent_mixed , influent_var
