[cig-commits] r13212 - seismo/2D/SPECFEM2D/trunk

[dkomati1 at geodynamics.org]

Author: dkomati1
Date: 2008-11-01 11:13:44 -0700 (Sat, 01 Nov 2008)
New Revision: 13212

Modified:
   seismo/2D/SPECFEM2D/trunk/specfem2D.F90
Log:
added printing of the total number of time steps and also of the percentage of time steps already computed


Modified: seismo/2D/SPECFEM2D/trunk/specfem2D.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/specfem2D.F90	2008-11-01 18:01:17 UTC (rev 13211)
+++ seismo/2D/SPECFEM2D/trunk/specfem2D.F90	2008-11-01 18:13:44 UTC (rev 13212)
@@ -2855,16 +2855,17 @@
          potential_dot_acoustic,potential_dot_dot_acoustic,TURN_ATTENUATION_ON,TURN_ANISOTROPY_ON,Mu_nu1,Mu_nu2,N_SLS)
 
 !----  display time step and max of norm of displacement
-  if(mod(it,NTSTEP_BETWEEN_OUTPUT_INFO) == 0 .or. it == 5) then
+  if(mod(it,NTSTEP_BETWEEN_OUTPUT_INFO) == 0 .or. it == 5 .or. it == NSTEP) then
 
     if (myrank == 0) then
-    write(IOUT,*)
-    if(time >= 1.d-3 .and. time < 1000.d0) then
-      write(IOUT,"('Time step number ',i7,'   t = ',f9.4,' s')") it,time
-    else
-      write(IOUT,"('Time step number ',i7,'   t = ',1pe12.6,' s')") it,time
+      write(IOUT,*)
+      if(time >= 1.d-3 .and. time < 1000.d0) then
+        write(IOUT,"('Time step number ',i7,'   t = ',f9.4,' s out of ',i7)") it,time,NSTEP
+      else
+        write(IOUT,"('Time step number ',i7,'   t = ',1pe12.6,' s out of ',i7)") it,time,NSTEP
+      endif
+      write(IOUT,*) 'We have done ',sngl(100.d0*dble(it-1)/dble(NSTEP-1)),'% of the total'
     endif
-    endif
 
     if(any_elastic_glob) then
       if(any_elastic) then