fix issue when getting subregion

NOAA-GFDL · Aug 23, 2023 · e0cc678 · e0cc678
1 parent b8a83aa
commit e0cc678
Showing 1 changed file with 12 additions and 6 deletions.
diff --git a/diag_manager/fms_diag_axis_object.F90 b/diag_manager/fms_diag_axis_object.F90
@@ -1128,19 +1128,25 @@ subroutine define_subaxis_latlon(diag_axis, axis_ids, naxis, subRegion, is_cube_
           !< Get the starting and ending indices of the subregion relative to the global grid
           if (parent_axis%cart_name .eq. "X") then
             select type(adata=>parent_axis%axis_data)
-            type is (real)
-              lon_indices(1) = nearest_index(lon(1), adata)
-              lon_indices(2) = nearest_index(lon(2), adata) + 1
+            type is (real(kind=r8_kind))
+              lon_indices(1) = nearest_index(real(lon(1), kind=r8_kind), adata)
+              lon_indices(2) = nearest_index(real(lon(2), kind=r8_kind), adata) + 1
+            type is (real(kind=r4_kind))
+              lon_indices(1) = nearest_index(real(lon(1), kind=r4_kind), adata)
+              lon_indices(2) = nearest_index(real(lon(2), kind=r4_kind), adata) + 1
             end select
             call parent_axis%get_indices(compute_idx, lon_indices, starting_index(1), ending_index(1), &
               need_to_define_axis(1))
               parent_axis_ids(1) = axis_ids(i)
               compute_idx_2(1,:) = compute_idx
           else if (parent_axis%cart_name .eq. "Y") then
             select type(adata=>parent_axis%axis_data)
-            type is (real)
-              lat_indices(1) = nearest_index(lat(1), adata)
-              lat_indices(2) = nearest_index(lat(2), adata) + 1
+            type is (real(kind=r8_kind))
+              lat_indices(1) = nearest_index(real(lat(1), kind=r8_kind), adata)
+              lat_indices(2) = nearest_index(real(lat(2), kind=r8_kind), adata) + 1
+            type is (real(kind=r4_kind))
+              lat_indices(1) = nearest_index(real(lat(1), kind=r4_kind), adata)
+              lat_indices(2) = nearest_index(real(lat(2), kind=r4_kind), adata) + 1
             end select
             call parent_axis%get_indices(compute_idx, lat_indices, starting_index(2), ending_index(2), &
               need_to_define_axis(2))