Commit 403e6135 authored by brunner's avatar brunner
Browse files

No commit message

No commit message
parent c3e5e418
......@@ -160,8 +160,8 @@ class ObservationOperator(object):
with Pool(self.forecast_nmembers) as pool:
pool.starmap(self.extract_model_data, args)
pool.close()
pool.join()
# pool.close()
# pool.join()
for i in range(0,self.forecast_nmembers):
idx = str(i).zfill(3)
......
......@@ -296,91 +296,93 @@ class StateVector(object):
rands_bg = np.random.uniform(low=-0.05, high=0.05, size=1)
# if member == 1 and lag == 0:
# rands=np.array([-0.2236048,0.1492591,-0.4134172,0.2030501,-0.3056341,0.6603037,0.3247675,0.6550926,0.8024981,0.1593491,-1,0.08171216,-0.4006799,-0.3123662, 0.1022556, 0.1977177, -0.4678199, 0.8751256, 0.9531838,0.9820679,0.3594339,-1,0.03129737])
# rands=np.array([-0.5402309, 0.3021888, 0.346095, 0.4942706, 0.6276228, 0.7450897, 0.5155229, -0.3006688, 0.8386819, -0.1197997, -1, -0.9025682, -1, 0.5018581, 0.923963, -0.2253573, 0.4716174, 0.4637589, 0.7411034, 0.6526616, -0.3188075, -1, 0.01435228])
# if member == 2 and lag == 0:
# rands=np.array([-0.1477133,-0.6238207,-0.8987545,-1,-1,0.1647741,0.6627979,0.4307449,-0.4548403,-0.1154654,-1,-0.05427124,-0.9482519,0.7116647,0.03384373,-0.09207293,0.4631868,-0.5646629,0.7750968,-0.8704156,0.1906742,-1,0.004112379])
# rands=np.array([0.1876527, -0.4323478, 0.2562822, -0.7032343, 0.83417, 0.1378682, 0.3969375, 0.2101818, 0.5018877, 0.1390382, -1, 0.9267841, 0.1679351, 0.1019922, 0.4781876, -0.08182608, 0.9295725, 0.05694205, -0.119032, -0.8025782, -0.8610154, -1, -0.03559723])
# if member == 3 and lag == 0:
# rands=np.array([-0.5120361,-0.6635633,0.4805151,0.3135487,0.02116868,0.0115789,-0.2708648, 0.6396611,-0.4579386,-0.8980371,-1,0.4571423,0.2513106,0.3329439,0.5895487,-0.2265133,-0.3762723,0.9784772,0.6362864,-0.0876907,-0.1079749,-1,0.00497188])
# rands=np.array([-0.8030154, 0.08172598, -0.7622297, -0.4516532, 0.4969012, 0.8023733, -0.08237687, -0.8089882, -0.8048904, -0.4082271, -1, -0.1105984, -0.4859152, 0.3398715, -0.09020101, 0.462075, 0.5386102, -0.865896, 0.5063729, -0.01304542, -0.2682804, -1, 0.04435889])
# if member == 4 and lag == 0:
# rands=np.array([0.5773118,0.8732759,1.126802,-0.3782958,0.662484,1.199195,-0.5346411,-0.0550764,-0.6245533,0.6759356,-1,-0.8961901,0.3352093,-0.5532659,-0.7597447,0.5648068,-0.3713143,0.124998,-0.6321908,0.2655768,-0.02420997,-1,-0.01510416])
# rands=np.array([-0.555479, -0.6089986, -0.664399, -0.03338498, -1, -1, 0.1820136, 0.1071932, -0.9634671, -1, -1, -0.8547595, 0.273797, -0.802806, 0.1928162, -0.5642408, 0.6471559, 0.2174408, -0.05968733, -0.004075899, 0.9075829, -1, -0.02335107])
# if member == 5 and lag == 0:
# rands=np.array([0.1289302,-0.08508042,-0.3619103,-0.5414423,-0.07692706,-0.9641049,-0.268249,-0.7405945,-0.2146955,-0.1976647,-1,-0.1085423,-0.946457,-0.008138934,-0.530806,0.315332,-0.2215554,0.7034806,0.527541,-0.4504287,-0.3679401,-1,0.02846291])
# rands=np.array([-0.674361, 0.3454686, -0.4433188, -0.8486685, -0.7747395, 0.203481, -1, -0.09257206, -0.3892125, -0.02775251, -1, -0.6137948, -0.8518382, -0.8779979, -0.6405048, -0.8992591, -0.6312428, 0.4434872, -0.05931631, -0.2671082, -0.7183358, -1, -0.02273439])
# if member == 6 and lag == 0:
# rands=np.array([-0.2574656,-0.3259587,0.8530074,0.04455876,0.5508205,0.553733,0.5282572,-0.7782328,0.7970823,-0.4580753,-1,-0.1998393,0.04793931,0.01853116,0.6679623,-0.4434216,-0.1505937,-0.1531717,-0.1313033,0.3774659,-0.8304669,-1,0.0197916,])
# rands=np.array([-0.4753621, -1, 0.6610132, -0.3525974, -0.987463, -0.2421917, -0.8256455, 0.1535834, 0.6603956, 0.2776085, -1, -0.6808311, -0.6428144, -1, -0.518185, -1, -1, 0.4222396, 0.7343495, 0.6562976, -0.2627398, -1, 0.00420737])
#if member == 7 and lag == 0:
# rands=np.array([0.6881989,0.3826058,-0.5495766,0.1132997,-0.57166,0.3214818,-0.9898017,0.6837509,-0.3501504,0.7928728,-1,-0.2640067,0.5886974,0.262185,0.2781147,0.7389106,0.6068769,0.1230468,-0.4471386,0.09954312,-0.4604041,-1,-0.04947798])
# rands=np.array([0.9235092, -0.08976277, -0.693751, 0.7479168, -0.5282985, -0.325265, 0.5994204, 0.1087474, -0.8495662, 0.5450585, -1, -0.4888758, -0.6259907, 0.5597795, 0.7957715, 0.04843958, 0.2360165, -0.7671149, 0.2514055, -0.872794, -0.5365307, -1, -0.02834922])
# if member == 8 and lag == 0:
# rands=np.array([-0.8718347,-0.03896444,-0.384823,-0.5786896,0.2542168,-0.5035606,0.7728764,-0.6899376,-0.5994669,0.1059537,-1,-0.859535,0.2382655,-0.3679568,0.450976,0.6301627,-0.6471004,-0.8724778,0.7950776,-0.1000319,-0.1237619,-1,-0.04825296])
# rands=np.array([-0.2851942, 0.2424788, 0.6036716, 0.1861704, 0.2724828, -0.1842902, 0.4256651, 0.9268349, -0.5386767, -0.282816, -1, 0.1387041, 0.6110855, -0.5958148, -0.9803032, -0.5043366, -0.6089351, 0.2659722, -0.7406561, 0.409389, -0.07663137, -1, -0.04040922])
# if member == 9 and lag == 0:
# rands=np.array([0.06249616,0.78708,-0.5142601,-0.1251998,-0.5136412,-0.9106983,0.7841995,-0.8327041,-0.03150041,-0.7065354,-1,-0.1541739,0.7130719,0.6997306,0.5323737,-0.6959949,0.1497954,0.2615451,-0.451424,-0.8742322,-0.4250237,-1,0.03196495])
# rands=np.array([0.7698311, 0.2813426, 0.4236566, -0.5569826, 0.06671345, 0.1696645, -0.5119113, 0.5269542, 0.4109541, 0.5580103, -1, -0.4272517, -0.556591, 0.1682963, 0.02460092, -0.2566992, -0.9167355, -0.7835014, 0.7101904, -0.5695945, -0.1746321, -1, 0.04223957])
# if member == 10 and lag == 0:
# rands=np.array([-0.1827562,0.3512708,0.3367097,-0.6161442,0.2948794,0.3662989,0.7885535,-0.6752307,0.4505339,0.9947404,-1,-0.1427709,-0.4374289,-0.4191162,-0.9822996,-0.2934023,-0.2808521,-0.8965437, 0.7602426,-0.1778743,-0.000989437,-1,-0.01274748])
# rands=np.array([-0.6941028, -0.9384637, -0.9250847, -1, -0.959329, -1, 0.04787683, -0.896835, -0.4911715, -0.01541466, -1, 0.9626108, 0.6731427, 0.5944214, 1.05593, 1.0805, -0.1627321, -0.154229, 0.6146008, 0.420991, 0.7912573, -1, -0.02473733])
# if member == 11 and lag == 0:
# rands=np.array([0.3784615,-0.7083653,-0.821948,0.06502473,-0.6035533,-0.7092476,0.9460906,0.2481164,0.772471,0.6445355,-1,-0.1532128,-0.5447619,-0.5863768,0.2508805,-0.2562239,-0.9050175,-0.6196863,-1,-0.3799991,0.5817245,-1,-0.02425112])
# rands=np.array([0.96286, 0.9415601, -0.7406601, 0.119507, 0.5506695, 0.5813184, 0.9465097, -0.658381, -0.5169885, 0.3601404, -1, -0.1221153, -0.911568, 0.4350488, 0.983259, 0.4603367, 0.4221282, -0.3391019, -0.192896, 0.3793332, 0.8009401, -1, 0.03354285])
# if member == 12 and lag == 0:
# rands=np.array([-0.3772987,-0.09633347,-0.9725795,-0.9913203,-0.1097328,0.3683293,-0.3654553,0.3282615,0.5734389,-0.03295924,-1,0.6370476,-0.4982804,0.174771,-0.740661,-0.02542321,-0.7887223,0.7312356,-0.3222232,-0.1027924,0.02633788,-1,-0.040399])
# rands=np.array([0.4491677, -0.5364176, 0.1777764, 0.3457868, 0.4886276, 0.6140389, -0.9525508, 0.7758317, 0.7040095, 0.8019212, -1, 0.2732197, 0.994552, -0.199925, 0.7970242, 0.5645052, -0.6018554, 0.01864898, 0.0417577, 0.3765279, -0.2793771, -1, -0.008472145])
# if member == 13 and lag == 0:
# rands=np.array([0.7854867,0.3317518,0.6469278,0.5409229,-0.7125678,-0.5082978,0.4885957,-0.490857,-0.1638525,-0.3768145,-1,-0.8093954,0.5902002,0.03597473,-0.7890476,0.5711794,-0.3381689,-0.67668,-1,0.7070413,-0.1368402,-1,0.01318651])
# rands=np.array([-0.8257746, 0.5533802, -0.7910783, 0.4864411, -0.5154006, -0.5978034, -0.202793, -0.9381504, 0.2701632, 0.3198348, -1, 0.3341525, 1.001176, -0.2380723, 0.1933595, -0.6892359, -0.4271372, 0.6341484, -0.5026917, -0.3053397, 0.750175, -1, -0.02964199])
#if member == 14 and lag == 0:
# rands=np.array([0.7588015,0.4496917,-0.1509738,0.5066012,0.7388234,0.310389,-0.2899449,-0.6862295,-0.5775338,0.02781044,-1,-0.3429568,-0.7947903,-0.6043629,-0.8869596,-0.5914347,-0.2657907,-0.6107077,-0.2889067,-0.283401,0.8258069,-1,-0.03961588])
# rands=np.array([0.701148, -0.1619855, -0.063532, -0.4197367, -0.6209376, -1, -0.8872383, -0.4309417, -1, -1, -1, 0.1452325, -0.2525761, 0.1687515, -0.3141938, -0.1218985, -0.4126042, 0.1641715, 0.3630136, -0.5947977, 0.7974821, -1, -0.004756214])
# if member == 15 and lag == 0:
# rands=np.array([-0.1068678,-0.6150778,0.2749358,-0.006487241,0.3801498,-0.2595585,-0.1739224,0.4673905,0.2893905,-0.7730854,-1,-0.6906102,0.3519762,0.6020219,-0.1512664,0.4392967,0.4572283,-0.109156,0.1201844,-0.8120148,-1,-1,0.008621466])
# rands=np.array([-0.2590746, -0.5643597, -0.2787321, -0.2195021, -0.7023828, 0.5482118, 0.6232124, 0.5892509, 0.6535864, -0.5035124, -1, -0.6891443, 0.3709087, -0.1523427, -0.9473059, 0.1686965, 0.04820674, 0.07875252, 0.5189905, 0.2113448, -0.708707, -1, -0.01257968])
# if member == 16 and lag == 0:
# rands=np.array([-0.4203033,0.6120592,-0.07144696,0.07225665,-0.7220178,-0.5390107,-0.3260586,0.7489389,0.4817865,0.04893194,-1,0.3428221,0.3335567,1.011872,0.3968997,1.008683,0.6700837,0.886711,0.5548246,0.4422186,-0.3110386,-1,0.04391268])
# rands=np.array([0.4056181, 0.04536337, -0.6866099, -0.6401204, -0.2879847, -0.7785423, 0.6233089, 0.2740796, -0.8484439, 0.9647028, -1, 0.789863, 1.163193, -0.2164443, 0.08588136, 0.5577631, -0.3278782, -0.1751566, -0.7633625, 0.7689492, -0.819483, -1, -0.03903025])
# if member == 17 and lag == 0:
# rands=np.array([0.6939667,-0.02326318,-0.6011316,0.5709327,-0.6587268,-0.2600601,0.09742739,-0.6605019,0.8611858,-0.3578327,-1,0.8069043,1.003004,0.5246789,-0.2310518,-0.3621024,-0.1325169,-0.2894324,0.5083101,-0.04920054,0.2199543,-1,0.009945367])
# rands=np.array([0.2155399, 0.1660601, 0.4649725, 0.5862293, -0.4452171, -0.7389292, -0.07185336, 0.9012561, -0.2824533, -0.7226023, -1, 0.4194106, -0.2461459, 0.109305, 0.2449136, -0.09785557, 0.7447503, -0.1461458, 0.3223822, 0.4597907, 0.9065241, -1, -0.0007600483])
# if member == 18 and lag == 0:
# rands=np.array([0.3914793,-0.06666035,0.7248741,-0.4003048,-0.6606267,-0.807146,-0.5127164,-0.3519507,0.3180182,0.1996123,-1,-0.9029648,0.1214889,0.02123361,0.6533288,-0.3655816,0.7919229,-0.02215478,0.246939,0.4821936,-0.5654917,-1,-0.003499507])
# rands=np.array([0.04116038, 0.7824114, 0.4636176, 0.1557975, 0.4527292, -0.4133459, 0.994531, -0.3580682, 0.5178631, 0.2745949, -1, 0.5335826, 0.265704, 0.1912572, 0.536295, -0.7093778, 0.07351156, -0.2118817, -0.7678121, 0.43079, 0.512942, -1, 0.03622137])
# if member == 19 and lag == 0:
# rands=np.array([-0.8513823,-1,-0.9796664,0.3713398,-0.4331474,0.3044456,-0.761591,-0.5368128,0.2107361,0.3060099,-1,-0.2859453,-0.5175856,-0.4946607,0.3281591,-0.2430111,0.05990436,0.01103085,-0.6508834,-0.002251296,0.3083046,-1,0.005067795])
# rands=np.array([0.3184969, 0.2478539, 0.1378153, -0.6565253, 0.6140021, -0.2936001, 0.6872, 0.8872706, 0.4869822, -0.2930091, -1, 0.2853186, 0.9346947, 0.7716588, -0.1080804, 0.3278734, -0.2541114, -0.7900179, -0.259804, 0.6708161, -0.8029982, -1, 0.000798168])
# if member == 20 and lag == 0:
# rands=np.array([0.1491016,-0.5305333,-0.4203315,-0.08089796,0.4151005,-0.7401157,0.3687969,-0.5454835,0.714487,1.023472,-1,0.3852504,0.2557997,0.1142526,-0.3280363,0.6893721,-0.4994218,0.7723718,1.014927,0.1936809,0.4526471,-1,0.01181981])
# rands=np.array([-0.7732743, -0.3873762, -0.8547782, -0.07859742, -1, -0.4397394, 0.2048925, -0.3029643, 0.7081758, -0.06017716, -1, -0.466458, -0.2800255, 0.1250267, 0.632107, -0.7519485, 0.4637258, 0.5653031, 0.4806009, -0.01624972, 0.05281375, -1, 0.01537185])
#if member == 1 and lag == 1:
# rands=np.array([-0.8534008,-1,-0.01288013,-0.9141065,0.4783214,-0.1723187,-0.8843083,-0.1162029,-0.009133766,-0.6961404,-1,-0.992416,-0.2270698,-0.795909,-0.9905862,-1,-1,0.5997077,-0.08789548,0.9465016,-0.3178105,-1,0.02080303])
# rands=np.array([0.4061777, -0.5013675, -0.2862317, 0.6116837, 0.05803703, -0.3041813, 0.4205766, -0.04346657, -0.4199566, -0.2088348, -1, 0.8685411, -0.1089062, 0.8366363, -0.278753, -0.4832861, -0.147602, 0.3227859, 0.2312384, 0.357722, -0.03119155, -1, -0.03759673])
# if member == 2 and lag == 1:
# rands=np.array([0.2119585,-0.06496529,-0.113645,-0.3219482,0.3354161,-0.2561546,-0.8578026,0.488537,0.4687303,0.9537283,-1,0.2903501,-0.2037528,0.2844993,-0.2640329,0.6626308,1.016451,-0.2284834,-0.9198995,-0.772432,-0.9456592,-1,-0.00193539])
# rands=np.array([-0.3528758, 0.2302615, 0.8945007, -0.2294943, -0.7223361, 0.4254291, 0.2018437, 0.555488, -0.03991819, 0.2198012, -1, 0.918696, 0.3752462, 1.003467, -0.3229041, 0.6925161, 1.036227, 0.7551609, 1.008353, 0.6003665, -0.6694191, -1, 0.03411322])
# if member == 3 and lag == 1:
# rands=np.array([-0.8951237,-0.03746896,0.6338864,-0.7385425,-0.5215907,-0.2300481,-1,-0.5272307,-0.3447395,-0.161121,-1,-0.3376194,-0.5966941,0.5816291,0.4801027,0.8620167,0.02183707,-0.5763896,0.3380404,-0.4263453,0.5804913,-1,0.04326432])
# rands=np.array([0.732138, 0.0511689, 0.5308715, -0.3620436, 0.1138427, 0.3123921, 0.7225392, 0.8311706, -0.7413127, 0.02754115, -1, 0.2934451, -0.7651439, -0.574313, -0.9932647, -0.5833654, -0.5632083, 0.2126374, -0.4117518, -0.6983544, -0.5544686, -1, 0.03195827])
# if member == 4 and lag == 1:
# rands=np.array([-0.3111334,0.1525611,-0.5140871,-0.9283414,0.08702877,-0.8027726,0.7558329,-0.3327861,-0.3074336,0.4950007,-1,-0.3270844,0.6640525,-0.2391291,-0.08126307,-0.01146206,0.01725924,0.7988746,1.104426,0.1875366,0.5716885,-1,-0.01126271])
# rands=np.array([-0.5030815, -0.02779967, -0.9290833, 0.1386678, -0.510511, 0.297296, -0.6832802, -0.1534745, -0.7899511, 0.4862387, -1, -0.991854, -0.2682375, -0.7712117, -0.8905275, 0.3489431, 0.6761383, 0.1792146, 0.8675478, -0.7689523, -0.8875247, -1, 0.01273509])
# if member == 5 and lag == 1:
# rands=np.array([0.7825249,-0.4134479,-0.4143417,-0.9416251,0.497882,-0.7672758,0.4174937,0.8023338,0.9017133,0.4341536,-1,0.1975201,0.7627142,0.2342599,-0.5487284,-0.8285513,-0.1792318,0.7551883,0.7257428,0.6304003,0.4330811,-1,-0.03357022])
# rands=np.array([0.3466186, 0.6708195, -0.4234087, 0.791083, 0.5170874, -0.4746583, 0.9374641, 0.9230136, -0.5804763, -0.02489379, -1, 0.9292546, 0.3582832, -0.2736852, -0.4632118, 1.018919, 0.09223775, 0.3311332, 0.2550537, 0.7294372, 0.4026192, -1, 0.03795173])
# if member == 6 and lag == 1:
# rands=np.array([-0.3209096,-0.5156574,0.4877877,-0.235484,-0.4256767,0.6646661,0.2886034,0.9879024,-0.06839389,0.6481704,-1,-0.5804082,-0.6514803,0.7308455,0.4101343,-0.8554199,-1,-0.254794,-0.1424363,-0.1296123,0.1170656,-1,0.007876628])
# rands=np.array([0.3490212, 0.2809639, 0.08980897, -0.6552433, 0.3108646, 0.9600278, -0.7445735, -0.929427, 0.6871993, -0.1122855, -1, 0.2775505, 0.7744503, 0.5353921, 0.8993617, 0.5068132, 0.6897388, -0.5868602, -0.5443915, 0.2293534, -0.880945, -1, -0.04769203])
# if member == 7 and lag == 1:
# rands=np.array([-0.4625655,0.6606569,-0.1525902,-0.627054,-0.5481565,-0.7688485,-0.6856943,-0.3053931,-0.03701457,-0.3725091,-1,0.155129,0.136141,-0.7843549,-0.0782493,-0.5983688,0.4498052,0.7143801,-0.1611456,1.007675,0.9826942,-1,0.04915765])
# rands=np.array([0.7321754, 0.06772549, 0.2179533, 0.2919183, 1.062581, 0.7554626, 0.5785612, 1.011077, 0.6109563, -0.350915, -1, -0.579572, -0.07610841, -0.3409401, 0.7316601, -0.904253, -1, -0.01014853, 0.5529011, 0.534495, -0.3398449, -1, -0.007357216])
#if member == 8 and lag == 1:
# rands=np.array([-0.3430965,-0.4227155,-0.07735894,0.1337408,0.3767924,-0.3956357,-0.2084762,0.2861364,0.9596422,0.718588,-1,0.7803586,-0.2091809,0.02144417,-0.269459,0.5019733,0.4792759,0.4641314,1.016747,1.19465,0.3354242,-1,0.02544297])
# rands=np.array([-0.4717625, -0.881941, -0.1655097, 0.2257761, 0.8264135, -0.7091568, 0.2949457, -0.296845, 0.794414, -0.743104, -1, 0.9034545, -0.2101449, 0.3021894, -0.4184734, -0.4638716, 0.6389658, 0.9005892, -0.5002754, 0.6826805, -0.3747995, -1, 0.03224342])
# if member == 9 and lag == 1:
# rands=np.array([-0.7412553,-0.5466392,0.6302599,0.6477938,0.3541165,-0.3394596,0.346593,0.7322681,-0.7179851,-0.731463,-1,0.5184832,-0.1218895,0.7145109,1.139154,0.1859014,-0.6142078,-0.6477512,0.799903,-0.2647903,0.264973,-1,0.01599488])
# rands=np.array([-0.5808461, -0.8636444, 0.7052112, 0.7909082, -0.8140187, -0.02347049, 0.5638158, -0.7550521, 0.5693654, 0.07545193, -1, -0.2417251, 0.6826159, 0.9508564, -0.3956808, 0.8499226, 0.09175592, 0.5388746, 0.8177778, 0.4029195, 0.4111119, -1, -0.01329621])
# if member == 10 and lag == 1:
# rands=np.array([-0.03720611,-0.4162885,0.3607477,0.223355,0.1092008,-0.3264467,-0.721734,-0.6520889,-0.6977112,-0.9723021,-1,-0.693766,-0.4895594,-0.6218721,-0.9569306,-0.009422544,-0.09222607,0.6750132,-0.7697481,-0.6925188,0.6922073,-1,0.04024281])
# rands=np.array([0.2166863, -0.437274, -0.7773396, 0.3759225, -0.7887743, -0.5783228, 0.2769741, -0.8951628, 0.5877202, -0.7473742, -1, 0.8568208, -0.3967594, -0.02381929, 0.911634, 0.1854728, 0.347165, -0.7240146, 0.6684847, -0.05459603, -0.2197302, -1, -0.04627957])
# if member == 11 and lag == 1:
# rands=np.array([-0.5116026,0.343556,-0.9518977,0.4320086,0.7536819,-0.6906458,-0.8697743,-0.5209014,-0.0006806303,-0.7983823,-1,-0.6035937,-0.6045027,0.2035089,-0.01446889,-0.4622703,0.04925987,-0.1290628,-0.8433856,-0.5047902,-0.4363325,-1,0.04609365])
# rands=np.array([-0.5553841, 0.1981317, -0.01706193, 0.3010842, 0.3364336, 0.6311802, -0.4822653, -1, -0.2353798, -0.4258329, -1, -0.3820901, -0.1350132, 0.1199797, -0.4477318, -0.03276956, 0.2006144, 0.7612231, -0.4382411, 0.3140147, -0.7996647, -1, -0.0351413])
# if member == 12 and lag == 1:
# rands=np.array([-0.03203867,0.1992727,0.01177486,-0.8682783,-0.6371214,-0.8578926,-0.5281639,0.4448915,0.05798108,0.4604711,-1,-0.6847421,-1,-0.2765961,-0.4123928,-0.1756036,0.459426,0.303864,-0.6993556,0.1463735,0.06817089,-1,-0.02779672])
# rands=np.array([0.1795571, 0.7117739, 0.3948795, 0.01603693, 0.4342803, 0.3915263, -0.1187998, 0.3347789, -0.1160258, -0.6199519, -1, 0.9331298, 0.2128856, 0.1029978, 0.5622161, -0.03326744, 0.1113134, 0.718336, -0.7506785, 0.8058456, -0.165978, -1, 0.03335277])
# if member == 13 and lag == 1:
# rands=np.array([-0.1434858,0.6530933,0.4572614,1.004458,-0.1688031,0.2899839,-0.4890676,-0.08017457,0.2399375,0.8467381,-1,-0.9818414,-0.7904143,0.6110988,-0.2375934,0.008658743,-0.6596891,0.804915,0.8108714,1.026871,-0.632525,-1,0.008359899])
# rands=np.array([0.1133219, 0.1868301, 0.6434304, 0.06698409, 0.7134196, 0.6549824, -0.3687601, 0.2397162, 0.3933871, -0.8959517, -1, -0.9083285, 0.1290343, -0.8501528, -0.1640191, -0.04551715, 0.2021651, 0.1527653, 0.4235286, 0.7385602, -0.1721378, -1, 0.008510991])
# if member == 14 and lag == 1:
# rands=np.array([0.1951734,0.1356803,-0.7984297,0.6222939,-0.3783839,-0.2109271,-0.5047454,0.7663568,-0.2467663,0.7969002,-1,-0.1039866,-0.0799861,-0.9506351,-0.2622426,-0.5555511,-0.6443623,-0.3867293,0.8174118,-0.2623511,-0.6504129,-1,0.02077365])
# rands=np.array([0.5128248, -0.025835, -0.3690266, 0.5315114, -0.3562845, -0.8088912, 0.6070702, 0.8863306, 0.127363, 0.4851948, -1, 0.2463859, 0.5403649, 0.1009353, -0.3769973, 0.8082039, 0.642301, -0.3246575, -0.8921521, 0.485802, 0.06472221, -1, 0.042198])
#if member == 15 and lag == 1:
# rands=np.array([0.5399131,0.9709167,0.1353998,0.9513429,0.1061087,-0.1906423,0.2546594,0.6244627,0.4253108,-0.5599692,-1,-0.6188933,-0.009267714,0.5136899,0.7686136,-0.7976044,0.6507949,0.2095091,-0.8000021,0.7896537,-0.6757631,-1,-0.04675981])
# rands=np.array([-0.6925536, -0.1888134, 0.7955419, -0.1718084, -0.4925962, -0.63966, 0.712366, 0.7542216, 0.6485603, 0.9584286, -1, 0.2036391, -0.547542, -0.1605409, 0.5345823, -0.8189099, -0.4620013, 0.2116281, 0.5273049, -0.5836464, -0.7646005, -1, 0.03012722])
# if member == 16 and lag == 1:
# rands=np.array([0.2039032,0.4437117,-0.8730499,-0.8436514,-0.3197937,0.1300565,0.2079234,0.07507716,-0.1368318,-0.8760222,-1,0.7775679,0.1503317,-0.6096121,-0.4738224,-0.4769458,0.01603237,-0.5856133,0.1603515,0.639724,0.8044817,-1,-0.03076992])
# rands=np.array([-0.9031112, -0.9747161, 0.6262707, 0.1299344, -0.656873, -1, 0.7236382, 0.9888599, 0.599555, -0.3344736, -1, 0.1587437, 0.9665449, 0.7327783, 0.2025174, 1.045059, 0.7680195, 0.202516, 0.3561991, -0.08263882, -0.7926366, -1, -0.0340635])
# if member == 17 and lag == 1:
# rands=np.array([0.6912341,0.03829452,0.0645803,-0.1688701,-0.06769498,0.5935946,0.794823,-0.1764786,-0.5644639,0.1006729,-1,0.3818233,0.8611205,0.8617951,1.227797,0.06026031,-0.3915816,-0.2597351,0.03994489,0.04319654,0.4627483,-1,0.02341671])
# rands=np.array([-0.4468101, -0.6611752, 0.02045458, -0.818266, -0.06696166, -0.2870016, -0.8372102, -1, 0.6159038, 0.1100874, -1, -0.2941307, -0.8388062, 0.7502818, -0.02831579, -0.4171109, 0.3452238, 0.3286322, -0.1746102, 0.741174, -0.1962416, -1, -0.03235182])
# if member == 18 and lag == 1:
# rands=np.array([0.1631074,-0.7144772,0.5909197,0.8623705,-0.157497,0.5982791,-0.2239446,0.5731422,0.6889938,-0.04528522,-1,0.07149924,-0.6744667,0.4025916,0.9270787,0.1585069,0.8808218,0.5500813,1.050181,0.6608209,-0.4115532,-1,0.0001819044])
# rands=np.array([-0.117545, 0.6358077, -0.8399339, -0.4401485, -0.4653139, -0.1412946, 0.368727, 0.5386308, 0.07112807, 0.2995664, -1, -0.3923609, -0.8697481, 0.283159, -0.1124206, -0.2029478, -1, -0.9128151, 0.2968773, -0.4699179, -0.1645328, -1, -0.04468708])
# if member == 19 and lag == 1:
# rands=np.array([-0.4331279,-0.2030248,-0.2923387,-0.1422065,0.04285587,-0.6101277,-0.5285338,0.820045,0.4597438,-0.7786282,-1,0.3376869,0.1703605,-0.5970429,-0.2520662,-0.6760406,0.3587383,0.1665516,-0.351487,-0.3638605,0.4790523,-1,-0.01560173])
# rands=np.array([-0.3586554, -0.8109571, -0.4370517, 0.1238612, -0.996357, -0.6696278, 0.5044641, 0.6546122, -0.0149131, 0.8158504, -1, 0.8035823, 0.565383, 0.2036764, -0.6208936, 0.4221191, 0.4371626, -0.6936532, -0.6731901, -1, 0.1557615, -1, -0.02435511])
# if member == 20 and lag == 1:
# rands=np.array([0.4845566,0.6727741,0.9107583,-0.2242058,0.8334284,-0.0900023,-0.1698219,-0.1709382,-0.1608661,-0.9621966,-1,0.9589482,-0.5061898,0.02963674,-0.5164318,-0.05424574,-0.8521295,0.7594227,-0.5421007,0.8706058,0.7356571,-1,-0.03166761])
#newmember.param_values = rands + newmean
# rands=np.array([-0.6186742, 0.1736774, 0.3401965, 0.8706512, -0.6920347, -0.2222075, 0.1594849, 0.3773506, 1.039722, 0.3430154, -1, 0.8352696, 1.13694, -0.20392, 0.1833831, 1.04152, 0.1299644, 0.07491443, -0.6444398, 0.2261413, 0.01262438, -1, -0.03116111])
newmember = EnsembleMember(member)
# newmember.param_values = rands + newmean
newmember.param_values = (np.hstack((np.dot(C, rands[0:11]),np.dot(C, rands[11:22]), rands_bg)) + newmean).ravel()
newmember.param_values[10] = 0.
newmember.param_values[21] = 0.
newmember.param_values[newmember.param_values<0.] = 0.
newmember.param_values[newmember.param_values>2.] = 2.
self.ensemble_members[lag].append(newmember)
logging.debug('%d new ensemble members were added to the state vector # %d' % (self.nmembers, (lag + 1)))
......@@ -579,6 +581,7 @@ class StateVector(object):
#data = mem.param_values[0:11]
data = 2.*np.ones(11)-mem.param_values[0:11] # unblock this maybe
print(str(lag),data)
#data = mem.param_values[0:(self.nparams-1)//2]
savedict = io.std_savedict.copy()
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment