Fixed mistake in passing of ov arrays from r580

src/move_surface.f90

@@ -129,12 +129,15 @@ do i=iproc+1,surface%nsurface,nproc
       zi=(z0+z)/2.d0
       if (flag1==1) then
         ! Only vertical velocity from isostasy
-        call octree_interpolate_three (3,ov%octree,ov%noctree,ov%icon,ov%nleaves, &
-                                      ov%nnode,xi,yi,zi,0.d0,u,0.d0,v,ov%wnodeiso,w)
+        call octree_interpolate_three (3,ov%octree,ov%noctree,ov%icon,         &
+                                      ov%nleaves,ov%nnode,xi,yi,zi,            &
+                                      0.d0*ov%unode,u,0.d0*ov%vnode,v,         &
+                                      ov%wnodeiso,w)
       else
         ! Full velocity field
-        call octree_interpolate_three (3,ov%octree,ov%noctree,ov%icon,ov%nleaves, &
-                                      ov%nnode,xi,yi,zi,ov%unode,u,ov%vnode,v,ov%wnode/params%vex,w)
+        call octree_interpolate_three (3,ov%octree,ov%noctree,ov%icon,         &
+                                      ov%nleaves,ov%nnode,xi,yi,zi,ov%unode,u, &
+                                      ov%vnode,v,ov%wnode/params%vex,w)
       endif
 
       ! This section is by Cedric and is very restrictive in the sense that it
@@ -196,9 +199,9 @@ do i=iproc+1,surface%nsurface,nproc
         ! Only vertical velocity from isostasy
         call octree_interpolate_three_derivative (3,ov%octree,ov%noctree,      &
                                                  ov%icon,ov%nleaves,ov%nnode,  &
-                                                 xi,yi,zi,0.d0,u,ux,uy,uz,0.d0,&
-                                                 v,vx,vy,vz,ov%wnodeiso,w,wx,  &
-                                                 wy,wz)
+                                                 xi,yi,zi,0.d0*ov%unode,u,ux,  &
+                                                 uy,uz,0.d0*ov%vnode,v,vx,vy,  &
+                                                 vz,ov%wnodeiso,w,wx,wy,wz)
       else
         ! Full velocity field
         call octree_interpolate_three_derivative (3,ov%octree,ov%noctree,      &
@@ -207,7 +210,6 @@ do i=iproc+1,surface%nsurface,nproc
                                                  ov%vnode,v,vx,vy,vz, &
                                                  ov%wnode/params%vex,w,wx,wy,wz)
       endif
-
       xn0=surface%xn(i)
       yn0=surface%yn(i)
       zn0=surface%zn(i)