chgres_cube_grib2_atmos: This commit references #66268

Fixed a new indexing error in creating the staggered grid coordinates.
Also removed the call to iso2sig when data is in isobaric coordinates.
3d pressure array is now generated from pressure levels read from file.

sorc/chgres_cube.fd/input_data.F90

@@ -2617,7 +2617,15 @@ subroutine read_input_atm_grib2_file(localpet)
     call error_handler("IN FieldGet", rc) 
   end do
 
-  call iso2sig(rlevs,vcoord,lev_input,levp1_input,psptr,atm,clb,cub,5+num_tracers, iret)
+  !call iso2sig(rlevs,vcoord,lev_input,levp1_input,psptr,atm,clb,cub,5+num_tracers, iret)
+  do i = clb(1),cub(1)
+    do j = clb(2),cub(2)
+      atm(1)%var(i,j,:) = rlevs(lev_input:1:-1)
+      do n = 2,num_tracers+5
+        atm(n)%var(i,j,:) = atm(n)%var(i,j,lev_input:1:-1)
+      end do
+    end do
+  end do  
   deallocate(vcoord)
 
  else

sorc/chgres_cube.fd/model_grid.F90

@@ -663,6 +663,7 @@ subroutine define_input_grid_gfs_grib2(localpet, npets)
  enddo
 
  deltalon = abs(longitude(2,1)-longitude(1,1))
+ if(localpet==0) print*, "deltalon = ", deltalon
 
  print*,"- CALL FieldScatter FOR INPUT GRID LONGITUDE."
  call ESMF_FieldScatter(longitude_input_grid, longitude, rootpet=0, rc=rc)
@@ -707,7 +708,9 @@ subroutine define_input_grid_gfs_grib2(localpet, npets)
      lat_src_ptr(i,j) = latitude(i,j)
    enddo
  enddo
-
+ if(localpet==0) print*, "lon first = ", lon_src_ptr(1:10,1)
+ if(localpet==0) print*, "lat first = ", lat_src_ptr(1,1:10)
+
  print*,"- CALL GridAddCoord FOR INPUT GRID."
  call ESMF_GridAddCoord(input_grid, &
                         staggerloc=ESMF_STAGGERLOC_CORNER, rc=rc)
@@ -782,7 +785,7 @@ subroutine define_input_grid_grib2(localpet, npets)
                                           lon_target(i_target,j_target)
                                           !lat_corners(ip1_input,jp1_input), &
                                           !lon_corners(ip1_input,jp1_input)
- real(esmf_kind_r8)                    :: deltalon, dx
+ real(esmf_kind_r8)                    :: deltalon, dx, deltalat
  integer                               :: ncid,id_var, id_dim 
  real(esmf_kind_r8), pointer           :: lat_src_ptr(:,:), lon_src_ptr(:,:)
  character(len=10000)            :: cmdline_msg, temp_msg, temp_msg2
@@ -1010,10 +1013,23 @@ subroutine define_input_grid_grib2(localpet, npets)
      if (trim(input_grid_type) == "latlon") then
 
       deltalon = abs(longitude_one_tile(2,1)-longitude_one_tile(1,1))
+      deltalat = abs(latitude_one_tile(1,2) - latitude_one_tile(1,1))
        do j = clb(2), cub(2)
          do i = clb(1), cub(1)
-           lon_src_ptr(i,j) = longitude_one_tile(i,1) - (0.5_esmf_kind_r8*deltalon)
-           lat_src_ptr(i,j) = 0.5_esmf_kind_r8 * (latitude_one_tile(i,j-1)+ latitude_one_tile(i,j))
+
+           if (i == ip1_input) then
+             lon_src_ptr(i,j) = longitude_one_tile(i-1,1)+(0.5_esmf_kind_r8*deltalon)
+           else 
+             lon_src_ptr(i,j) = longitude_one_tile(i,1) - (0.5_esmf_kind_r8*deltalon)
+           endif
+
+           if (j == jp1_input) then 
+             lat_src_ptr(i,j) = latitude_one_tile(1,j-1)+(0.5_esmf_kind_r8*deltalat)
+           else
+             lat_src_ptr(i,j) = latitude_one_tile(1,j) - (0.5_esmf_kind_r8*deltalat)
+           endif
+
+           !lat_src_ptr(i,j) = 0.5_esmf_kind_r8 * (latitude_one_tile(i,j-1)+ latitude_one_tile(i,j))
          enddo
        enddo
      else