[cig-commits] [commit] devel: put display of NSTEP back in src/shared/save_header_file.F90 (6faf52b)

[cig_noreply at geodynamics.org]

Repository : https://github.com/geodynamics/specfem3d_globe

On branch  : devel
Link       : https://github.com/geodynamics/specfem3d_globe/compare/9fae3b1bba1a593823a8a794065969a26d3d7c5e...3b4e39336ee6cdf596346281c3ae85e3019fef73

>---------------------------------------------------------------

commit 6faf52bd93bae378b023d8d906f088ca32c1b3a8
Author: Dimitri Komatitsch <komatitsch at lma.cnrs-mrs.fr>
Date:   Sun Jul 13 19:24:35 2014 +0200

    put display of NSTEP back in src/shared/save_header_file.F90


>---------------------------------------------------------------

6faf52bd93bae378b023d8d906f088ca32c1b3a8
 src/create_header_file/create_header_file.f90 | 5 ++++-
 src/shared/save_header_file.F90               | 7 ++++++-
 2 files changed, 10 insertions(+), 2 deletions(-)

diff --git a/src/create_header_file/create_header_file.f90 b/src/create_header_file/create_header_file.f90
index c65d3c0..7b0b149 100644
--- a/src/create_header_file/create_header_file.f90
+++ b/src/create_header_file/create_header_file.f90
@@ -95,7 +95,10 @@
   print *,'total elements per slice = ',sum(NSPEC)
   print *,'total points per slice = ',sum(nglob)
   print *
-  print *,'time-stepping of the solver will be: ',DT
+  print *,'the total number of time steps will be NSTEP = ',NSTEP
+  print *,'the time step of the solver will be DT = ',sngl(DT)
+  print *,'the total duration will thus be ',sngl(DT*(NSTEP-1)),' seconds'
+  print *,'                    i.e. ',sngl(DT*(NSTEP-1)/60.d0),' minutes'
   print *
   if(MOVIE_SURFACE .or. MOVIE_VOLUME) then
     print *,'MOVIE_VOLUME :',MOVIE_VOLUME
diff --git a/src/shared/save_header_file.F90 b/src/shared/save_header_file.F90
index 65247cc..7113173 100644
--- a/src/shared/save_header_file.F90
+++ b/src/shared/save_header_file.F90
@@ -71,7 +71,7 @@
     ANGULAR_WIDTH_XI_IN_DEGREES,ANGULAR_WIDTH_ETA_IN_DEGREES,NCHUNKS, &
     INCLUDE_CENTRAL_CUBE,CENTER_LONGITUDE_IN_DEGREES, &
     CENTER_LATITUDE_IN_DEGREES,GAMMA_ROTATION_AZIMUTH, &
-    NSTEP,NEX_XI,NEX_ETA, &
+    NSTEP,DT,NEX_XI,NEX_ETA, &
     NPROC_XI,NPROC_ETA, &
     SAVE_REGULAR_KL, &
     PARTIAL_PHYS_DISPERSION_ONLY, &
@@ -161,6 +161,11 @@
   write(IOUT,*) '! total elements per slice = ',sum(NSPEC)
   write(IOUT,*) '! total points per slice = ',sum(NGLOB)
   write(IOUT,*) '!'
+  write(IOUT,*) '! the total number of time steps will be NSTEP = ',NSTEP
+  write(IOUT,*) '! the time step of the solver will be DT = ',sngl(DT)
+  write(IOUT,*) '! the total duration will thus be ',sngl(DT*(NSTEP-1)),' seconds'
+  write(IOUT,*) '!        i.e. ',sngl(DT*(NSTEP-1)/60.d0),' minutes'
+  write(IOUT,*) '!'
 
   write(IOUT,'(1x,a,i1,a)') '! total for full ',NCHUNKS,'-chunk mesh:'
   write(IOUT,*) '! ---------------------------'