diff --git a/CMakeLists.txt b/CMakeLists.txt
index a7472af6..ee33723e 100644
--- a/CMakeLists.txt
+++ b/CMakeLists.txt
@@ -17,7 +17,7 @@ option(USE_HDF5 "Use HDF5?" off)
 option(USE_REGEX "Use Regex?" on)
 option(USE_TIGGE "Use tigge?" on)
 option(USE_MYSQL "Use MySQL?" off)
-option(USE_IPOLATES "Use Ipolates" on)
+option(USE_IPOLATES "Use Ipolates" off)
 option(USE_UDF "Use UDF?" off)
 option(USE_OPENMP "Use OpenMP?" on)
 option(USE_PROJ4 "Use Proj4?" off)
diff --git a/tests/CMakeLists.txt b/tests/CMakeLists.txt
index 0f25a28e..d069ec38 100644
--- a/tests/CMakeLists.txt
+++ b/tests/CMakeLists.txt
@@ -106,8 +106,6 @@ copy_test_data(ref_sec0.gdas.t12z.pgrb2.1p00.anl.75r.grib2.txt)
 copy_test_data(ref_sec_len.gdaswave.t00z.wcoast.0p16.f000.grib2.txt)
 copy_test_data(ref_sec_len.simple_packing.grib2.txt)
 copy_test_data(ref_simple_packing.grib2.spread.txt)
-copy_test_data(ref_new_grid_gdas.t12z.pgrb2.1p00.anl.75r.grib2.txt)
-copy_test_data(ref_new_grid_earth_gdas.t12z.pgrb2.1p00.anl.75r.grib2.txt)
 
 # Run these shell tests.
 shell_test(run_wgrib2_tests)
@@ -124,6 +122,35 @@ if (USE_NETCDF3)
 endif()
 if (USE_IPOLATES)
   shell_test(run_ipolates_tests)
+  copy_test_data(ref_new_grid_gdas.t12z.pgrb2.1p00.anl.75r.grib2.txt)
+  copy_test_data(ref_new_grid_earth_gdas.t12z.pgrb2.1p00.anl.75r.grib2.txt)
+  copy_test_data(ref_new_grid_ncep_2.txt)
+  copy_test_data(ref_new_grid_ncep_3.txt)
+  copy_test_data(ref_new_grid_ncep_4.txt)
+  copy_test_data(ref_new_grid_ncep_45.txt)
+  copy_test_data(ref_new_grid_ncep_98.txt)
+  copy_test_data(ref_new_grid_ncep_126.txt)
+  copy_test_data(ref_new_grid_ncep_127.txt)
+  copy_test_data(ref_new_grid_ncep_128.txt)
+  copy_test_data(ref_new_grid_ncep_129.txt)
+  copy_test_data(ref_new_grid_ncep_170.txt)
+  copy_test_data(ref_new_grid_ncep_173.txt)
+  copy_test_data(ref_new_grid_ncep_184.txt)
+  copy_test_data(ref_new_grid_ncep_194.txt)
+  copy_test_data(ref_new_grid_ncep_221.txt)
+  copy_test_data(ref_new_grid_ncep_230.txt)
+  copy_test_data(ref_new_grid_ncep_242.txt)
+  copy_test_data(ref_new_grid_ncep_249.txt)
+  copy_test_data(ref_new_grid_ncep_t62.txt)
+  copy_test_data(ref_new_grid_ncep_t126.txt)
+  copy_test_data(ref_new_grid_ncep_t170.txt)
+  copy_test_data(ref_new_grid_ncep_t190.txt)
+  copy_test_data(ref_new_grid_ncep_t254.txt)
+  copy_test_data(ref_new_grid_ncep_t382.txt)
+  copy_test_data(ref_new_grid_ncep_t574.txt)
+  copy_test_data(ref_new_grid_ncep_t1148.txt)
+  copy_test_data(ref_new_grid_ncep_t1534.txt)
+  copy_test_data(ref_mercator_grid.txt)
 endif()
 if (USE_AEC)
   shell_test(run_wgrib2_aec_tests)
