Skip to content

Commit e90d40e

Browse files
nshankarnshankar
committed
add pandas import
1 parent 062d5cc commit e90d40e

1 file changed

Lines changed: 76 additions & 3 deletions

File tree

‎check_pygrib_vs_herbie_crs_extraction.ipynb‎

Lines changed: 76 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -15,21 +15,94 @@
1515
},
1616
{
1717
"cell_type": "code",
18-
"execution_count": 2,
18+
"execution_count": 1,
1919
"metadata": {},
2020
"outputs": [],
2121
"source": [
2222
"from herbie import Herbie\n",
2323
"import pygrib\n",
24+
"import pandas as pd\n",
2425
"from pyproj import CRS\n",
2526
"from herbie.crs import get_cf_crs"
2627
]
2728
},
2829
{
2930
"cell_type": "code",
30-
"execution_count": null,
31+
"execution_count": 2,
3132
"metadata": {},
32-
"outputs": [],
33+
"outputs": [
34+
{
35+
"name": "stdout",
36+
"output_type": "stream",
37+
"text": [
38+
"\n",
39+
"Model: HRRR\n",
40+
" pygrib {'a': 6371229, 'b': 6371229, 'lat_0': 38.5, 'lat_1': 38.5, 'lat_2': 38.5, 'lon_0': 262.5, 'proj': 'lcc'}\n",
41+
" Herbie {'a': 6371229, 'b': 6371229, 'lat_0': 38.5, 'lat_1': 38.5, 'lat_2': 38.5, 'lon_0': 262.5, 'proj': 'lcc'}\n",
42+
" equal= True\n",
43+
"\n",
44+
"Model: HRRRAK\n",
45+
" pygrib {'a': 6371229, 'b': 6371229, 'lat_0': 90.0, 'lat_ts': 60.0, 'lon_0': 225.0, 'proj': 'stere'}\n",
46+
" Herbie {'a': 6371229, 'b': 6371229, 'lat_0': 90, 'lat_ts': 60.0, 'lon_0': 225.0, 'proj': 'stere'}\n",
47+
" equal= True\n",
48+
"\n",
49+
"Model: GFS\n",
50+
" pygrib {'a': 6371229, 'b': 6371229, 'proj': 'longlat'}\n",
51+
" Herbie {'a': 6371229, 'b': 6371229, 'proj': 'longlat'}\n",
52+
" equal= True\n",
53+
"\n",
54+
"Model: GRAPHCAST\n",
55+
" pygrib {'a': 4326.0, 'b': 4326.0, 'proj': 'longlat'}\n",
56+
" Herbie {'a': 4326.0, 'b': 4326.0, 'proj': 'longlat'}\n",
57+
" equal= True\n",
58+
"\n",
59+
"Model: CFS\n",
60+
" pygrib {'a': 6371229, 'b': 6371229, 'proj': 'longlat'}\n",
61+
" Herbie {'a': 6371229, 'b': 6371229, 'proj': 'longlat'}\n",
62+
" equal= True\n",
63+
"\n",
64+
"Model: IFS\n",
65+
" pygrib {'a': 6371229, 'b': 6371229, 'proj': 'longlat'}\n",
66+
" Herbie {'a': 6371229, 'b': 6371229, 'proj': 'longlat'}\n",
67+
" equal= True\n",
68+
"\n",
69+
"Model: AIFS\n",
70+
" pygrib {'a': 6371229, 'b': 6371229, 'proj': 'longlat'}\n",
71+
" Herbie {'a': 6371229, 'b': 6371229, 'proj': 'longlat'}\n",
72+
" equal= True\n",
73+
"\n",
74+
"Model: GEFS\n",
75+
" pygrib {'a': 6371229, 'b': 6371229, 'proj': 'longlat'}\n",
76+
" Herbie {'a': 6371229, 'b': 6371229, 'proj': 'longlat'}\n",
77+
" equal= True\n",
78+
"\n",
79+
"Model: GEFS\n",
80+
" pygrib {'a': 6371229, 'b': 6371229, 'proj': 'longlat'}\n",
81+
" Herbie {'a': 6371229, 'b': 6371229, 'proj': 'longlat'}\n",
82+
" equal= True\n",
83+
"\n",
84+
"Model: NAM\n",
85+
" pygrib {'a': 6371229, 'b': 6371229, 'lat_0': 38.5, 'lat_1': 38.5, 'lat_2': 38.5, 'lon_0': 262.5, 'proj': 'lcc'}\n",
86+
" Herbie {'a': 6371229, 'b': 6371229, 'lat_0': 38.5, 'lat_1': 38.5, 'lat_2': 38.5, 'lon_0': 262.5, 'proj': 'lcc'}\n",
87+
" equal= True\n",
88+
"\n",
89+
"Model: URMA\n",
90+
" pygrib {'a': 6371200.0, 'b': 6371200.0, 'lat_0': 25.0, 'lat_1': 25.0, 'lat_2': 25.0, 'lon_0': 265.0, 'proj': 'lcc'}\n",
91+
" Herbie {'a': 6371200.0, 'b': 6371200.0, 'lat_0': 25.0, 'lat_1': 25.0, 'lat_2': 25.0, 'lon_0': 265.0, 'proj': 'lcc'}\n",
92+
" equal= True\n",
93+
"\n",
94+
"Model: RTMA\n",
95+
" pygrib {'a': 6371200.0, 'b': 6371200.0, 'lat_0': 25.0, 'lat_1': 25.0, 'lat_2': 25.0, 'lon_0': 265.0, 'proj': 'lcc'}\n",
96+
" Herbie {'a': 6371200.0, 'b': 6371200.0, 'lat_0': 25.0, 'lat_1': 25.0, 'lat_2': 25.0, 'lon_0': 265.0, 'proj': 'lcc'}\n",
97+
" equal= True\n",
98+
"\n",
99+
"Model: HRDPS\n",
100+
" pygrib {'a': 6371229, 'b': 6371229, 'lon_0': 245.305142, 'o_lat_p': 36.08852, 'o_lon_p': 0.0, 'o_proj': 'longlat', 'proj': 'ob_tran'}\n",
101+
" Herbie {'a': 6371229, 'b': 6371229, 'lon_0': 245.305142, 'o_lat_p': 36.08852, 'o_lon_p': 0.0, 'o_proj': 'longlat', 'proj': 'ob_tran'}\n",
102+
" equal= True\n"
103+
]
104+
}
105+
],
33106
"source": [
34107
"for model, search in [\n",
35108
" (dict(model=\"hrrr\"), \":TMP:2 m above\"),\n",

0 commit comments

Comments
 (0)