diff --git a/src/initialise_mod.f90 b/src/initialise_mod.f90
index d497a1930b1ab9917f10db330908437de48a4517..798f2acf07830c72462c25207347950759eb19c1 100644
--- a/src/initialise_mod.f90
+++ b/src/initialise_mod.f90
@@ -933,9 +933,9 @@ subroutine init_domainfill
   call alloc_domainfill
 
   nx_we(1)=max(int(xpoint1(1)),0)
-  nx_we(2)=min((int(xpoint2(1))+1),nxmin1)
+  nx_we(2)=min((ceiling(xpoint2(1))),nxmin1)
   ny_sn(1)=max(int(ypoint1(1)),0)
-  ny_sn(2)=min((int(ypoint2(1))+1),nymin1)
+  ny_sn(2)=min((ceiling(ypoint2(1))),nymin1)
 
   ! For global simulations (both global wind data and global domain-filling),
   ! set a switch, such that no boundary conditions are used