add test

Author: mlee03

test_fms/horiz_interp/test_horiz_interp_save_weights_as_fregrid.F90

@@ -0,0 +1,37 @@
+program main
+
+  use constants_mod, only: DEG_TO_RAD
+  use horiz_interp_mod, only: horiz_interp_type, horiz_interp_new
+  use mpp_mod, only: FATAL, mpp_error
+  use platform_mod, only: r8_kind
+
+  implicit none
+
+  integer, parameter :: nlon = 12
+  integer, parameter :: nlat = 8
+
+  type(horiz_interp_type) :: interp
+  real(r8_kind),allocatable :: lon(:,:), lat(:,:), answer_area(:,:)
+  integer :: i, j
+
+  allocate(lon(nlon+1, nlat+1))
+  allocate(lat(nlon+1, nlat+1))
+  allocate(answer_area(nlon, nlat))
+
+  do j=1, nlat+1
+    do i=1, nlon+1
+      lon(i,j) = real(i, r8_kind)*DEG_TO_RAD
+      lat(i,j) = real(j, r8_kind)*DEG_TO_RAD
+    end do
+  end do
+
+  call get_grid_area(nlon, nlat,lon, lat, answer_area)
+
+  call horiz_interp_new(interp, lon, lat, lon, lat, interp_method="conservative", &
+    is_latlon_in=.false., is_latlon_out=.false., save_weights_as_fregrid=.true.)
+
+  if(any(interp%xgrid_area /= pack(answer_area, .true.))) then
+    call mpp_error(FATAL, "saved xgrid_area does not match answers")
+  end if
+
+end program main