Skip to content

Commit

Permalink
Dump STF at time 0
Browse files Browse the repository at this point in the history
 - STF dump previously started after first timestep, which led to first sample (at t=0) missing.
 - Also removes the zero-valued sample at the end of the STF
  • Loading branch information
sstaehler committed Aug 27, 2014
1 parent 42a3e78 commit a068ba8
Showing 1 changed file with 2 additions and 2 deletions.
4 changes: 2 additions & 2 deletions SOLVER/nc_routines.F90
Original file line number Diff line number Diff line change
Expand Up @@ -534,17 +534,17 @@ subroutine nc_dump_stf(stf)
do i = 1, niter
! Dumping the STF in the fine time stepping of the seismogram output
if ( mod(i,seis_it) == 0) then
it_s = it_s + 1
stf_seis_dumpvar(it_s) = stf(i)
stf_d_seis_dumpvar(it_s) = stf_d_dumpvar(i)
it_s = it_s + 1
end if

if (dump_wavefields) then
! Dumping the STF in the coarse time stepping of the strain (KERNER) output
if ( mod(i,strain_it) == 0) then
it_d = it_d + 1
stf_dump_dumpvar(it_d) = stf(i)
stf_d_dump_dumpvar(it_d) = stf_d_dumpvar(i)
it_d = it_d + 1
end if
end if
end do
Expand Down

0 comments on commit a068ba8

Please sign in to comment.