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:         integer,intent(in) :: nx,ny
          integer :: i, j
          integer :: xdeltas(4), ydeltas(4), dx, dy
          real :: best_err, current_err, tri1_err, tri2_err
          logical :: tri_1, tri_2
          type(fourpos) :: search_point
  
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
  