Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Add several NCEP grids #239

Merged
merged 12 commits into from
Oct 28, 2024
Merged

Add several NCEP grids #239

Changes from 6 commits
Commits
Show all changes
12 commits
Select commit Hold shift + click to select a range
c3262f8
Add NCEP grids
arcticskew Aug 22, 2024
aa01294
Mixed tab indent fix
arcticskew Aug 22, 2024
62bddc3
Merge branch 'develop' into grids
edwardhartnett Aug 23, 2024
06f1c0c
Merge branch 'develop' into grids
austincrosswx Sep 9, 2024
879b7bf
Merge branch 'develop' into grids
AlysonStahl-NOAA Sep 26, 2024
83a66ed
Merge branch 'develop' into grids
AlysonStahl-NOAA Sep 30, 2024
a298321
Merge branch 'develop' into grids
AlysonStahl-NOAA Oct 4, 2024
16f3c3c
Merge branch 'develop' into grids
AlysonStahl-NOAA Oct 18, 2024
10ffdf5
tests
AlysonStahl-NOAA Oct 18, 2024
64c7550
test data
AlysonStahl-NOAA Oct 18, 2024
636e656
rename files
AlysonStahl-NOAA Oct 18, 2024
881b623
Merge branch 'develop' into grids
AlysonStahl-NOAA Oct 28, 2024
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
82 changes: 65 additions & 17 deletions wgrib2/ncep_grids.c
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -51,6 +51,18 @@ void ncep_grids(const char **arg1, const char **arg2, const char **arg3) {
*arg2 = "244.723:2345:2540";
*arg3 = "19.299:1597:2540";
return;
} /* HRRR 3 km */
if (strcmp(*arg3,"hrrr") == 0) {
*arg1 = "lambert:262.5:38.5";
*arg2 = "237.28:1799:3000";
*arg3 = "21.138:1059:3000";
return;
} /* NAM 12 km */
if (strcmp(*arg3,"nam12") == 0) {
*arg1 = "nps:260:60";
*arg2 = "217.5:1196:12679";
*arg3 = "-2.5:817:12679";
return;
}
if (strcmp(*arg3,"core") == 0) {

Expand All @@ -74,37 +86,37 @@ void ncep_grids(const char **arg1, const char **arg2, const char **arg3) {
else grid_no = atoi(*arg3);

switch(grid_no) {
case 2:
case 2: /* Global 2.5 deg */
*arg1="latlon";
*arg2="0:144:2.5";
*arg3="90:73:-2.5";
break;

case 3:
case 3: /* Global 1 deg */
*arg1="latlon";
*arg2="0:360:1";
*arg3="90:181:-1";
break;

case 4:
case 4: /* Global 0.5 deg */
*arg1="latlon";
*arg2="0:720:0.5";
*arg3="90:361:-0.5";
break;

case 45:
case 45: /* Global 1.25 deg */
*arg1="latlon";
*arg2="0:288:1.25";
*arg3="90:145:-1.25";
break;

case 98:
case 98: /* t62 gaussian grid */
*arg1 = "gaussian";
*arg2 = "0:192:1.875";
*arg3 = "88.5419501372975:94";
break;

case 126:
case 126: /* t126 gaussian grid */
*arg1 = "gaussian";
*arg2 = "0:384:0.9375";
*arg3 = "89.2767128781058:190";
Expand All @@ -128,49 +140,85 @@ void ncep_grids(const char **arg1, const char **arg2, const char **arg3) {
*arg3 = "89.8435135178685:880";
break;

case 130: /* NCEP 130 (RAP) */
*arg1="lambert:265:25";
*arg2="233.862:451:13545.087";
*arg3="16.281:337:13545.087";
break;

case 132: /* NCEP 132 (SREF) */
*arg1="lambert:265:25";
*arg2="214.500:697:16231.74";
*arg3="1.000:553:16231.74";
break;

case 170: /* t170 gaussian grid */
*arg1 = "gaussian";
*arg2 = "0:512:0.703125";
*arg3 = "89.4628215685774:256";
break;

case 173:
case 173: /* Global 0.5 degree */
*arg1 = "latlon";
*arg2 = "0.041666666667:4320:0.083333333333";
*arg3 = "89.95833333333:2160:-0.083333333333";
break;

case 184:
*arg1 = "lambert:265:25";
*arg2 = "238.446:2145:2540";
*arg3 = "20.192:1377:2540";
break;
case 184: /* CONUS NDFD 2.5 km */
*arg1 = "lambert:265:25";
*arg2 = "238.446:2145:2540";
*arg3 = "20.192:1377:2540";
break;

case 187: /* CONUS RTMA 2.5 km */
*arg1 = "lambert:265:25";
*arg2 = "238.446:2145:2540";
*arg3 = "20.192:1597:2540";
break;

case 194:
case 194: /* Puerto Rico NDFD 3 km Mercator */
*arg1 = "mercator:20";
*arg2 = "284.5:544:2500:297.491";
*arg3 = "15:310:2500:22.005";
break;

case 221:
case 196: /* Hawaii NDFD 2.5 km Mercator */
*arg1 = "mercator:20";
*arg2 = "198.475:321:2500:206.131";
*arg3 = "18.073:225:2500:23.088";
break;

case 197: /* CONUS NDFD 5 km */
*arg1="lambert:265:25";
*arg2="238.445999:1073:5079";
*arg3="20.191999:689:5079";
break;

case 221: /* North America 32 km Lambert (NAM, SREF, RAP) */
*arg1="lambert:253:50";
*arg2="214.5:349:32463.41";
*arg3="1:277:32463.41";
break;

case 230:
case 230: /* Global 0.5 deg */
*arg1="latlon";
*arg2="0:720:0.5";
*arg3="90:361:-0.5";
break;

case 242:
case 236: /* CONUS 40 km */
*arg1="lambert:265:25";
*arg2="233.862:151:40635.25";
*arg3="16.281:113:40635.25";
break;

case 242: /* Alaska NDFD 11.25 km */
*arg1="nps:225:60";
*arg2="187:553:11250";
*arg3="30:425:11250";
break;

case 249:
case 249: /* Alaska 10 km polar stereographic */
*arg1="nps:210.0:60.0";
*arg2="188.4:367:9867.89";
*arg3="45.4:343:9867.89";
Expand Down
Loading