[cig-commits] r21094 - seismo/3D/ADJOINT_TOMO/measure_adj
liuqy at geodynamics.org
liuqy at geodynamics.org
Fri Nov 30 14:51:01 PST 2012
Author: liuqy
Date: 2012-11-30 14:51:01 -0800 (Fri, 30 Nov 2012)
New Revision: 21094
Modified:
seismo/3D/ADJOINT_TOMO/measure_adj/ma_sub.f90
Log:
Now istart is indeed the first index of selected window [tstart, tend]
Modified: seismo/3D/ADJOINT_TOMO/measure_adj/ma_sub.f90
===================================================================
--- seismo/3D/ADJOINT_TOMO/measure_adj/ma_sub.f90 2012-11-30 22:09:56 UTC (rev 21093)
+++ seismo/3D/ADJOINT_TOMO/measure_adj/ma_sub.f90 2012-11-30 22:51:01 UTC (rev 21094)
@@ -800,7 +800,7 @@
! compute traveltime and amplitude adjoint sources for imeas
do i = 1,nlen
- i1 = istart + i ! start index in the full adjoint source array(1:npts)
+ i1 = istart + i -1 ! start index in the full adjoint source array(1:npts)
! waveform
if(imeas==1 .or. imeas==2) then
@@ -1010,9 +1010,9 @@
!call setlhv('leven',.false.,nerr) ! sets un-even sampling
!call setihv('iftype','itime',nerr) ! sets file type: time file
- ! write file with headers
- call wsac0(datafile,ti_sngl,dat_sngl,nerr)
-
+ ! write file with headers (LQY: bug with b in wsac0())
+ ! call wsac0(datafile,ti_sngl,dat_sngl,nerr)
+ call wsac1(datafile,dat_sngl,npt1,b1_sngl,dt1_sngl,nerr)
else
call wsac1(datafile,dat_sngl,npt1,b1_sngl,dt1_sngl,nerr)
endif
@@ -1213,7 +1213,8 @@
implicit none
double precision, dimension(NPT), intent(in) :: data, syn
- double precision, intent(in) :: tstart, tend, t0, dt
+ double precision, intent(in) :: tstart
+ double precision, intent(in) :: tend, t0, dt
integer, intent(in) :: NPT
double precision, dimension(NPT), intent(out) :: dat_win, syn_win
@@ -1223,7 +1224,10 @@
double precision :: time, t1
nlen = floor((tend-tstart)/dt) + 1
- istart = floor((tstart-t0)/dt)
+
+ istart = floor((tstart-t0)/dt) + 1
+ ! tstart = t0+(istart-1)*dt ! minor adjustments
+ !print *, '*** diff tstart = ', t0+(istart-1)*dt - tstart
! limits array bounds
if( nlen > NPT ) nlen = NPT
More information about the CIG-COMMITS
mailing list