anchor.diff:  !!------------------------------------------------------------
   module geo
       use data_structures
  +    use io_routines, only:io_write
       implicit none
   
       private
  
positive.code:                         lon(k) = lo%lon(xy%x(k), xy%y(k))
                      enddo
                      ! and calculate the weights to apply to each gridcell
                      lo%geolut%w(:,i,j) = bilin_weights(hi%lat(i,j), lat, hi%lon(i,j), lon)
                  endif
              enddo
          enddo
  
      end subroutine geo_LUT
  
      !>------------------------------------------------------------
  
negative.code: 
              do i=1, nx
                  curpos = find_location(lo, hi%lat(i,j), hi%lon(i,j), lastpos)
  
                  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
  