relaxed tolerance for grid_is_ok (#4906)

Author: weiyuan-jiang

CHANGELOG.md

@@ -11,6 +11,8 @@ and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0
 
 ### Fixed
 
+- Relaxed the comparison standard for grid_is_ok in case a grid is r4
+
 ### Added
 
 ### Changed

geom/CubedSphere/CubedSphereGetGlobalHorzIJIndex.F90

@@ -295,8 +295,8 @@ logical function grid_is_ok_(grid, stretched, IM_WORLD, dalpha, alpha, shift) re
       accurate_lon = 1.750d0*MAPL_PI_R8 - shift0
       accurate_lat = [(-alpha + (j-1)*dalpha, j = J1, J2)]
 
-      if (any(abs(accurate_lon - lonRe) > 2.0*tolerance) .or. &
-          any(abs(accurate_lat - latRe) > 2.0*tolerance)) then
+      if (any(abs(accurate_lon - lonRe) > 20.0*tolerance) .or. &
+          any(abs(accurate_lat - latRe) > 20.0*tolerance)) then
         print*, "Error in grid_is_ok_: grid may not have pi/18 Japan mountain shift,"
         print*, "  may not be gnomonic_ed, or corners are not accurate enough."
         ok = .false.