modified to use new zi (zisodisp) structure

src/write_global_output.f90

@@ -17,7 +17,7 @@
 !------------------------------------------------------------------------------|
 
 subroutine write_global_output (params,istep,iter,current_time,osolve,ov, &
-                                vo,surface,cl,zisodisp,outputtype)
+                                vo,surface,cl,zi,outputtype)
 
 !------------------------------------------------------------------------------|
 !(((((((((((((((( Purpose of the routine  ))))))))))))))))))))))))))))))))))))))
@@ -55,7 +55,8 @@ 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)
+type (ziso) zi
+!double precision zisodisp(2**params%levelmax_oct+1,2**params%levelmax_oct+1)
 character*5 outputtype
 
 !------------------------------------------------------------------------------|
@@ -184,7 +185,7 @@ if (iproc.eq.0) then
 
    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,& 
+       write (9) ((zi%zisodisp(i,j),j=1,2**params%levelmax_oct+1),i=1,& 
                  2**params%levelmax_oct+1) 
    end if