-
Notifications
You must be signed in to change notification settings - Fork 2
/
Copy pathprofiles.py
99 lines (93 loc) · 9.24 KB
/
profiles.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
import numpy as np
profile_data = {}
# CGILS CUMULUS
profile_data['cgils_s6'] = {}
profile_data['cgils_s6']['pressure'] = np.array([ 100000., 98994.07 , 94814.59 , 88429.96 , 80311.41 ,
71023.28 , 61186.85 , 51984.22 ,44161.84 , 37512.69 ,31860.8 , 27056.6 , 22972.96 ,
19501.8 , 16551.25 , 14043.24 ,11911.39 , 10099.28 , 8558.965 , 7005.914 ,5311.46 , 3723.029 ,
2394.463 , 1396.721 ,738.8813, 354.4638])
profile_data['cgils_s6']['temperature'] = np.array([ 296.7135, 296.3231, 294.7007, 292.2223, 288.489 , 283.3983,
276.8964, 269.1902, 261.0161, 252.5977, 244.0094, 235.4117,
226.9552, 218.8376, 211.1382, 203.8611, 197.9945, 196.3311,
199.3113, 204.4244, 209.3439, 215.0171, 220.7181, 227.077 ,
233.7758, 237.6904])
profile_data['cgils_s6']['vapor_mixing_ratio'] = np.array([ 1.44507000e-02, 1.47947400e-02, 1.58887600e-02,
1.05742000e-02, 4.72888800e-03, 2.73167200e-03, 2.02892600e-03, 1.35543900e-03, 7.92933200e-04,
4.09263300e-04, 2.04198800e-04, 1.23010300e-04, 7.22973200e-05, 3.76380300e-05, 1.82118600e-05,
8.45097900e-06, 4.26736600e-06, 3.02022000e-06, 2.48533000e-06, 2.48533000e-06, 2.48533000e-06,
2.48533000e-06, 2.48533000e-06, 2.48533000e-06, 2.48533000e-06, 2.48533000e-06])
#CGILS STRATUS
profile_data['cgils_s12'] = {}
profile_data['cgils_s12']['pressure'] = np.array( [100000.0, 98994.07, 94814.59, 88429.96,
80311.41, 71023.28, 61186.85, 51984.22, 44161.84, 37512.69, 31860.8, 27056.6,
22972.96, 19501.8, 16551.25, 14043.24, 11911.39, 10099.28, 8558.965, 7005.914,
5311.46, 3723.029, 2394.463, 1396.721 , 738.8813, 354.4638])
profile_data['cgils_s12']['temperature'] = np.array([ 289.682 , 289.8017, 290.2993, 291.0594, 288.004 , 282.5101,
275.4304, 266.9344, 258.6857, 250.605 , 242.474 , 234.4359, 226.596 , 219.0945, 212.7398, 208.0131,
205.6474, 205.648 , 207.0625, 209.5425, 213.3837, 217.5103, 222.5818, 229.6124, 237.069 , 241.4264])
profile_data['cgils_s12']['vapor_mixing_ratio'] = np.array([ 9.34719900e-03, 9.87366400e-03, 3.20395900e-03,
3.60434200e-03, 3.26695700e-03, 2.57409700e-03, 1.82892500e-03, 1.13080000e-03, 6.36610800e-04,
3.35390600e-04, 1.73921200e-04, 7.88453700e-05, 3.32961700e-05, 1.37285700e-05, 6.20203700e-06,
3.38909300e-06, 2.55585200e-06, 2.48533000e-06, 2.48533000e-06, 2.48533000e-06, 2.48533000e-06,
2.48533000e-06, 2.48533000e-06, 2.48533000e-06, 2.48533000e-06, 2.48533000e-06])
# Arctic profiles
profile_data['arctic'] = {}
profile_data['arctic']['pressure'] = np.array([1.00575873e5, 1.00278284e5, 9.98399536e4, 9.92184082e4,
9.83804260e4, 9.73009582e4, 9.59621155e4, 9.43535889e4,
9.24715210e4, 9.03178955e4, 8.79002075e4, 8.52306702e4,
8.23262024e4, 7.92071045e4, 7.58969238e4, 7.24222778e4,
6.88113770e4, 6.50942261e4, 6.13018860e4, 5.74658386e4,
5.36173767e4, 4.97872192e4, 4.60047852e4, 4.22979370e4,
3.86919434e4, 3.52094727e4, 3.18696625e4, 2.86879181e4,
2.56749146e4, 2.28361450e4, 2.01715500e4, 1.76832245e4,
1.53775101e4, 1.32585281e4, 1.13282722e4, 9.58641815e3,
8.03020325e3, 6.65455856e3, 5.45598755e3, 4.42831612e3,
3.57421494e3, 2.88478088e3, 2.32831478e3, 1.87929726e3,
1.51679745e3, 1.22420835e3, 9.88086700e2, 7.97505760e2,
6.43681290e2, 5.19531150e2, 4.18798970e2, 3.35839870e2,
2.66320660e2, 2.07438110e2, 1.57348970e2, 1.14929540e2,
7.95491400e1, 5.09767600e1, 2.91778800e1, 9.98834000e0])
profile_data['arctic']['temperature'] = np.array([266.8943176, 267.0693359, 267.3087158, 267.5513916, 267.6365662,
267.4140015, 266.8858337, 266.1845093, 265.3601990, 264.5728455,
263.9873962, 263.5353088, 262.6596069, 261.3655701, 259.7712097,
257.8913879, 255.7699127, 253.3448029, 250.6314545, 247.6959229,
244.4879303, 240.9742279, 237.2752991, 233.4430847, 229.6099701,
225.9502258, 222.7632751, 220.5044250, 219.6378784, 220.4424591,
222.1228943, 223.3230896, 223.6871643, 223.5032501, 223.3511658,
223.3572693, 222.9019165, 221.5261230, 219.9896393, 218.5794678,
217.5787201, 217.0293579, 216.6041870, 216.1675110, 215.6414032,
215.1298370, 214.7327423, 215.0572052, 216.5581055, 219.0173950,
221.4697571, 223.0003967, 223.6576843, 224.3420258, 226.3999939,
230.8148956, 237.6799164, 244.1721497, 245.6276703, 238.1793976])
profile_data['arctic']['specific_humidity'] = np.array([2.03999292e-3, 2.07812231e-3, 2.13174864e-3, 2.20352874e-3,
2.26384269e-3, 2.27144509e-3, 2.23518157e-3, 2.15544183e-3,
2.05782512e-3, 1.92552052e-3, 1.77560851e-3, 1.63220366e-3,
1.48773784e-3, 1.32572761e-3, 1.16327132e-3, 1.00650553e-3,
8.66593564e-4, 7.30944929e-4, 6.02148798e-4, 4.82887407e-4,
3.72127570e-4, 2.69140044e-4, 1.90601564e-4, 1.34933990e-4,
9.35404494e-5, 6.32730963e-5, 4.19374412e-5, 2.68939767e-5,
1.80862729e-5, 1.19694567e-5, 6.16606198e-6, 3.18518985e-6,
2.30199470e-6, 2.18249524e-6, 2.18069524e-6, 2.18809521e-6,
2.18759521e-6, 2.18679522e-6, 2.19889516e-6, 2.24759495e-6,
2.32919457e-6, 2.42719411e-6, 2.53789356e-6, 2.63449306e-6,
2.73629251e-6, 2.84559190e-6, 2.98549109e-6, 3.12929021e-6,
3.26818932e-6, 3.40558840e-6, 3.53848748e-6, 3.64988668e-6,
3.75878587e-6, 3.88308492e-6, 3.96388429e-6, 4.00488396e-6,
4.01368389e-6, 3.99578403e-6, 3.96328429e-6, 3.93488452e-6])
profile_data['arctic']['vapor_mixing_ratio'] = profile_data['arctic']['specific_humidity'] / (1.0 -
profile_data['arctic']['specific_humidity'])
profile_data['arctic']['o3_mr'] = np.array([2.29865000e-5, 2.30413000e-5, 2.31245000e-5, 2.32422000e-5,
2.33976000e-5, 2.36059000e-5, 2.38659000e-5, 2.41841000e-5,
2.45592000e-5, 2.49994000e-5, 2.55411000e-5, 2.61804000e-5,
2.68984000e-5, 2.76856000e-5, 2.80097000e-5, 2.82597000e-5,
2.85266000e-5, 2.87814000e-5, 2.90502000e-5, 3.03922000e-5,
3.20686000e-5, 3.36867000e-5, 3.55035000e-5, 3.99821000e-5,
4.48590000e-5, 5.01222000e-5, 6.56444000e-5, 8.59159000e-5,
1.14569600e-4, 1.61802800e-4, 2.21155900e-4, 2.83780500e-4,
3.56959900e-4, 4.57441900e-4, 6.14007900e-4, 8.55658500e-4,
1.17563900e-3, 1.61028430e-3, 2.15759020e-3, 2.84742160e-3,
3.69521780e-3, 4.24513960e-3, 4.68018010e-3, 5.03616810e-3,
5.41168980e-3, 5.82303240e-3, 6.20196440e-3, 6.48212480e-3,
6.64787910e-3, 6.66918370e-3, 6.50225020e-3, 6.11083460e-3,
5.51963280e-3, 4.80870440e-3, 3.99442720e-3, 3.06811450e-3,
2.24217030e-3, 1.52163050e-3, 9.70994100e-4, 4.19090500e-4])