anchor.diff:         lat = GMTR_lat(n,l)
          lon = GMTR_lon(n,l)
   
  -       z(ADM_kmin-1) = GRD_vz(n,2,l,GRD_ZH)
  -       do k = ADM_kmin, ADM_kmax+1
  +!       z(ADM_kmin-1) = GRD_vz(n,2,l,GRD_ZH)
  +!       do k = ADM_kmin, ADM_kmax+1
  +!          z(k) = GRD_vz(n,k,l,GRD_Z)
  +!       enddo
  +
  +       do k = 1, kdim
             z(k) = GRD_vz(n,k,l,GRD_Z)
  -       enddo
  -
  -       do k = 1, kdim
             call supercell_test( lon,       &  ! [IN ]
                                  lat,       &  ! [IN ]
                                  p,         &  ! [INOUT]
  
positive.code:        lat = GMTR_lat(n,l)
         lon = GMTR_lon(n,l)
  
         z(ADM_kmin-1) = GRD_vz(n,2,l,GRD_ZH)
         do k = ADM_kmin, ADM_kmax+1
            z(k) = GRD_vz(n,k,l,GRD_Z)
         enddo
  
         do k = 1, kdim
            call baroclinic_wave_test( deep,       &  ! [IN ]
                                       moist,      &  ! [IN ]
                                       pertt,      &  ! [IN ]
  
negative.code:        lat = GMTR_lat(n,l)
         lon = GMTR_lon(n,l)
  
  !       z(ADM_kmin-1) = GRD_vz(n,2,l,GRD_ZH)
  !       do k = ADM_kmin, ADM_kmax+1
  !          z(k) = GRD_vz(n,k,l,GRD_Z)
  !       enddo
  
         do k = 1, kdim
            z(k) = GRD_vz(n,k,l,GRD_Z)
            call tropical_cyclone_test( lon,       &  ! [IN ]
                                        lat,       &  ! [IN ]
                                        p,         &  ! [INOUT]
  