In [1]:
from pygchem import grid
In [2]:
grid.get_supported_models()
Out[2]:
('GEOSFP_REDUCED',
'MERRA_NATIVE',
'GEOS5_REDUCED',
'GEOSFP_47L',
'GEOSFP_NATIVE',
'GENERIC',
'GEOS5_NATIVE',
'GEOS_STRAT_46L',
'GEOS4_30L',
'GEOS',
'GEOS_STRAT',
'MERRA_47L',
'GEOSFP',
'GEOS5',
'GEOS4',
'GEOS3',
'GEOS2',
'GEOS1',
'GEOS57_REDUCED',
'GEOS5_47L',
'GEOS57',
'GEOS2_70L',
'FVDAS',
'MERRA',
'GEOS57_47L',
'GEOS3_30L',
'GEOS3_REDUCED',
'GEOS4_REDUCED',
'MERRA_REDUCED',
'GEOS57_NATIVE')
In [3]:
grid.get_model_info('GEOSFP_47L')
Out[3]:
{'Ap': array([ 0.00000000e+00, 4.80482600e-02, 6.59375200e+00,
1.31348000e+01, 1.96131100e+01, 2.60920100e+01,
3.25708100e+01, 3.89820100e+01, 4.53390100e+01,
5.16961100e+01, 5.80532100e+01, 6.43626400e+01,
7.06219800e+01, 7.88342200e+01, 8.90999200e+01,
9.93652100e+01, 1.09181700e+02, 1.18958600e+02,
1.28695900e+02, 1.42910000e+02, 1.56260000e+02,
1.69609000e+02, 1.81619000e+02, 1.93097000e+02,
2.03259000e+02, 2.12150000e+02, 2.18776000e+02,
2.23898000e+02, 2.24363000e+02, 2.16865000e+02,
2.01192000e+02, 1.76930000e+02, 1.50393000e+02,
1.27837000e+02, 1.08663000e+02, 9.23657200e+01,
7.85123100e+01, 5.63879100e+01, 4.01754100e+01,
2.83678100e+01, 1.97916000e+01, 9.29294200e+00,
4.07657100e+00, 1.65079000e+00, 6.16779100e-01,
2.11349000e-01, 6.60000100e-02, 1.00000000e-02]),
'Bp': array([ 1.00000000e+00, 9.84952000e-01, 9.63406000e-01,
9.41865000e-01, 9.20387000e-01, 8.98908000e-01,
8.77429000e-01, 8.56018000e-01, 8.34660900e-01,
8.13303900e-01, 7.91946900e-01, 7.70637500e-01,
7.49378200e-01, 7.21166000e-01, 6.85899900e-01,
6.50634900e-01, 6.15818400e-01, 5.81041500e-01,
5.46304200e-01, 4.94590200e-01, 4.43740200e-01,
3.92891100e-01, 3.43381100e-01, 2.94403100e-01,
2.46741100e-01, 2.00350100e-01, 1.56224100e-01,
1.13602100e-01, 6.37200600e-02, 2.80100400e-02,
6.96002500e-03, 8.17541300e-09, 0.00000000e+00,
0.00000000e+00, 0.00000000e+00, 0.00000000e+00,
0.00000000e+00, 0.00000000e+00, 0.00000000e+00,
0.00000000e+00, 0.00000000e+00, 0.00000000e+00,
0.00000000e+00, 0.00000000e+00, 0.00000000e+00,
0.00000000e+00, 0.00000000e+00, 0.00000000e+00]),
'Nlayers': 47,
'Ntrop': 38,
'Ptop': 0.01,
'center180': True,
'description': 'GEOS-5.11.x (FP) hybrid reduced',
'halfpolar': True,
'hybrid': True,
'model_family': 'GEOS',
'model_name': 'GEOSFP_47L',
'resolution': (5, 4)}
In [4]:
geos5_47l = grid.grid_from_model('GEOS5_47L')
In [5]:
geos5_47l.resolution
Out[5]:
(5, 4)
In [6]:
geos5_47l.lonlat_centers
Out[6]:
(array([-180., -175., -170., -165., -160., -155., -150., -145., -140.,
-135., -130., -125., -120., -115., -110., -105., -100., -95.,
-90., -85., -80., -75., -70., -65., -60., -55., -50.,
-45., -40., -35., -30., -25., -20., -15., -10., -5.,
0., 5., 10., 15., 20., 25., 30., 35., 40.,
45., 50., 55., 60., 65., 70., 75., 80., 85.,
90., 95., 100., 105., 110., 115., 120., 125., 130.,
135., 140., 145., 150., 155., 160., 165., 170., 175.]),
array([-89., -86., -82., -78., -74., -70., -66., -62., -58., -54., -50.,
-46., -42., -38., -34., -30., -26., -22., -18., -14., -10., -6.,
-2., 2., 6., 10., 14., 18., 22., 26., 30., 34., 38.,
42., 46., 50., 54., 58., 62., 66., 70., 74., 78., 82.,
86., 89.]))
In [7]:
geos5_47l.lonlat_edges
Out[7]:
(array([-182.5, -177.5, -172.5, -167.5, -162.5, -157.5, -152.5, -147.5,
-142.5, -137.5, -132.5, -127.5, -122.5, -117.5, -112.5, -107.5,
-102.5, -97.5, -92.5, -87.5, -82.5, -77.5, -72.5, -67.5,
-62.5, -57.5, -52.5, -47.5, -42.5, -37.5, -32.5, -27.5,
-22.5, -17.5, -12.5, -7.5, -2.5, 2.5, 7.5, 12.5,
17.5, 22.5, 27.5, 32.5, 37.5, 42.5, 47.5, 52.5,
57.5, 62.5, 67.5, 72.5, 77.5, 82.5, 87.5, 92.5,
97.5, 102.5, 107.5, 112.5, 117.5, 122.5, 127.5, 132.5,
137.5, 142.5, 147.5, 152.5, 157.5, 162.5, 167.5, 172.5,
177.5]),
array([-90., -88., -84., -80., -76., -72., -68., -64., -60., -56., -52.,
-48., -44., -40., -36., -32., -28., -24., -20., -16., -12., -8.,
-4., 0., 4., 8., 12., 16., 20., 24., 28., 32., 36.,
40., 44., 48., 52., 56., 60., 64., 68., 72., 76., 80.,
84., 88., 90.]))
In [8]:
geos5_47l.altitude_centers
Out[8]:
array([ 5.84809864e-02, 1.88511025e-01, 3.20422704e-01,
4.53889340e-01, 5.88961308e-01, 7.25693038e-01,
8.64128312e-01, 1.00431865e+00, 1.14632573e+00,
1.29021114e+00, 1.43603600e+00, 1.58386313e+00,
1.75894694e+00, 1.98839405e+00, 2.24937317e+00,
2.51706707e+00, 2.79191944e+00, 3.07444912e+00,
3.43922041e+00, 3.89554534e+00, 4.37451566e+00,
4.87902036e+00, 5.41256618e+00, 5.97953320e+00,
6.58532597e+00, 7.23681814e+00, 7.94307038e+00,
8.84583568e+00, 9.93645629e+00, 1.10208259e+01,
1.20855254e+01, 1.31334999e+01, 1.41696044e+01,
1.51981805e+01, 1.62217433e+01, 1.72425829e+01,
1.87272466e+01, 2.08357014e+01, 2.30204591e+01,
2.53067610e+01, 2.86542829e+01, 3.40243859e+01,
4.01660947e+01, 4.71353891e+01, 5.48335979e+01,
6.30527627e+01, 7.21800704e+01])
In [9]:
geos5_47l.altitude_edges
Out[9]:
array([ -5.79950393e-03, 1.23125821e-01, 2.54275401e-01,
3.86960358e-01, 5.21220339e-01, 6.57116603e-01,
7.94696691e-01, 9.34000598e-01, 1.07509155e+00,
1.21802967e+00, 1.36287809e+00, 1.50969590e+00,
1.65854976e+00, 1.86030574e+00, 2.11806800e+00,
2.38236479e+00, 2.65356628e+00, 2.93219253e+00,
3.21876227e+00, 3.66474047e+00, 4.13202035e+00,
4.62341489e+00, 5.14191462e+00, 5.69160126e+00,
6.27721410e+00, 6.90492866e+00, 7.58245981e+00,
8.32046789e+00, 9.40901582e+00, 1.05036909e+01,
1.15776984e+01, 1.26329969e+01, 1.36736878e+01,
1.47055246e+01, 1.57312199e+01, 1.67530675e+01,
1.77733638e+01, 1.98548954e+01, 2.20035941e+01,
2.42399292e+01, 2.65964320e+01, 3.17162997e+01,
3.75744782e+01, 4.42857808e+01, 5.17877845e+01,
5.99238598e+01, 6.83920984e+01, 8.05806112e+01])
In [10]:
geos5_47l.pressure_centers
Out[10]:
array([ 1.00565033e+03, 9.90407772e+02, 9.75122196e+02,
9.59837374e+02, 9.44552889e+02, 9.29268143e+02,
9.13983996e+02, 8.98700708e+02, 8.83417727e+02,
8.68134847e+02, 8.52852247e+02, 8.37570264e+02,
8.19742555e+02, 7.96821832e+02, 7.71354508e+02,
7.45890358e+02, 7.20429297e+02, 6.94968765e+02,
6.63146075e+02, 6.24966639e+02, 5.86792832e+02,
5.48627903e+02, 5.10475420e+02, 4.72335180e+02,
4.34212079e+02, 3.96112404e+02, 3.58037699e+02,
3.13966339e+02, 2.67086762e+02, 2.26745209e+02,
1.92587127e+02, 1.63661504e+02, 1.39115000e+02,
1.18250000e+02, 1.00514360e+02, 8.54390150e+01,
6.74501100e+01, 4.82816600e+01, 3.42716100e+01,
2.40797050e+01, 1.45422710e+01, 6.68475650e+00,
2.86368050e+00, 1.13378455e+00, 4.14064050e-01,
1.38674505e-01, 3.80000050e-02])
In [11]:
geos5_47l.pressure_edges
Out[11]:
array([ 1.01325000e+03, 9.98050662e+02, 9.82764882e+02,
9.67479511e+02, 9.52195238e+02, 9.36910541e+02,
9.21625744e+02, 9.06342248e+02, 8.91059167e+02,
8.75776287e+02, 8.60493406e+02, 8.45211087e+02,
8.29929441e+02, 8.09555669e+02, 7.84087994e+02,
7.58621022e+02, 7.33159694e+02, 7.07698900e+02,
6.82238631e+02, 6.44053520e+02, 6.05879758e+02,
5.67705907e+02, 5.29549900e+02, 4.91400941e+02,
4.53269420e+02, 4.15154739e+02, 3.77070069e+02,
3.39005328e+02, 2.88927351e+02, 2.45246173e+02,
2.08244245e+02, 1.76930008e+02, 1.50393000e+02,
1.27837000e+02, 1.08663000e+02, 9.23657200e+01,
7.85123100e+01, 5.63879100e+01, 4.01754100e+01,
2.83678100e+01, 1.97916000e+01, 9.29294200e+00,
4.07657100e+00, 1.65079000e+00, 6.16779100e-01,
2.11349000e-01, 6.60000100e-02, 1.00000000e-02])
In [12]:
geos5_47l.eta_centers
Out[12]:
array([ 9.92499636e-01, 9.77456251e-01, 9.62370412e-01,
9.47285317e-01, 9.32200554e-01, 9.17115533e-01,
9.02031105e-01, 8.86947523e-01, 8.71864244e-01,
8.56781065e-01, 8.41698163e-01, 8.26615870e-01,
8.09021116e-01, 7.86399897e-01, 7.61265355e-01,
7.36133945e-01, 7.11005583e-01, 6.85877744e-01,
6.54470881e-01, 6.16790335e-01, 5.79115345e-01,
5.41449117e-01, 5.03795172e-01, 4.66153310e-01,
4.28528364e-01, 3.90926537e-01, 3.53349353e-01,
3.09853874e-01, 2.63586872e-01, 2.23772462e-01,
1.90060723e-01, 1.61513071e-01, 1.37287316e-01,
1.16694959e-01, 9.91910702e-02, 8.43127147e-02,
6.65588706e-02, 4.76408946e-02, 3.38139138e-02,
2.37551863e-02, 1.43423779e-02, 6.58753750e-03,
2.81639148e-03, 1.10910007e-03, 3.98784148e-04,
1.26993116e-04, 2.76341291e-05])
In [13]:
geos5_47l.eta_edges
Out[13]:
array([ 1.00000000e+00, 9.84999272e-01, 9.69913230e-01,
9.54827594e-01, 9.39743040e-01, 9.24658068e-01,
9.09572998e-01, 8.94489211e-01, 8.79405834e-01,
8.64322655e-01, 8.49239476e-01, 8.34156850e-01,
8.19074890e-01, 7.98967342e-01, 7.73832452e-01,
7.48698257e-01, 7.23569632e-01, 6.98441534e-01,
6.73313954e-01, 6.35627808e-01, 5.97952862e-01,
5.60277829e-01, 5.22620405e-01, 4.84969939e-01,
4.47336682e-01, 4.09720045e-01, 3.72133028e-01,
3.34565678e-01, 2.85142070e-01, 2.42031674e-01,
2.05513250e-01, 1.74608196e-01, 1.48417946e-01,
1.26156685e-01, 1.07233232e-01, 9.11489085e-02,
7.74765209e-02, 5.56412202e-02, 3.96405689e-02,
2.79872587e-02, 1.95231140e-02, 9.16164186e-03,
4.01343315e-03, 1.61934981e-03, 5.98850322e-04,
1.98717974e-04, 5.52682583e-05, 0.00000000e+00])
In [14]:
geos5_47l.sigma_centers
In [15]:
geos5_47l.sigma_edges
In [16]:
geosfp_2x25 = grid.grid_from_model('GEOSFP_47L', resolution=(2.5, 2))
In [17]:
geosfp_2x25.lonlat_centers
Out[17]:
(array([-180. , -177.5, -175. , -172.5, -170. , -167.5, -165. , -162.5,
-160. , -157.5, -155. , -152.5, -150. , -147.5, -145. , -142.5,
-140. , -137.5, -135. , -132.5, -130. , -127.5, -125. , -122.5,
-120. , -117.5, -115. , -112.5, -110. , -107.5, -105. , -102.5,
-100. , -97.5, -95. , -92.5, -90. , -87.5, -85. , -82.5,
-80. , -77.5, -75. , -72.5, -70. , -67.5, -65. , -62.5,
-60. , -57.5, -55. , -52.5, -50. , -47.5, -45. , -42.5,
-40. , -37.5, -35. , -32.5, -30. , -27.5, -25. , -22.5,
-20. , -17.5, -15. , -12.5, -10. , -7.5, -5. , -2.5,
0. , 2.5, 5. , 7.5, 10. , 12.5, 15. , 17.5,
20. , 22.5, 25. , 27.5, 30. , 32.5, 35. , 37.5,
40. , 42.5, 45. , 47.5, 50. , 52.5, 55. , 57.5,
60. , 62.5, 65. , 67.5, 70. , 72.5, 75. , 77.5,
80. , 82.5, 85. , 87.5, 90. , 92.5, 95. , 97.5,
100. , 102.5, 105. , 107.5, 110. , 112.5, 115. , 117.5,
120. , 122.5, 125. , 127.5, 130. , 132.5, 135. , 137.5,
140. , 142.5, 145. , 147.5, 150. , 152.5, 155. , 157.5,
160. , 162.5, 165. , 167.5, 170. , 172.5, 175. , 177.5]),
array([-89.5, -88. , -86. , -84. , -82. , -80. , -78. , -76. , -74. ,
-72. , -70. , -68. , -66. , -64. , -62. , -60. , -58. , -56. ,
-54. , -52. , -50. , -48. , -46. , -44. , -42. , -40. , -38. ,
-36. , -34. , -32. , -30. , -28. , -26. , -24. , -22. , -20. ,
-18. , -16. , -14. , -12. , -10. , -8. , -6. , -4. , -2. ,
0. , 2. , 4. , 6. , 8. , 10. , 12. , 14. , 16. ,
18. , 20. , 22. , 24. , 26. , 28. , 30. , 32. , 34. ,
36. , 38. , 40. , 42. , 44. , 46. , 48. , 50. , 52. ,
54. , 56. , 58. , 60. , 62. , 64. , 66. , 68. , 70. ,
72. , 74. , 76. , 78. , 80. , 82. , 84. , 86. , 88. , 89.5]))
In [18]:
geosfp_2x25.lonlat_edges
Out[18]:
(array([-181.25, -178.75, -176.25, -173.75, -171.25, -168.75, -166.25,
-163.75, -161.25, -158.75, -156.25, -153.75, -151.25, -148.75,
-146.25, -143.75, -141.25, -138.75, -136.25, -133.75, -131.25,
-128.75, -126.25, -123.75, -121.25, -118.75, -116.25, -113.75,
-111.25, -108.75, -106.25, -103.75, -101.25, -98.75, -96.25,
-93.75, -91.25, -88.75, -86.25, -83.75, -81.25, -78.75,
-76.25, -73.75, -71.25, -68.75, -66.25, -63.75, -61.25,
-58.75, -56.25, -53.75, -51.25, -48.75, -46.25, -43.75,
-41.25, -38.75, -36.25, -33.75, -31.25, -28.75, -26.25,
-23.75, -21.25, -18.75, -16.25, -13.75, -11.25, -8.75,
-6.25, -3.75, -1.25, 1.25, 3.75, 6.25, 8.75,
11.25, 13.75, 16.25, 18.75, 21.25, 23.75, 26.25,
28.75, 31.25, 33.75, 36.25, 38.75, 41.25, 43.75,
46.25, 48.75, 51.25, 53.75, 56.25, 58.75, 61.25,
63.75, 66.25, 68.75, 71.25, 73.75, 76.25, 78.75,
81.25, 83.75, 86.25, 88.75, 91.25, 93.75, 96.25,
98.75, 101.25, 103.75, 106.25, 108.75, 111.25, 113.75,
116.25, 118.75, 121.25, 123.75, 126.25, 128.75, 131.25,
133.75, 136.25, 138.75, 141.25, 143.75, 146.25, 148.75,
151.25, 153.75, 156.25, 158.75, 161.25, 163.75, 166.25,
168.75, 171.25, 173.75, 176.25, 178.75]),
array([-90., -89., -87., -85., -83., -81., -79., -77., -75., -73., -71.,
-69., -67., -65., -63., -61., -59., -57., -55., -53., -51., -49.,
-47., -45., -43., -41., -39., -37., -35., -33., -31., -29., -27.,
-25., -23., -21., -19., -17., -15., -13., -11., -9., -7., -5.,
-3., -1., 1., 3., 5., 7., 9., 11., 13., 15., 17.,
19., 21., 23., 25., 27., 29., 31., 33., 35., 37., 39.,
41., 43., 45., 47., 49., 51., 53., 55., 57., 59., 61.,
63., 65., 67., 69., 71., 73., 75., 77., 79., 81., 83.,
85., 87., 89., 90.]))
In [ ]:
Content source: benbovy/PyGChem_examples
Similar notebooks: