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

[danielpeter at geodynamics.org]

Author: danielpeter
Date: 2009-09-27 11:02:04 -0700 (Sun, 27 Sep 2009)
New Revision: 15694

Modified:
   seismo/3D/ADJOINT_TOMO/measure_adj/mt_measure_adj.f90
   seismo/3D/ADJOINT_TOMO/measure_adj/mt_sub.f90
Log:
minor change for wsac0 routine call

Modified: seismo/3D/ADJOINT_TOMO/measure_adj/mt_measure_adj.f90
===================================================================
--- seismo/3D/ADJOINT_TOMO/measure_adj/mt_measure_adj.f90	2009-09-26 23:30:54 UTC (rev 15693)
+++ seismo/3D/ADJOINT_TOMO/measure_adj/mt_measure_adj.f90	2009-09-27 18:02:04 UTC (rev 15694)
@@ -14,7 +14,7 @@
   character(len=150) :: out_dir,datafile,synfile,file_prefix,file_prefix0,file_prefix2,measure_file_prefix,adj_file_prefix
 
   integer :: num_files, num_meas, i, j, k, iker, iker0, &
-     ios, is_mtm, is_mtm0, npt1, npt2, npts, nn, nerr
+     ios, is_mtm, is_mtm0, npt1, npt2, npts, nn
   double precision, dimension(NDIM) :: data, syn, adj_syn_all, tr_adj_src, am_adj_src, recon_cc_all, tseis_recon_cc
   double precision :: t01, dt1, t02, dt2, t0, dt, tstart, tend, chi, tt, dtt, df
   double precision, dimension(NCHI) :: window_chi
@@ -402,7 +402,7 @@
 
     ! write SAC file
     call dwsac1(trim(file_prefix2)//'.recon.cc.sac',recon_cc_all,npts,t0,dt)
-    if (nerr > 0) stop 'Error writing reconstructed CC file'
+    
     !call dwrite_sacfile_f(trim(datafile),trim(file_prefix2)//'.recon.cc',t0,npts,recon_cc_all)
 
   enddo ! npairs

Modified: seismo/3D/ADJOINT_TOMO/measure_adj/mt_sub.f90
===================================================================
--- seismo/3D/ADJOINT_TOMO/measure_adj/mt_sub.f90	2009-09-26 23:30:54 UTC (rev 15693)
+++ seismo/3D/ADJOINT_TOMO/measure_adj/mt_sub.f90	2009-09-27 18:02:04 UTC (rev 15694)
@@ -980,6 +980,8 @@
        call setfhv('depmin',xmin_sngl,nerr)
        call setfhv('depmax',xmax_sngl,nerr)
 
+       call setnhv('npts',npt1,nerr)          ! sets number of points
+
        ! write file with headers
        call wsac0(datafile,ti_sngl,dat_sngl,nerr)