Add testing for Arawaka grib2 cases

NOAA-EMC · Mar 2, 2024 · e3da1f3 · e3da1f3
1 parent 70841b8
commit e3da1f3
Show file tree

Hide file tree

Showing 6 changed files with 16 additions and 0 deletions.
diff --git a/tests/CMakeLists.txt b/tests/CMakeLists.txt
@@ -104,6 +104,8 @@ foreach(kind ${kinds})
   add_test(test_polar_stereo_neighbor_budget_vector_grib2_${kind} test_vector_grib2_${kind} 212 6)
   add_test(test_rotatedB_spectral_vector_grib2_${kind} test_vector_grib2_${kind} 205 4)
   add_test(test_rotatedE_budget_vector_grib2_${kind} test_vector_grib2_${kind} 203 3)
+  add_test(test_rotatedB_direct_spectral_vector_grib2_${kind} test_vector_grib2_${kind} 32769 4)
+  add_test(test_rotatedE_direct_budget_vector_grib2_${kind} test_vector_grib2_${kind} 32768 3)
 
   # vector station point tests
   add_test(test_station_points_bilinear_vector_grib2_${kind} test_vector_grib2_${kind} -1 0)

diff --git a/tests/data/baseline_data/grib2/vector/4_byte_bin/grid32768.opt3.bin_4 b/tests/data/baseline_data/grib2/vector/4_byte_bin/grid32768.opt3.bin_4
diff --git a/tests/data/baseline_data/grib2/vector/4_byte_bin/grid32769.opt4.bin_4 b/tests/data/baseline_data/grib2/vector/4_byte_bin/grid32769.opt4.bin_4
diff --git a/tests/data/baseline_data/grib2/vector/8_byte_bin/grid32768.opt3.bin_8 b/tests/data/baseline_data/grib2/vector/8_byte_bin/grid32768.opt3.bin_8
diff --git a/tests/data/baseline_data/grib2/vector/8_byte_bin/grid32769.opt4.bin_8 b/tests/data/baseline_data/grib2/vector/8_byte_bin/grid32769.opt4.bin_8
diff --git a/tests/interp_mod_grib2.F90 b/tests/interp_mod_grib2.F90
@@ -540,6 +540,20 @@ subroutine interp_vector(grid, interp_opt)
        output_gdtmpl=gdtmpl218
        i_output = output_gdtmpl(8)
        j_output = output_gdtmpl(9)
+    case ('32768')
+       output_gdtnum = 32768
+       output_gdtlen = gdtlen203
+       allocate(output_gdtmpl(output_gdtlen))
+       output_gdtmpl=gdtmpl203
+       i_output = output_gdtmpl(8)
+       j_output = output_gdtmpl(9)
+    case ('32769')
+       output_gdtnum = 32769
+       output_gdtlen = gdtlen205
+       allocate(output_gdtmpl(output_gdtlen))
+       output_gdtmpl=gdtmpl205
+       i_output = output_gdtmpl(8)
+       j_output = output_gdtmpl(9)
     case default
        print*,"ERROR: ENTER VALID GRID NUMBER."
        stop 55