diff --git a/tests/data/ref_mercator_grid.txt b/tests/data/ref_mercator_grid.txt
new file mode 100644
index 00000000..8029eed3
--- /dev/null
+++ b/tests/data/ref_mercator_grid.txt
@@ -0,0 +1,70 @@
+1:0:grid_template=10:winds(N/S):
+	Mercator grid: (361 x 91) LatD 0.000000 input WE:SN output WE:SN res 48
+	lat 0.000000 to 90.000000 by 1.000000 m
+	lon 0.000000 to 360.000000 by 1.000000 m
+	orientation 0.000000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl:
+2:28924:grid_template=10:winds(N/S):
+	Mercator grid: (361 x 91) LatD 0.000000 input WE:SN output WE:SN res 48
+	lat 0.000000 to 90.000000 by 1.000000 m
+	lon 0.000000 to 360.000000 by 1.000000 m
+	orientation 0.000000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl:
+3:57848:grid_template=10:winds(N/S):
+	Mercator grid: (361 x 91) LatD 0.000000 input WE:SN output WE:SN res 48
+	lat 0.000000 to 90.000000 by 1.000000 m
+	lon 0.000000 to 360.000000 by 1.000000 m
+	orientation 0.000000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl:
+4:82666:grid_template=10:winds(N/S):
+	Mercator grid: (361 x 91) LatD 0.000000 input WE:SN output WE:SN res 48
+	lat 0.000000 to 90.000000 by 1.000000 m
+	lon 0.000000 to 360.000000 by 1.000000 m
+	orientation 0.000000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl:
+5:107484:grid_template=10:winds(N/S):
+	Mercator grid: (361 x 91) LatD 0.000000 input WE:SN output WE:SN res 48
+	lat 0.000000 to 90.000000 by 1.000000 m
+	lon 0.000000 to 360.000000 by 1.000000 m
+	orientation 0.000000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl:
+6:132302:grid_template=10:winds(N/S):
+	Mercator grid: (361 x 91) LatD 0.000000 input WE:SN output WE:SN res 48
+	lat 0.000000 to 90.000000 by 1.000000 m
+	lon 0.000000 to 360.000000 by 1.000000 m
+	orientation 0.000000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl:
+7:157120:grid_template=10:winds(N/S):
+	Mercator grid: (361 x 91) LatD 0.000000 input WE:SN output WE:SN res 48
+	lat 0.000000 to 90.000000 by 1.000000 m
+	lon 0.000000 to 360.000000 by 1.000000 m
+	orientation 0.000000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl:
+8:181938:grid_template=10:winds(N/S):
+	Mercator grid: (361 x 91) LatD 0.000000 input WE:SN output WE:SN res 48
+	lat 0.000000 to 90.000000 by 1.000000 m
+	lon 0.000000 to 360.000000 by 1.000000 m
+	orientation 0.000000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl:
+9:206756:grid_template=10:winds(N/S):
+	Mercator grid: (361 x 91) LatD 0.000000 input WE:SN output WE:SN res 48
+	lat 0.000000 to 90.000000 by 1.000000 m
+	lon 0.000000 to 360.000000 by 1.000000 m
+	orientation 0.000000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl:
+10:231574:grid_template=10:winds(N/S):
+	Mercator grid: (361 x 91) LatD 0.000000 input WE:SN output WE:SN res 48
+	lat 0.000000 to 90.000000 by 1.000000 m
+	lon 0.000000 to 360.000000 by 1.000000 m
+	orientation 0.000000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl:
+11:256392:grid_template=10:winds(N/S):
+	Mercator grid: (361 x 91) LatD 0.000000 input WE:SN output WE:SN res 48
+	lat 0.000000 to 90.000000 by 1.000000 m
+	lon 0.000000 to 360.000000 by 1.000000 m
+	orientation 0.000000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl:
+12:281210:grid_template=10:winds(N/S):
+	Mercator grid: (361 x 91) LatD 0.000000 input WE:SN output WE:SN res 48
+	lat 0.000000 to 90.000000 by 1.000000 m
+	lon 0.000000 to 360.000000 by 1.000000 m
+	orientation 0.000000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl:
+13:306028:grid_template=10:winds(N/S):
+	Mercator grid: (361 x 91) LatD 0.000000 input WE:SN output WE:SN res 48
+	lat 0.000000 to 90.000000 by 1.000000 m
+	lon 0.000000 to 360.000000 by 1.000000 m
+	orientation 0.000000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl:
+14:330846:grid_template=10:winds(N/S):
+	Mercator grid: (361 x 91) LatD 0.000000 input WE:SN output WE:SN res 48
+	lat 0.000000 to 90.000000 by 1.000000 m
+	lon 0.000000 to 360.000000 by 1.000000 m
+	orientation 0.000000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl:
diff --git a/tests/data/ref_new_grid_ncep_126.txt b/tests/data/ref_new_grid_ncep_126.txt
new file mode 100644
index 00000000..95a604f8
--- /dev/null
+++ b/tests/data/ref_new_grid_ncep_126.txt
@@ -0,0 +1,98 @@
+1:0:grid_template=40:winds(N/S):
+	Gaussian grid: (384 x 190) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=95 #points=72960
+	lat 89.276713 to -89.276713
+	lon 0.000000 to 359.062500 by 0.937500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89276713 to -89276713 lon=0 to 359062500:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl::lon=10.312500,lat=11.810982,i=41100,ix=12,iy=108,val=-42.503:lon=19.687500,lat=79.841419,i=68758,ix=22,iy=180,val=-17.503
+2:73139:grid_template=40:winds(N/S):
+	Gaussian grid: (384 x 190) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=95 #points=72960
+	lat 89.276713 to -89.276713
+	lon 0.000000 to 359.062500 by 0.937500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89276713 to -89276713 lon=0 to 359062500:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl::lon=10.312500,lat=11.810982,i=41100,ix=12,iy=108,val=-68.2211:lon=19.687500,lat=79.841419,i=68758,ix=22,iy=180,val=-17.2211
+3:146278:grid_template=40:winds(N/S):
+	Gaussian grid: (384 x 190) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=95 #points=72960
+	lat 89.276713 to -89.276713
+	lon 0.000000 to 359.062500 by 0.937500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89276713 to -89276713 lon=0 to 359062500:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl::lon=10.312500,lat=11.810982,i=41100,ix=12,iy=108,val=-27.6889:lon=19.687500,lat=79.841419,i=68758,ix=22,iy=180,val=-13.6889
+4:219417:grid_template=40:winds(N/S):
+	Gaussian grid: (384 x 190) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=95 #points=72960
+	lat 89.276713 to -89.276713
+	lon 0.000000 to 359.062500 by 0.937500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89276713 to -89276713 lon=0 to 359062500:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl::lon=10.312500,lat=11.810982,i=41100,ix=12,iy=108,val=-58.6219:lon=19.687500,lat=79.841419,i=68758,ix=22,iy=180,val=-22.6219
+5:283436:grid_template=40:winds(N/S):
+	Gaussian grid: (384 x 190) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=95 #points=72960
+	lat 89.276713 to -89.276713
+	lon 0.000000 to 359.062500 by 0.937500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89276713 to -89276713 lon=0 to 359062500:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl::lon=10.312500,lat=11.810982,i=41100,ix=12,iy=108,val=-6.00193:lon=19.687500,lat=79.841419,i=68758,ix=22,iy=180,val=0.998074
+6:356575:grid_template=40:winds(N/S):
+	Gaussian grid: (384 x 190) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=95 #points=72960
+	lat 89.276713 to -89.276713
+	lon 0.000000 to 359.062500 by 0.937500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89276713 to -89276713 lon=0 to 359062500:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl::lon=10.312500,lat=11.810982,i=41100,ix=12,iy=108,val=-26.8221:lon=19.687500,lat=79.841419,i=68758,ix=22,iy=180,val=-15.8221
+7:420594:grid_template=40:winds(N/S):
+	Gaussian grid: (384 x 190) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=95 #points=72960
+	lat 89.276713 to -89.276713
+	lon 0.000000 to 359.062500 by 0.937500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89276713 to -89276713 lon=0 to 359062500:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl::lon=10.312500,lat=11.810982,i=41100,ix=12,iy=108,val=-2.67741:lon=19.687500,lat=79.841419,i=68758,ix=22,iy=180,val=10.3226
+8:493733:grid_template=40:winds(N/S):
+	Gaussian grid: (384 x 190) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=95 #points=72960
+	lat 89.276713 to -89.276713
+	lon 0.000000 to 359.062500 by 0.937500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89276713 to -89276713 lon=0 to 359062500:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl::lon=10.312500,lat=11.810982,i=41100,ix=12,iy=108,val=7.26558:lon=19.687500,lat=79.841419,i=68758,ix=22,iy=180,val=-2.73442
+9:557752:grid_template=40:winds(N/S):
+	Gaussian grid: (384 x 190) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=95 #points=72960
+	lat 89.276713 to -89.276713
+	lon 0.000000 to 359.062500 by 0.937500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89276713 to -89276713 lon=0 to 359062500:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl::lon=10.312500,lat=11.810982,i=41100,ix=12,iy=108,val=-10.8744:lon=19.687500,lat=79.841419,i=68758,ix=22,iy=180,val=8.12563
+10:630891:grid_template=40:winds(N/S):
+	Gaussian grid: (384 x 190) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=95 #points=72960
+	lat 89.276713 to -89.276713
+	lon 0.000000 to 359.062500 by 0.937500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89276713 to -89276713 lon=0 to 359062500:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl::lon=10.312500,lat=11.810982,i=41100,ix=12,iy=108,val=11.4455:lon=19.687500,lat=79.841419,i=68758,ix=22,iy=180,val=-0.554543
+11:694910:grid_template=40:winds(N/S):
+	Gaussian grid: (384 x 190) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=95 #points=72960
+	lat 89.276713 to -89.276713
+	lon 0.000000 to 359.062500 by 0.937500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89276713 to -89276713 lon=0 to 359062500:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl::lon=10.312500,lat=11.810982,i=41100,ix=12,iy=108,val=-7.9589:lon=19.687500,lat=79.841419,i=68758,ix=22,iy=180,val=9.0411
+12:768049:grid_template=40:winds(N/S):
+	Gaussian grid: (384 x 190) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=95 #points=72960
+	lat 89.276713 to -89.276713
+	lon 0.000000 to 359.062500 by 0.937500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89276713 to -89276713 lon=0 to 359062500:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl::lon=10.312500,lat=11.810982,i=41100,ix=12,iy=108,val=13.0509:lon=19.687500,lat=79.841419,i=68758,ix=22,iy=180,val=8.05088
+13:832068:grid_template=40:winds(N/S):
+	Gaussian grid: (384 x 190) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=95 #points=72960
+	lat 89.276713 to -89.276713
+	lon 0.000000 to 359.062500 by 0.937500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89276713 to -89276713 lon=0 to 359062500:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl::lon=10.312500,lat=11.810982,i=41100,ix=12,iy=108,val=-2.99871:lon=19.687500,lat=79.841419,i=68758,ix=22,iy=180,val=-7.99871
+14:905207:grid_template=40:winds(N/S):
+	Gaussian grid: (384 x 190) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=95 #points=72960
+	lat 89.276713 to -89.276713
+	lon 0.000000 to 359.062500 by 0.937500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89276713 to -89276713 lon=0 to 359062500:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl::lon=10.312500,lat=11.810982,i=41100,ix=12,iy=108,val=14.8462:lon=19.687500,lat=79.841419,i=68758,ix=22,iy=180,val=12.8462
diff --git a/tests/data/ref_new_grid_ncep_127.txt b/tests/data/ref_new_grid_ncep_127.txt
new file mode 100644
index 00000000..f16a5ffe
--- /dev/null
+++ b/tests/data/ref_new_grid_ncep_127.txt
@@ -0,0 +1,98 @@
+1:0:grid_template=40:winds(N/S):
+	Gaussian grid: (768 x 384) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=192 #points=294912
+	lat 89.641648 to -89.641648
+	lon 0.000000 to 359.531250 by 0.468750
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89641648 to -89641648 lon=0 to 359531250:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl::lon=9.843750,lat=11.937571,i=166678,ix=22,iy=218,val=-41.8533:lon=20.156250,lat=79.817676,i=278060,ix=44,iy=363,val=-17.8533
+2:295091:grid_template=40:winds(N/S):
+	Gaussian grid: (768 x 384) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=192 #points=294912
+	lat 89.641648 to -89.641648
+	lon 0.000000 to 359.531250 by 0.468750
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89641648 to -89641648 lon=0 to 359531250:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl::lon=9.843750,lat=11.937571,i=166678,ix=22,iy=218,val=-65.7432:lon=20.156250,lat=79.817676,i=278060,ix=44,iy=363,val=-16.7432
+3:590182:grid_template=40:winds(N/S):
+	Gaussian grid: (768 x 384) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=192 #points=294912
+	lat 89.641648 to -89.641648
+	lon 0.000000 to 359.531250 by 0.468750
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89641648 to -89641648 lon=0 to 359531250:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl::lon=9.843750,lat=11.937571,i=166678,ix=22,iy=218,val=-27.6892:lon=20.156250,lat=79.817676,i=278060,ix=44,iy=363,val=-14.6892
+4:885273:grid_template=40:winds(N/S):
+	Gaussian grid: (768 x 384) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=192 #points=294912
+	lat 89.641648 to -89.641648
+	lon 0.000000 to 359.531250 by 0.468750
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89641648 to -89641648 lon=0 to 359531250:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl::lon=9.843750,lat=11.937571,i=166678,ix=22,iy=218,val=-57.8751:lon=20.156250,lat=79.817676,i=278060,ix=44,iy=363,val=-21.8751
+5:1180364:grid_template=40:winds(N/S):
+	Gaussian grid: (768 x 384) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=192 #points=294912
+	lat 89.641648 to -89.641648
+	lon 0.000000 to 359.531250 by 0.468750
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89641648 to -89641648 lon=0 to 359531250:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl::lon=9.843750,lat=11.937571,i=166678,ix=22,iy=218,val=-6.00474:lon=20.156250,lat=79.817676,i=278060,ix=44,iy=363,val=0.995262
+6:1475455:grid_template=40:winds(N/S):
+	Gaussian grid: (768 x 384) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=192 #points=294912
+	lat 89.641648 to -89.641648
+	lon 0.000000 to 359.531250 by 0.468750
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89641648 to -89641648 lon=0 to 359531250:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl::lon=9.843750,lat=11.937571,i=166678,ix=22,iy=218,val=-26.9561:lon=20.156250,lat=79.817676,i=278060,ix=44,iy=363,val=-15.9561
+7:1733682:grid_template=40:winds(N/S):
+	Gaussian grid: (768 x 384) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=192 #points=294912
+	lat 89.641648 to -89.641648
+	lon 0.000000 to 359.531250 by 0.468750
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89641648 to -89641648 lon=0 to 359531250:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl::lon=9.843750,lat=11.937571,i=166678,ix=22,iy=218,val=-0.7798:lon=20.156250,lat=79.817676,i=278060,ix=44,iy=363,val=10.2202
+8:2028773:grid_template=40:winds(N/S):
+	Gaussian grid: (768 x 384) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=192 #points=294912
+	lat 89.641648 to -89.641648
+	lon 0.000000 to 359.531250 by 0.468750
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89641648 to -89641648 lon=0 to 359531250:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl::lon=9.843750,lat=11.937571,i=166678,ix=22,iy=218,val=6.15591:lon=20.156250,lat=79.817676,i=278060,ix=44,iy=363,val=-2.84409
+9:2287000:grid_template=40:winds(N/S):
+	Gaussian grid: (768 x 384) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=192 #points=294912
+	lat 89.641648 to -89.641648
+	lon 0.000000 to 359.531250 by 0.468750
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89641648 to -89641648 lon=0 to 359531250:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl::lon=9.843750,lat=11.937571,i=166678,ix=22,iy=218,val=-9.88248:lon=20.156250,lat=79.817676,i=278060,ix=44,iy=363,val=8.11752
+10:2582091:grid_template=40:winds(N/S):
+	Gaussian grid: (768 x 384) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=192 #points=294912
+	lat 89.641648 to -89.641648
+	lon 0.000000 to 359.531250 by 0.468750
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89641648 to -89641648 lon=0 to 359531250:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl::lon=9.843750,lat=11.937571,i=166678,ix=22,iy=218,val=11.0449:lon=20.156250,lat=79.817676,i=278060,ix=44,iy=363,val=-0.955051
+11:2840318:grid_template=40:winds(N/S):
+	Gaussian grid: (768 x 384) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=192 #points=294912
+	lat 89.641648 to -89.641648
+	lon 0.000000 to 359.531250 by 0.468750
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89641648 to -89641648 lon=0 to 359531250:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl::lon=9.843750,lat=11.937571,i=166678,ix=22,iy=218,val=-7.9535:lon=20.156250,lat=79.817676,i=278060,ix=44,iy=363,val=9.0465
+12:3135409:grid_template=40:winds(N/S):
+	Gaussian grid: (768 x 384) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=192 #points=294912
+	lat 89.641648 to -89.641648
+	lon 0.000000 to 359.531250 by 0.468750
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89641648 to -89641648 lon=0 to 359531250:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl::lon=9.843750,lat=11.937571,i=166678,ix=22,iy=218,val=13.6152:lon=20.156250,lat=79.817676,i=278060,ix=44,iy=363,val=8.6152
+13:3393636:grid_template=40:winds(N/S):
+	Gaussian grid: (768 x 384) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=192 #points=294912
+	lat 89.641648 to -89.641648
+	lon 0.000000 to 359.531250 by 0.468750
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89641648 to -89641648 lon=0 to 359531250:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl::lon=9.843750,lat=11.937571,i=166678,ix=22,iy=218,val=-2.99985:lon=20.156250,lat=79.817676,i=278060,ix=44,iy=363,val=-7.99985
+14:3688727:grid_template=40:winds(N/S):
+	Gaussian grid: (768 x 384) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=192 #points=294912
+	lat 89.641648 to -89.641648
+	lon 0.000000 to 359.531250 by 0.468750
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89641648 to -89641648 lon=0 to 359531250:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl::lon=9.843750,lat=11.937571,i=166678,ix=22,iy=218,val=14.4508:lon=20.156250,lat=79.817676,i=278060,ix=44,iy=363,val=13.4508
diff --git a/tests/data/ref_new_grid_ncep_128.txt b/tests/data/ref_new_grid_ncep_128.txt
new file mode 100644
index 00000000..4841c8b0
--- /dev/null
+++ b/tests/data/ref_new_grid_ncep_128.txt
@@ -0,0 +1,98 @@
+1:0:grid_template=40:winds(N/S):
+	Gaussian grid: (1152 x 576) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=288 #points=663552
+	lat 89.760995 to -89.760995
+	lon 0.000000 to 359.687500 by 0.312500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89760995 to -89760995 lon=0 to 359687500:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl::lon=10.000000,lat=12.020811,i=375585,ix=33,iy=327,val=-42.8616:lon=20.000000,lat=80.086607,i=626753,ix=65,iy=545,val=-16.8616
+2:663731:grid_template=40:winds(N/S):
+	Gaussian grid: (1152 x 576) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=288 #points=663552
+	lat 89.760995 to -89.760995
+	lon 0.000000 to 359.687500 by 0.312500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89760995 to -89760995 lon=0 to 359687500:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl::lon=10.000000,lat=12.020811,i=375585,ix=33,iy=327,val=-66.9752:lon=20.000000,lat=80.086607,i=626753,ix=65,iy=545,val=-18.9752
+3:1327462:grid_template=40:winds(N/S):
+	Gaussian grid: (1152 x 576) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=288 #points=663552
+	lat 89.760995 to -89.760995
+	lon 0.000000 to 359.687500 by 0.312500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89760995 to -89760995 lon=0 to 359687500:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl::lon=10.000000,lat=12.020811,i=375585,ix=33,iy=327,val=-28:lon=20.000000,lat=80.086607,i=626753,ix=65,iy=545,val=-14
+4:1991193:grid_template=40:winds(N/S):
+	Gaussian grid: (1152 x 576) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=288 #points=663552
+	lat 89.760995 to -89.760995
+	lon 0.000000 to 359.687500 by 0.312500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89760995 to -89760995 lon=0 to 359687500:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl::lon=10.000000,lat=12.020811,i=375585,ix=33,iy=327,val=-57.9584:lon=20.000000,lat=80.086607,i=626753,ix=65,iy=545,val=-22.9584
+5:2654924:grid_template=40:winds(N/S):
+	Gaussian grid: (1152 x 576) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=288 #points=663552
+	lat 89.760995 to -89.760995
+	lon 0.000000 to 359.687500 by 0.312500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89760995 to -89760995 lon=0 to 359687500:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl::lon=10.000000,lat=12.020811,i=375585,ix=33,iy=327,val=-6.00572:lon=20.000000,lat=80.086607,i=626753,ix=65,iy=545,val=0.994282
+6:3318655:grid_template=40:winds(N/S):
+	Gaussian grid: (1152 x 576) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=288 #points=663552
+	lat 89.760995 to -89.760995
+	lon 0.000000 to 359.687500 by 0.312500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89760995 to -89760995 lon=0 to 359687500:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl::lon=10.000000,lat=12.020811,i=375585,ix=33,iy=327,val=-26.7165:lon=20.000000,lat=80.086607,i=626753,ix=65,iy=545,val=-15.7165
+7:3899442:grid_template=40:winds(N/S):
+	Gaussian grid: (1152 x 576) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=288 #points=663552
+	lat 89.760995 to -89.760995
+	lon 0.000000 to 359.687500 by 0.312500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89760995 to -89760995 lon=0 to 359687500:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl::lon=10.000000,lat=12.020811,i=375585,ix=33,iy=327,val=-0.86166:lon=20.000000,lat=80.086607,i=626753,ix=65,iy=545,val=10.1383
+8:4563173:grid_template=40:winds(N/S):
+	Gaussian grid: (1152 x 576) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=288 #points=663552
+	lat 89.760995 to -89.760995
+	lon 0.000000 to 359.687500 by 0.312500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89760995 to -89760995 lon=0 to 359687500:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl::lon=10.000000,lat=12.020811,i=375585,ix=33,iy=327,val=6.03753:lon=20.000000,lat=80.086607,i=626753,ix=65,iy=545,val=-2.96247
+9:5143960:grid_template=40:winds(N/S):
+	Gaussian grid: (1152 x 576) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=288 #points=663552
+	lat 89.760995 to -89.760995
+	lon 0.000000 to 359.687500 by 0.312500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89760995 to -89760995 lon=0 to 359687500:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl::lon=10.000000,lat=12.020811,i=375585,ix=33,iy=327,val=-9.99627:lon=20.000000,lat=80.086607,i=626753,ix=65,iy=545,val=8.00373
+10:5807691:grid_template=40:winds(N/S):
+	Gaussian grid: (1152 x 576) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=288 #points=663552
+	lat 89.760995 to -89.760995
+	lon 0.000000 to 359.687500 by 0.312500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89760995 to -89760995 lon=0 to 359687500:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl::lon=10.000000,lat=12.020811,i=375585,ix=33,iy=327,val=11.1528:lon=20.000000,lat=80.086607,i=626753,ix=65,iy=545,val=-0.847157
+11:6388478:grid_template=40:winds(N/S):
+	Gaussian grid: (1152 x 576) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=288 #points=663552
+	lat 89.760995 to -89.760995
+	lon 0.000000 to 359.687500 by 0.312500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89760995 to -89760995 lon=0 to 359687500:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl::lon=10.000000,lat=12.020811,i=375585,ix=33,iy=327,val=-7.97073:lon=20.000000,lat=80.086607,i=626753,ix=65,iy=545,val=9.02927
+12:7052209:grid_template=40:winds(N/S):
+	Gaussian grid: (1152 x 576) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=288 #points=663552
+	lat 89.760995 to -89.760995
+	lon 0.000000 to 359.687500 by 0.312500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89760995 to -89760995 lon=0 to 359687500:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl::lon=10.000000,lat=12.020811,i=375585,ix=33,iy=327,val=14.1994:lon=20.000000,lat=80.086607,i=626753,ix=65,iy=545,val=8.19939
+13:7632996:grid_template=40:winds(N/S):
+	Gaussian grid: (1152 x 576) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=288 #points=663552
+	lat 89.760995 to -89.760995
+	lon 0.000000 to 359.687500 by 0.312500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89760995 to -89760995 lon=0 to 359687500:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl::lon=10.000000,lat=12.020811,i=375585,ix=33,iy=327,val=-2.99987:lon=20.000000,lat=80.086607,i=626753,ix=65,iy=545,val=-7.99987
+14:8296727:grid_template=40:winds(N/S):
+	Gaussian grid: (1152 x 576) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=288 #points=663552
+	lat 89.760995 to -89.760995
+	lon 0.000000 to 359.687500 by 0.312500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89760995 to -89760995 lon=0 to 359687500:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl::lon=10.000000,lat=12.020811,i=375585,ix=33,iy=327,val=15.1791:lon=20.000000,lat=80.086607,i=626753,ix=65,iy=545,val=13.1791
diff --git a/tests/data/ref_new_grid_ncep_129.txt b/tests/data/ref_new_grid_ncep_129.txt
new file mode 100644
index 00000000..fb0efc30
--- /dev/null
+++ b/tests/data/ref_new_grid_ncep_129.txt
@@ -0,0 +1,98 @@
+1:0:grid_template=40:winds(N/S):
+	Gaussian grid: (1760 x 880) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=440 #points=1548800
+	lat 89.843514 to -89.843514
+	lon 0.000000 to 359.795455 by 0.204545
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89843514 to -89843514 lon=0 to 359795455:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl::lon=10.022727,lat=11.959112,i=876530,ix=50,iy=499,val=-42.8673:lon=20.045455,lat=80.034019,i=1462659,ix=99,iy=832,val=-16.8673
+2:1548979:grid_template=40:winds(N/S):
+	Gaussian grid: (1760 x 880) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=440 #points=1548800
+	lat 89.843514 to -89.843514
+	lon 0.000000 to 359.795455 by 0.204545
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89843514 to -89843514 lon=0 to 359795455:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl::lon=10.022727,lat=11.959112,i=876530,ix=50,iy=499,val=-66.8333:lon=20.045455,lat=80.034019,i=1462659,ix=99,iy=832,val=-17.8333
+3:3097958:grid_template=40:winds(N/S):
+	Gaussian grid: (1760 x 880) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=440 #points=1548800
+	lat 89.843514 to -89.843514
+	lon 0.000000 to 359.795455 by 0.204545
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89843514 to -89843514 lon=0 to 359795455:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl::lon=10.022727,lat=11.959112,i=876530,ix=50,iy=499,val=-27.9102:lon=20.045455,lat=80.034019,i=1462659,ix=99,iy=832,val=-13.9102
+4:4646937:grid_template=40:winds(N/S):
+	Gaussian grid: (1760 x 880) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=440 #points=1548800
+	lat 89.843514 to -89.843514
+	lon 0.000000 to 359.795455 by 0.204545
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89843514 to -89843514 lon=0 to 359795455:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl::lon=10.022727,lat=11.959112,i=876530,ix=50,iy=499,val=-57.9182:lon=20.045455,lat=80.034019,i=1462659,ix=99,iy=832,val=-22.9182
+5:6195916:grid_template=40:winds(N/S):
+	Gaussian grid: (1760 x 880) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=440 #points=1548800
+	lat 89.843514 to -89.843514
+	lon 0.000000 to 359.795455 by 0.204545
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89843514 to -89843514 lon=0 to 359795455:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl::lon=10.022727,lat=11.959112,i=876530,ix=50,iy=499,val=-6.00612:lon=20.045455,lat=80.034019,i=1462659,ix=99,iy=832,val=0.993881
+6:7744895:grid_template=40:winds(N/S):
+	Gaussian grid: (1760 x 880) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=440 #points=1548800
+	lat 89.843514 to -89.843514
+	lon 0.000000 to 359.795455 by 0.204545
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89843514 to -89843514 lon=0 to 359795455:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl::lon=10.022727,lat=11.959112,i=876530,ix=50,iy=499,val=-26.9514:lon=20.045455,lat=80.034019,i=1462659,ix=99,iy=832,val=-15.9514
+7:9100274:grid_template=40:winds(N/S):
+	Gaussian grid: (1760 x 880) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=440 #points=1548800
+	lat 89.843514 to -89.843514
+	lon 0.000000 to 359.795455 by 0.204545
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89843514 to -89843514 lon=0 to 359795455:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl::lon=10.022727,lat=11.959112,i=876530,ix=50,iy=499,val=-0.930981:lon=20.045455,lat=80.034019,i=1462659,ix=99,iy=832,val=10.069
+8:10649253:grid_template=40:winds(N/S):
+	Gaussian grid: (1760 x 880) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=440 #points=1548800
+	lat 89.843514 to -89.843514
+	lon 0.000000 to 359.795455 by 0.204545
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89843514 to -89843514 lon=0 to 359795455:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl::lon=10.022727,lat=11.959112,i=876530,ix=50,iy=499,val=6.07594:lon=20.045455,lat=80.034019,i=1462659,ix=99,iy=832,val=-2.92406
+9:12004632:grid_template=40:winds(N/S):
+	Gaussian grid: (1760 x 880) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=440 #points=1548800
+	lat 89.843514 to -89.843514
+	lon 0.000000 to 359.795455 by 0.204545
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89843514 to -89843514 lon=0 to 359795455:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl::lon=10.022727,lat=11.959112,i=876530,ix=50,iy=499,val=-9.98491:lon=20.045455,lat=80.034019,i=1462659,ix=99,iy=832,val=8.01509
+10:13553611:grid_template=40:winds(N/S):
+	Gaussian grid: (1760 x 880) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=440 #points=1548800
+	lat 89.843514 to -89.843514
+	lon 0.000000 to 359.795455 by 0.204545
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89843514 to -89843514 lon=0 to 359795455:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl::lon=10.022727,lat=11.959112,i=876530,ix=50,iy=499,val=11.0038:lon=20.045455,lat=80.034019,i=1462659,ix=99,iy=832,val=-0.996151
+11:14908990:grid_template=40:winds(N/S):
+	Gaussian grid: (1760 x 880) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=440 #points=1548800
+	lat 89.843514 to -89.843514
+	lon 0.000000 to 359.795455 by 0.204545
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89843514 to -89843514 lon=0 to 359795455:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl::lon=10.022727,lat=11.959112,i=876530,ix=50,iy=499,val=-7.96378:lon=20.045455,lat=80.034019,i=1462659,ix=99,iy=832,val=9.03622
+12:16457969:grid_template=40:winds(N/S):
+	Gaussian grid: (1760 x 880) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=440 #points=1548800
+	lat 89.843514 to -89.843514
+	lon 0.000000 to 359.795455 by 0.204545
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89843514 to -89843514 lon=0 to 359795455:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl::lon=10.022727,lat=11.959112,i=876530,ix=50,iy=499,val=14.2607:lon=20.045455,lat=80.034019,i=1462659,ix=99,iy=832,val=8.26074
+13:17813348:grid_template=40:winds(N/S):
+	Gaussian grid: (1760 x 880) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=440 #points=1548800
+	lat 89.843514 to -89.843514
+	lon 0.000000 to 359.795455 by 0.204545
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89843514 to -89843514 lon=0 to 359795455:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl::lon=10.022727,lat=11.959112,i=876530,ix=50,iy=499,val=-3:lon=20.045455,lat=80.034019,i=1462659,ix=99,iy=832,val=-8
+14:19362327:grid_template=40:winds(N/S):
+	Gaussian grid: (1760 x 880) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=440 #points=1548800
+	lat 89.843514 to -89.843514
+	lon 0.000000 to 359.795455 by 0.204545
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89843514 to -89843514 lon=0 to 359795455:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl::lon=10.022727,lat=11.959112,i=876530,ix=50,iy=499,val=15.2262:lon=20.045455,lat=80.034019,i=1462659,ix=99,iy=832,val=13.2262
diff --git a/tests/data/ref_new_grid_ncep_170.txt b/tests/data/ref_new_grid_ncep_170.txt
new file mode 100644
index 00000000..29724d08
--- /dev/null
+++ b/tests/data/ref_new_grid_ncep_170.txt
@@ -0,0 +1,98 @@
+1:0:grid_template=40:winds(N/S):
+	Gaussian grid: (512 x 256) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=128 #points=131072
+	lat 89.462822 to -89.462822
+	lon 0.000000 to 359.296875 by 0.703125
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89462822 to -89462822 lon=0 to 359296875:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl::lon=9.843750,lat=12.280678,i=74255,ix=15,iy=146,val=-44.7159:lon=19.687500,lat=80.350237,i=123933,ix=29,iy=243,val=-15.7159
+2:131251:grid_template=40:winds(N/S):
+	Gaussian grid: (512 x 256) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=128 #points=131072
+	lat 89.462822 to -89.462822
+	lon 0.000000 to 359.296875 by 0.703125
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89462822 to -89462822 lon=0 to 359296875:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl::lon=9.843750,lat=12.280678,i=74255,ix=15,iy=146,val=-65.657:lon=19.687500,lat=80.350237,i=123933,ix=29,iy=243,val=-20.657
+3:262502:grid_template=40:winds(N/S):
+	Gaussian grid: (512 x 256) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=128 #points=131072
+	lat 89.462822 to -89.462822
+	lon 0.000000 to 359.296875 by 0.703125
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89462822 to -89462822 lon=0 to 359296875:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl::lon=9.843750,lat=12.280678,i=74255,ix=15,iy=146,val=-28.9228:lon=19.687500,lat=80.350237,i=123933,ix=29,iy=243,val=-11.9228
+4:393753:grid_template=40:winds(N/S):
+	Gaussian grid: (512 x 256) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=128 #points=131072
+	lat 89.462822 to -89.462822
+	lon 0.000000 to 359.296875 by 0.703125
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89462822 to -89462822 lon=0 to 359296875:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl::lon=9.843750,lat=12.280678,i=74255,ix=15,iy=146,val=-56.4385:lon=19.687500,lat=80.350237,i=123933,ix=29,iy=243,val=-24.4385
+5:508620:grid_template=40:winds(N/S):
+	Gaussian grid: (512 x 256) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=128 #points=131072
+	lat 89.462822 to -89.462822
+	lon 0.000000 to 359.296875 by 0.703125
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89462822 to -89462822 lon=0 to 359296875:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl::lon=9.843750,lat=12.280678,i=74255,ix=15,iy=146,val=-6.00541:lon=19.687500,lat=80.350237,i=123933,ix=29,iy=243,val=-0.00540543
+6:639871:grid_template=40:winds(N/S):
+	Gaussian grid: (512 x 256) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=128 #points=131072
+	lat 89.462822 to -89.462822
+	lon 0.000000 to 359.296875 by 0.703125
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89462822 to -89462822 lon=0 to 359296875:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl::lon=9.843750,lat=12.280678,i=74255,ix=15,iy=146,val=-26.4699:lon=19.687500,lat=80.350237,i=123933,ix=29,iy=243,val=-15.4699
+7:754738:grid_template=40:winds(N/S):
+	Gaussian grid: (512 x 256) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=128 #points=131072
+	lat 89.462822 to -89.462822
+	lon 0.000000 to 359.296875 by 0.703125
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89462822 to -89462822 lon=0 to 359296875:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl::lon=9.843750,lat=12.280678,i=74255,ix=15,iy=146,val=0.508083:lon=19.687500,lat=80.350237,i=123933,ix=29,iy=243,val=9.50808
+8:885989:grid_template=40:winds(N/S):
+	Gaussian grid: (512 x 256) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=128 #points=131072
+	lat 89.462822 to -89.462822
+	lon 0.000000 to 359.296875 by 0.703125
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89462822 to -89462822 lon=0 to 359296875:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl::lon=9.843750,lat=12.280678,i=74255,ix=15,iy=146,val=4.24332:lon=19.687500,lat=80.350237,i=123933,ix=29,iy=243,val=-2.75668
+9:1000856:grid_template=40:winds(N/S):
+	Gaussian grid: (512 x 256) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=128 #points=131072
+	lat 89.462822 to -89.462822
+	lon 0.000000 to 359.296875 by 0.703125
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89462822 to -89462822 lon=0 to 359296875:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl::lon=9.843750,lat=12.280678,i=74255,ix=15,iy=146,val=-7.99564:lon=19.687500,lat=80.350237,i=123933,ix=29,iy=243,val=7.00436
+10:1132107:grid_template=40:winds(N/S):
+	Gaussian grid: (512 x 256) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=128 #points=131072
+	lat 89.462822 to -89.462822
+	lon 0.000000 to 359.296875 by 0.703125
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89462822 to -89462822 lon=0 to 359296875:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl::lon=9.843750,lat=12.280678,i=74255,ix=15,iy=146,val=10.1171:lon=19.687500,lat=80.350237,i=123933,ix=29,iy=243,val=-0.882854
+11:1246974:grid_template=40:winds(N/S):
+	Gaussian grid: (512 x 256) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=128 #points=131072
+	lat 89.462822 to -89.462822
+	lon 0.000000 to 359.296875 by 0.703125
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89462822 to -89462822 lon=0 to 359296875:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl::lon=9.843750,lat=12.280678,i=74255,ix=15,iy=146,val=-7.91026:lon=19.687500,lat=80.350237,i=123933,ix=29,iy=243,val=8.08974
+12:1378225:grid_template=40:winds(N/S):
+	Gaussian grid: (512 x 256) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=128 #points=131072
+	lat 89.462822 to -89.462822
+	lon 0.000000 to 359.296875 by 0.703125
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89462822 to -89462822 lon=0 to 359296875:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl::lon=9.843750,lat=12.280678,i=74255,ix=15,iy=146,val=14.5461:lon=19.687500,lat=80.350237,i=123933,ix=29,iy=243,val=7.54612
+13:1493092:grid_template=40:winds(N/S):
+	Gaussian grid: (512 x 256) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=128 #points=131072
+	lat 89.462822 to -89.462822
+	lon 0.000000 to 359.296875 by 0.703125
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89462822 to -89462822 lon=0 to 359296875:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl::lon=9.843750,lat=12.280678,i=74255,ix=15,iy=146,val=-2.99963:lon=19.687500,lat=80.350237,i=123933,ix=29,iy=243,val=-7.99963
+14:1624343:grid_template=40:winds(N/S):
+	Gaussian grid: (512 x 256) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=128 #points=131072
+	lat 89.462822 to -89.462822
+	lon 0.000000 to 359.296875 by 0.703125
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89462822 to -89462822 lon=0 to 359296875:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl::lon=9.843750,lat=12.280678,i=74255,ix=15,iy=146,val=14.3613:lon=19.687500,lat=80.350237,i=123933,ix=29,iy=243,val=12.3613
diff --git a/tests/data/ref_new_grid_ncep_173.txt b/tests/data/ref_new_grid_ncep_173.txt
new file mode 100644
index 00000000..5de7d9ab
--- /dev/null
+++ b/tests/data/ref_new_grid_ncep_173.txt
@@ -0,0 +1,84 @@
+1:0:grid_template=0:winds(N/S):
+	lat-lon grid:(4320 x 2160) units 1e-06 input WE:NS output WE:SN res 48
+	lat 89.958333 to -89.958333 by 0.083333
+	lon 0.041667 to 359.958333 by 0.083333 #points=9331200
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89958333 to -89958333 lon=41667 to 359958333:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl::lon=9.958334,lat=12.041667,i=5287800,ix=120,iy=1225,val=-42.9188:lon=19.958334,lat=80.041666,i=8813040,ix=240,iy=2041,val=-16.9188
+2:9331379:grid_template=0:winds(N/S):
+	lat-lon grid:(4320 x 2160) units 1e-06 input WE:NS output WE:SN res 48
+	lat 89.958333 to -89.958333 by 0.083333
+	lon 0.041667 to 359.958333 by 0.083333 #points=9331200
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89958333 to -89958333 lon=41667 to 359958333:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl::lon=9.958334,lat=12.041667,i=5287800,ix=120,iy=1225,val=-66.9569:lon=19.958334,lat=80.041666,i=8813040,ix=240,iy=2041,val=-17.9569
+3:18662758:grid_template=0:winds(N/S):
+	lat-lon grid:(4320 x 2160) units 1e-06 input WE:NS output WE:SN res 48
+	lat 89.958333 to -89.958333 by 0.083333
+	lon 0.041667 to 359.958333 by 0.083333 #points=9331200
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89958333 to -89958333 lon=41667 to 359958333:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl::lon=9.958334,lat=12.041667,i=5287800,ix=120,iy=1225,val=-27.9589:lon=19.958334,lat=80.041666,i=8813040,ix=240,iy=2041,val=-13.9589
+4:27994137:grid_template=0:winds(N/S):
+	lat-lon grid:(4320 x 2160) units 1e-06 input WE:NS output WE:SN res 48
+	lat 89.958333 to -89.958333 by 0.083333
+	lon 0.041667 to 359.958333 by 0.083333 #points=9331200
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89958333 to -89958333 lon=41667 to 359958333:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl::lon=9.958334,lat=12.041667,i=5287800,ix=120,iy=1225,val=-57.9166:lon=19.958334,lat=80.041666,i=8813040,ix=240,iy=2041,val=-22.9166
+5:37325516:grid_template=0:winds(N/S):
+	lat-lon grid:(4320 x 2160) units 1e-06 input WE:NS output WE:SN res 48
+	lat 89.958333 to -89.958333 by 0.083333
+	lon 0.041667 to 359.958333 by 0.083333 #points=9331200
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89958333 to -89958333 lon=41667 to 359958333:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl::lon=9.958334,lat=12.041667,i=5287800,ix=120,iy=1225,val=-6.00643:lon=19.958334,lat=80.041666,i=8813040,ix=240,iy=2041,val=0.993568
+6:46656895:grid_template=0:winds(N/S):
+	lat-lon grid:(4320 x 2160) units 1e-06 input WE:NS output WE:SN res 48
+	lat 89.958333 to -89.958333 by 0.083333
+	lon 0.041667 to 359.958333 by 0.083333 #points=9331200
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89958333 to -89958333 lon=41667 to 359958333:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl::lon=9.958334,lat=12.041667,i=5287800,ix=120,iy=1225,val=-26.8775:lon=19.958334,lat=80.041666,i=8813040,ix=240,iy=2041,val=-15.8775
+7:54821874:grid_template=0:winds(N/S):
+	lat-lon grid:(4320 x 2160) units 1e-06 input WE:NS output WE:SN res 48
+	lat 89.958333 to -89.958333 by 0.083333
+	lon 0.041667 to 359.958333 by 0.083333 #points=9331200
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89958333 to -89958333 lon=41667 to 359958333:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl::lon=9.958334,lat=12.041667,i=5287800,ix=120,iy=1225,val=-0.918118:lon=19.958334,lat=80.041666,i=8813040,ix=240,iy=2041,val=10.0819
+8:64153253:grid_template=0:winds(N/S):
+	lat-lon grid:(4320 x 2160) units 1e-06 input WE:NS output WE:SN res 48
+	lat 89.958333 to -89.958333 by 0.083333
+	lon 0.041667 to 359.958333 by 0.083333 #points=9331200
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89958333 to -89958333 lon=41667 to 359958333:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl::lon=9.958334,lat=12.041667,i=5287800,ix=120,iy=1225,val=6.03912:lon=19.958334,lat=80.041666,i=8813040,ix=240,iy=2041,val=-2.96088
+9:72318232:grid_template=0:winds(N/S):
+	lat-lon grid:(4320 x 2160) units 1e-06 input WE:NS output WE:SN res 48
+	lat 89.958333 to -89.958333 by 0.083333
+	lon 0.041667 to 359.958333 by 0.083333 #points=9331200
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89958333 to -89958333 lon=41667 to 359958333:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl::lon=9.958334,lat=12.041667,i=5287800,ix=120,iy=1225,val=-9.99788:lon=19.958334,lat=80.041666,i=8813040,ix=240,iy=2041,val=8.00212
+10:81649611:grid_template=0:winds(N/S):
+	lat-lon grid:(4320 x 2160) units 1e-06 input WE:NS output WE:SN res 48
+	lat 89.958333 to -89.958333 by 0.083333
+	lon 0.041667 to 359.958333 by 0.083333 #points=9331200
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89958333 to -89958333 lon=41667 to 359958333:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl::lon=9.958334,lat=12.041667,i=5287800,ix=120,iy=1225,val=11.0354:lon=19.958334,lat=80.041666,i=8813040,ix=240,iy=2041,val=-0.964577
+11:89814590:grid_template=0:winds(N/S):
+	lat-lon grid:(4320 x 2160) units 1e-06 input WE:NS output WE:SN res 48
+	lat 89.958333 to -89.958333 by 0.083333
+	lon 0.041667 to 359.958333 by 0.083333 #points=9331200
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89958333 to -89958333 lon=41667 to 359958333:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl::lon=9.958334,lat=12.041667,i=5287800,ix=120,iy=1225,val=-7.96125:lon=19.958334,lat=80.041666,i=8813040,ix=240,iy=2041,val=9.03875
+12:99145969:grid_template=0:winds(N/S):
+	lat-lon grid:(4320 x 2160) units 1e-06 input WE:NS output WE:SN res 48
+	lat 89.958333 to -89.958333 by 0.083333
+	lon 0.041667 to 359.958333 by 0.083333 #points=9331200
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89958333 to -89958333 lon=41667 to 359958333:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl::lon=9.958334,lat=12.041667,i=5287800,ix=120,iy=1225,val=14.1222:lon=19.958334,lat=80.041666,i=8813040,ix=240,iy=2041,val=8.1222
+13:107310948:grid_template=0:winds(N/S):
+	lat-lon grid:(4320 x 2160) units 1e-06 input WE:NS output WE:SN res 48
+	lat 89.958333 to -89.958333 by 0.083333
+	lon 0.041667 to 359.958333 by 0.083333 #points=9331200
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89958333 to -89958333 lon=41667 to 359958333:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl::lon=9.958334,lat=12.041667,i=5287800,ix=120,iy=1225,val=-2.99997:lon=19.958334,lat=80.041666,i=8813040,ix=240,iy=2041,val=-7.99997
+14:116642327:grid_template=0:winds(N/S):
+	lat-lon grid:(4320 x 2160) units 1e-06 input WE:NS output WE:SN res 48
+	lat 89.958333 to -89.958333 by 0.083333
+	lon 0.041667 to 359.958333 by 0.083333 #points=9331200
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89958333 to -89958333 lon=41667 to 359958333:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl::lon=9.958334,lat=12.041667,i=5287800,ix=120,iy=1225,val=15.1212:lon=19.958334,lat=80.041666,i=8813040,ix=240,iy=2041,val=13.1212
diff --git a/tests/data/ref_new_grid_ncep_184.txt b/tests/data/ref_new_grid_ncep_184.txt
new file mode 100644
index 00000000..640e2664
--- /dev/null
+++ b/tests/data/ref_new_grid_ncep_184.txt
@@ -0,0 +1,84 @@
+1:0:grid_template=30:winds(N/S):
+	Lambert Conformal: (2145 x 1377) input WE:SN output WE:SN res 48
+	Lat1 20.192000 Lon1 238.446000 LoV 265.000000
+	LatD 25.000000 Latin1 25.000000 Latin2 25.000000
+	LatSP -90.000000 LonSP 0.000000
+	North Pole (2145 x 1377) Dx 2540.000000 m Dy 2540.000000 m :12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl:
+2:2584645:grid_template=30:winds(N/S):
+	Lambert Conformal: (2145 x 1377) input WE:SN output WE:SN res 48
+	Lat1 20.192000 Lon1 238.446000 LoV 265.000000
+	LatD 25.000000 Latin1 25.000000 Latin2 25.000000
+	LatSP -90.000000 LonSP 0.000000
+	North Pole (2145 x 1377) Dx 2540.000000 m Dy 2540.000000 m :12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl:
+3:5169290:grid_template=30:winds(N/S):
+	Lambert Conformal: (2145 x 1377) input WE:SN output WE:SN res 48
+	Lat1 20.192000 Lon1 238.446000 LoV 265.000000
+	LatD 25.000000 Latin1 25.000000 Latin2 25.000000
+	LatSP -90.000000 LonSP 0.000000
+	North Pole (2145 x 1377) Dx 2540.000000 m Dy 2540.000000 m :12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl:
+4:7753935:grid_template=30:winds(N/S):
+	Lambert Conformal: (2145 x 1377) input WE:SN output WE:SN res 48
+	Lat1 20.192000 Lon1 238.446000 LoV 265.000000
+	LatD 25.000000 Latin1 25.000000 Latin2 25.000000
+	LatSP -90.000000 LonSP 0.000000
+	North Pole (2145 x 1377) Dx 2540.000000 m Dy 2540.000000 m :12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl:
+5:10338580:grid_template=30:winds(N/S):
+	Lambert Conformal: (2145 x 1377) input WE:SN output WE:SN res 48
+	Lat1 20.192000 Lon1 238.446000 LoV 265.000000
+	LatD 25.000000 Latin1 25.000000 Latin2 25.000000
+	LatSP -90.000000 LonSP 0.000000
+	North Pole (2145 x 1377) Dx 2540.000000 m Dy 2540.000000 m :12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl:
+6:12554017:grid_template=30:winds(N/S):
+	Lambert Conformal: (2145 x 1377) input WE:SN output WE:SN res 48
+	Lat1 20.192000 Lon1 238.446000 LoV 265.000000
+	LatD 25.000000 Latin1 25.000000 Latin2 25.000000
+	LatSP -90.000000 LonSP 0.000000
+	North Pole (2145 x 1377) Dx 2540.000000 m Dy 2540.000000 m :12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl:
+7:15138662:grid_template=30:winds(N/S):
+	Lambert Conformal: (2145 x 1377) input WE:SN output WE:SN res 48
+	Lat1 20.192000 Lon1 238.446000 LoV 265.000000
+	LatD 25.000000 Latin1 25.000000 Latin2 25.000000
+	LatSP -90.000000 LonSP 0.000000
+	North Pole (2145 x 1377) Dx 2540.000000 m Dy 2540.000000 m :12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl:
+8:17723307:grid_template=30:winds(N/S):
+	Lambert Conformal: (2145 x 1377) input WE:SN output WE:SN res 48
+	Lat1 20.192000 Lon1 238.446000 LoV 265.000000
+	LatD 25.000000 Latin1 25.000000 Latin2 25.000000
+	LatSP -90.000000 LonSP 0.000000
+	North Pole (2145 x 1377) Dx 2540.000000 m Dy 2540.000000 m :12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl:
+9:20307952:grid_template=30:winds(N/S):
+	Lambert Conformal: (2145 x 1377) input WE:SN output WE:SN res 48
+	Lat1 20.192000 Lon1 238.446000 LoV 265.000000
+	LatD 25.000000 Latin1 25.000000 Latin2 25.000000
+	LatSP -90.000000 LonSP 0.000000
+	North Pole (2145 x 1377) Dx 2540.000000 m Dy 2540.000000 m :12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl:
+10:22523389:grid_template=30:winds(N/S):
+	Lambert Conformal: (2145 x 1377) input WE:SN output WE:SN res 48
+	Lat1 20.192000 Lon1 238.446000 LoV 265.000000
+	LatD 25.000000 Latin1 25.000000 Latin2 25.000000
+	LatSP -90.000000 LonSP 0.000000
+	North Pole (2145 x 1377) Dx 2540.000000 m Dy 2540.000000 m :12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl:
+11:24738826:grid_template=30:winds(N/S):
+	Lambert Conformal: (2145 x 1377) input WE:SN output WE:SN res 48
+	Lat1 20.192000 Lon1 238.446000 LoV 265.000000
+	LatD 25.000000 Latin1 25.000000 Latin2 25.000000
+	LatSP -90.000000 LonSP 0.000000
+	North Pole (2145 x 1377) Dx 2540.000000 m Dy 2540.000000 m :12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl:
+12:26954263:grid_template=30:winds(N/S):
+	Lambert Conformal: (2145 x 1377) input WE:SN output WE:SN res 48
+	Lat1 20.192000 Lon1 238.446000 LoV 265.000000
+	LatD 25.000000 Latin1 25.000000 Latin2 25.000000
+	LatSP -90.000000 LonSP 0.000000
+	North Pole (2145 x 1377) Dx 2540.000000 m Dy 2540.000000 m :12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl:
+13:29169700:grid_template=30:winds(N/S):
+	Lambert Conformal: (2145 x 1377) input WE:SN output WE:SN res 48
+	Lat1 20.192000 Lon1 238.446000 LoV 265.000000
+	LatD 25.000000 Latin1 25.000000 Latin2 25.000000
+	LatSP -90.000000 LonSP 0.000000
+	North Pole (2145 x 1377) Dx 2540.000000 m Dy 2540.000000 m :12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl:
+14:31015929:grid_template=30:winds(N/S):
+	Lambert Conformal: (2145 x 1377) input WE:SN output WE:SN res 48
+	Lat1 20.192000 Lon1 238.446000 LoV 265.000000
+	LatD 25.000000 Latin1 25.000000 Latin2 25.000000
+	LatSP -90.000000 LonSP 0.000000
+	North Pole (2145 x 1377) Dx 2540.000000 m Dy 2540.000000 m :12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl:
diff --git a/tests/data/ref_new_grid_ncep_194.txt b/tests/data/ref_new_grid_ncep_194.txt
new file mode 100644
index 00000000..6c7a2c51
--- /dev/null
+++ b/tests/data/ref_new_grid_ncep_194.txt
@@ -0,0 +1,70 @@
+1:0:grid_template=10:winds(N/S):
+	Mercator grid: (544 x 310) LatD 20.000000 input WE:SN output WE:SN res 48
+	lat 15.000000 to 22.005000 by 2500.000000 m
+	lon 284.500000 to 297.491000 by 2500.000000 m
+	orientation 0.000000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl:
+2:105579:grid_template=10:winds(N/S):
+	Mercator grid: (544 x 310) LatD 20.000000 input WE:SN output WE:SN res 48
+	lat 15.000000 to 22.005000 by 2500.000000 m
+	lon 284.500000 to 297.491000 by 2500.000000 m
+	orientation 0.000000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl:
+3:232238:grid_template=10:winds(N/S):
+	Mercator grid: (544 x 310) LatD 20.000000 input WE:SN output WE:SN res 48
+	lat 15.000000 to 22.005000 by 2500.000000 m
+	lon 284.500000 to 297.491000 by 2500.000000 m
+	orientation 0.000000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl:
+4:337817:grid_template=10:winds(N/S):
+	Mercator grid: (544 x 310) LatD 20.000000 input WE:SN output WE:SN res 48
+	lat 15.000000 to 22.005000 by 2500.000000 m
+	lon 284.500000 to 297.491000 by 2500.000000 m
+	orientation 0.000000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl:
+5:443396:grid_template=10:winds(N/S):
+	Mercator grid: (544 x 310) LatD 20.000000 input WE:SN output WE:SN res 48
+	lat 15.000000 to 22.005000 by 2500.000000 m
+	lon 284.500000 to 297.491000 by 2500.000000 m
+	orientation 0.000000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl:
+6:548975:grid_template=10:winds(N/S):
+	Mercator grid: (544 x 310) LatD 20.000000 input WE:SN output WE:SN res 48
+	lat 15.000000 to 22.005000 by 2500.000000 m
+	lon 284.500000 to 297.491000 by 2500.000000 m
+	orientation 0.000000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl:
+7:654554:grid_template=10:winds(N/S):
+	Mercator grid: (544 x 310) LatD 20.000000 input WE:SN output WE:SN res 48
+	lat 15.000000 to 22.005000 by 2500.000000 m
+	lon 284.500000 to 297.491000 by 2500.000000 m
+	orientation 0.000000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl:
+8:760133:grid_template=10:winds(N/S):
+	Mercator grid: (544 x 310) LatD 20.000000 input WE:SN output WE:SN res 48
+	lat 15.000000 to 22.005000 by 2500.000000 m
+	lon 284.500000 to 297.491000 by 2500.000000 m
+	orientation 0.000000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl:
+9:865712:grid_template=10:winds(N/S):
+	Mercator grid: (544 x 310) LatD 20.000000 input WE:SN output WE:SN res 48
+	lat 15.000000 to 22.005000 by 2500.000000 m
+	lon 284.500000 to 297.491000 by 2500.000000 m
+	orientation 0.000000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl:
+10:950211:grid_template=10:winds(N/S):
+	Mercator grid: (544 x 310) LatD 20.000000 input WE:SN output WE:SN res 48
+	lat 15.000000 to 22.005000 by 2500.000000 m
+	lon 284.500000 to 297.491000 by 2500.000000 m
+	orientation 0.000000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl:
+11:1034710:grid_template=10:winds(N/S):
+	Mercator grid: (544 x 310) LatD 20.000000 input WE:SN output WE:SN res 48
+	lat 15.000000 to 22.005000 by 2500.000000 m
+	lon 284.500000 to 297.491000 by 2500.000000 m
+	orientation 0.000000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl:
+12:1119209:grid_template=10:winds(N/S):
+	Mercator grid: (544 x 310) LatD 20.000000 input WE:SN output WE:SN res 48
+	lat 15.000000 to 22.005000 by 2500.000000 m
+	lon 284.500000 to 297.491000 by 2500.000000 m
+	orientation 0.000000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl:
+13:1203708:grid_template=10:winds(N/S):
+	Mercator grid: (544 x 310) LatD 20.000000 input WE:SN output WE:SN res 48
+	lat 15.000000 to 22.005000 by 2500.000000 m
+	lon 284.500000 to 297.491000 by 2500.000000 m
+	orientation 0.000000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl:
+14:1267127:grid_template=10:winds(N/S):
+	Mercator grid: (544 x 310) LatD 20.000000 input WE:SN output WE:SN res 48
+	lat 15.000000 to 22.005000 by 2500.000000 m
+	lon 284.500000 to 297.491000 by 2500.000000 m
+	orientation 0.000000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl:
diff --git a/tests/data/ref_new_grid_ncep_2.txt b/tests/data/ref_new_grid_ncep_2.txt
new file mode 100644
index 00000000..09427901
--- /dev/null
+++ b/tests/data/ref_new_grid_ncep_2.txt
@@ -0,0 +1,84 @@
+1:0:grid_template=0:winds(N/S):
+	lat-lon grid:(144 x 73) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 2.500000
+	lon 0.000000 to 357.500000 by 2.500000 #points=10512
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 357500000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl::lon=10.000000,lat=12.500000,i=5909,ix=5,iy=42,val=-46.2517:lon=20.000000,lat=80.000000,i=9801,ix=9,iy=69,val=-17.2517
+2:10691:grid_template=0:winds(N/S):
+	lat-lon grid:(144 x 73) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 2.500000
+	lon 0.000000 to 357.500000 by 2.500000 #points=10512
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 357500000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl::lon=10.000000,lat=12.500000,i=5909,ix=5,iy=42,val=-66.5:lon=20.000000,lat=80.000000,i=9801,ix=9,iy=69,val=-17.5
+3:21382:grid_template=0:winds(N/S):
+	lat-lon grid:(144 x 73) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 2.500000
+	lon 0.000000 to 357.500000 by 2.500000 #points=10512
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 357500000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl::lon=10.000000,lat=12.500000,i=5909,ix=5,iy=42,val=-29:lon=20.000000,lat=80.000000,i=9801,ix=9,iy=69,val=-14
+4:32073:grid_template=0:winds(N/S):
+	lat-lon grid:(144 x 73) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 2.500000
+	lon 0.000000 to 357.500000 by 2.500000 #points=10512
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 357500000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl::lon=10.000000,lat=12.500000,i=5909,ix=5,iy=42,val=-56:lon=20.000000,lat=80.000000,i=9801,ix=9,iy=69,val=-23
+5:41450:grid_template=0:winds(N/S):
+	lat-lon grid:(144 x 73) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 2.500000
+	lon 0.000000 to 357.500000 by 2.500000 #points=10512
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 357500000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl::lon=10.000000,lat=12.500000,i=5909,ix=5,iy=42,val=-6.00665:lon=20.000000,lat=80.000000,i=9801,ix=9,iy=69,val=0.993347
+6:52141:grid_template=0:winds(N/S):
+	lat-lon grid:(144 x 73) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 2.500000
+	lon 0.000000 to 357.500000 by 2.500000 #points=10512
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 357500000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl::lon=10.000000,lat=12.500000,i=5909,ix=5,iy=42,val=-27.0025:lon=20.000000,lat=80.000000,i=9801,ix=9,iy=69,val=-16.0025
+7:61518:grid_template=0:winds(N/S):
+	lat-lon grid:(144 x 73) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 2.500000
+	lon 0.000000 to 357.500000 by 2.500000 #points=10512
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 357500000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl::lon=10.000000,lat=12.500000,i=5909,ix=5,iy=42,val=1:lon=20.000000,lat=80.000000,i=9801,ix=9,iy=69,val=10
+8:72209:grid_template=0:winds(N/S):
+	lat-lon grid:(144 x 73) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 2.500000
+	lon 0.000000 to 357.500000 by 2.500000 #points=10512
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 357500000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl::lon=10.000000,lat=12.500000,i=5909,ix=5,iy=42,val=3.4976:lon=20.000000,lat=80.000000,i=9801,ix=9,iy=69,val=-2.5024
+9:81586:grid_template=0:winds(N/S):
+	lat-lon grid:(144 x 73) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 2.500000
+	lon 0.000000 to 357.500000 by 2.500000 #points=10512
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 357500000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl::lon=10.000000,lat=12.500000,i=5909,ix=5,iy=42,val=-8:lon=20.000000,lat=80.000000,i=9801,ix=9,iy=69,val=8
+10:92277:grid_template=0:winds(N/S):
+	lat-lon grid:(144 x 73) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 2.500000
+	lon 0.000000 to 357.500000 by 2.500000 #points=10512
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 357500000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl::lon=10.000000,lat=12.500000,i=5909,ix=5,iy=42,val=9.99758:lon=20.000000,lat=80.000000,i=9801,ix=9,iy=69,val=-1.00242
+11:101654:grid_template=0:winds(N/S):
+	lat-lon grid:(144 x 73) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 2.500000
+	lon 0.000000 to 357.500000 by 2.500000 #points=10512
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 357500000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl::lon=10.000000,lat=12.500000,i=5909,ix=5,iy=42,val=-8:lon=20.000000,lat=80.000000,i=9801,ix=9,iy=69,val=9
+12:112345:grid_template=0:winds(N/S):
+	lat-lon grid:(144 x 73) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 2.500000
+	lon 0.000000 to 357.500000 by 2.500000 #points=10512
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 357500000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl::lon=10.000000,lat=12.500000,i=5909,ix=5,iy=42,val=15:lon=20.000000,lat=80.000000,i=9801,ix=9,iy=69,val=8
+13:121722:grid_template=0:winds(N/S):
+	lat-lon grid:(144 x 73) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 2.500000
+	lon 0.000000 to 357.500000 by 2.500000 #points=10512
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 357500000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl::lon=10.000000,lat=12.500000,i=5909,ix=5,iy=42,val=-3.99862:lon=20.000000,lat=80.000000,i=9801,ix=9,iy=69,val=-7.99862
+14:132413:grid_template=0:winds(N/S):
+	lat-lon grid:(144 x 73) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 2.500000
+	lon 0.000000 to 357.500000 by 2.500000 #points=10512
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 357500000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl::lon=10.000000,lat=12.500000,i=5909,ix=5,iy=42,val=14.9954:lon=20.000000,lat=80.000000,i=9801,ix=9,iy=69,val=12.9954
diff --git a/tests/data/ref_new_grid_ncep_221.txt b/tests/data/ref_new_grid_ncep_221.txt
new file mode 100644
index 00000000..21e8648c
--- /dev/null
+++ b/tests/data/ref_new_grid_ncep_221.txt
@@ -0,0 +1,84 @@
+1:0:grid_template=30:winds(N/S):
+	Lambert Conformal: (349 x 277) input WE:SN output WE:SN res 48
+	Lat1 1.000000 Lon1 214.500000 LoV 253.000000
+	LatD 50.000000 Latin1 50.000000 Latin2 50.000000
+	LatSP -90.000000 LonSP 0.000000
+	North Pole (349 x 277) Dx 32463.410000 m Dy 32463.410000 m :12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl:
+2:84777:grid_template=30:winds(N/S):
+	Lambert Conformal: (349 x 277) input WE:SN output WE:SN res 48
+	Lat1 1.000000 Lon1 214.500000 LoV 253.000000
+	LatD 50.000000 Latin1 50.000000 Latin2 50.000000
+	LatSP -90.000000 LonSP 0.000000
+	North Pole (349 x 277) Dx 32463.410000 m Dy 32463.410000 m :12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl:
+3:169554:grid_template=30:winds(N/S):
+	Lambert Conformal: (349 x 277) input WE:SN output WE:SN res 48
+	Lat1 1.000000 Lon1 214.500000 LoV 253.000000
+	LatD 50.000000 Latin1 50.000000 Latin2 50.000000
+	LatSP -90.000000 LonSP 0.000000
+	North Pole (349 x 277) Dx 32463.410000 m Dy 32463.410000 m :12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl:
+4:254331:grid_template=30:winds(N/S):
+	Lambert Conformal: (349 x 277) input WE:SN output WE:SN res 48
+	Lat1 1.000000 Lon1 214.500000 LoV 253.000000
+	LatD 50.000000 Latin1 50.000000 Latin2 50.000000
+	LatSP -90.000000 LonSP 0.000000
+	North Pole (349 x 277) Dx 32463.410000 m Dy 32463.410000 m :12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl:
+5:339108:grid_template=30:winds(N/S):
+	Lambert Conformal: (349 x 277) input WE:SN output WE:SN res 48
+	Lat1 1.000000 Lon1 214.500000 LoV 253.000000
+	LatD 50.000000 Latin1 50.000000 Latin2 50.000000
+	LatSP -90.000000 LonSP 0.000000
+	North Pole (349 x 277) Dx 32463.410000 m Dy 32463.410000 m :12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl:
+6:423885:grid_template=30:winds(N/S):
+	Lambert Conformal: (349 x 277) input WE:SN output WE:SN res 48
+	Lat1 1.000000 Lon1 214.500000 LoV 253.000000
+	LatD 50.000000 Latin1 50.000000 Latin2 50.000000
+	LatSP -90.000000 LonSP 0.000000
+	North Pole (349 x 277) Dx 32463.410000 m Dy 32463.410000 m :12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl:
+7:508662:grid_template=30:winds(N/S):
+	Lambert Conformal: (349 x 277) input WE:SN output WE:SN res 48
+	Lat1 1.000000 Lon1 214.500000 LoV 253.000000
+	LatD 50.000000 Latin1 50.000000 Latin2 50.000000
+	LatSP -90.000000 LonSP 0.000000
+	North Pole (349 x 277) Dx 32463.410000 m Dy 32463.410000 m :12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl:
+8:593439:grid_template=30:winds(N/S):
+	Lambert Conformal: (349 x 277) input WE:SN output WE:SN res 48
+	Lat1 1.000000 Lon1 214.500000 LoV 253.000000
+	LatD 50.000000 Latin1 50.000000 Latin2 50.000000
+	LatSP -90.000000 LonSP 0.000000
+	North Pole (349 x 277) Dx 32463.410000 m Dy 32463.410000 m :12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl:
+9:678216:grid_template=30:winds(N/S):
+	Lambert Conformal: (349 x 277) input WE:SN output WE:SN res 48
+	Lat1 1.000000 Lon1 214.500000 LoV 253.000000
+	LatD 50.000000 Latin1 50.000000 Latin2 50.000000
+	LatSP -90.000000 LonSP 0.000000
+	North Pole (349 x 277) Dx 32463.410000 m Dy 32463.410000 m :12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl:
+10:762993:grid_template=30:winds(N/S):
+	Lambert Conformal: (349 x 277) input WE:SN output WE:SN res 48
+	Lat1 1.000000 Lon1 214.500000 LoV 253.000000
+	LatD 50.000000 Latin1 50.000000 Latin2 50.000000
+	LatSP -90.000000 LonSP 0.000000
+	North Pole (349 x 277) Dx 32463.410000 m Dy 32463.410000 m :12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl:
+11:847770:grid_template=30:winds(N/S):
+	Lambert Conformal: (349 x 277) input WE:SN output WE:SN res 48
+	Lat1 1.000000 Lon1 214.500000 LoV 253.000000
+	LatD 50.000000 Latin1 50.000000 Latin2 50.000000
+	LatSP -90.000000 LonSP 0.000000
+	North Pole (349 x 277) Dx 32463.410000 m Dy 32463.410000 m :12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl:
+12:920463:grid_template=30:winds(N/S):
+	Lambert Conformal: (349 x 277) input WE:SN output WE:SN res 48
+	Lat1 1.000000 Lon1 214.500000 LoV 253.000000
+	LatD 50.000000 Latin1 50.000000 Latin2 50.000000
+	LatSP -90.000000 LonSP 0.000000
+	North Pole (349 x 277) Dx 32463.410000 m Dy 32463.410000 m :12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl:
+13:993156:grid_template=30:winds(N/S):
+	Lambert Conformal: (349 x 277) input WE:SN output WE:SN res 48
+	Lat1 1.000000 Lon1 214.500000 LoV 253.000000
+	LatD 50.000000 Latin1 50.000000 Latin2 50.000000
+	LatSP -90.000000 LonSP 0.000000
+	North Pole (349 x 277) Dx 32463.410000 m Dy 32463.410000 m :12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl:
+14:1065849:grid_template=30:winds(N/S):
+	Lambert Conformal: (349 x 277) input WE:SN output WE:SN res 48
+	Lat1 1.000000 Lon1 214.500000 LoV 253.000000
+	LatD 50.000000 Latin1 50.000000 Latin2 50.000000
+	LatSP -90.000000 LonSP 0.000000
+	North Pole (349 x 277) Dx 32463.410000 m Dy 32463.410000 m :12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl:
diff --git a/tests/data/ref_new_grid_ncep_230.txt b/tests/data/ref_new_grid_ncep_230.txt
new file mode 100644
index 00000000..a3a1c128
--- /dev/null
+++ b/tests/data/ref_new_grid_ncep_230.txt
@@ -0,0 +1,84 @@
+1:0:grid_template=0:winds(N/S):
+	lat-lon grid:(720 x 361) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 0.500000
+	lon 0.000000 to 359.500000 by 0.500000 #points=259920
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 359500000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl:
+2:260099:grid_template=0:winds(N/S):
+	lat-lon grid:(720 x 361) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 0.500000
+	lon 0.000000 to 359.500000 by 0.500000 #points=259920
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 359500000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl:
+3:520198:grid_template=0:winds(N/S):
+	lat-lon grid:(720 x 361) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 0.500000
+	lon 0.000000 to 359.500000 by 0.500000 #points=259920
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 359500000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl:
+4:780297:grid_template=0:winds(N/S):
+	lat-lon grid:(720 x 361) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 0.500000
+	lon 0.000000 to 359.500000 by 0.500000 #points=259920
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 359500000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl:
+5:1040396:grid_template=0:winds(N/S):
+	lat-lon grid:(720 x 361) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 0.500000
+	lon 0.000000 to 359.500000 by 0.500000 #points=259920
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 359500000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl:
+6:1300495:grid_template=0:winds(N/S):
+	lat-lon grid:(720 x 361) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 0.500000
+	lon 0.000000 to 359.500000 by 0.500000 #points=259920
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 359500000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl:
+7:1528104:grid_template=0:winds(N/S):
+	lat-lon grid:(720 x 361) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 0.500000
+	lon 0.000000 to 359.500000 by 0.500000 #points=259920
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 359500000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl:
+8:1788203:grid_template=0:winds(N/S):
+	lat-lon grid:(720 x 361) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 0.500000
+	lon 0.000000 to 359.500000 by 0.500000 #points=259920
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 359500000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl:
+9:2015812:grid_template=0:winds(N/S):
+	lat-lon grid:(720 x 361) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 0.500000
+	lon 0.000000 to 359.500000 by 0.500000 #points=259920
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 359500000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl:
+10:2275911:grid_template=0:winds(N/S):
+	lat-lon grid:(720 x 361) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 0.500000
+	lon 0.000000 to 359.500000 by 0.500000 #points=259920
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 359500000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl:
+11:2503520:grid_template=0:winds(N/S):
+	lat-lon grid:(720 x 361) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 0.500000
+	lon 0.000000 to 359.500000 by 0.500000 #points=259920
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 359500000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl:
+12:2763619:grid_template=0:winds(N/S):
+	lat-lon grid:(720 x 361) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 0.500000
+	lon 0.000000 to 359.500000 by 0.500000 #points=259920
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 359500000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl:
+13:2991228:grid_template=0:winds(N/S):
+	lat-lon grid:(720 x 361) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 0.500000
+	lon 0.000000 to 359.500000 by 0.500000 #points=259920
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 359500000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl:
+14:3251327:grid_template=0:winds(N/S):
+	lat-lon grid:(720 x 361) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 0.500000
+	lon 0.000000 to 359.500000 by 0.500000 #points=259920
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 359500000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl:
diff --git a/tests/data/ref_new_grid_ncep_242.txt b/tests/data/ref_new_grid_ncep_242.txt
new file mode 100644
index 00000000..8ac8bbbb
--- /dev/null
+++ b/tests/data/ref_new_grid_ncep_242.txt
@@ -0,0 +1,42 @@
+1:0:grid_template=20:winds(N/S):
+	polar stereographic grid: (553 x 425) input WE:SN output WE:SN res 48
+	North pole lat1 30.000000 lon1 187.000000 latD 60.000000 lonV 225.000000 dx 11250.000000 m dy 11250.000000 m:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl:
+2:205819:grid_template=20:winds(N/S):
+	polar stereographic grid: (553 x 425) input WE:SN output WE:SN res 48
+	North pole lat1 30.000000 lon1 187.000000 latD 60.000000 lonV 225.000000 dx 11250.000000 m dy 11250.000000 m:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl:
+3:411638:grid_template=20:winds(N/S):
+	polar stereographic grid: (553 x 425) input WE:SN output WE:SN res 48
+	North pole lat1 30.000000 lon1 187.000000 latD 60.000000 lonV 225.000000 dx 11250.000000 m dy 11250.000000 m:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl:
+4:617457:grid_template=20:winds(N/S):
+	polar stereographic grid: (553 x 425) input WE:SN output WE:SN res 48
+	North pole lat1 30.000000 lon1 187.000000 latD 60.000000 lonV 225.000000 dx 11250.000000 m dy 11250.000000 m:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl:
+5:823276:grid_template=20:winds(N/S):
+	polar stereographic grid: (553 x 425) input WE:SN output WE:SN res 48
+	North pole lat1 30.000000 lon1 187.000000 latD 60.000000 lonV 225.000000 dx 11250.000000 m dy 11250.000000 m:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl:
+6:999717:grid_template=20:winds(N/S):
+	polar stereographic grid: (553 x 425) input WE:SN output WE:SN res 48
+	North pole lat1 30.000000 lon1 187.000000 latD 60.000000 lonV 225.000000 dx 11250.000000 m dy 11250.000000 m:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl:
+7:1176158:grid_template=20:winds(N/S):
+	polar stereographic grid: (553 x 425) input WE:SN output WE:SN res 48
+	North pole lat1 30.000000 lon1 187.000000 latD 60.000000 lonV 225.000000 dx 11250.000000 m dy 11250.000000 m:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl:
+8:1352599:grid_template=20:winds(N/S):
+	polar stereographic grid: (553 x 425) input WE:SN output WE:SN res 48
+	North pole lat1 30.000000 lon1 187.000000 latD 60.000000 lonV 225.000000 dx 11250.000000 m dy 11250.000000 m:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl:
+9:1558418:grid_template=20:winds(N/S):
+	polar stereographic grid: (553 x 425) input WE:SN output WE:SN res 48
+	North pole lat1 30.000000 lon1 187.000000 latD 60.000000 lonV 225.000000 dx 11250.000000 m dy 11250.000000 m:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl:
+10:1734859:grid_template=20:winds(N/S):
+	polar stereographic grid: (553 x 425) input WE:SN output WE:SN res 48
+	North pole lat1 30.000000 lon1 187.000000 latD 60.000000 lonV 225.000000 dx 11250.000000 m dy 11250.000000 m:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl:
+11:1940678:grid_template=20:winds(N/S):
+	polar stereographic grid: (553 x 425) input WE:SN output WE:SN res 48
+	North pole lat1 30.000000 lon1 187.000000 latD 60.000000 lonV 225.000000 dx 11250.000000 m dy 11250.000000 m:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl:
+12:2117119:grid_template=20:winds(N/S):
+	polar stereographic grid: (553 x 425) input WE:SN output WE:SN res 48
+	North pole lat1 30.000000 lon1 187.000000 latD 60.000000 lonV 225.000000 dx 11250.000000 m dy 11250.000000 m:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl:
+13:2293560:grid_template=20:winds(N/S):
+	polar stereographic grid: (553 x 425) input WE:SN output WE:SN res 48
+	North pole lat1 30.000000 lon1 187.000000 latD 60.000000 lonV 225.000000 dx 11250.000000 m dy 11250.000000 m:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl:
+14:2470001:grid_template=20:winds(N/S):
+	polar stereographic grid: (553 x 425) input WE:SN output WE:SN res 48
+	North pole lat1 30.000000 lon1 187.000000 latD 60.000000 lonV 225.000000 dx 11250.000000 m dy 11250.000000 m:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl:
diff --git a/tests/data/ref_new_grid_ncep_249.txt b/tests/data/ref_new_grid_ncep_249.txt
new file mode 100644
index 00000000..23bda81a
--- /dev/null
+++ b/tests/data/ref_new_grid_ncep_249.txt
@@ -0,0 +1,42 @@
+1:0:grid_template=20:winds(N/S):
+	polar stereographic grid: (367 x 343) input WE:SN output WE:SN res 48
+	North pole lat1 45.400000 lon1 188.400000 latD 60.000000 lonV 210.000000 dx 9867.890000 m dy 9867.890000 m:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl:
+2:110318:grid_template=20:winds(N/S):
+	polar stereographic grid: (367 x 343) input WE:SN output WE:SN res 48
+	North pole lat1 45.400000 lon1 188.400000 latD 60.000000 lonV 210.000000 dx 9867.890000 m dy 9867.890000 m:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl:
+3:220636:grid_template=20:winds(N/S):
+	polar stereographic grid: (367 x 343) input WE:SN output WE:SN res 48
+	North pole lat1 45.400000 lon1 188.400000 latD 60.000000 lonV 210.000000 dx 9867.890000 m dy 9867.890000 m:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl:
+4:330954:grid_template=20:winds(N/S):
+	polar stereographic grid: (367 x 343) input WE:SN output WE:SN res 48
+	North pole lat1 45.400000 lon1 188.400000 latD 60.000000 lonV 210.000000 dx 9867.890000 m dy 9867.890000 m:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl:
+5:441272:grid_template=20:winds(N/S):
+	polar stereographic grid: (367 x 343) input WE:SN output WE:SN res 48
+	North pole lat1 45.400000 lon1 188.400000 latD 60.000000 lonV 210.000000 dx 9867.890000 m dy 9867.890000 m:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl:
+6:535855:grid_template=20:winds(N/S):
+	polar stereographic grid: (367 x 343) input WE:SN output WE:SN res 48
+	North pole lat1 45.400000 lon1 188.400000 latD 60.000000 lonV 210.000000 dx 9867.890000 m dy 9867.890000 m:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl:
+7:630438:grid_template=20:winds(N/S):
+	polar stereographic grid: (367 x 343) input WE:SN output WE:SN res 48
+	North pole lat1 45.400000 lon1 188.400000 latD 60.000000 lonV 210.000000 dx 9867.890000 m dy 9867.890000 m:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl:
+8:725021:grid_template=20:winds(N/S):
+	polar stereographic grid: (367 x 343) input WE:SN output WE:SN res 48
+	North pole lat1 45.400000 lon1 188.400000 latD 60.000000 lonV 210.000000 dx 9867.890000 m dy 9867.890000 m:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl:
+9:819604:grid_template=20:winds(N/S):
+	polar stereographic grid: (367 x 343) input WE:SN output WE:SN res 48
+	North pole lat1 45.400000 lon1 188.400000 latD 60.000000 lonV 210.000000 dx 9867.890000 m dy 9867.890000 m:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl:
+10:914187:grid_template=20:winds(N/S):
+	polar stereographic grid: (367 x 343) input WE:SN output WE:SN res 48
+	North pole lat1 45.400000 lon1 188.400000 latD 60.000000 lonV 210.000000 dx 9867.890000 m dy 9867.890000 m:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl:
+11:1008770:grid_template=20:winds(N/S):
+	polar stereographic grid: (367 x 343) input WE:SN output WE:SN res 48
+	North pole lat1 45.400000 lon1 188.400000 latD 60.000000 lonV 210.000000 dx 9867.890000 m dy 9867.890000 m:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl:
+12:1103353:grid_template=20:winds(N/S):
+	polar stereographic grid: (367 x 343) input WE:SN output WE:SN res 48
+	North pole lat1 45.400000 lon1 188.400000 latD 60.000000 lonV 210.000000 dx 9867.890000 m dy 9867.890000 m:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl:
+13:1197936:grid_template=20:winds(N/S):
+	polar stereographic grid: (367 x 343) input WE:SN output WE:SN res 48
+	North pole lat1 45.400000 lon1 188.400000 latD 60.000000 lonV 210.000000 dx 9867.890000 m dy 9867.890000 m:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl:
+14:1276784:grid_template=20:winds(N/S):
+	polar stereographic grid: (367 x 343) input WE:SN output WE:SN res 48
+	North pole lat1 45.400000 lon1 188.400000 latD 60.000000 lonV 210.000000 dx 9867.890000 m dy 9867.890000 m:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl:
diff --git a/tests/data/ref_new_grid_ncep_3.txt b/tests/data/ref_new_grid_ncep_3.txt
new file mode 100644
index 00000000..d97bba80
--- /dev/null
+++ b/tests/data/ref_new_grid_ncep_3.txt
@@ -0,0 +1,84 @@
+1:0:grid_template=0:winds(N/S):
+	lat-lon grid:(360 x 181) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 1.000000
+	lon 0.000000 to 359.000000 by 1.000000 #points=65160
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 359000000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl::lon=10.000000,lat=12.000000,i=36731,ix=11,iy=103,val=-43:lon=20.000000,lat=80.000000,i=61221,ix=21,iy=171,val=-17
+2:65339:grid_template=0:winds(N/S):
+	lat-lon grid:(360 x 181) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 1.000000
+	lon 0.000000 to 359.000000 by 1.000000 #points=65160
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 359000000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl::lon=10.000000,lat=12.000000,i=36731,ix=11,iy=103,val=-67:lon=20.000000,lat=80.000000,i=61221,ix=21,iy=171,val=-18
+3:130678:grid_template=0:winds(N/S):
+	lat-lon grid:(360 x 181) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 1.000000
+	lon 0.000000 to 359.000000 by 1.000000 #points=65160
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 359000000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl::lon=10.000000,lat=12.000000,i=36731,ix=11,iy=103,val=-28:lon=20.000000,lat=80.000000,i=61221,ix=21,iy=171,val=-14
+4:196017:grid_template=0:winds(N/S):
+	lat-lon grid:(360 x 181) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 1.000000
+	lon 0.000000 to 359.000000 by 1.000000 #points=65160
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 359000000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl::lon=10.000000,lat=12.000000,i=36731,ix=11,iy=103,val=-58:lon=20.000000,lat=80.000000,i=61221,ix=21,iy=171,val=-23
+5:261356:grid_template=0:winds(N/S):
+	lat-lon grid:(360 x 181) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 1.000000
+	lon 0.000000 to 359.000000 by 1.000000 #points=65160
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 359000000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl::lon=10.000000,lat=12.000000,i=36731,ix=11,iy=103,val=-6:lon=20.000000,lat=80.000000,i=61221,ix=21,iy=171,val=1
+6:326695:grid_template=0:winds(N/S):
+	lat-lon grid:(360 x 181) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 1.000000
+	lon 0.000000 to 359.000000 by 1.000000 #points=65160
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 359000000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl::lon=10.000000,lat=12.000000,i=36731,ix=11,iy=103,val=-27:lon=20.000000,lat=80.000000,i=61221,ix=21,iy=171,val=-16
+7:383889:grid_template=0:winds(N/S):
+	lat-lon grid:(360 x 181) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 1.000000
+	lon 0.000000 to 359.000000 by 1.000000 #points=65160
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 359000000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl::lon=10.000000,lat=12.000000,i=36731,ix=11,iy=103,val=-1:lon=20.000000,lat=80.000000,i=61221,ix=21,iy=171,val=10
+8:449228:grid_template=0:winds(N/S):
+	lat-lon grid:(360 x 181) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 1.000000
+	lon 0.000000 to 359.000000 by 1.000000 #points=65160
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 359000000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl::lon=10.000000,lat=12.000000,i=36731,ix=11,iy=103,val=6:lon=20.000000,lat=80.000000,i=61221,ix=21,iy=171,val=-3
+9:506422:grid_template=0:winds(N/S):
+	lat-lon grid:(360 x 181) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 1.000000
+	lon 0.000000 to 359.000000 by 1.000000 #points=65160
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 359000000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl::lon=10.000000,lat=12.000000,i=36731,ix=11,iy=103,val=-10:lon=20.000000,lat=80.000000,i=61221,ix=21,iy=171,val=8
+10:571761:grid_template=0:winds(N/S):
+	lat-lon grid:(360 x 181) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 1.000000
+	lon 0.000000 to 359.000000 by 1.000000 #points=65160
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 359000000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl::lon=10.000000,lat=12.000000,i=36731,ix=11,iy=103,val=11:lon=20.000000,lat=80.000000,i=61221,ix=21,iy=171,val=-1
+11:628955:grid_template=0:winds(N/S):
+	lat-lon grid:(360 x 181) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 1.000000
+	lon 0.000000 to 359.000000 by 1.000000 #points=65160
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 359000000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl::lon=10.000000,lat=12.000000,i=36731,ix=11,iy=103,val=-8:lon=20.000000,lat=80.000000,i=61221,ix=21,iy=171,val=9
+12:694294:grid_template=0:winds(N/S):
+	lat-lon grid:(360 x 181) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 1.000000
+	lon 0.000000 to 359.000000 by 1.000000 #points=65160
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 359000000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl::lon=10.000000,lat=12.000000,i=36731,ix=11,iy=103,val=14:lon=20.000000,lat=80.000000,i=61221,ix=21,iy=171,val=8
+13:751488:grid_template=0:winds(N/S):
+	lat-lon grid:(360 x 181) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 1.000000
+	lon 0.000000 to 359.000000 by 1.000000 #points=65160
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 359000000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl::lon=10.000000,lat=12.000000,i=36731,ix=11,iy=103,val=-3:lon=20.000000,lat=80.000000,i=61221,ix=21,iy=171,val=-8
+14:816827:grid_template=0:winds(N/S):
+	lat-lon grid:(360 x 181) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 1.000000
+	lon 0.000000 to 359.000000 by 1.000000 #points=65160
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 359000000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl::lon=10.000000,lat=12.000000,i=36731,ix=11,iy=103,val=15:lon=20.000000,lat=80.000000,i=61221,ix=21,iy=171,val=13
diff --git a/tests/data/ref_new_grid_ncep_4.txt b/tests/data/ref_new_grid_ncep_4.txt
new file mode 100644
index 00000000..9f07779c
--- /dev/null
+++ b/tests/data/ref_new_grid_ncep_4.txt
@@ -0,0 +1,84 @@
+1:0:grid_template=0:winds(N/S):
+	lat-lon grid:(720 x 361) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 0.500000
+	lon 0.000000 to 359.500000 by 0.500000 #points=259920
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 359500000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl::lon=10.000000,lat=12.000000,i=146901,ix=21,iy=205,val=-43:lon=20.000000,lat=80.000000,i=244841,ix=41,iy=341,val=-17
+2:260099:grid_template=0:winds(N/S):
+	lat-lon grid:(720 x 361) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 0.500000
+	lon 0.000000 to 359.500000 by 0.500000 #points=259920
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 359500000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl::lon=10.000000,lat=12.000000,i=146901,ix=21,iy=205,val=-67.0017:lon=20.000000,lat=80.000000,i=244841,ix=41,iy=341,val=-18.0017
+3:520198:grid_template=0:winds(N/S):
+	lat-lon grid:(720 x 361) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 0.500000
+	lon 0.000000 to 359.500000 by 0.500000 #points=259920
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 359500000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl::lon=10.000000,lat=12.000000,i=146901,ix=21,iy=205,val=-28:lon=20.000000,lat=80.000000,i=244841,ix=41,iy=341,val=-14
+4:780297:grid_template=0:winds(N/S):
+	lat-lon grid:(720 x 361) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 0.500000
+	lon 0.000000 to 359.500000 by 0.500000 #points=259920
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 359500000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl::lon=10.000000,lat=12.000000,i=146901,ix=21,iy=205,val=-58:lon=20.000000,lat=80.000000,i=244841,ix=41,iy=341,val=-23
+5:1040396:grid_template=0:winds(N/S):
+	lat-lon grid:(720 x 361) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 0.500000
+	lon 0.000000 to 359.500000 by 0.500000 #points=259920
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 359500000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl::lon=10.000000,lat=12.000000,i=146901,ix=21,iy=205,val=-6.00665:lon=20.000000,lat=80.000000,i=244841,ix=41,iy=341,val=0.993347
+6:1300495:grid_template=0:winds(N/S):
+	lat-lon grid:(720 x 361) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 0.500000
+	lon 0.000000 to 359.500000 by 0.500000 #points=259920
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 359500000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl::lon=10.000000,lat=12.000000,i=146901,ix=21,iy=205,val=-27.0025:lon=20.000000,lat=80.000000,i=244841,ix=41,iy=341,val=-16.0025
+7:1528104:grid_template=0:winds(N/S):
+	lat-lon grid:(720 x 361) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 0.500000
+	lon 0.000000 to 359.500000 by 0.500000 #points=259920
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 359500000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl::lon=10.000000,lat=12.000000,i=146901,ix=21,iy=205,val=-1:lon=20.000000,lat=80.000000,i=244841,ix=41,iy=341,val=10
+8:1788203:grid_template=0:winds(N/S):
+	lat-lon grid:(720 x 361) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 0.500000
+	lon 0.000000 to 359.500000 by 0.500000 #points=259920
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 359500000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl::lon=10.000000,lat=12.000000,i=146901,ix=21,iy=205,val=5.99762:lon=20.000000,lat=80.000000,i=244841,ix=41,iy=341,val=-3.00238
+9:2015812:grid_template=0:winds(N/S):
+	lat-lon grid:(720 x 361) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 0.500000
+	lon 0.000000 to 359.500000 by 0.500000 #points=259920
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 359500000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl::lon=10.000000,lat=12.000000,i=146901,ix=21,iy=205,val=-10:lon=20.000000,lat=80.000000,i=244841,ix=41,iy=341,val=8
+10:2275911:grid_template=0:winds(N/S):
+	lat-lon grid:(720 x 361) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 0.500000
+	lon 0.000000 to 359.500000 by 0.500000 #points=259920
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 359500000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl::lon=10.000000,lat=12.000000,i=146901,ix=21,iy=205,val=10.9935:lon=20.000000,lat=80.000000,i=244841,ix=41,iy=341,val=-1.00645
+11:2503520:grid_template=0:winds(N/S):
+	lat-lon grid:(720 x 361) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 0.500000
+	lon 0.000000 to 359.500000 by 0.500000 #points=259920
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 359500000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl::lon=10.000000,lat=12.000000,i=146901,ix=21,iy=205,val=-8.00293:lon=20.000000,lat=80.000000,i=244841,ix=41,iy=341,val=8.99707
+12:2763619:grid_template=0:winds(N/S):
+	lat-lon grid:(720 x 361) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 0.500000
+	lon 0.000000 to 359.500000 by 0.500000 #points=259920
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 359500000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl::lon=10.000000,lat=12.000000,i=146901,ix=21,iy=205,val=14:lon=20.000000,lat=80.000000,i=244841,ix=41,iy=341,val=8
+13:2991228:grid_template=0:winds(N/S):
+	lat-lon grid:(720 x 361) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 0.500000
+	lon 0.000000 to 359.500000 by 0.500000 #points=259920
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 359500000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl::lon=10.000000,lat=12.000000,i=146901,ix=21,iy=205,val=-3:lon=20.000000,lat=80.000000,i=244841,ix=41,iy=341,val=-8
+14:3251327:grid_template=0:winds(N/S):
+	lat-lon grid:(720 x 361) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 0.500000
+	lon 0.000000 to 359.500000 by 0.500000 #points=259920
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 359500000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl::lon=10.000000,lat=12.000000,i=146901,ix=21,iy=205,val=15:lon=20.000000,lat=80.000000,i=244841,ix=41,iy=341,val=13
diff --git a/tests/data/ref_new_grid_ncep_45.txt b/tests/data/ref_new_grid_ncep_45.txt
new file mode 100644
index 00000000..a3e7e61d
--- /dev/null
+++ b/tests/data/ref_new_grid_ncep_45.txt
@@ -0,0 +1,84 @@
+1:0:grid_template=0:winds(N/S):
+	lat-lon grid:(288 x 145) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 1.250000
+	lon 0.000000 to 358.750000 by 1.250000 #points=41760
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 358750000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl::lon=10.000000,lat=12.500000,i=23625,ix=9,iy=83,val=-46.2517:lon=20.000000,lat=80.000000,i=39185,ix=17,iy=137,val=-17.2517
+2:41939:grid_template=0:winds(N/S):
+	lat-lon grid:(288 x 145) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 1.250000
+	lon 0.000000 to 358.750000 by 1.250000 #points=41760
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 358750000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl::lon=10.000000,lat=12.500000,i=23625,ix=9,iy=83,val=-66.3756:lon=20.000000,lat=80.000000,i=39185,ix=17,iy=137,val=-18.3756
+3:83878:grid_template=0:winds(N/S):
+	lat-lon grid:(288 x 145) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 1.250000
+	lon 0.000000 to 358.750000 by 1.250000 #points=41760
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 358750000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl::lon=10.000000,lat=12.500000,i=23625,ix=9,iy=83,val=-29:lon=20.000000,lat=80.000000,i=39185,ix=17,iy=137,val=-14
+4:125817:grid_template=0:winds(N/S):
+	lat-lon grid:(288 x 145) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 1.250000
+	lon 0.000000 to 358.750000 by 1.250000 #points=41760
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 358750000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl::lon=10.000000,lat=12.500000,i=23625,ix=9,iy=83,val=-56:lon=20.000000,lat=80.000000,i=39185,ix=17,iy=137,val=-23
+5:162536:grid_template=0:winds(N/S):
+	lat-lon grid:(288 x 145) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 1.250000
+	lon 0.000000 to 358.750000 by 1.250000 #points=41760
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 358750000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl::lon=10.000000,lat=12.500000,i=23625,ix=9,iy=83,val=-6.00665:lon=20.000000,lat=80.000000,i=39185,ix=17,iy=137,val=0.993347
+6:204475:grid_template=0:winds(N/S):
+	lat-lon grid:(288 x 145) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 1.250000
+	lon 0.000000 to 358.750000 by 1.250000 #points=41760
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 358750000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl::lon=10.000000,lat=12.500000,i=23625,ix=9,iy=83,val=-27.0025:lon=20.000000,lat=80.000000,i=39185,ix=17,iy=137,val=-16.0025
+7:241194:grid_template=0:winds(N/S):
+	lat-lon grid:(288 x 145) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 1.250000
+	lon 0.000000 to 358.750000 by 1.250000 #points=41760
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 358750000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl::lon=10.000000,lat=12.500000,i=23625,ix=9,iy=83,val=0.439747:lon=20.000000,lat=80.000000,i=39185,ix=17,iy=137,val=10.4397
+8:283133:grid_template=0:winds(N/S):
+	lat-lon grid:(288 x 145) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 1.250000
+	lon 0.000000 to 358.750000 by 1.250000 #points=41760
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 358750000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl::lon=10.000000,lat=12.500000,i=23625,ix=9,iy=83,val=3.4976:lon=20.000000,lat=80.000000,i=39185,ix=17,iy=137,val=-2.5024
+9:319852:grid_template=0:winds(N/S):
+	lat-lon grid:(288 x 145) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 1.250000
+	lon 0.000000 to 358.750000 by 1.250000 #points=41760
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 358750000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl::lon=10.000000,lat=12.500000,i=23625,ix=9,iy=83,val=-7.87319:lon=20.000000,lat=80.000000,i=39185,ix=17,iy=137,val=8.12681
+10:361791:grid_template=0:winds(N/S):
+	lat-lon grid:(288 x 145) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 1.250000
+	lon 0.000000 to 358.750000 by 1.250000 #points=41760
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 358750000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl::lon=10.000000,lat=12.500000,i=23625,ix=9,iy=83,val=10.2451:lon=20.000000,lat=80.000000,i=39185,ix=17,iy=137,val=-0.754887
+11:398510:grid_template=0:winds(N/S):
+	lat-lon grid:(288 x 145) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 1.250000
+	lon 0.000000 to 358.750000 by 1.250000 #points=41760
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 358750000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl::lon=10.000000,lat=12.500000,i=23625,ix=9,iy=83,val=-8:lon=20.000000,lat=80.000000,i=39185,ix=17,iy=137,val=9
+12:440449:grid_template=0:winds(N/S):
+	lat-lon grid:(288 x 145) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 1.250000
+	lon 0.000000 to 358.750000 by 1.250000 #points=41760
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 358750000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl::lon=10.000000,lat=12.500000,i=23625,ix=9,iy=83,val=15:lon=20.000000,lat=80.000000,i=39185,ix=17,iy=137,val=8
+13:477168:grid_template=0:winds(N/S):
+	lat-lon grid:(288 x 145) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 1.250000
+	lon 0.000000 to 358.750000 by 1.250000 #points=41760
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 358750000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl::lon=10.000000,lat=12.500000,i=23625,ix=9,iy=83,val=-3.99862:lon=20.000000,lat=80.000000,i=39185,ix=17,iy=137,val=-7.99862
+14:519107:grid_template=0:winds(N/S):
+	lat-lon grid:(288 x 145) units 1e-06 input WE:NS output WE:SN res 48
+	lat 90.000000 to -90.000000 by 1.250000
+	lon 0.000000 to 358.750000 by 1.250000 #points=41760
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=90000000 to -90000000 lon=0 to 358750000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl::lon=10.000000,lat=12.500000,i=23625,ix=9,iy=83,val=14.6878:lon=20.000000,lat=80.000000,i=39185,ix=17,iy=137,val=12.6878
diff --git a/tests/data/ref_new_grid_ncep_98.txt b/tests/data/ref_new_grid_ncep_98.txt
new file mode 100644
index 00000000..0fc35bab
--- /dev/null
+++ b/tests/data/ref_new_grid_ncep_98.txt
@@ -0,0 +1,98 @@
+1:0:grid_template=40:winds(N/S):
+	Gaussian grid: (192 x 94) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=47 #points=18048
+	lat 88.541950 to -88.541950
+	lon 0.000000 to 358.125000 by 1.875000
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=88541950 to -88541950 lon=0 to 358125000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl::lon=9.375000,lat=12.380776,i=10182,ix=6,iy=54,val=-43.7909:lon=20.625000,lat=80.947359,i=17100,ix=12,iy=90,val=-13.7909
+2:18227:grid_template=40:winds(N/S):
+	Gaussian grid: (192 x 94) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=47 #points=18048
+	lat 88.541950 to -88.541950
+	lon 0.000000 to 358.125000 by 1.875000
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=88541950 to -88541950 lon=0 to 358125000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl::lon=9.375000,lat=12.380776,i=10182,ix=6,iy=54,val=-63.2865:lon=20.625000,lat=80.947359,i=17100,ix=12,iy=90,val=-24.2865
+3:36454:grid_template=40:winds(N/S):
+	Gaussian grid: (192 x 94) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=47 #points=18048
+	lat 88.541950 to -88.541950
+	lon 0.000000 to 358.125000 by 1.875000
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=88541950 to -88541950 lon=0 to 358125000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl::lon=9.375000,lat=12.380776,i=10182,ix=6,iy=54,val=-29:lon=20.625000,lat=80.947359,i=17100,ix=12,iy=90,val=-11
+4:54681:grid_template=40:winds(N/S):
+	Gaussian grid: (192 x 94) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=47 #points=18048
+	lat 88.541950 to -88.541950
+	lon 0.000000 to 358.125000 by 1.875000
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=88541950 to -88541950 lon=0 to 358125000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl::lon=9.375000,lat=12.380776,i=10182,ix=6,iy=54,val=-55.2366:lon=20.625000,lat=80.947359,i=17100,ix=12,iy=90,val=-26.2366
+5:70652:grid_template=40:winds(N/S):
+	Gaussian grid: (192 x 94) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=47 #points=18048
+	lat 88.541950 to -88.541950
+	lon 0.000000 to 358.125000 by 1.875000
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=88541950 to -88541950 lon=0 to 358125000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl::lon=9.375000,lat=12.380776,i=10182,ix=6,iy=54,val=-6.76678:lon=20.625000,lat=80.947359,i=17100,ix=12,iy=90,val=-0.766781
+6:88879:grid_template=40:winds(N/S):
+	Gaussian grid: (192 x 94) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=47 #points=18048
+	lat 88.541950 to -88.541950
+	lon 0.000000 to 358.125000 by 1.875000
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=88541950 to -88541950 lon=0 to 358125000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl::lon=9.375000,lat=12.380776,i=10182,ix=6,iy=54,val=-26.4428:lon=20.625000,lat=80.947359,i=17100,ix=12,iy=90,val=-16.4428
+7:104850:grid_template=40:winds(N/S):
+	Gaussian grid: (192 x 94) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=47 #points=18048
+	lat 88.541950 to -88.541950
+	lon 0.000000 to 358.125000 by 1.875000
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=88541950 to -88541950 lon=0 to 358125000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl::lon=9.375000,lat=12.380776,i=10182,ix=6,iy=54,val=0.988941:lon=20.625000,lat=80.947359,i=17100,ix=12,iy=90,val=9.98894
+8:123077:grid_template=40:winds(N/S):
+	Gaussian grid: (192 x 94) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=47 #points=18048
+	lat 88.541950 to -88.541950
+	lon 0.000000 to 358.125000 by 1.875000
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=88541950 to -88541950 lon=0 to 358125000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl::lon=9.375000,lat=12.380776,i=10182,ix=6,iy=54,val=3.84836:lon=20.625000,lat=80.947359,i=17100,ix=12,iy=90,val=-2.15164
+9:139048:grid_template=40:winds(N/S):
+	Gaussian grid: (192 x 94) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=47 #points=18048
+	lat 88.541950 to -88.541950
+	lon 0.000000 to 358.125000 by 1.875000
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=88541950 to -88541950 lon=0 to 358125000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl::lon=9.375000,lat=12.380776,i=10182,ix=6,iy=54,val=-6.45764:lon=20.625000,lat=80.947359,i=17100,ix=12,iy=90,val=7.54236
+10:157275:grid_template=40:winds(N/S):
+	Gaussian grid: (192 x 94) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=47 #points=18048
+	lat 88.541950 to -88.541950
+	lon 0.000000 to 358.125000 by 1.875000
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=88541950 to -88541950 lon=0 to 358125000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl::lon=9.375000,lat=12.380776,i=10182,ix=6,iy=54,val=9.4843:lon=20.625000,lat=80.947359,i=17100,ix=12,iy=90,val=-0.515701
+11:173246:grid_template=40:winds(N/S):
+	Gaussian grid: (192 x 94) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=47 #points=18048
+	lat 88.541950 to -88.541950
+	lon 0.000000 to 358.125000 by 1.875000
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=88541950 to -88541950 lon=0 to 358125000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl::lon=9.375000,lat=12.380776,i=10182,ix=6,iy=54,val=-7.57658:lon=20.625000,lat=80.947359,i=17100,ix=12,iy=90,val=8.42342
+12:191473:grid_template=40:winds(N/S):
+	Gaussian grid: (192 x 94) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=47 #points=18048
+	lat 88.541950 to -88.541950
+	lon 0.000000 to 358.125000 by 1.875000
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=88541950 to -88541950 lon=0 to 358125000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl::lon=9.375000,lat=12.380776,i=10182,ix=6,iy=54,val=14.6009:lon=20.625000,lat=80.947359,i=17100,ix=12,iy=90,val=7.60093
+13:207444:grid_template=40:winds(N/S):
+	Gaussian grid: (192 x 94) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=47 #points=18048
+	lat 88.541950 to -88.541950
+	lon 0.000000 to 358.125000 by 1.875000
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=88541950 to -88541950 lon=0 to 358125000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl::lon=9.375000,lat=12.380776,i=10182,ix=6,iy=54,val=-3.99902:lon=20.625000,lat=80.947359,i=17100,ix=12,iy=90,val=-5.99902
+14:225671:grid_template=40:winds(N/S):
+	Gaussian grid: (192 x 94) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=47 #points=18048
+	lat 88.541950 to -88.541950
+	lon 0.000000 to 358.125000 by 1.875000
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=88541950 to -88541950 lon=0 to 358125000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl::lon=9.375000,lat=12.380776,i=10182,ix=6,iy=54,val=14.3316:lon=20.625000,lat=80.947359,i=17100,ix=12,iy=90,val=12.3316
diff --git a/tests/data/ref_new_grid_ncep_t1148.txt b/tests/data/ref_new_grid_ncep_t1148.txt
new file mode 100644
index 00000000..9785ad96
--- /dev/null
+++ b/tests/data/ref_new_grid_ncep_t1148.txt
@@ -0,0 +1,98 @@
+1:0:grid_template=40:winds(N/S):
+	Gaussian grid: (2304 x 1152) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=576 #points=2654208
+	lat 89.880446 to -89.880446
+	lon 0.000000 to 359.843750 by 0.156250
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89880446 to -89880446 lon=0 to 359843750:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl:
+2:2654387:grid_template=40:winds(N/S):
+	Gaussian grid: (2304 x 1152) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=576 #points=2654208
+	lat 89.880446 to -89.880446
+	lon 0.000000 to 359.843750 by 0.156250
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89880446 to -89880446 lon=0 to 359843750:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl:
+3:5308774:grid_template=40:winds(N/S):
+	Gaussian grid: (2304 x 1152) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=576 #points=2654208
+	lat 89.880446 to -89.880446
+	lon 0.000000 to 359.843750 by 0.156250
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89880446 to -89880446 lon=0 to 359843750:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl:
+4:7963161:grid_template=40:winds(N/S):
+	Gaussian grid: (2304 x 1152) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=576 #points=2654208
+	lat 89.880446 to -89.880446
+	lon 0.000000 to 359.843750 by 0.156250
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89880446 to -89880446 lon=0 to 359843750:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl:
+5:10617548:grid_template=40:winds(N/S):
+	Gaussian grid: (2304 x 1152) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=576 #points=2654208
+	lat 89.880446 to -89.880446
+	lon 0.000000 to 359.843750 by 0.156250
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89880446 to -89880446 lon=0 to 359843750:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl:
+6:13271935:grid_template=40:winds(N/S):
+	Gaussian grid: (2304 x 1152) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=576 #points=2654208
+	lat 89.880446 to -89.880446
+	lon 0.000000 to 359.843750 by 0.156250
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89880446 to -89880446 lon=0 to 359843750:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl:
+7:15594546:grid_template=40:winds(N/S):
+	Gaussian grid: (2304 x 1152) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=576 #points=2654208
+	lat 89.880446 to -89.880446
+	lon 0.000000 to 359.843750 by 0.156250
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89880446 to -89880446 lon=0 to 359843750:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl:
+8:18248933:grid_template=40:winds(N/S):
+	Gaussian grid: (2304 x 1152) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=576 #points=2654208
+	lat 89.880446 to -89.880446
+	lon 0.000000 to 359.843750 by 0.156250
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89880446 to -89880446 lon=0 to 359843750:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl:
+9:20571544:grid_template=40:winds(N/S):
+	Gaussian grid: (2304 x 1152) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=576 #points=2654208
+	lat 89.880446 to -89.880446
+	lon 0.000000 to 359.843750 by 0.156250
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89880446 to -89880446 lon=0 to 359843750:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl:
+10:23225931:grid_template=40:winds(N/S):
+	Gaussian grid: (2304 x 1152) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=576 #points=2654208
+	lat 89.880446 to -89.880446
+	lon 0.000000 to 359.843750 by 0.156250
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89880446 to -89880446 lon=0 to 359843750:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl:
+11:25548542:grid_template=40:winds(N/S):
+	Gaussian grid: (2304 x 1152) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=576 #points=2654208
+	lat 89.880446 to -89.880446
+	lon 0.000000 to 359.843750 by 0.156250
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89880446 to -89880446 lon=0 to 359843750:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl:
+12:28202929:grid_template=40:winds(N/S):
+	Gaussian grid: (2304 x 1152) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=576 #points=2654208
+	lat 89.880446 to -89.880446
+	lon 0.000000 to 359.843750 by 0.156250
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89880446 to -89880446 lon=0 to 359843750:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl:
+13:30525540:grid_template=40:winds(N/S):
+	Gaussian grid: (2304 x 1152) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=576 #points=2654208
+	lat 89.880446 to -89.880446
+	lon 0.000000 to 359.843750 by 0.156250
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89880446 to -89880446 lon=0 to 359843750:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl:
+14:33179927:grid_template=40:winds(N/S):
+	Gaussian grid: (2304 x 1152) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=576 #points=2654208
+	lat 89.880446 to -89.880446
+	lon 0.000000 to 359.843750 by 0.156250
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89880446 to -89880446 lon=0 to 359843750:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl:
diff --git a/tests/data/ref_new_grid_ncep_t126.txt b/tests/data/ref_new_grid_ncep_t126.txt
new file mode 100644
index 00000000..2ad356da
--- /dev/null
+++ b/tests/data/ref_new_grid_ncep_t126.txt
@@ -0,0 +1,98 @@
+1:0:grid_template=40:winds(N/S):
+	Gaussian grid: (384 x 190) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=95 #points=72960
+	lat 89.276713 to -89.276713
+	lon 0.000000 to 359.062500 by 0.937500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89276713 to -89276713 lon=0 to 359062500:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl:
+2:73139:grid_template=40:winds(N/S):
+	Gaussian grid: (384 x 190) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=95 #points=72960
+	lat 89.276713 to -89.276713
+	lon 0.000000 to 359.062500 by 0.937500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89276713 to -89276713 lon=0 to 359062500:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl:
+3:146278:grid_template=40:winds(N/S):
+	Gaussian grid: (384 x 190) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=95 #points=72960
+	lat 89.276713 to -89.276713
+	lon 0.000000 to 359.062500 by 0.937500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89276713 to -89276713 lon=0 to 359062500:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl:
+4:219417:grid_template=40:winds(N/S):
+	Gaussian grid: (384 x 190) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=95 #points=72960
+	lat 89.276713 to -89.276713
+	lon 0.000000 to 359.062500 by 0.937500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89276713 to -89276713 lon=0 to 359062500:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl:
+5:283436:grid_template=40:winds(N/S):
+	Gaussian grid: (384 x 190) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=95 #points=72960
+	lat 89.276713 to -89.276713
+	lon 0.000000 to 359.062500 by 0.937500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89276713 to -89276713 lon=0 to 359062500:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl:
+6:356575:grid_template=40:winds(N/S):
+	Gaussian grid: (384 x 190) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=95 #points=72960
+	lat 89.276713 to -89.276713
+	lon 0.000000 to 359.062500 by 0.937500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89276713 to -89276713 lon=0 to 359062500:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl:
+7:420594:grid_template=40:winds(N/S):
+	Gaussian grid: (384 x 190) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=95 #points=72960
+	lat 89.276713 to -89.276713
+	lon 0.000000 to 359.062500 by 0.937500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89276713 to -89276713 lon=0 to 359062500:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl:
+8:493733:grid_template=40:winds(N/S):
+	Gaussian grid: (384 x 190) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=95 #points=72960
+	lat 89.276713 to -89.276713
+	lon 0.000000 to 359.062500 by 0.937500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89276713 to -89276713 lon=0 to 359062500:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl:
+9:557752:grid_template=40:winds(N/S):
+	Gaussian grid: (384 x 190) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=95 #points=72960
+	lat 89.276713 to -89.276713
+	lon 0.000000 to 359.062500 by 0.937500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89276713 to -89276713 lon=0 to 359062500:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl:
+10:630891:grid_template=40:winds(N/S):
+	Gaussian grid: (384 x 190) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=95 #points=72960
+	lat 89.276713 to -89.276713
+	lon 0.000000 to 359.062500 by 0.937500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89276713 to -89276713 lon=0 to 359062500:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl:
+11:694910:grid_template=40:winds(N/S):
+	Gaussian grid: (384 x 190) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=95 #points=72960
+	lat 89.276713 to -89.276713
+	lon 0.000000 to 359.062500 by 0.937500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89276713 to -89276713 lon=0 to 359062500:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl:
+12:768049:grid_template=40:winds(N/S):
+	Gaussian grid: (384 x 190) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=95 #points=72960
+	lat 89.276713 to -89.276713
+	lon 0.000000 to 359.062500 by 0.937500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89276713 to -89276713 lon=0 to 359062500:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl:
+13:832068:grid_template=40:winds(N/S):
+	Gaussian grid: (384 x 190) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=95 #points=72960
+	lat 89.276713 to -89.276713
+	lon 0.000000 to 359.062500 by 0.937500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89276713 to -89276713 lon=0 to 359062500:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl:
+14:905207:grid_template=40:winds(N/S):
+	Gaussian grid: (384 x 190) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=95 #points=72960
+	lat 89.276713 to -89.276713
+	lon 0.000000 to 359.062500 by 0.937500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89276713 to -89276713 lon=0 to 359062500:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl:
diff --git a/tests/data/ref_new_grid_ncep_t1534.txt b/tests/data/ref_new_grid_ncep_t1534.txt
new file mode 100644
index 00000000..f0f9e29c
--- /dev/null
+++ b/tests/data/ref_new_grid_ncep_t1534.txt
@@ -0,0 +1,98 @@
+1:0:grid_template=40:winds(N/S):
+	Gaussian grid: (3072 x 1536) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=768 #points=4718592
+	lat 89.910325 to -89.910325
+	lon 0.000000 to 359.882813 by 0.117188
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89910325 to -89910325 lon=0 to 359882813:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl:
+2:4718771:grid_template=40:winds(N/S):
+	Gaussian grid: (3072 x 1536) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=768 #points=4718592
+	lat 89.910325 to -89.910325
+	lon 0.000000 to 359.882813 by 0.117188
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89910325 to -89910325 lon=0 to 359882813:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl:
+3:9437542:grid_template=40:winds(N/S):
+	Gaussian grid: (3072 x 1536) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=768 #points=4718592
+	lat 89.910325 to -89.910325
+	lon 0.000000 to 359.882813 by 0.117188
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89910325 to -89910325 lon=0 to 359882813:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl:
+4:14156313:grid_template=40:winds(N/S):
+	Gaussian grid: (3072 x 1536) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=768 #points=4718592
+	lat 89.910325 to -89.910325
+	lon 0.000000 to 359.882813 by 0.117188
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89910325 to -89910325 lon=0 to 359882813:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl:
+5:18875084:grid_template=40:winds(N/S):
+	Gaussian grid: (3072 x 1536) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=768 #points=4718592
+	lat 89.910325 to -89.910325
+	lon 0.000000 to 359.882813 by 0.117188
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89910325 to -89910325 lon=0 to 359882813:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl:
+6:23593855:grid_template=40:winds(N/S):
+	Gaussian grid: (3072 x 1536) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=768 #points=4718592
+	lat 89.910325 to -89.910325
+	lon 0.000000 to 359.882813 by 0.117188
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89910325 to -89910325 lon=0 to 359882813:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl:
+7:27722802:grid_template=40:winds(N/S):
+	Gaussian grid: (3072 x 1536) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=768 #points=4718592
+	lat 89.910325 to -89.910325
+	lon 0.000000 to 359.882813 by 0.117188
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89910325 to -89910325 lon=0 to 359882813:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl:
+8:32441573:grid_template=40:winds(N/S):
+	Gaussian grid: (3072 x 1536) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=768 #points=4718592
+	lat 89.910325 to -89.910325
+	lon 0.000000 to 359.882813 by 0.117188
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89910325 to -89910325 lon=0 to 359882813:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl:
+9:36570520:grid_template=40:winds(N/S):
+	Gaussian grid: (3072 x 1536) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=768 #points=4718592
+	lat 89.910325 to -89.910325
+	lon 0.000000 to 359.882813 by 0.117188
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89910325 to -89910325 lon=0 to 359882813:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl:
+10:41289291:grid_template=40:winds(N/S):
+	Gaussian grid: (3072 x 1536) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=768 #points=4718592
+	lat 89.910325 to -89.910325
+	lon 0.000000 to 359.882813 by 0.117188
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89910325 to -89910325 lon=0 to 359882813:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl:
+11:45418238:grid_template=40:winds(N/S):
+	Gaussian grid: (3072 x 1536) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=768 #points=4718592
+	lat 89.910325 to -89.910325
+	lon 0.000000 to 359.882813 by 0.117188
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89910325 to -89910325 lon=0 to 359882813:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl:
+12:50137009:grid_template=40:winds(N/S):
+	Gaussian grid: (3072 x 1536) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=768 #points=4718592
+	lat 89.910325 to -89.910325
+	lon 0.000000 to 359.882813 by 0.117188
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89910325 to -89910325 lon=0 to 359882813:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl:
+13:54265956:grid_template=40:winds(N/S):
+	Gaussian grid: (3072 x 1536) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=768 #points=4718592
+	lat 89.910325 to -89.910325
+	lon 0.000000 to 359.882813 by 0.117188
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89910325 to -89910325 lon=0 to 359882813:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl:
+14:58984727:grid_template=40:winds(N/S):
+	Gaussian grid: (3072 x 1536) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=768 #points=4718592
+	lat 89.910325 to -89.910325
+	lon 0.000000 to 359.882813 by 0.117188
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89910325 to -89910325 lon=0 to 359882813:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl:
diff --git a/tests/data/ref_new_grid_ncep_t170.txt b/tests/data/ref_new_grid_ncep_t170.txt
new file mode 100644
index 00000000..4bb22222
--- /dev/null
+++ b/tests/data/ref_new_grid_ncep_t170.txt
@@ -0,0 +1,98 @@
+1:0:grid_template=40:winds(N/S):
+	Gaussian grid: (512 x 256) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=128 #points=131072
+	lat 89.462822 to -89.462822
+	lon 0.000000 to 359.296875 by 0.703125
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89462822 to -89462822 lon=0 to 359296875:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl:
+2:131251:grid_template=40:winds(N/S):
+	Gaussian grid: (512 x 256) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=128 #points=131072
+	lat 89.462822 to -89.462822
+	lon 0.000000 to 359.296875 by 0.703125
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89462822 to -89462822 lon=0 to 359296875:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl:
+3:262502:grid_template=40:winds(N/S):
+	Gaussian grid: (512 x 256) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=128 #points=131072
+	lat 89.462822 to -89.462822
+	lon 0.000000 to 359.296875 by 0.703125
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89462822 to -89462822 lon=0 to 359296875:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl:
+4:393753:grid_template=40:winds(N/S):
+	Gaussian grid: (512 x 256) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=128 #points=131072
+	lat 89.462822 to -89.462822
+	lon 0.000000 to 359.296875 by 0.703125
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89462822 to -89462822 lon=0 to 359296875:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl:
+5:508620:grid_template=40:winds(N/S):
+	Gaussian grid: (512 x 256) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=128 #points=131072
+	lat 89.462822 to -89.462822
+	lon 0.000000 to 359.296875 by 0.703125
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89462822 to -89462822 lon=0 to 359296875:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl:
+6:639871:grid_template=40:winds(N/S):
+	Gaussian grid: (512 x 256) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=128 #points=131072
+	lat 89.462822 to -89.462822
+	lon 0.000000 to 359.296875 by 0.703125
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89462822 to -89462822 lon=0 to 359296875:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl:
+7:754738:grid_template=40:winds(N/S):
+	Gaussian grid: (512 x 256) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=128 #points=131072
+	lat 89.462822 to -89.462822
+	lon 0.000000 to 359.296875 by 0.703125
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89462822 to -89462822 lon=0 to 359296875:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl:
+8:885989:grid_template=40:winds(N/S):
+	Gaussian grid: (512 x 256) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=128 #points=131072
+	lat 89.462822 to -89.462822
+	lon 0.000000 to 359.296875 by 0.703125
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89462822 to -89462822 lon=0 to 359296875:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl:
+9:1000856:grid_template=40:winds(N/S):
+	Gaussian grid: (512 x 256) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=128 #points=131072
+	lat 89.462822 to -89.462822
+	lon 0.000000 to 359.296875 by 0.703125
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89462822 to -89462822 lon=0 to 359296875:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl:
+10:1132107:grid_template=40:winds(N/S):
+	Gaussian grid: (512 x 256) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=128 #points=131072
+	lat 89.462822 to -89.462822
+	lon 0.000000 to 359.296875 by 0.703125
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89462822 to -89462822 lon=0 to 359296875:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl:
+11:1246974:grid_template=40:winds(N/S):
+	Gaussian grid: (512 x 256) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=128 #points=131072
+	lat 89.462822 to -89.462822
+	lon 0.000000 to 359.296875 by 0.703125
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89462822 to -89462822 lon=0 to 359296875:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl:
+12:1378225:grid_template=40:winds(N/S):
+	Gaussian grid: (512 x 256) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=128 #points=131072
+	lat 89.462822 to -89.462822
+	lon 0.000000 to 359.296875 by 0.703125
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89462822 to -89462822 lon=0 to 359296875:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl:
+13:1493092:grid_template=40:winds(N/S):
+	Gaussian grid: (512 x 256) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=128 #points=131072
+	lat 89.462822 to -89.462822
+	lon 0.000000 to 359.296875 by 0.703125
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89462822 to -89462822 lon=0 to 359296875:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl:
+14:1624343:grid_template=40:winds(N/S):
+	Gaussian grid: (512 x 256) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=128 #points=131072
+	lat 89.462822 to -89.462822
+	lon 0.000000 to 359.296875 by 0.703125
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89462822 to -89462822 lon=0 to 359296875:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl:
diff --git a/tests/data/ref_new_grid_ncep_t190.txt b/tests/data/ref_new_grid_ncep_t190.txt
new file mode 100644
index 00000000..b68cd9b4
--- /dev/null
+++ b/tests/data/ref_new_grid_ncep_t190.txt
@@ -0,0 +1,98 @@
+1:0:grid_template=40:winds(N/S):
+	Gaussian grid: (576 x 288) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=144 #points=165888
+	lat 89.522405 to -89.522405
+	lon 0.000000 to 359.375000 by 0.625000
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89522405 to -89522405 lon=0 to 359375000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl:
+2:166067:grid_template=40:winds(N/S):
+	Gaussian grid: (576 x 288) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=144 #points=165888
+	lat 89.522405 to -89.522405
+	lon 0.000000 to 359.375000 by 0.625000
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89522405 to -89522405 lon=0 to 359375000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl:
+3:332134:grid_template=40:winds(N/S):
+	Gaussian grid: (576 x 288) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=144 #points=165888
+	lat 89.522405 to -89.522405
+	lon 0.000000 to 359.375000 by 0.625000
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89522405 to -89522405 lon=0 to 359375000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl:
+4:498201:grid_template=40:winds(N/S):
+	Gaussian grid: (576 x 288) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=144 #points=165888
+	lat 89.522405 to -89.522405
+	lon 0.000000 to 359.375000 by 0.625000
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89522405 to -89522405 lon=0 to 359375000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl:
+5:643532:grid_template=40:winds(N/S):
+	Gaussian grid: (576 x 288) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=144 #points=165888
+	lat 89.522405 to -89.522405
+	lon 0.000000 to 359.375000 by 0.625000
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89522405 to -89522405 lon=0 to 359375000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl:
+6:809599:grid_template=40:winds(N/S):
+	Gaussian grid: (576 x 288) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=144 #points=165888
+	lat 89.522405 to -89.522405
+	lon 0.000000 to 359.375000 by 0.625000
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89522405 to -89522405 lon=0 to 359375000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl:
+7:954930:grid_template=40:winds(N/S):
+	Gaussian grid: (576 x 288) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=144 #points=165888
+	lat 89.522405 to -89.522405
+	lon 0.000000 to 359.375000 by 0.625000
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89522405 to -89522405 lon=0 to 359375000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl:
+8:1120997:grid_template=40:winds(N/S):
+	Gaussian grid: (576 x 288) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=144 #points=165888
+	lat 89.522405 to -89.522405
+	lon 0.000000 to 359.375000 by 0.625000
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89522405 to -89522405 lon=0 to 359375000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl:
+9:1266328:grid_template=40:winds(N/S):
+	Gaussian grid: (576 x 288) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=144 #points=165888
+	lat 89.522405 to -89.522405
+	lon 0.000000 to 359.375000 by 0.625000
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89522405 to -89522405 lon=0 to 359375000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl:
+10:1432395:grid_template=40:winds(N/S):
+	Gaussian grid: (576 x 288) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=144 #points=165888
+	lat 89.522405 to -89.522405
+	lon 0.000000 to 359.375000 by 0.625000
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89522405 to -89522405 lon=0 to 359375000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl:
+11:1577726:grid_template=40:winds(N/S):
+	Gaussian grid: (576 x 288) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=144 #points=165888
+	lat 89.522405 to -89.522405
+	lon 0.000000 to 359.375000 by 0.625000
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89522405 to -89522405 lon=0 to 359375000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl:
+12:1743793:grid_template=40:winds(N/S):
+	Gaussian grid: (576 x 288) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=144 #points=165888
+	lat 89.522405 to -89.522405
+	lon 0.000000 to 359.375000 by 0.625000
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89522405 to -89522405 lon=0 to 359375000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl:
+13:1889124:grid_template=40:winds(N/S):
+	Gaussian grid: (576 x 288) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=144 #points=165888
+	lat 89.522405 to -89.522405
+	lon 0.000000 to 359.375000 by 0.625000
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89522405 to -89522405 lon=0 to 359375000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl:
+14:2055191:grid_template=40:winds(N/S):
+	Gaussian grid: (576 x 288) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=144 #points=165888
+	lat 89.522405 to -89.522405
+	lon 0.000000 to 359.375000 by 0.625000
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89522405 to -89522405 lon=0 to 359375000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl:
diff --git a/tests/data/ref_new_grid_ncep_t254.txt b/tests/data/ref_new_grid_ncep_t254.txt
new file mode 100644
index 00000000..7b49fc11
--- /dev/null
+++ b/tests/data/ref_new_grid_ncep_t254.txt
@@ -0,0 +1,98 @@
+1:0:grid_template=40:winds(N/S):
+	Gaussian grid: (768 x 384) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=192 #points=294912
+	lat 89.641648 to -89.641648
+	lon 0.000000 to 359.531250 by 0.468750
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89641648 to -89641648 lon=0 to 359531250:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl:
+2:295091:grid_template=40:winds(N/S):
+	Gaussian grid: (768 x 384) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=192 #points=294912
+	lat 89.641648 to -89.641648
+	lon 0.000000 to 359.531250 by 0.468750
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89641648 to -89641648 lon=0 to 359531250:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl:
+3:590182:grid_template=40:winds(N/S):
+	Gaussian grid: (768 x 384) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=192 #points=294912
+	lat 89.641648 to -89.641648
+	lon 0.000000 to 359.531250 by 0.468750
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89641648 to -89641648 lon=0 to 359531250:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl:
+4:885273:grid_template=40:winds(N/S):
+	Gaussian grid: (768 x 384) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=192 #points=294912
+	lat 89.641648 to -89.641648
+	lon 0.000000 to 359.531250 by 0.468750
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89641648 to -89641648 lon=0 to 359531250:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl:
+5:1180364:grid_template=40:winds(N/S):
+	Gaussian grid: (768 x 384) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=192 #points=294912
+	lat 89.641648 to -89.641648
+	lon 0.000000 to 359.531250 by 0.468750
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89641648 to -89641648 lon=0 to 359531250:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl:
+6:1475455:grid_template=40:winds(N/S):
+	Gaussian grid: (768 x 384) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=192 #points=294912
+	lat 89.641648 to -89.641648
+	lon 0.000000 to 359.531250 by 0.468750
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89641648 to -89641648 lon=0 to 359531250:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl:
+7:1733682:grid_template=40:winds(N/S):
+	Gaussian grid: (768 x 384) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=192 #points=294912
+	lat 89.641648 to -89.641648
+	lon 0.000000 to 359.531250 by 0.468750
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89641648 to -89641648 lon=0 to 359531250:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl:
+8:2028773:grid_template=40:winds(N/S):
+	Gaussian grid: (768 x 384) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=192 #points=294912
+	lat 89.641648 to -89.641648
+	lon 0.000000 to 359.531250 by 0.468750
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89641648 to -89641648 lon=0 to 359531250:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl:
+9:2287000:grid_template=40:winds(N/S):
+	Gaussian grid: (768 x 384) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=192 #points=294912
+	lat 89.641648 to -89.641648
+	lon 0.000000 to 359.531250 by 0.468750
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89641648 to -89641648 lon=0 to 359531250:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl:
+10:2582091:grid_template=40:winds(N/S):
+	Gaussian grid: (768 x 384) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=192 #points=294912
+	lat 89.641648 to -89.641648
+	lon 0.000000 to 359.531250 by 0.468750
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89641648 to -89641648 lon=0 to 359531250:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl:
+11:2840318:grid_template=40:winds(N/S):
+	Gaussian grid: (768 x 384) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=192 #points=294912
+	lat 89.641648 to -89.641648
+	lon 0.000000 to 359.531250 by 0.468750
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89641648 to -89641648 lon=0 to 359531250:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl:
+12:3135409:grid_template=40:winds(N/S):
+	Gaussian grid: (768 x 384) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=192 #points=294912
+	lat 89.641648 to -89.641648
+	lon 0.000000 to 359.531250 by 0.468750
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89641648 to -89641648 lon=0 to 359531250:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl:
+13:3393636:grid_template=40:winds(N/S):
+	Gaussian grid: (768 x 384) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=192 #points=294912
+	lat 89.641648 to -89.641648
+	lon 0.000000 to 359.531250 by 0.468750
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89641648 to -89641648 lon=0 to 359531250:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl:
+14:3688727:grid_template=40:winds(N/S):
+	Gaussian grid: (768 x 384) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=192 #points=294912
+	lat 89.641648 to -89.641648
+	lon 0.000000 to 359.531250 by 0.468750
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89641648 to -89641648 lon=0 to 359531250:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl:
diff --git a/tests/data/ref_new_grid_ncep_t382.txt b/tests/data/ref_new_grid_ncep_t382.txt
new file mode 100644
index 00000000..4be18344
--- /dev/null
+++ b/tests/data/ref_new_grid_ncep_t382.txt
@@ -0,0 +1,98 @@
+1:0:grid_template=40:winds(N/S):
+	Gaussian grid: (1152 x 576) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=288 #points=663552
+	lat 89.760995 to -89.760995
+	lon 0.000000 to 359.687500 by 0.312500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89760995 to -89760995 lon=0 to 359687500:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl:
+2:663731:grid_template=40:winds(N/S):
+	Gaussian grid: (1152 x 576) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=288 #points=663552
+	lat 89.760995 to -89.760995
+	lon 0.000000 to 359.687500 by 0.312500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89760995 to -89760995 lon=0 to 359687500:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl:
+3:1327462:grid_template=40:winds(N/S):
+	Gaussian grid: (1152 x 576) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=288 #points=663552
+	lat 89.760995 to -89.760995
+	lon 0.000000 to 359.687500 by 0.312500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89760995 to -89760995 lon=0 to 359687500:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl:
+4:1991193:grid_template=40:winds(N/S):
+	Gaussian grid: (1152 x 576) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=288 #points=663552
+	lat 89.760995 to -89.760995
+	lon 0.000000 to 359.687500 by 0.312500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89760995 to -89760995 lon=0 to 359687500:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl:
+5:2654924:grid_template=40:winds(N/S):
+	Gaussian grid: (1152 x 576) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=288 #points=663552
+	lat 89.760995 to -89.760995
+	lon 0.000000 to 359.687500 by 0.312500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89760995 to -89760995 lon=0 to 359687500:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl:
+6:3318655:grid_template=40:winds(N/S):
+	Gaussian grid: (1152 x 576) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=288 #points=663552
+	lat 89.760995 to -89.760995
+	lon 0.000000 to 359.687500 by 0.312500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89760995 to -89760995 lon=0 to 359687500:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl:
+7:3899442:grid_template=40:winds(N/S):
+	Gaussian grid: (1152 x 576) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=288 #points=663552
+	lat 89.760995 to -89.760995
+	lon 0.000000 to 359.687500 by 0.312500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89760995 to -89760995 lon=0 to 359687500:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl:
+8:4563173:grid_template=40:winds(N/S):
+	Gaussian grid: (1152 x 576) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=288 #points=663552
+	lat 89.760995 to -89.760995
+	lon 0.000000 to 359.687500 by 0.312500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89760995 to -89760995 lon=0 to 359687500:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl:
+9:5143960:grid_template=40:winds(N/S):
+	Gaussian grid: (1152 x 576) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=288 #points=663552
+	lat 89.760995 to -89.760995
+	lon 0.000000 to 359.687500 by 0.312500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89760995 to -89760995 lon=0 to 359687500:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl:
+10:5807691:grid_template=40:winds(N/S):
+	Gaussian grid: (1152 x 576) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=288 #points=663552
+	lat 89.760995 to -89.760995
+	lon 0.000000 to 359.687500 by 0.312500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89760995 to -89760995 lon=0 to 359687500:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl:
+11:6388478:grid_template=40:winds(N/S):
+	Gaussian grid: (1152 x 576) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=288 #points=663552
+	lat 89.760995 to -89.760995
+	lon 0.000000 to 359.687500 by 0.312500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89760995 to -89760995 lon=0 to 359687500:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl:
+12:7052209:grid_template=40:winds(N/S):
+	Gaussian grid: (1152 x 576) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=288 #points=663552
+	lat 89.760995 to -89.760995
+	lon 0.000000 to 359.687500 by 0.312500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89760995 to -89760995 lon=0 to 359687500:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl:
+13:7632996:grid_template=40:winds(N/S):
+	Gaussian grid: (1152 x 576) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=288 #points=663552
+	lat 89.760995 to -89.760995
+	lon 0.000000 to 359.687500 by 0.312500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89760995 to -89760995 lon=0 to 359687500:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl:
+14:8296727:grid_template=40:winds(N/S):
+	Gaussian grid: (1152 x 576) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=288 #points=663552
+	lat 89.760995 to -89.760995
+	lon 0.000000 to 359.687500 by 0.312500
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89760995 to -89760995 lon=0 to 359687500:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl:
diff --git a/tests/data/ref_new_grid_ncep_t574.txt b/tests/data/ref_new_grid_ncep_t574.txt
new file mode 100644
index 00000000..b1ffe0ee
--- /dev/null
+++ b/tests/data/ref_new_grid_ncep_t574.txt
@@ -0,0 +1,98 @@
+1:0:grid_template=40:winds(N/S):
+	Gaussian grid: (1760 x 880) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=440 #points=1548800
+	lat 89.843514 to -89.843514
+	lon 0.000000 to 359.795455 by 0.204545
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89843514 to -89843514 lon=0 to 359795455:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl:
+2:1548979:grid_template=40:winds(N/S):
+	Gaussian grid: (1760 x 880) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=440 #points=1548800
+	lat 89.843514 to -89.843514
+	lon 0.000000 to 359.795455 by 0.204545
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89843514 to -89843514 lon=0 to 359795455:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl:
+3:3097958:grid_template=40:winds(N/S):
+	Gaussian grid: (1760 x 880) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=440 #points=1548800
+	lat 89.843514 to -89.843514
+	lon 0.000000 to 359.795455 by 0.204545
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89843514 to -89843514 lon=0 to 359795455:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl:
+4:4646937:grid_template=40:winds(N/S):
+	Gaussian grid: (1760 x 880) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=440 #points=1548800
+	lat 89.843514 to -89.843514
+	lon 0.000000 to 359.795455 by 0.204545
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89843514 to -89843514 lon=0 to 359795455:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl:
+5:6195916:grid_template=40:winds(N/S):
+	Gaussian grid: (1760 x 880) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=440 #points=1548800
+	lat 89.843514 to -89.843514
+	lon 0.000000 to 359.795455 by 0.204545
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89843514 to -89843514 lon=0 to 359795455:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl:
+6:7744895:grid_template=40:winds(N/S):
+	Gaussian grid: (1760 x 880) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=440 #points=1548800
+	lat 89.843514 to -89.843514
+	lon 0.000000 to 359.795455 by 0.204545
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89843514 to -89843514 lon=0 to 359795455:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl:
+7:9100274:grid_template=40:winds(N/S):
+	Gaussian grid: (1760 x 880) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=440 #points=1548800
+	lat 89.843514 to -89.843514
+	lon 0.000000 to 359.795455 by 0.204545
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89843514 to -89843514 lon=0 to 359795455:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl:
+8:10649253:grid_template=40:winds(N/S):
+	Gaussian grid: (1760 x 880) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=440 #points=1548800
+	lat 89.843514 to -89.843514
+	lon 0.000000 to 359.795455 by 0.204545
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89843514 to -89843514 lon=0 to 359795455:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl:
+9:12004632:grid_template=40:winds(N/S):
+	Gaussian grid: (1760 x 880) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=440 #points=1548800
+	lat 89.843514 to -89.843514
+	lon 0.000000 to 359.795455 by 0.204545
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89843514 to -89843514 lon=0 to 359795455:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl:
+10:13553611:grid_template=40:winds(N/S):
+	Gaussian grid: (1760 x 880) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=440 #points=1548800
+	lat 89.843514 to -89.843514
+	lon 0.000000 to 359.795455 by 0.204545
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89843514 to -89843514 lon=0 to 359795455:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl:
+11:14908990:grid_template=40:winds(N/S):
+	Gaussian grid: (1760 x 880) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=440 #points=1548800
+	lat 89.843514 to -89.843514
+	lon 0.000000 to 359.795455 by 0.204545
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89843514 to -89843514 lon=0 to 359795455:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl:
+12:16457969:grid_template=40:winds(N/S):
+	Gaussian grid: (1760 x 880) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=440 #points=1548800
+	lat 89.843514 to -89.843514
+	lon 0.000000 to 359.795455 by 0.204545
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89843514 to -89843514 lon=0 to 359795455:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl:
+13:17813348:grid_template=40:winds(N/S):
+	Gaussian grid: (1760 x 880) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=440 #points=1548800
+	lat 89.843514 to -89.843514
+	lon 0.000000 to 359.795455 by 0.204545
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89843514 to -89843514 lon=0 to 359795455:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl:
+14:19362327:grid_template=40:winds(N/S):
+	Gaussian grid: (1760 x 880) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=440 #points=1548800
+	lat 89.843514 to -89.843514
+	lon 0.000000 to 359.795455 by 0.204545
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=89843514 to -89843514 lon=0 to 359795455:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl:
diff --git a/tests/data/ref_new_grid_ncep_t62.txt b/tests/data/ref_new_grid_ncep_t62.txt
new file mode 100644
index 00000000..19eae8e5
--- /dev/null
+++ b/tests/data/ref_new_grid_ncep_t62.txt
@@ -0,0 +1,98 @@
+1:0:grid_template=40:winds(N/S):
+	Gaussian grid: (192 x 94) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=47 #points=18048
+	lat 88.541950 to -88.541950
+	lon 0.000000 to 358.125000 by 1.875000
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=88541950 to -88541950 lon=0 to 358125000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl:
+2:18227:grid_template=40:winds(N/S):
+	Gaussian grid: (192 x 94) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=47 #points=18048
+	lat 88.541950 to -88.541950
+	lon 0.000000 to 358.125000 by 1.875000
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=88541950 to -88541950 lon=0 to 358125000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,1) lvl2=(255,missing):0.01 mb:anl:
+3:36454:grid_template=40:winds(N/S):
+	Gaussian grid: (192 x 94) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=47 #points=18048
+	lat 88.541950 to -88.541950
+	lon 0.000000 to 358.125000 by 1.875000
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=88541950 to -88541950 lon=0 to 358125000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl:
+4:54681:grid_template=40:winds(N/S):
+	Gaussian grid: (192 x 94) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=47 #points=18048
+	lat 88.541950 to -88.541950
+	lon 0.000000 to 358.125000 by 1.875000
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=88541950 to -88541950 lon=0 to 358125000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,2) lvl2=(255,missing):0.02 mb:anl:
+5:70652:grid_template=40:winds(N/S):
+	Gaussian grid: (192 x 94) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=47 #points=18048
+	lat 88.541950 to -88.541950
+	lon 0.000000 to 358.125000 by 1.875000
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=88541950 to -88541950 lon=0 to 358125000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl:
+6:88879:grid_template=40:winds(N/S):
+	Gaussian grid: (192 x 94) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=47 #points=18048
+	lat 88.541950 to -88.541950
+	lon 0.000000 to 358.125000 by 1.875000
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=88541950 to -88541950 lon=0 to 358125000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,4) lvl2=(255,missing):0.04 mb:anl:
+7:104850:grid_template=40:winds(N/S):
+	Gaussian grid: (192 x 94) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=47 #points=18048
+	lat 88.541950 to -88.541950
+	lon 0.000000 to 358.125000 by 1.875000
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=88541950 to -88541950 lon=0 to 358125000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl:
+8:123077:grid_template=40:winds(N/S):
+	Gaussian grid: (192 x 94) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=47 #points=18048
+	lat 88.541950 to -88.541950
+	lon 0.000000 to 358.125000 by 1.875000
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=88541950 to -88541950 lon=0 to 358125000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,7) lvl2=(255,missing):0.07 mb:anl:
+9:139048:grid_template=40:winds(N/S):
+	Gaussian grid: (192 x 94) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=47 #points=18048
+	lat 88.541950 to -88.541950
+	lon 0.000000 to 358.125000 by 1.875000
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=88541950 to -88541950 lon=0 to 358125000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl:
+10:157275:grid_template=40:winds(N/S):
+	Gaussian grid: (192 x 94) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=47 #points=18048
+	lat 88.541950 to -88.541950
+	lon 0.000000 to 358.125000 by 1.875000
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=88541950 to -88541950 lon=0 to 358125000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,10) lvl2=(255,missing):0.1 mb:anl:
+11:173246:grid_template=40:winds(N/S):
+	Gaussian grid: (192 x 94) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=47 #points=18048
+	lat 88.541950 to -88.541950
+	lon 0.000000 to 358.125000 by 1.875000
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=88541950 to -88541950 lon=0 to 358125000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl:
+12:191473:grid_template=40:winds(N/S):
+	Gaussian grid: (192 x 94) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=47 #points=18048
+	lat 88.541950 to -88.541950
+	lon 0.000000 to 358.125000 by 1.875000
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=88541950 to -88541950 lon=0 to 358125000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,20) lvl2=(255,missing):0.2 mb:anl:
+13:207444:grid_template=40:winds(N/S):
+	Gaussian grid: (192 x 94) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=47 #points=18048
+	lat 88.541950 to -88.541950
+	lon 0.000000 to 358.125000 by 1.875000
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=88541950 to -88541950 lon=0 to 358125000:12Z26apr2024:UGRD U-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl:
+14:225671:grid_template=40:winds(N/S):
+	Gaussian grid: (192 x 94) units 1e-06 input WE:NS output WE:SN
+	number of latitudes between pole-equator=47 #points=18048
+	lat 88.541950 to -88.541950
+	lon 0.000000 to 358.125000 by 1.875000
+	basic_ang=0 sub_angle=1000000 units=0.000001
+	unscaled lat=88541950 to -88541950 lon=0 to 358125000:12Z26apr2024:VGRD V-Component of Wind [m/s]:lvl1=(100,40) lvl2=(255,missing):0.4 mb:anl:
diff --git a/tests/run_ipolates_tests.sh b/tests/run_ipolates_tests.sh
index 19a43e21..dd1694ab 100644
--- a/tests/run_ipolates_tests.sh
+++ b/tests/run_ipolates_tests.sh
@@ -5,11 +5,11 @@
 
 set -e
 
