Call io_write_porotrack here.

MintonGroup · Dec 29, 2016 · c135716 · c135716
1 parent ab45ebe
commit c135716
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/src/io/io_write_surf.f90 b/src/io/io_write_surf.f90
@@ -65,6 +65,10 @@ subroutine io_write_surf(user,surf)
       !call io_write_regodist(user,surf)
       call io_write_regotrack(user,surf) 
    end if
+
+   ! write data from porosity analysis
+   if (user%doporosity) call io_write_porotrack(user, surf)
+
 !   if (user%docrustal_thinning) then
 !      recsize = sizeof(itmp) * user%gridsize * user%gridsize
 !      open(LUN,file=THICKFILE,status='replace',form='unformatted',recl=recsize,access='direct')