diff --git a/CMakeLists.txt b/CMakeLists.txt index ea49fdd8..7bed756a 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" off) +option(USE_IPOLATES "Use Ipolates" on) option(USE_UDF "Use UDF?" off) option(USE_OPENMP "Use OpenMP?" on) option(USE_PROJ4 "Use Proj4?" off) @@ -94,7 +94,7 @@ endif() message(STATUS "Checking if the user want to use NCEPLIBS-ip...") if(USE_IPOLATES) find_package(ip CONFIG REQUIRED) - list(APPEND definitions_list -DIPOLATES_LIB="ipolates_lib_4") + list(APPEND definitions_list -DIPOLATES_LIB="ipolates_lib_d") list(APPEND definitions_list -DUSE_IPOLATES) endif() diff --git a/tests/data/ref_new_grid_earth_gdas.t12z.pgrb2.1p00.anl.75r.grib2.txt b/tests/data/ref_new_grid_earth_gdas.t12z.pgrb2.1p00.anl.75r.grib2.txt index 0c1ab74d..07ccd65c 100644 --- a/tests/data/ref_new_grid_earth_gdas.t12z.pgrb2.1p00.anl.75r.grib2.txt +++ b/tests/data/ref_new_grid_earth_gdas.t12z.pgrb2.1p00.anl.75r.grib2.txt @@ -1,84 +1,84 @@ -1:0:grid_template=0:winds(N/S): - lat-lon grid:(360 x 181) units 1e-06 input WE:SN output WE:SN res 48 - lat 0.000000 to 180.000000 by 1.000000 - lon 0.000000 to 359.000000 by 1.000000 #points=65160 +1:0:grid_template=0:winds(grid): + lat-lon grid:(360 x 91) units 1e-06 input WE:SN output WE:SN res 56 + lat 0.000000 to 90.000000 by 1.000000 + lon 0.000000 to 359.000000 by 1.000000 #points=32760 basic_ang=0 sub_angle=1000000 units=0.000001 - unscaled lat=0 to 180000000 lon=0 to 359000000 -2:8324:grid_template=0:winds(N/S): - lat-lon grid:(360 x 181) units 1e-06 input WE:SN output WE:SN res 48 - lat 0.000000 to 180.000000 by 1.000000 - lon 0.000000 to 359.000000 by 1.000000 #points=65160 + unscaled lat=0 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=4331,ix=11,iy=13,val=-43:lon=20.000000,lat=80.000000,i=28821,ix=21,iy=81,val=-17 +2:28844:grid_template=0:winds(grid): + lat-lon grid:(360 x 91) units 1e-06 input WE:SN output WE:SN res 56 + lat 0.000000 to 90.000000 by 1.000000 + lon 0.000000 to 359.000000 by 1.000000 #points=32760 basic_ang=0 sub_angle=1000000 units=0.000001 - unscaled lat=0 to 180000000 lon=0 to 359000000 -3:16648:grid_template=0:winds(N/S): - lat-lon grid:(360 x 181) units 1e-06 input WE:SN output WE:SN res 48 - lat 0.000000 to 180.000000 by 1.000000 - lon 0.000000 to 359.000000 by 1.000000 #points=65160 + unscaled lat=0 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=4331,ix=11,iy=13,val=-67:lon=20.000000,lat=80.000000,i=28821,ix=21,iy=81,val=-18 +3:57688:grid_template=0:winds(grid): + lat-lon grid:(360 x 91) units 1e-06 input WE:SN output WE:SN res 56 + lat 0.000000 to 90.000000 by 1.000000 + lon 0.000000 to 359.000000 by 1.000000 #points=32760 basic_ang=0 sub_angle=1000000 units=0.000001 - unscaled lat=0 to 180000000 lon=0 to 359000000 -4:24972:grid_template=0:winds(N/S): - lat-lon grid:(360 x 181) units 1e-06 input WE:SN output WE:SN res 48 - lat 0.000000 to 180.000000 by 1.000000 - lon 0.000000 to 359.000000 by 1.000000 #points=65160 + unscaled lat=0 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=4331,ix=11,iy=13,val=-28:lon=20.000000,lat=80.000000,i=28821,ix=21,iy=81,val=-14 +4:86532:grid_template=0:winds(grid): + lat-lon grid:(360 x 91) units 1e-06 input WE:SN output WE:SN res 56 + lat 0.000000 to 90.000000 by 1.000000 + lon 0.000000 to 359.000000 by 1.000000 #points=32760 basic_ang=0 sub_angle=1000000 units=0.000001 - unscaled lat=0 to 180000000 lon=0 to 359000000 -5:33296:grid_template=0:winds(N/S): - lat-lon grid:(360 x 181) units 1e-06 input WE:SN output WE:SN res 48 - lat 0.000000 to 180.000000 by 1.000000 - lon 0.000000 to 359.000000 by 1.000000 #points=65160 + unscaled lat=0 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=4331,ix=11,iy=13,val=-58:lon=20.000000,lat=80.000000,i=28821,ix=21,iy=81,val=-23 +5:115376:grid_template=0:winds(grid): + lat-lon grid:(360 x 91) units 1e-06 input WE:SN output WE:SN res 56 + lat 0.000000 to 90.000000 by 1.000000 + lon 0.000000 to 359.000000 by 1.000000 #points=32760 basic_ang=0 sub_angle=1000000 units=0.000001 - unscaled lat=0 to 180000000 lon=0 to 359000000 -6:41620:grid_template=0:winds(N/S): - lat-lon grid:(360 x 181) units 1e-06 input WE:SN output WE:SN res 48 - lat 0.000000 to 180.000000 by 1.000000 - lon 0.000000 to 359.000000 by 1.000000 #points=65160 + unscaled lat=0 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=4331,ix=11,iy=13,val=-6:lon=20.000000,lat=80.000000,i=28821,ix=21,iy=81,val=1 +6:144220:grid_template=0:winds(grid): + lat-lon grid:(360 x 91) units 1e-06 input WE:SN output WE:SN res 56 + lat 0.000000 to 90.000000 by 1.000000 + lon 0.000000 to 359.000000 by 1.000000 #points=32760 basic_ang=0 sub_angle=1000000 units=0.000001 - unscaled lat=0 to 180000000 lon=0 to 359000000 -7:49944:grid_template=0:winds(N/S): - lat-lon grid:(360 x 181) units 1e-06 input WE:SN output WE:SN res 48 - lat 0.000000 to 180.000000 by 1.000000 - lon 0.000000 to 359.000000 by 1.000000 #points=65160 + unscaled lat=0 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=4331,ix=11,iy=13,val=-27:lon=20.000000,lat=80.000000,i=28821,ix=21,iy=81,val=-16 +7:173064:grid_template=0:winds(grid): + lat-lon grid:(360 x 91) units 1e-06 input WE:SN output WE:SN res 56 + lat 0.000000 to 90.000000 by 1.000000 + lon 0.000000 to 359.000000 by 1.000000 #points=32760 basic_ang=0 sub_angle=1000000 units=0.000001 - unscaled lat=0 to 180000000 lon=0 to 359000000 -8:58268:grid_template=0:winds(N/S): - lat-lon grid:(360 x 181) units 1e-06 input WE:SN output WE:SN res 48 - lat 0.000000 to 180.000000 by 1.000000 - lon 0.000000 to 359.000000 by 1.000000 #points=65160 + unscaled lat=0 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=4331,ix=11,iy=13,val=-1:lon=20.000000,lat=80.000000,i=28821,ix=21,iy=81,val=10 +8:201908:grid_template=0:winds(grid): + lat-lon grid:(360 x 91) units 1e-06 input WE:SN output WE:SN res 56 + lat 0.000000 to 90.000000 by 1.000000 + lon 0.000000 to 359.000000 by 1.000000 #points=32760 basic_ang=0 sub_angle=1000000 units=0.000001 - unscaled lat=0 to 180000000 lon=0 to 359000000 -9:66592:grid_template=0:winds(N/S): - lat-lon grid:(360 x 181) units 1e-06 input WE:SN output WE:SN res 48 - lat 0.000000 to 180.000000 by 1.000000 - lon 0.000000 to 359.000000 by 1.000000 #points=65160 + unscaled lat=0 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=4331,ix=11,iy=13,val=6:lon=20.000000,lat=80.000000,i=28821,ix=21,iy=81,val=-3 +9:230752:grid_template=0:winds(grid): + lat-lon grid:(360 x 91) units 1e-06 input WE:SN output WE:SN res 56 + lat 0.000000 to 90.000000 by 1.000000 + lon 0.000000 to 359.000000 by 1.000000 #points=32760 basic_ang=0 sub_angle=1000000 units=0.000001 - unscaled lat=0 to 180000000 lon=0 to 359000000 -10:74916:grid_template=0:winds(N/S): - lat-lon grid:(360 x 181) units 1e-06 input WE:SN output WE:SN res 48 - lat 0.000000 to 180.000000 by 1.000000 - lon 0.000000 to 359.000000 by 1.000000 #points=65160 + unscaled lat=0 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=4331,ix=11,iy=13,val=-10:lon=20.000000,lat=80.000000,i=28821,ix=21,iy=81,val=8 +10:259596:grid_template=0:winds(grid): + lat-lon grid:(360 x 91) units 1e-06 input WE:SN output WE:SN res 56 + lat 0.000000 to 90.000000 by 1.000000 + lon 0.000000 to 359.000000 by 1.000000 #points=32760 basic_ang=0 sub_angle=1000000 units=0.000001 - unscaled lat=0 to 180000000 lon=0 to 359000000 -11:83240:grid_template=0:winds(N/S): - lat-lon grid:(360 x 181) units 1e-06 input WE:SN output WE:SN res 48 - lat 0.000000 to 180.000000 by 1.000000 - lon 0.000000 to 359.000000 by 1.000000 #points=65160 + unscaled lat=0 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=4331,ix=11,iy=13,val=11:lon=20.000000,lat=80.000000,i=28821,ix=21,iy=81,val=-1 +11:288440:grid_template=0:winds(grid): + lat-lon grid:(360 x 91) units 1e-06 input WE:SN output WE:SN res 56 + lat 0.000000 to 90.000000 by 1.000000 + lon 0.000000 to 359.000000 by 1.000000 #points=32760 basic_ang=0 sub_angle=1000000 units=0.000001 - unscaled lat=0 to 180000000 lon=0 to 359000000 -12:91564:grid_template=0:winds(N/S): - lat-lon grid:(360 x 181) units 1e-06 input WE:SN output WE:SN res 48 - lat 0.000000 to 180.000000 by 1.000000 - lon 0.000000 to 359.000000 by 1.000000 #points=65160 + unscaled lat=0 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=4331,ix=11,iy=13,val=-8:lon=20.000000,lat=80.000000,i=28821,ix=21,iy=81,val=9 +12:313189:grid_template=0:winds(grid): + lat-lon grid:(360 x 91) units 1e-06 input WE:SN output WE:SN res 56 + lat 0.000000 to 90.000000 by 1.000000 + lon 0.000000 to 359.000000 by 1.000000 #points=32760 basic_ang=0 sub_angle=1000000 units=0.000001 - unscaled lat=0 to 180000000 lon=0 to 359000000 -13:99888:grid_template=0:winds(N/S): - lat-lon grid:(360 x 181) units 1e-06 input WE:SN output WE:SN res 48 - lat 0.000000 to 180.000000 by 1.000000 - lon 0.000000 to 359.000000 by 1.000000 #points=65160 + unscaled lat=0 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=4331,ix=11,iy=13,val=14:lon=20.000000,lat=80.000000,i=28821,ix=21,iy=81,val=8 +13:337938:grid_template=0:winds(grid): + lat-lon grid:(360 x 91) units 1e-06 input WE:SN output WE:SN res 56 + lat 0.000000 to 90.000000 by 1.000000 + lon 0.000000 to 359.000000 by 1.000000 #points=32760 basic_ang=0 sub_angle=1000000 units=0.000001 - unscaled lat=0 to 180000000 lon=0 to 359000000 -14:108212:grid_template=0:winds(N/S): - lat-lon grid:(360 x 181) units 1e-06 input WE:SN output WE:SN res 48 - lat 0.000000 to 180.000000 by 1.000000 - lon 0.000000 to 359.000000 by 1.000000 #points=65160 + unscaled lat=0 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=4331,ix=11,iy=13,val=-3:lon=20.000000,lat=80.000000,i=28821,ix=21,iy=81,val=-8 +14:362687:grid_template=0:winds(grid): + lat-lon grid:(360 x 91) units 1e-06 input WE:SN output WE:SN res 56 + lat 0.000000 to 90.000000 by 1.000000 + lon 0.000000 to 359.000000 by 1.000000 #points=32760 basic_ang=0 sub_angle=1000000 units=0.000001 - unscaled lat=0 to 180000000 lon=0 to 359000000 + unscaled lat=0 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=4331,ix=11,iy=13,val=15:lon=20.000000,lat=80.000000,i=28821,ix=21,iy=81,val=13 diff --git a/tests/data/ref_new_grid_gdas.t12z.pgrb2.1p00.anl.75r.grib2.txt b/tests/data/ref_new_grid_gdas.t12z.pgrb2.1p00.anl.75r.grib2.txt index b8a43f49..07ccd65c 100644 --- a/tests/data/ref_new_grid_gdas.t12z.pgrb2.1p00.anl.75r.grib2.txt +++ b/tests/data/ref_new_grid_gdas.t12z.pgrb2.1p00.anl.75r.grib2.txt @@ -1,84 +1,84 @@ 1:0:grid_template=0:winds(grid): - lat-lon grid:(360 x 181) units 1e-06 input WE:SN output WE:SN res 56 - lat 0.000000 to 180.000000 by 1.000000 - lon 0.000000 to 359.000000 by 1.000000 #points=65160 + lat-lon grid:(360 x 91) units 1e-06 input WE:SN output WE:SN res 56 + lat 0.000000 to 90.000000 by 1.000000 + lon 0.000000 to 359.000000 by 1.000000 #points=32760 basic_ang=0 sub_angle=1000000 units=0.000001 - unscaled lat=0 to 180000000 lon=0 to 359000000 -2:179:grid_template=0:winds(grid): - lat-lon grid:(360 x 181) units 1e-06 input WE:SN output WE:SN res 56 - lat 0.000000 to 180.000000 by 1.000000 - lon 0.000000 to 359.000000 by 1.000000 #points=65160 + unscaled lat=0 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=4331,ix=11,iy=13,val=-43:lon=20.000000,lat=80.000000,i=28821,ix=21,iy=81,val=-17 +2:28844:grid_template=0:winds(grid): + lat-lon grid:(360 x 91) units 1e-06 input WE:SN output WE:SN res 56 + lat 0.000000 to 90.000000 by 1.000000 + lon 0.000000 to 359.000000 by 1.000000 #points=32760 basic_ang=0 sub_angle=1000000 units=0.000001 - unscaled lat=0 to 180000000 lon=0 to 359000000 -3:358:grid_template=0:winds(grid): - lat-lon grid:(360 x 181) units 1e-06 input WE:SN output WE:SN res 56 - lat 0.000000 to 180.000000 by 1.000000 - lon 0.000000 to 359.000000 by 1.000000 #points=65160 + unscaled lat=0 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=4331,ix=11,iy=13,val=-67:lon=20.000000,lat=80.000000,i=28821,ix=21,iy=81,val=-18 +3:57688:grid_template=0:winds(grid): + lat-lon grid:(360 x 91) units 1e-06 input WE:SN output WE:SN res 56 + lat 0.000000 to 90.000000 by 1.000000 + lon 0.000000 to 359.000000 by 1.000000 #points=32760 basic_ang=0 sub_angle=1000000 units=0.000001 - unscaled lat=0 to 180000000 lon=0 to 359000000 -4:537:grid_template=0:winds(grid): - lat-lon grid:(360 x 181) units 1e-06 input WE:SN output WE:SN res 56 - lat 0.000000 to 180.000000 by 1.000000 - lon 0.000000 to 359.000000 by 1.000000 #points=65160 + unscaled lat=0 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=4331,ix=11,iy=13,val=-28:lon=20.000000,lat=80.000000,i=28821,ix=21,iy=81,val=-14 +4:86532:grid_template=0:winds(grid): + lat-lon grid:(360 x 91) units 1e-06 input WE:SN output WE:SN res 56 + lat 0.000000 to 90.000000 by 1.000000 + lon 0.000000 to 359.000000 by 1.000000 #points=32760 basic_ang=0 sub_angle=1000000 units=0.000001 - unscaled lat=0 to 180000000 lon=0 to 359000000 -5:716:grid_template=0:winds(grid): - lat-lon grid:(360 x 181) units 1e-06 input WE:SN output WE:SN res 56 - lat 0.000000 to 180.000000 by 1.000000 - lon 0.000000 to 359.000000 by 1.000000 #points=65160 + unscaled lat=0 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=4331,ix=11,iy=13,val=-58:lon=20.000000,lat=80.000000,i=28821,ix=21,iy=81,val=-23 +5:115376:grid_template=0:winds(grid): + lat-lon grid:(360 x 91) units 1e-06 input WE:SN output WE:SN res 56 + lat 0.000000 to 90.000000 by 1.000000 + lon 0.000000 to 359.000000 by 1.000000 #points=32760 basic_ang=0 sub_angle=1000000 units=0.000001 - unscaled lat=0 to 180000000 lon=0 to 359000000 -6:895:grid_template=0:winds(grid): - lat-lon grid:(360 x 181) units 1e-06 input WE:SN output WE:SN res 56 - lat 0.000000 to 180.000000 by 1.000000 - lon 0.000000 to 359.000000 by 1.000000 #points=65160 + unscaled lat=0 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=4331,ix=11,iy=13,val=-6:lon=20.000000,lat=80.000000,i=28821,ix=21,iy=81,val=1 +6:144220:grid_template=0:winds(grid): + lat-lon grid:(360 x 91) units 1e-06 input WE:SN output WE:SN res 56 + lat 0.000000 to 90.000000 by 1.000000 + lon 0.000000 to 359.000000 by 1.000000 #points=32760 basic_ang=0 sub_angle=1000000 units=0.000001 - unscaled lat=0 to 180000000 lon=0 to 359000000 -7:1074:grid_template=0:winds(grid): - lat-lon grid:(360 x 181) units 1e-06 input WE:SN output WE:SN res 56 - lat 0.000000 to 180.000000 by 1.000000 - lon 0.000000 to 359.000000 by 1.000000 #points=65160 + unscaled lat=0 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=4331,ix=11,iy=13,val=-27:lon=20.000000,lat=80.000000,i=28821,ix=21,iy=81,val=-16 +7:173064:grid_template=0:winds(grid): + lat-lon grid:(360 x 91) units 1e-06 input WE:SN output WE:SN res 56 + lat 0.000000 to 90.000000 by 1.000000 + lon 0.000000 to 359.000000 by 1.000000 #points=32760 basic_ang=0 sub_angle=1000000 units=0.000001 - unscaled lat=0 to 180000000 lon=0 to 359000000 -8:1253:grid_template=0:winds(grid): - lat-lon grid:(360 x 181) units 1e-06 input WE:SN output WE:SN res 56 - lat 0.000000 to 180.000000 by 1.000000 - lon 0.000000 to 359.000000 by 1.000000 #points=65160 + unscaled lat=0 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=4331,ix=11,iy=13,val=-1:lon=20.000000,lat=80.000000,i=28821,ix=21,iy=81,val=10 +8:201908:grid_template=0:winds(grid): + lat-lon grid:(360 x 91) units 1e-06 input WE:SN output WE:SN res 56 + lat 0.000000 to 90.000000 by 1.000000 + lon 0.000000 to 359.000000 by 1.000000 #points=32760 basic_ang=0 sub_angle=1000000 units=0.000001 - unscaled lat=0 to 180000000 lon=0 to 359000000 -9:1432:grid_template=0:winds(grid): - lat-lon grid:(360 x 181) units 1e-06 input WE:SN output WE:SN res 56 - lat 0.000000 to 180.000000 by 1.000000 - lon 0.000000 to 359.000000 by 1.000000 #points=65160 + unscaled lat=0 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=4331,ix=11,iy=13,val=6:lon=20.000000,lat=80.000000,i=28821,ix=21,iy=81,val=-3 +9:230752:grid_template=0:winds(grid): + lat-lon grid:(360 x 91) units 1e-06 input WE:SN output WE:SN res 56 + lat 0.000000 to 90.000000 by 1.000000 + lon 0.000000 to 359.000000 by 1.000000 #points=32760 basic_ang=0 sub_angle=1000000 units=0.000001 - unscaled lat=0 to 180000000 lon=0 to 359000000 -10:1611:grid_template=0:winds(grid): - lat-lon grid:(360 x 181) units 1e-06 input WE:SN output WE:SN res 56 - lat 0.000000 to 180.000000 by 1.000000 - lon 0.000000 to 359.000000 by 1.000000 #points=65160 + unscaled lat=0 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=4331,ix=11,iy=13,val=-10:lon=20.000000,lat=80.000000,i=28821,ix=21,iy=81,val=8 +10:259596:grid_template=0:winds(grid): + lat-lon grid:(360 x 91) units 1e-06 input WE:SN output WE:SN res 56 + lat 0.000000 to 90.000000 by 1.000000 + lon 0.000000 to 359.000000 by 1.000000 #points=32760 basic_ang=0 sub_angle=1000000 units=0.000001 - unscaled lat=0 to 180000000 lon=0 to 359000000 -11:1790:grid_template=0:winds(grid): - lat-lon grid:(360 x 181) units 1e-06 input WE:SN output WE:SN res 56 - lat 0.000000 to 180.000000 by 1.000000 - lon 0.000000 to 359.000000 by 1.000000 #points=65160 + unscaled lat=0 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=4331,ix=11,iy=13,val=11:lon=20.000000,lat=80.000000,i=28821,ix=21,iy=81,val=-1 +11:288440:grid_template=0:winds(grid): + lat-lon grid:(360 x 91) units 1e-06 input WE:SN output WE:SN res 56 + lat 0.000000 to 90.000000 by 1.000000 + lon 0.000000 to 359.000000 by 1.000000 #points=32760 basic_ang=0 sub_angle=1000000 units=0.000001 - unscaled lat=0 to 180000000 lon=0 to 359000000 -12:1969:grid_template=0:winds(grid): - lat-lon grid:(360 x 181) units 1e-06 input WE:SN output WE:SN res 56 - lat 0.000000 to 180.000000 by 1.000000 - lon 0.000000 to 359.000000 by 1.000000 #points=65160 + unscaled lat=0 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=4331,ix=11,iy=13,val=-8:lon=20.000000,lat=80.000000,i=28821,ix=21,iy=81,val=9 +12:313189:grid_template=0:winds(grid): + lat-lon grid:(360 x 91) units 1e-06 input WE:SN output WE:SN res 56 + lat 0.000000 to 90.000000 by 1.000000 + lon 0.000000 to 359.000000 by 1.000000 #points=32760 basic_ang=0 sub_angle=1000000 units=0.000001 - unscaled lat=0 to 180000000 lon=0 to 359000000 -13:2148:grid_template=0:winds(grid): - lat-lon grid:(360 x 181) units 1e-06 input WE:SN output WE:SN res 56 - lat 0.000000 to 180.000000 by 1.000000 - lon 0.000000 to 359.000000 by 1.000000 #points=65160 + unscaled lat=0 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=4331,ix=11,iy=13,val=14:lon=20.000000,lat=80.000000,i=28821,ix=21,iy=81,val=8 +13:337938:grid_template=0:winds(grid): + lat-lon grid:(360 x 91) units 1e-06 input WE:SN output WE:SN res 56 + lat 0.000000 to 90.000000 by 1.000000 + lon 0.000000 to 359.000000 by 1.000000 #points=32760 basic_ang=0 sub_angle=1000000 units=0.000001 - unscaled lat=0 to 180000000 lon=0 to 359000000 -14:2327:grid_template=0:winds(grid): - lat-lon grid:(360 x 181) units 1e-06 input WE:SN output WE:SN res 56 - lat 0.000000 to 180.000000 by 1.000000 - lon 0.000000 to 359.000000 by 1.000000 #points=65160 + unscaled lat=0 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=4331,ix=11,iy=13,val=-3:lon=20.000000,lat=80.000000,i=28821,ix=21,iy=81,val=-8 +14:362687:grid_template=0:winds(grid): + lat-lon grid:(360 x 91) units 1e-06 input WE:SN output WE:SN res 56 + lat 0.000000 to 90.000000 by 1.000000 + lon 0.000000 to 359.000000 by 1.000000 #points=32760 basic_ang=0 sub_angle=1000000 units=0.000001 - unscaled lat=0 to 180000000 lon=0 to 359000000 + unscaled lat=0 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=4331,ix=11,iy=13,val=15:lon=20.000000,lat=80.000000,i=28821,ix=21,iy=81,val=13 diff --git a/tests/run_ipolates_tests.sh b/tests/run_ipolates_tests.sh index b15574f0..19a43e21 100644 --- a/tests/run_ipolates_tests.sh +++ b/tests/run_ipolates_tests.sh @@ -5,19 +5,22 @@ 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 "*** Testing conversion from earth to grid" -../wgrib2/wgrib2 data/gdas.t12z.pgrb2.1p00.anl.75r.grib2 -match "(UGRD|VGRD)" -new_grid_winds grid \ - -new_grid latlon 0:360:1 00:181:1 new_grid.grb -../wgrib2/wgrib2 new_grid.grb -grid -v2 > new_grid.txt +../wgrib2/wgrib2 test.grb -match "(UGRD|VGRD)" -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 diff -w new_grid.txt data/ref_new_grid_gdas.t12z.pgrb2.1p00.anl.75r.grib2.txt echo "*** Testing conversion from grid to earth" ../wgrib2/wgrib2 new_grid.grb -new_grid_winds earth \ - -new_grid latlon 0:360:1 00:181:1 new_grid_earth.grb -../wgrib2/wgrib2 new_grid_earth.grb -grid -v2 > new_grid_earth.txt + -new_grid latlon 0:360:1 00:91:1 new_grid_earth.grb +../wgrib2/wgrib2 new_grid.grb -grid -v2 -s -lon 10 12 -lon 20 80 > new_grid_earth.txt 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 "*** SUCCESS!" -exit 0 \ No newline at end of file +exit 0 diff --git a/wgrib2/CMakeLists.txt b/wgrib2/CMakeLists.txt index 597d6959..7794dd14 100644 --- a/wgrib2/CMakeLists.txt +++ b/wgrib2/CMakeLists.txt @@ -106,7 +106,7 @@ if(USE_G2CLIB) endif() if(USE_IPOLATES) - target_link_libraries(obj_lib PUBLIC ip::ip_4) + target_link_libraries(obj_lib PUBLIC ip::ip_d) # Link to the Fortran runtime library for each compiler if using ip2. # The wgrib2 exectuable is created using the C compiler and @@ -133,9 +133,6 @@ target_link_libraries(wgrib2_exe PRIVATE obj_lib) if(BUILD_LIB) set(headers wgrib2_api.h wgrib2.h ${CMAKE_BINARY_DIR}/wgrib2/wgrib2_meta.h) target_link_libraries(wgrib2_lib PUBLIC gctpc) - # if(USE_IPOLATES EQUAL 1) - # target_link_libraries(wgrib2_lib PRIVATE ip:ip_4 sp:sp_4) - # endif() set_target_properties(wgrib2_lib PROPERTIES PUBLIC_HEADER "${headers}") target_compile_definitions(wgrib2_lib PUBLIC ${definitions_list})