Revert back to original computation of input grid lat/lons

to facilitate testing. Fixes ufs-community#591.
GeorgeGayno-NOAA · Nov 3, 2021 · 4ab7ca1 · 4ab7ca1
1 parent a4072e7
commit 4ab7ca1
Showing 1 changed file with 10 additions and 10 deletions.
diff --git a/sorc/chgres_cube.fd/model_grid.F90 b/sorc/chgres_cube.fd/model_grid.F90
@@ -1138,10 +1138,10 @@ subroutine define_input_grid_grib2(localpet, npets)
 
     do j = clb(2),cub(2)
       do i = clb(1), cub(1)
-!       lon_src_ptr(i,j)=real(longitude_one_tile(i,j),esmf_kind_r8)
-!       lat_src_ptr(i,j)=real(latitude_one_tile(i,j),esmf_kind_r8)
-        lon_src_ptr(i,j)=rlon(i,j)
-        lat_src_ptr(i,j)=rlat(i,j)
+        lon_src_ptr(i,j)=real(longitude_one_tile(i,j),esmf_kind_r8)
+        lat_src_ptr(i,j)=real(latitude_one_tile(i,j),esmf_kind_r8)
+!       lon_src_ptr(i,j)=rlon(i,j)
+!       lat_src_ptr(i,j)=rlat(i,j)
         if (abs(rlon(i,j)-longitude_one_tile(i,j)) > 0.1 .or. &
             abs(rlat(i,j)-latitude_one_tile(i,j)) > 0.1) then
           print*,'lat/lon diff ', rlon(i,j), longitude_one_tile(i,j), &
@@ -1227,8 +1227,8 @@ subroutine define_input_grid_grib2(localpet, npets)
           print*,'lat/lon corn diff ', rlonc(i,j), lon_src_ptr(i,j), &
                    rlatc(i,j),lat_src_ptr(i,j)
         endif
-          lon_src_ptr(i,j)=rlonc(i,j)
-          lat_src_ptr(i,j)=rlatc(i,j)
+!         lon_src_ptr(i,j)=rlonc(i,j)
+!         lat_src_ptr(i,j)=rlatc(i,j)
         enddo
        enddo
 
@@ -1250,10 +1250,10 @@ subroutine define_input_grid_grib2(localpet, npets)
 
     do j = clb(2),cub(2)
       do i = clb(1), cub(1)
-!       lon_src_ptr(i,j)=real(lon_corners(i,j),esmf_kind_r8)
-!       lat_src_ptr(i,j)=real(lat_corners(i,j),esmf_kind_r8)
-        lon_src_ptr(i,j)=rlonc(i,j)
-        lat_src_ptr(i,j)=rlatc(i,j)
+        lon_src_ptr(i,j)=real(lon_corners(i,j),esmf_kind_r8)
+        lat_src_ptr(i,j)=real(lat_corners(i,j),esmf_kind_r8)
+!       lon_src_ptr(i,j)=rlonc(i,j)
+!       lat_src_ptr(i,j)=rlatc(i,j)
         if (abs(rlonc(i,j)-lon_corners(i,j)) > 0.1 .or. &
             abs(rlatc(i,j)-lat_corners(i,j)) > 0.1) then
           print*,'lat/lon corn diff ', rlonc(i,j), lon_corners(i,j), &