[cig-commits] [commit] master: deltat is needed earlier (c7f8675)

[cig_noreply at geodynamics.org]

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

On branch  : master
Link       : https://github.com/geodynamics/axisem/compare/088416bb7afde58cdb85f91f5a5589149afc1141...c7f867511ca351fbf339112efbdc0143d1f725e8

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

commit c7f867511ca351fbf339112efbdc0143d1f725e8
Author: martinvandriel <vandriel at erdw.ethz.ch>
Date:   Mon Dec 29 17:06:03 2014 +0100

    deltat is needed earlier


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

c7f867511ca351fbf339112efbdc0143d1f725e8
 SOLVER/parameters.F90     | 1 -
 SOLVER/time_evol_wave.F90 | 6 +++---
 2 files changed, 3 insertions(+), 4 deletions(-)

diff --git a/SOLVER/parameters.F90 b/SOLVER/parameters.F90
index e4af391..4a81189 100644
--- a/SOLVER/parameters.F90
+++ b/SOLVER/parameters.F90
@@ -779,7 +779,6 @@ subroutine compute_numerical_parameters
   use attenuation, only: dump_memory_vars
   use data_mesh
 
-
   real(kind=dp)         :: s,z,r,theta,s_max,dshift
   real(kind=dp)         :: dsaxis(0:npol-1,0:npol), dzaxis(0:npol-1) 
   real(kind=dp)         :: minds(nelem),maxds(nelem),mindz(nelem),maxdz(nelem)
diff --git a/SOLVER/time_evol_wave.F90 b/SOLVER/time_evol_wave.F90
index 67ab643..9a347d4 100644
--- a/SOLVER/time_evol_wave.F90
+++ b/SOLVER/time_evol_wave.F90
@@ -79,14 +79,14 @@ subroutine prepare_waves
                         .or. trim(dump_type) == 'strain_only')) &
      call build_kwf_grid()
 
+  ! compute/output some more parameters
+  call compute_numerical_parameters
+
   ! Define velocity/density model (velocities in m/s, density in kg/m^3 ) AND 
   ! compute all global matrices (Jacobian, mapping, mass matrix, S/F boundary)
   ! for solid and fluid domains respectively
   call read_model_compute_terms
 
-  ! compute/output some more parameters 
-  call compute_numerical_parameters
-
   ! compute source time function
   call compute_stf