[cig-commits] r21094 - seismo/3D/ADJOINT_TOMO/measure_adj

[liuqy at geodynamics.org]

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