Bug fix to LSM resolution

prep_regions/read_lsm.f90

@@ -89,8 +89,8 @@ subroutine read_lsm(files)
   endif
     
   ! indices to variable
-  ix = minloc(lon, dim=1, mask=abs(lon-dx/2.-rllx).eq.minval(abs(lon-dx/2.-rllx)))
-  jy = minloc(lat, dim=1, mask=abs(lat-dy/2.-rlly).eq.minval(abs(lat-dy/2.-rlly)))
+  ix = minloc(lon, dim=1, mask=abs(lon-rdx/2.-rllx).eq.minval(abs(lon-rdx/2.-rllx)))
+  jy = minloc(lat, dim=1, mask=abs(lat-rdy/2.-rlly).eq.minval(abs(lat-rdy/2.-rlly)))
 
   print*, 'read_lsm: lsm indices lon(ix), lat(jy) = ',lon(ix),lat(jy)
 
@@ -154,8 +154,8 @@ subroutine read_lsm(files)
     call check( nf90_get_var(ncid,latid,lat) )
 
     ! indices to variable
-    ix = minloc(lon, dim=1, mask=abs(lon-dx/2.-rllx).eq.minval(abs(lon-dx/2.-rllx)))
-    jy = minloc(lat, dim=1, mask=abs(lat-dy/2.-rlly).eq.minval(abs(lat-dy/2.-rlly)))
+    ix = minloc(lon, dim=1, mask=abs(lon-rdx/2.-rllx).eq.minval(abs(lon-rdx/2.-rllx)))
+    jy = minloc(lat, dim=1, mask=abs(lat-rdy/2.-rlly).eq.minval(abs(lat-rdy/2.-rlly)))
 
     print*, 'read_lsm: ice indices lon(ix), lat(jy) = ',lon(ix),lat(jy)