Add some hooks to call the search routine.

Fixes #279

sorc/chgres_cube.fd/search_util.F90

@@ -118,7 +118,7 @@ subroutine search (field, mask, idim, jdim, tile, field_num, latitude, terrain_l
    case (230) ! soil type on the input grid
      default_value = 11.0
    case default
-     print*,'- FATAL ERROR.  UNIDENTIFIED FIELD NUMBER : ', field
+     print*,'- FATAL ERROR.  UNIDENTIFIED FIELD NUMBER : ', field_num
      call mpi_abort(mpi_comm_world, 77, ierr)
  end select
 

tests/chres_cube/CMakeLists.txt

@@ -70,5 +70,13 @@ if(OpenMP_Fortran_FOUND)
 endif()
 
 add_executable(ftst_search_util
+  "${CMAKE_SOURCE_DIR}/sorc/chgres_cube.fd/search_util.F90"
   ftst_search_util.F90)
 add_test(NAME ftst_search_util COMMAND ftst_search_util)
+target_link_libraries(
+  ftst_search_util
+  esmf
+  MPI::MPI_Fortran)
+if(OpenMP_Fortran_FOUND)
+  target_link_libraries(ftst_search_util OpenMP::OpenMP_Fortran)
+endif()

tests/chres_cube/ftst_search_util.F90

@@ -1,7 +1,42 @@
  program test_search_util
 
+ use mpi
+ use esmf
+ use search_util
+
  implicit none
 
+ integer, parameter :: idim = 3
+ integer, parameter :: jdim = 3
+ integer, parameter :: tile = 1
+
+ integer :: field_num, ierr
+ integer(esmf_kind_i8) :: mask(idim,jdim)
+
+ real(esmf_kind_r8) :: field(idim,jdim)
+ real(esmf_kind_r8) :: latitude(idim,jdim)
+ real(esmf_kind_r8) :: terrain_land(idim,jdim)
+ real(esmf_kind_r8) :: soilt_climo(idim,jdim)
+
+ data mask /0, 0, 0,  1, 1, 1, 0, 0, 0/
+
+ data latitude /-30.0, -30.0, -30.0, 0., 0., 0., 25.0, 25.0, 25.0/
+
+ data terrain_land /0., 0., 0., 50.0, 75.0, 80.0, 0., 0., 0./
+
+ data soilt_climo /285.0, 285.0, 285.0, 0., 0., 0., 278.0, 278.0, 278.0/
+
+ data field/0., 0., 0.,  0., 0., 0.,  0., 0., 0./
+
+ call mpi_init(ierr)
+
  print*,'hello world'
 
+ field_num = 0
+
+ call search (field, mask, idim, jdim, tile, field_num, latitude, terrain_land, &
+       soilt_climo)
+
+ call mpi_finalize(ierr)
+
  end program test_search_util