Fix petsc finaliser

Author: HarryMclean

src/heatflow.f90

@@ -65,6 +65,8 @@ program HEATFLOW_V0_3
    !-------------------------------------------------------------!
    ! run simulation for 'ntime' time steps                       !
    !-------------------------------------------------------------!
+   CALL petsc_init()                                 
+
    do itime=1,ntime 
 
       if (iverb.eq.0) then
@@ -81,17 +83,17 @@ program HEATFLOW_V0_3
       if (itime .eq. 1) CALL initial_evolve                      
 
       ! run the time evolution  
-      CALL petsc_init()                                 
       CALL simulate(itime)
-      CALL petsc_finalize()
 
 
 
       ! Write results                           
       CALL data_write(itime) 
       if (IVERB.ge.3) CALL final_print                           
 
-   end do                                                        
+   end do  
+   CALL petsc_finalize()
+                                                      
    !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^!
 
    !-------------------------------------------------------------!

src/heatflow/mod_evolve.f90

@@ -152,7 +152,6 @@ subroutine simulate(itime)
     iter= 0
     err=E
 
-   print *, "Calling solver..."
    !  call bicgstab(acsr, ia, ja, S, itmax, Temp_p, x, iter)
   if (.not. allocated(ia32)) then
      allocate(ia32(size(ia)), ja32(size(ja)))
@@ -164,7 +163,6 @@ subroutine simulate(itime)
    call solve_petsc_csr(NA32, ia32, ja32, acsr, S, x, tol, itmax)
   if (allocated(ia32)) deallocate(ia32, ja32)
 
-   print *, "Solver finished."
 
    ! CALL solve_pardiso(acsr, S, ia, ja, x)
    !  CALL linbcg(S,x,itol=int(itol,I4B),tol=tol, itmax=int(itmax,I4B), iter=iter, &