[cig-commits] r12381 - seismo/3D/ADJOINT_TOMO/flexwin

[alessia at geodynamics.org]

Author: alessia
Date: 2008-07-07 02:42:14 -0700 (Mon, 07 Jul 2008)
New Revision: 12381

Added:
   seismo/3D/ADJOINT_TOMO/flexwin/flexwin.f90
Removed:
   seismo/3D/ADJOINT_TOMO/flexwin/test2.f90
Modified:
   seismo/3D/ADJOINT_TOMO/flexwin/make_gfortran
Log:
renamed test2 to flexwin and fixed make_gfortran makefile

Copied: seismo/3D/ADJOINT_TOMO/flexwin/flexwin.f90 (from rev 12377, seismo/3D/ADJOINT_TOMO/flexwin/test2.f90)
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/flexwin.f90	                        (rev 0)
+++ seismo/3D/ADJOINT_TOMO/flexwin/flexwin.f90	2008-07-07 09:42:14 UTC (rev 12381)
@@ -0,0 +1,48 @@
+!----------------------------------------------------
+program flexwin
+
+use seismo_variables
+
+character*120, dimension (1000):: basename
+character*240, dimension (1000):: obs_name, syn_name
+integer :: n_seis
+
+! read the parameter file
+call read_parameter_file()
+
+write(*,*) 'Number of seismograms to measure : '
+read(*,*) n_seis
+write(*,*) n_seis
+
+write(*,*) 'For each seismogram: '
+do i = 1, n_seis
+  write(*,*) 'Observed seismogram'
+  read(*,'(a)') obs_name(i)
+  write(*,'(a)') trim(obs_name(i))
+  write(*,*) 'Synthetic seismogram'
+  read(*,'(a)') syn_name(i)
+  write(*,'(a)') trim(syn_name(i))
+  write(*,*) 'Output basename'
+  read(*,'(a)') basename(i)
+  write(*,'(a)') trim(basename(i))
+
+  if (DEBUG) write(*,*) 'DEBUG : reading sac files'
+  call read_sac_files(syn_name(i),obs_name(i))
+
+  if (DEBUG) write(*,*) 'DEBUG : selecting windows'
+  call select_windows_stalta2()
+
+
+  if(MAKE_SEISMO_PLOTS) then
+     if (DEBUG) write(*,*) 'DEBUG : writing output seismos'
+     call write_seismos_gmt(basename(i))
+  endif
+
+  if(MAKE_WINDOW_FILES) then
+     if (DEBUG) write(*,*) 'DEBUG : writing mt input'
+     call write_mt_input_2(basename(i),obs_name(i),syn_name(i))
+  endif
+
+enddo
+
+end program flexwin

Modified: seismo/3D/ADJOINT_TOMO/flexwin/make_gfortran
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/make_gfortran	2008-07-07 09:35:58 UTC (rev 12380)
+++ seismo/3D/ADJOINT_TOMO/flexwin/make_gfortran	2008-07-07 09:42:14 UTC (rev 12381)
@@ -8,15 +8,15 @@
 A_CODE=/home/alessia/code
 
 SACLIBDIR=/opt/sac/lib
-QDSACLIBDIR=qd_sac
+QDSACLIBDIR=../SacLib
 TAULIBDIR=ttimes_mod
 #LIBDIR=${A_CODE}/lib
 #INCDIR=${A_CODE}/include
 #SHARE=${A_CODE}/share
 #LIBS= -lSacLib -lsac -ltau -lm 
-LIBS= -lsacio -lqd_sac -ltau -lm 
+LIBS= -lsacio -lSacLib -ltau -lm 
 
-PROGS= test2
+PROGS= flexwin
 default:${PROGS} 
 
 LIBSRC2=distaz.f user_parameters.f90  seismo_subs.f90  user_functions.f90 maxima.f90 travel_times.f90 measurement_module.f90 xcorr-measure.f90 measure_windows_xcorr.f90 io_subs.f90 select_windows_stalta2.f90
@@ -39,8 +39,8 @@
 	$(FC) ${FFLAGS} -c -o $@ $< 
 
 
-test2:test2.f90 ${LIB2}
-	$(FC) ${FFLAGS} test2.f90 -o test2 ${LIB2} -Lttimes_mod -Lqd_sac -L${SACLIBDIR} ${LIBS} 
+flexwin:flexwin.f90 ${LIB2}
+	$(FC) ${FFLAGS} flexwin.f90 -o flexwin ${LIB2} -Lttimes_mod -L${QDSACLIBDIR} -L${SACLIBDIR} ${LIBS} 
 
 clean:
 	rm -f  ${LIB2} ${PROGS} *.o *.mod

Deleted: seismo/3D/ADJOINT_TOMO/flexwin/test2.f90
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/test2.f90	2008-07-07 09:35:58 UTC (rev 12380)
+++ seismo/3D/ADJOINT_TOMO/flexwin/test2.f90	2008-07-07 09:42:14 UTC (rev 12381)
@@ -1,57 +0,0 @@
-!
-! $Id:$
-!
-
-!----------------------------------------------------
-program test2
-
-use seismo_variables
-
-character*120, dimension (1000):: basename
-character*240, dimension (1000):: obs_name, syn_name
-integer :: n_seis
-
-! read the parameter file
-call read_parameter_file()
-
-write(*,*) 'Number of seismograms to measure : '
-read(*,*) n_seis
-write(*,*) n_seis
-
-write(*,*) 'For each seismogram: '
-do i = 1, n_seis
-  write(*,*) 'Observed seismogram'
-  read(*,'(a)') obs_name(i)
-  write(*,'(a)') trim(obs_name(i))
-  write(*,*) 'Synthetic seismogram'
-  read(*,'(a)') syn_name(i)
-  write(*,'(a)') trim(syn_name(i))
-  write(*,*) 'Output basename'
-  read(*,'(a)') basename(i)
-  write(*,'(a)') trim(basename(i))
-
-  if (DEBUG) write(*,*) 'DEBUG : reading sac files'
-  call read_sac_files(syn_name(i),obs_name(i))
-
-  if (DEBUG) write(*,*) 'DEBUG : selecting windows'
-  call select_windows_stalta2()
-
-  !if (DEBUG) write(*,*) 'DEBUG : making xcorr measurements'
-  !call measure_windows_xcorr
-
-  if(MAKE_SEISMO_PLOTS) then
-     if (DEBUG) write(*,*) 'DEBUG : writing output seismos'
-     call write_seismos_gmt(basename(i))
-     !if (DEBUG) write(*,*) 'DEBUG : writing output measurements'
-     !call write_measurements_gmt(basename(i))
-  endif
-
-  if(MAKE_WINDOW_FILES) then
-     if (DEBUG) write(*,*) 'DEBUG : writing mt input'
-     !call write_mt_input(basename(i),obs_name(i),syn_name(i))
-     call write_mt_input_2(basename(i),obs_name(i),syn_name(i))
-  endif
-
-enddo
-
-end program test2