-echo "create a new grib file with intger values"
-../wgrib2/wgrib2 data/gdas.t12z.pgrb2.1p00.anl.75r.grib2 -rpn floor -set_scaling 0 0 -grib_out test.grb
+echo "create a new grib file with integer values"
+../wgrib2/wgrib2 data/gdas.t12z.pgrb2.1p00.anl.75r.grib2 -match "(UGRD|VGRD)" -rpn floor -set_scaling 0 0 -grib_out test.grb
 
 echo "*** Testing conversion from earth to grid"
-../wgrib2/wgrib2 test.grb -match "(UGRD|VGRD)" -new_grid_winds grid  \
+../wgrib2/wgrib2 test.grb -new_grid_winds grid  \
     -new_grid latlon 0:360:1 00:91:1 new_grid.grb
 ../wgrib2/wgrib2 new_grid.grb -grid -v2 -s -lon 10 12 -lon 20 80 > new_grid.txt
 touch new_grid.txt
@@ -22,5 +22,183 @@ echo "*** Testing conversion from grid to earth"
 touch new_grid_earth.txt
 diff -w new_grid_earth.txt data/ref_new_grid_earth_gdas.t12z.pgrb2.1p00.anl.75r.grib2.txt
 
+echo "*** Testing new_grid on file in incorrect order. This will return an incomplete output file and error message."
+../wgrib2/wgrib2 test.grb -match ":UGRD:" -grib_out test_badorder.grb
+../wgrib2/wgrib2 test.grb -match ":VGRD:" -append -grib_out test_badorder.grb
+{ out_err=$(../wgrib2/wgrib2 test_badorder.grb -new_grid_winds grid  \
+    -new_grid latlon 0:360:1 00:91:1 junk_badorder.grb 2>&1 1>&$out); } {out}>&1
+if [[ -z "$out_err" ]]; then 
+    exit 10
+fi
+
+echo "*** Testing new_grid_order on file in incorrect order."
+../wgrib2/wgrib2 test_badorder.grb -new_grid_order - junk | \
+    ../wgrib2/wgrib2 - -new_grid_winds grid -new_grid latlon 0:360:1 00:91:1 new_grid_reorder.grb
+../wgrib2/wgrib2 new_grid_reorder.grb -grid -v2 -s -lon 10 12 -lon 20 80 > new_grid_reorder.txt
+touch new_grid_reorder.txt
+diff -w new_grid_reorder.txt new_grid.txt
+
+echo "*** Testing conversion to NCEP grid definition 2"
+../wgrib2/wgrib2 test.grb -new_grid_winds earth -new_grid ncep grid 2 junk.grb
+../wgrib2/wgrib2 junk.grb -grid -v2 -s -lon 10 12 -lon 20 80 > ncep_grid_2.txt
+touch ncep_grid_2.txt
+diff -w ncep_grid_2.txt data/ref_new_grid_ncep_2.txt
+
+echo "*** Testing conversion to NCEP grid definition 3"
+../wgrib2/wgrib2 test.grb -new_grid_winds earth -new_grid ncep grid 3 junk.grb
+../wgrib2/wgrib2 junk.grb -grid -v2 -s -lon 10 12 -lon 20 80 > ncep_grid_3.txt
+touch ncep_grid_3.txt
+diff -w ncep_grid_3.txt data/ref_new_grid_ncep_3.txt
+
+echo "*** Testing conversion to NCEP grid definition 4"
+../wgrib2/wgrib2 test.grb -new_grid_winds earth -new_grid ncep grid 4 junk.grb
+../wgrib2/wgrib2 junk.grb -grid -v2 -s -lon 10 12 -lon 20 80 > ncep_grid_4.txt
+touch ncep_grid_4.txt
+diff -w ncep_grid_4.txt data/ref_new_grid_ncep_4.txt
+
+echo "*** Testing conversion to NCEP grid definition 45"
+../wgrib2/wgrib2 test.grb -new_grid_winds earth -new_grid ncep grid 45 junk.grb
+../wgrib2/wgrib2 junk.grb -grid -v2 -s -lon 10 12 -lon 20 80 > ncep_grid_45.txt
+touch ncep_grid_45.txt
+diff -w ncep_grid_45.txt data/ref_new_grid_ncep_45.txt
+
+echo "*** Testing conversion to NCEP grid definition 98"
+../wgrib2/wgrib2 test.grb -new_grid_winds earth -new_grid ncep grid 98 junk.grb
+../wgrib2/wgrib2 junk.grb -grid -v2 -s -lon 10 12 -lon 20 80 > ncep_grid_98.txt
+touch ncep_grid_98.txt
+diff -w ncep_grid_98.txt data/ref_new_grid_ncep_98.txt
+
+echo "*** Testing conversion to NCEP grid definition 126"
+../wgrib2/wgrib2 test.grb -new_grid_winds earth -new_grid ncep grid 126 junk.grb
+../wgrib2/wgrib2 junk.grb -grid -v2 -s -lon 10 12 -lon 20 80 > ncep_grid_126.txt
+touch ncep_grid_126.txt
+diff -w ncep_grid_126.txt data/ref_new_grid_ncep_126.txt
+
+echo "*** Testing conversion to NCEP grid definition 127"
+../wgrib2/wgrib2 test.grb -new_grid_winds earth -new_grid ncep grid 127 junk.grb
+../wgrib2/wgrib2 junk.grb -grid -v2 -s -lon 10 12 -lon 20 80 > ncep_grid_127.txt
+touch ncep_grid_127.txt
+diff -w ncep_grid_127.txt data/ref_new_grid_ncep_127.txt
+
+echo "*** Testing conversion to NCEP grid definition 128"
+../wgrib2/wgrib2 test.grb -new_grid_winds earth -new_grid ncep grid 128 junk.grb
+../wgrib2/wgrib2 junk.grb -grid -v2 -s -lon 10 12 -lon 20 80 > ncep_grid_128.txt
+touch ncep_grid_128.txt
+diff -w ncep_grid_128.txt data/ref_new_grid_ncep_128.txt
+
+echo "*** Testing conversion to NCEP grid definition 129"
+../wgrib2/wgrib2 test.grb -new_grid_winds earth -new_grid ncep grid 129 junk.grb
+../wgrib2/wgrib2 junk.grb -grid -v2 -s -lon 10 12 -lon 20 80 > ncep_grid_129.txt
+touch ncep_grid_129.txt
+diff -w ncep_grid_129.txt data/ref_new_grid_ncep_129.txt
+
+echo "*** Testing conversion to NCEP grid definition 170"
+../wgrib2/wgrib2 test.grb -new_grid_winds earth -new_grid ncep grid 170 junk.grb
+../wgrib2/wgrib2 junk.grb -grid -v2 -s -lon 10 12 -lon 20 80 > ncep_grid_170.txt
+touch ncep_grid_170.txt
+diff -w ncep_grid_170.txt data/ref_new_grid_ncep_170.txt
+
+echo "*** Testing conversion to NCEP grid definition 173"
+../wgrib2/wgrib2 test.grb -new_grid_winds earth -new_grid ncep grid 173 junk.grb
+../wgrib2/wgrib2 junk.grb -grid -v2 -s -lon 10 12 -lon 20 80 > ncep_grid_173.txt
+touch ncep_grid_173.txt
+diff -w ncep_grid_173.txt data/ref_new_grid_ncep_173.txt
+
+echo "*** Testing conversion to NCEP grid definition 184"
+../wgrib2/wgrib2 test.grb -new_grid_winds earth -new_grid ncep grid 184 junk.grb
+../wgrib2/wgrib2 junk.grb -grid -v2 -s > ncep_grid_184.txt
+touch ncep_grid_184.txt
+diff -w ncep_grid_184.txt data/ref_new_grid_ncep_184.txt
+
+echo "*** Testing conversion to NCEP grid definition 194"
+../wgrib2/wgrib2 test.grb -new_grid_winds earth -new_grid ncep grid 194 junk.grb
+../wgrib2/wgrib2 junk.grb -grid -v2 -s > ncep_grid_194.txt
+touch ncep_grid_194.txt
+diff -w ncep_grid_194.txt data/ref_new_grid_ncep_194.txt
+
+echo "*** Testing conversion to NCEP grid definition 221"
+../wgrib2/wgrib2 test.grb -new_grid_winds earth -new_grid ncep grid 221 junk.grb
+../wgrib2/wgrib2 junk.grb -grid -v2 -s > ncep_grid_221.txt
+touch ncep_grid_221.txt
+diff -w ncep_grid_221.txt data/ref_new_grid_ncep_221.txt
+
+echo "*** Testing conversion to NCEP grid definition 230"
+../wgrib2/wgrib2 test.grb -new_grid_winds earth -new_grid ncep grid 230 junk.grb
+../wgrib2/wgrib2 junk.grb -grid -v2 -s > ncep_grid_230.txt
+touch ncep_grid_230.txt
+diff -w ncep_grid_230.txt data/ref_new_grid_ncep_230.txt
+
+echo "*** Testing conversion to NCEP grid definition 242"
+../wgrib2/wgrib2 test.grb -new_grid_winds earth -new_grid ncep grid 242 junk.grb
+../wgrib2/wgrib2 junk.grb -grid -v2 -s > ncep_grid_242.txt
+touch ncep_grid_242.txt
+diff -w ncep_grid_242.txt data/ref_new_grid_ncep_242.txt
+
+echo "*** Testing conversion to NCEP grid definition 249"
+../wgrib2/wgrib2 test.grb -new_grid_winds earth -new_grid ncep grid 249 junk.grb
+../wgrib2/wgrib2 junk.grb -grid -v2 -s > ncep_grid_249.txt
+touch ncep_grid_249.txt
+diff -w ncep_grid_249.txt data/ref_new_grid_ncep_249.txt
+
+echo "*** Testing conversion to NCEP grid definition t62"
+../wgrib2/wgrib2 test.grb -new_grid_winds earth -new_grid ncep grid t62 junk.grb
+../wgrib2/wgrib2 junk.grb -grid -v2 -s > ncep_grid_t62.txt
+touch ncep_grid_t62.txt
+diff -w ncep_grid_t62.txt data/ref_new_grid_ncep_t62.txt
+
+echo "*** Testing conversion to NCEP grid definition t126"
+../wgrib2/wgrib2 test.grb -new_grid_winds earth -new_grid ncep grid t126 junk.grb
+../wgrib2/wgrib2 junk.grb -grid -v2 -s > ncep_grid_t126.txt
+touch ncep_grid_t126.txt
+diff -w ncep_grid_t126.txt data/ref_new_grid_ncep_t126.txt
+
+echo "*** Testing conversion to NCEP grid definition t170"
+../wgrib2/wgrib2 test.grb -new_grid_winds earth -new_grid ncep grid t170 junk.grb
+../wgrib2/wgrib2 junk.grb -grid -v2 -s > ncep_grid_t170.txt
+touch ncep_grid_t170.txt
+diff -w ncep_grid_t170.txt data/ref_new_grid_ncep_t170.txt
+
+echo "*** Testing conversion to NCEP grid definition t190"
+../wgrib2/wgrib2 test.grb -new_grid_winds earth -new_grid ncep grid t190 junk.grb
+../wgrib2/wgrib2 junk.grb -grid -v2 -s > ncep_grid_t190.txt
+touch ncep_grid_t190.txt
+diff -w ncep_grid_t190.txt data/ref_new_grid_ncep_t190.txt
+
+echo "*** Testing conversion to NCEP grid definition t254"
+../wgrib2/wgrib2 test.grb -new_grid_winds earth -new_grid ncep grid t254 junk.grb
+../wgrib2/wgrib2 junk.grb -grid -v2 -s > ncep_grid_t254.txt
+touch ncep_grid_t254.txt
+diff -w ncep_grid_t254.txt data/ref_new_grid_ncep_t254.txt
+
+echo "*** Testing conversion to NCEP grid definition t382"
+../wgrib2/wgrib2 test.grb -new_grid_winds earth -new_grid ncep grid t382 junk.grb
+../wgrib2/wgrib2 junk.grb -grid -v2 -s > ncep_grid_t382.txt
+touch ncep_grid_t382.txt
+diff -w ncep_grid_t382.txt data/ref_new_grid_ncep_t382.txt
+
+echo "*** Testing conversion to NCEP grid definition t574"
+../wgrib2/wgrib2 test.grb -new_grid_winds earth -new_grid ncep grid t574 junk.grb
+../wgrib2/wgrib2 junk.grb -grid -v2 -s > ncep_grid_t574.txt
+touch ncep_grid_t574.txt
+diff -w ncep_grid_t574.txt data/ref_new_grid_ncep_t574.txt
+
+echo "*** Testing conversion to NCEP grid definition t1148"
+../wgrib2/wgrib2 test.grb -new_grid_winds earth -new_grid ncep grid t1148 junk.grb
+../wgrib2/wgrib2 junk.grb -grid -v2 -s > ncep_grid_t1148.txt
+touch ncep_grid_t1148.txt
+diff -w ncep_grid_t1148.txt data/ref_new_grid_ncep_t1148.txt
+
+echo "*** Testing conversion to NCEP grid definition t1534"
+../wgrib2/wgrib2 test.grb -new_grid_winds earth -new_grid ncep grid t1534 junk.grb
+../wgrib2/wgrib2 junk.grb -grid -v2 -s > ncep_grid_t1534.txt
+touch ncep_grid_t1534.txt
+diff -w ncep_grid_t1534.txt data/ref_new_grid_ncep_t1534.txt
+
+echo "*** Testing conversion to Mercator grid"
+../wgrib2/wgrib2 test.grb -new_grid_winds earth -new_grid mercator:0 0:361:1:360 0:91:1:90 junk.grb
+../wgrib2/wgrib2 junk.grb -grid -v2 -s > mercator_grid.txt
+touch mercator_grid.txt
+diff -w mercator_grid.txt data/ref_mercator_grid.txt
+
 echo "*** SUCCESS!"
 exit 0