anchor.diff:          endif
           if (dy == 0) then
               if (.not.lo%dy_errors_printed) then
  -                print*, "ERROR : geo_find_location : DY = 0 !  Check your inputfile Latitude data"
  -                print*, "  Attempting to continue by assuming the grids are the same or can wrap arround"
  +                write(*,*) "ERROR : geo_find_location : DY = 0 !  Check your inputfile Latitude data"
  +                write(*,*) "  Attempting to continue by assuming the grids are the same or can wrap arround"
                   lo%dy_errors_printed=.True.
               endif
               find_location%x = -1
  
positive.code:                     endif
                  endif
              endif
              x0 = x1
              y0 = y1
          enddo
  
negative.code:         integer::nx,nz,ny
          integer:: i,j,k,l,localx,localy
          real::localw
          real :: local_center
  
          nx=size(fieldout,1)
          nz=size(fieldout,2)
  