diff --git a/GFDL_tools/fv_diag_column.F90 b/GFDL_tools/fv_diag_column.F90
index 66b58a99cc..7460a4a9f7 100644
--- a/GFDL_tools/fv_diag_column.F90
+++ b/GFDL_tools/fv_diag_column.F90
@@ -478,7 +478,7 @@ subroutine debug_column_dyn(pt, delp, delz, u, v, w, q, heat_source, cappa, akap
                 else
                    heats = 0.0
                 endif
-                write(unit,'(I4, F7.2, F8.3, I6, F8.3, F8.3, F8.3, F8.3, F9.5, F9.3, F9.3, G )') &
+                write(unit,'(I4, F7.2, F8.3, I6, F8.3, F8.3, F8.3, F8.3, F9.5, F9.3, F9.3, G9.3 )') &
                      k, temp, delp(i,j,k)*0.01, -int(delz(i,j,k)), u(i,j,k), v(i,j,k), w(i,j,k), &
                      q(i,j,k,sphum)*1000., cond*1000., pres*1.e-2, (pres-preshyd(k))*1.e-2, heats
              enddo
diff --git a/driver/SHiELD/atmosphere.F90 b/driver/SHiELD/atmosphere.F90
index 31aac118da..a815f6c30e 100644
--- a/driver/SHiELD/atmosphere.F90
+++ b/driver/SHiELD/atmosphere.F90
@@ -758,14 +758,22 @@ subroutine atmosphere_hgt (hgt, position, relative, flip)
    !--- if needed, flip the indexing during this step
    if (flip) then
      if (.not. relative) then
-       z(:,:,1) = Atm(mygrid)%phis(:,:)/grav
+       do j = jsc, jec
+       do i = isc, iec
+         z(i-isc+1,j-jsc+1,1) = Atm(mygrid)%phis(i,j)/grav
+       enddo
+       enddo
      endif
      do k = 2,npz+1
        z(:,:,k) = z(:,:,k-1) - dz(:,:,npz+2-k)
      enddo
    else
      if (.not. relative) then
-       z(:,:,npz+1) = Atm(mygrid)%phis(:,:)/grav
+       do j = jsc, jec
+       do i = isc, iec
+         z(i-isc+1,j-jsc+1,npz+1) = Atm(mygrid)%phis(i,j)/grav
+       enddo
+       enddo
      endif
      do k = npz,1,-1
        z(:,:,k) = z(:,:,k+1) - dz(:,:,k)
@@ -875,7 +883,7 @@ subroutine atmosphere_nggps_diag (Time, init)
    logical, optional, intent(in) :: init
 
    if (PRESENT(init)) then
-     if (init == .true.) then
+     if (init .eqv. .true.) then
        call fv_nggps_diag_init(Atm(mygrid:mygrid), Atm(mygrid)%atmos_axes, Time)
        return
      else
diff --git a/tools/test_cases.F90 b/tools/test_cases.F90
index 9fecc2e1cc..55c8e9e877 100644
--- a/tools/test_cases.F90
+++ b/tools/test_cases.F90
@@ -6808,6 +6808,7 @@ subroutine read_namelist_test_case_nml(nml_filename)
 
         character(*), intent(IN) :: nml_filename
         integer :: ierr, f_unit, unit, ios
+        namelist /test_case_nml/test_case, bubble_do, alpha, nsolitons, soliton_Umax, soliton_size
 
 #include<file_version.h>
 
@@ -6817,7 +6818,6 @@ subroutine read_namelist_test_case_nml(nml_filename)
         alpha = 0.
         bubble_do = .false.
         test_case = 11   ! (USGS terrain)
-        namelist /test_case_nml/test_case, bubble_do, alpha, nsolitons, soliton_Umax, soliton_size
 
 #ifdef INTERNAL_FILE_NML
         ! Read Test_Case namelist