anchor.diff:                  stop
               endif
   
  -            do i=ims, ime
  +            do i = ims, ime
                   curpos = find_location(lo, hi%lat(i,j), hi%lon(i,j), lastpos)
   
  -                if (curpos%x<1) then
  +                if (curpos%x < 1) then
                       ! something broke, try assuming that the grids are the same possibly wrapping (for ideal)
                       write(*,*) "Error in Geographic interpolation, check input lat / lon grids", i,j
                       curpos%x = mod(i-1, lo_nx) + 1
  
positive.code:         find_surrounding%x = -999
          find_surrounding%y = -999
  
          do i=1,4
              dx = xdeltas(i)
              dy = ydeltas(i)
              search_point%x = [pos%x, pos%x+dx, pos%x+dx, pos%x  ]
              search_point%y = [pos%y, pos%y,   pos%y+dy, pos%y+dy]
  
  
negative.code:                     lo%geolut%w(:,i,j) = tri_weights(hi%lat(i,j), lat, hi%lon(i,j), lon)
                      ! lo%geolut%w(:,i,j) = idw_weights(hi%lat(i,j), lat, hi%lon(i,j), lon)
                      ! lo%geolut%w(:,i,j) = bilin_weights(hi%lat(i,j), lat, hi%lon(i,j), lon)
  
                  endif
  
              enddo
  