Added output of zisodisp array

src/write_global_output.f90

@@ -17,7 +17,7 @@
 !------------------------------------------------------------------------------|
 
 subroutine write_global_output (params,istep,iter,current_time,osolve,ov, &
-                                vo,surface,cl,outputtype)
+                                vo,surface,cl,zisodisp,outputtype)
 
 !------------------------------------------------------------------------------|
 !(((((((((((((((( Purpose of the routine  ))))))))))))))))))))))))))))))))))))))
@@ -55,6 +55,7 @@ type (octreev) ov
 type (void) vo
 type (sheet) surface(params%ns)
 type (cloud) cl
+double precision zisodisp(2**params%levelmax_oct+1,2**params%levelmax_oct+1)
 character*5 outputtype
 
 !------------------------------------------------------------------------------|
@@ -180,6 +181,12 @@ if (iproc.eq.0) then
                 i=1,cl%np)
    close (9)
 
+   if (params%isobc) then 
+     ! write isostasy basal displacement array - dwhipp 11/09 
+       write (9) ((zisodisp(i,j),j=1,2**params%levelmax_oct+1),i=1,& 
+                 2**params%levelmax_oct+1) 
+   end if 
+
 !   do i=1,osolve%nnode
 !      if (abs(osolve%x(i)-0.5d0)<1.d-6 .and. osolve%z(i)<0.08d0) then
 !         write (1234,'(5f30.15)') osolve%y(i),osolve%z(i),osolve%v(i),osolve%w(i),ov%temporary_nodal_pressure(i)