anchor.diff:          type(geo_look_up_table),intent(in)::geolut
           integer::i,k,l,ny,nx,localx,localy
           real::localw
  -        
  +
           nx=size(fieldout,1)
           ny=size(fieldout,2)
           ! use the geographic lookup table generated earlier to
  
positive.code:             enddo
          endif
      end subroutine geo_interp
      
      !>------------------------------------------------------------
      !!  Interpolate fieldout to fieldin using geolut. 
      !!
      !!------------------------------------------------------------
      subroutine geo_interp2d(fieldout, fieldin, geolut)
  
negative.code:                 lastpos%y = j
                  write(*,*) "Error in Geographic interpolation, check input lat / lon grids", 1,j
              endif
  
              do i=1, nx
                  curpos = find_location(lo, hi%lat(i,j), hi%lon(i,j), lastpos)
                  if (curpos%x<1) then
  