[cig-commits] r16955 - seismo/3D/SPECFEM3D_GLOBE/branches/NOISE_TOMOGRAPHY

[yangl at geodynamics.org]

Author: yangl
Date: 2010-06-09 07:06:41 -0700 (Wed, 09 Jun 2010)
New Revision: 16955

Modified:
   seismo/3D/SPECFEM3D_GLOBE/branches/NOISE_TOMOGRAPHY/noise_tomography.f90
   seismo/3D/SPECFEM3D_GLOBE/branches/NOISE_TOMOGRAPHY/specfem3D.f90
Log:
minor changes on NOISE_TOMOGRAPHY

Modified: seismo/3D/SPECFEM3D_GLOBE/branches/NOISE_TOMOGRAPHY/noise_tomography.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/NOISE_TOMOGRAPHY/noise_tomography.f90	2010-06-09 09:42:11 UTC (rev 16954)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/NOISE_TOMOGRAPHY/noise_tomography.f90	2010-06-09 14:06:41 UTC (rev 16955)
@@ -215,7 +215,8 @@
 ! check for consistency of the parameters
   subroutine check_parameters_noise(myrank,NOISE_TOMOGRAPHY,SIMULATION_TYPE,SAVE_FORWARD, &
                                     NUMBER_OF_RUNS, NUMBER_OF_THIS_RUN,ROTATE_SEISMOGRAMS_RT, &
-                                    SAVE_ALL_SEISMOS_IN_ONE_FILE, USE_BINARY_FOR_LARGE_FILE)
+                                    SAVE_ALL_SEISMOS_IN_ONE_FILE, USE_BINARY_FOR_LARGE_FILE, &
+                                    MOVIE_COARSE)
   implicit none
   include 'mpif.h'
   include "precision.h"
@@ -223,7 +224,7 @@
   include "OUTPUT_FILES/values_from_mesher.h"
   ! input parameters
   integer :: myrank,NOISE_TOMOGRAPHY,SIMULATION_TYPE,NUMBER_OF_RUNS,NUMBER_OF_THIS_RUN
-  logical :: SAVE_FORWARD,ROTATE_SEISMOGRAMS_RT,SAVE_ALL_SEISMOS_IN_ONE_FILE, USE_BINARY_FOR_LARGE_FILE
+  logical :: SAVE_FORWARD,ROTATE_SEISMOGRAMS_RT,SAVE_ALL_SEISMOS_IN_ONE_FILE, USE_BINARY_FOR_LARGE_FILE,MOVIE_COARSE
   ! output parameters
   ! local parameters
 
@@ -254,7 +255,10 @@
      call exit_mpi(myrank,'Do NOT rotate seismograms in the code, change ROTATE_SEISMOGRAMS_RT in DATA/Par_file')
   if (SAVE_ALL_SEISMOS_IN_ONE_FILE .OR. USE_BINARY_FOR_LARGE_FILE) &
      call exit_mpi(myrank,'Please set SAVE_ALL_SEISMOS_IN_ONE_FILE and USE_BINARY_FOR_LARGE_FILE to be .false.')
+  if (MOVIE_COARSE) &
+     call exit_mpi(myrank,'Please set MOVIE_COARSE in DATA/Par_file to be .false.') 
 
+
   if (NOISE_TOMOGRAPHY==1) then
      if (SIMULATION_TYPE/=1) &
         call exit_mpi(myrank,'NOISE_TOMOGRAPHY=1 requires SIMULATION_TYPE=1! check DATA/Par_file')
@@ -376,14 +380,14 @@
   include "constants.h"
   include "OUTPUT_FILES/values_from_mesher.h"
   ! input parameters
-  integer :: myrank,nrec,NSTEP
+  integer :: myrank,nrec,NSTEP, irec_master_noise
   integer, dimension(nrec) :: islice_selected_rec,ispec_selected_rec
   integer, dimension(NGLLX,NGLLY,NGLLZ,NSPEC_CRUST_MANTLE) :: ibool_crust_mantle
   real(kind=CUSTOM_REAL), dimension(NDIM,NGLLX,NGLLY,NGLLZ,NSTEP) :: noise_sourcearray
   real(kind=CUSTOM_REAL),dimension(NDIM,NGLOB_CRUST_MANTLE) :: accel_crust_mantle  ! both input and output
   ! output parameters
   ! local parameters
-  integer :: irec_master_noise,i,j,k,iglob, it
+  integer :: i,j,k,iglob,it
 
 
   ! adds noise source (only if this proc carries the noise)

Modified: seismo/3D/SPECFEM3D_GLOBE/branches/NOISE_TOMOGRAPHY/specfem3D.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/NOISE_TOMOGRAPHY/specfem3D.f90	2010-06-09 09:42:11 UTC (rev 16954)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/NOISE_TOMOGRAPHY/specfem3D.f90	2010-06-09 14:06:41 UTC (rev 16955)
@@ -1600,7 +1600,8 @@
        if (myrank == 0) &
        call check_parameters_noise(myrank,NOISE_TOMOGRAPHY,SIMULATION_TYPE,SAVE_FORWARD, &
                                   NUMBER_OF_RUNS, NUMBER_OF_THIS_RUN,ROTATE_SEISMOGRAMS_RT, &
-                                  SAVE_ALL_SEISMOS_IN_ONE_FILE, USE_BINARY_FOR_LARGE_FILE)
+                                  SAVE_ALL_SEISMOS_IN_ONE_FILE, USE_BINARY_FOR_LARGE_FILE, &
+                                  MOVIE_COARSE)
     endif
 !>YANGL
 
@@ -2426,7 +2427,7 @@
                                 it,irec_master_noise)
     elseif ( NOISE_TOMOGRAPHY == 2 ) then
        ! second step of noise tomography, i.e., read the surface movie saved at every timestep
-       ! use the movie to drive the unsemble forward wavefield
+       ! use the movie to drive the ensemble forward wavefield
        call noise_read_add_surface_movie(myrank,nmovie_points,accel_crust_mantle, &
                               normal_x_noise,normal_y_noise,normal_z_noise,mask_noise, &
                               store_val_ux,store_val_uy,store_val_uz, &