Updating move_surface to move_cloud method

Using the same algorithm as in move_cloud.f90

src/move_surface.f90

@@ -116,9 +116,9 @@ do i=iproc+1,surface%nsurface,nproc
    u=0.d0
    v=0.d0
    w=0.d0
-   x0=surface%x(i)
-   y0=surface%y(i)
-   z0=surface%z(i)
+   x0=surface%x(i)*rat
+   y0=surface%y(i)*rat
+   z0=surface%z(i)*rat
    x=x0
    y=y0
    z=z0
@@ -154,23 +154,23 @@ do i=iproc+1,surface%nsurface,nproc
 
       if (surface%closed.eq.0) then
         if (i==1 .or. i==nside .or. i==2*nside-1 .or. i==3*nside-2) then !corner points
-          z=z0+dt*w
+          if (surface%z(i)*(surface%z(i)-1.d0).lt.0.d0) z=z0+dt*w
         elseif (i>1 .and. i<nside) then                                  !y=0 edge
-          x=x0+dt*u
-          z=z0+dt*w
+          if (surface%x(i)*(surface%x(i)-1.d0).lt.0.d0) x=x0+dt*u
+          if (surface%z(i)*(surface%z(i)-1.d0).lt.0.d0) z=z0+dt*w
         elseif (i>nside .and. i<2*nside-1) then                          !x=1 edge 
-          y=y0+dt*v
-          z=z0+dt*w
+          if (surface%y(i)*(surface%y(i)-1.d0).lt.0.d0) y=y0+dt*v
+          if (surface%z(i)*(surface%z(i)-1.d0).lt.0.d0) z=z0+dt*w
         elseif (i>2*nside-1 .and. i<3*nside-2) then                      !y=1 edge
-          x=x0+dt*u
-          z=z0+dt*w
+          if (surface%x(i)*(surface%x(i)-1.d0).lt.0.d0) x=x0+dt*u
+          if (surface%z(i)*(surface%z(i)-1.d0).lt.0.d0) z=z0+dt*w
         elseif (i>3*nside-2 .and. i<4*nside-3) then                      !x=0 edge
-          y=y0+dt*v
-          z=z0+dt*w
+          if (surface%y(i)*(surface%y(i)-1.d0).lt.0.d0) y=y0+dt*v
+          if (surface%z(i)*(surface%z(i)-1.d0).lt.0.d0) z=z0+dt*w
         else 
-          x=x0+dt*u
-          y=y0+dt*v
-          z=z0+dt*w
+          if (surface%x(i)*(surface%x(i)-1.d0).lt.0.d0) x=x0+dt*u
+          if (surface%y(i)*(surface%y(i)-1.d0).lt.0.d0) y=y0+dt*v
+          if (surface%z(i)*(surface%z(i)-1.d0).lt.0.d0) z=z0+dt*w
         end if
 
       ! So I have added this section such that the algorithm still works for

src/remove_point.f90

@@ -67,8 +67,11 @@ call mpi_comm_rank (mpi_comm_world,iproc,ierr)
 !nside=2**surface%levelt+1
 nside=2**surface%levelt+2   !opla
 
-if (iproc.eq.0) write(8,*) '=================================================='
-if (iproc.eq.0) write(8,*) 'removing point iii=',iii
+if (iproc.eq.0) then
+    write(8,*) '=================================================='
+    write(8,*) 'removing point iii=',iii
+    write(8,'(a,3F8.6)') 'iii= position',surface%x(iii),surface%y(iii),surface%z(iii)
+end if
 
 !===============================================================================
 !=====[find how many triangles have point iii as a vertex]======================