[cig-commits] r18161 - in seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples: . test_global test_regional
yangl at geodynamics.org
yangl at geodynamics.org
Fri Apr 1 08:57:38 PDT 2011
Author: yangl
Date: 2011-04-01 08:57:37 -0700 (Fri, 01 Apr 2011)
New Revision: 18161
Added:
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/CMTSOLUTION_NOISE
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/NOISE123.submit_atten
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/NOISE_clean_create
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/NOISE_collect
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/Par_file_NOISE_1_attenuation
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/Par_file_NOISE_1_noattenuation
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/Par_file_NOISE_2_attenuation
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/Par_file_NOISE_2_noattenuation
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/Par_file_NOISE_3_noattenuation
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/STATIONS_NOISE
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/STATIONS_NOISE_ADJOINT
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/S_squared
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/adj_traveltime_filter.f90
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/irec_master_noise
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/nu_master
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/CMTSOLUTION_NOISE
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/NOISE123.submit_atten
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/NOISE_clean_create
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/NOISE_collect
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_1_attenuation
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_1_noattenuation
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_2_attenuation
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_2_noattenuation
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_3_noattenuation
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/STATIONS_NOISE
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/STATIONS_NOISE_ADJOINT
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/S_squared
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/adj_traveltime_filter.f90
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/irec_master_noise
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/nu_master
Log:
two very short tests for NOISE_TOMOGRAPHY, just to make sure scripts work before submitting the real jobs
Added: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/CMTSOLUTION_NOISE
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/CMTSOLUTION_NOISE (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/CMTSOLUTION_NOISE 2011-04-01 15:57:37 UTC (rev 18161)
@@ -0,0 +1,13 @@
+PDE 1994 6 9 0 33 16.40 -13.8300 -67.5600 637.0 6.9 6.8 NORTHERN BOLIVIA
+event name: 060994A
+time shift: 29.0000
+half duration: 30.0000
+latitude: 0.0000
+longitude: 0.0000
+depth: 0.0000
+Mrr: 0.00e+27
+Mtt: 0.00e+27
+Mpp: 0.00e+27
+Mrt: 0.00e+28
+Mrp: 0.00e+26
+Mtp: 0.00e+27
Added: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/NOISE123.submit_atten
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/NOISE123.submit_atten (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/NOISE123.submit_atten 2011-04-01 15:57:37 UTC (rev 18161)
@@ -0,0 +1,136 @@
+#!/bin/bash
+#PBS -N NOISE_test_global
+#PBS -q tromp
+#PBS -o zzz_test_global/job_info/output
+#PBS -e zzz_test_global/job_info/error
+#PBS -l nodes=75:ppn=8,walltime=100:00:00
+
+DIR_RESULT="zzz_test_global"
+NSTEP=799
+
+cd $PBS_O_WORKDIR
+######################################## pre-simulation ###########################################
+## gather simulation infomation
+NPROC_XI=`grep NPROC_XI DATA/Par_file | cut -d = -f 2 `
+NPROC_ETA=`grep NPROC_ETA DATA/Par_file | cut -d = -f 2`
+NCHUNKS=`grep NCHUNKS DATA/Par_file | cut -d = -f 2 `
+numnodes=$(( $NCHUNKS * $NPROC_XI * $NPROC_ETA ))
+cat $PBS_NODEFILE > $DIR_RESULT/job_info/nodes
+echo $PBS_JOBID > $DIR_RESULT/job_info/id
+
+## check number of procs
+NPROC_1=`wc -l < $DIR_RESULT/job_info/nodes` # number of procs assigned to the job
+NPROC_2=$numnodes # number of procs needed in SPECFEM
+if [ $NPROC_1 -ne $NPROC_2 ]; then
+ echo "inconsistent number of procs, check your job submission script & DATA/Par_file" 1>&2
+ echo "NPROC from job submission: $NPROC_1" 1>&2
+ echo "NPROC from DATA/Par_file: $NPROC_2" 1>&2
+ exit
+fi
+
+## clean local nodes
+LOCAL_PATH=`grep LOCAL_PATH DATA/Par_file | cut -d = -f 2 | sed 's/ //g'`
+d=`date`
+echo "Start cleaning local disk on nodes, $d"
+pbsdsh $PBS_O_WORKDIR/NOISE_clean_create $LOCAL_PATH
+d=`date`
+echo "Finish cleaning local disk on nodes, $d"
+
+## clean OUTPUT_FILES
+rm -rf OUTPUT_FILES
+mkdir -p OUTPUT_FILES
+######################################## simulation ###############################################
+################################################################### step 1: ensemble forward source
+cp DATA/Par_file_NOISE_1_attenuation DATA/Par_file
+mpiexec $PBS_O_WORKDIR/xmeshfem3D_attenuation
+mpiexec $PBS_O_WORKDIR/xspecfem3D_attenuation
+echo "2" > temp_input_movie; echo "1" >> temp_input_movie; echo "$NSTEP" >> temp_input_movie
+$PBS_O_WORKDIR/xcreate_movie_AVS_DX_attenuation < temp_input_movie
+# backup step 1 results (with attenuation)
+DIR_TEMP=$DIR_RESULT/step_1
+rm -rf $DIR_TEMP
+mkdir -p $DIR_TEMP
+mv OUTPUT_FILES/* $DIR_TEMP/
+cp DATA/Par_file $DIR_TEMP/
+cp DATA/CMTSOLUTION $DIR_TEMP/
+cp DATA/STATIONS $DIR_TEMP/
+################################################################ step 2: ensemble forward wavefield
+cp DATA/Par_file_NOISE_2_attenuation DATA/Par_file
+mpiexec $PBS_O_WORKDIR/xmeshfem3D_attenuation
+mpiexec $PBS_O_WORKDIR/xspecfem3D_attenuation
+echo "2" > temp_input_movie; echo "1" >> temp_input_movie; echo "$NSTEP" >> temp_input_movie
+$PBS_O_WORKDIR/xcreate_movie_AVS_DX_attenuation < temp_input_movie
+# backup step 2 results (with attenuation)
+DIR_TEMP=$DIR_RESULT/step_2
+rm -rf $DIR_TEMP
+mkdir -p $DIR_TEMP
+mv OUTPUT_FILES/* $DIR_TEMP/
+cp DATA/Par_file $DIR_TEMP/
+cp DATA/CMTSOLUTION $DIR_TEMP/
+cp DATA/STATIONS $DIR_TEMP/
+############################################################## prepare adjoint sources ############
+cp $DIR_RESULT/step_2/A7.II.MXZ.sem.ascii ./SEM/
+./NOISE_adj
+################################################################### step 1: ensemble forward source
+cp DATA/Par_file_NOISE_1_noattenuation DATA/Par_file
+mpiexec $PBS_O_WORKDIR/xmeshfem3D_noattenuation
+mpiexec $PBS_O_WORKDIR/xspecfem3D_noattenuation
+echo "2" > temp_input_movie; echo "1" >> temp_input_movie; echo "$NSTEP" >> temp_input_movie
+$PBS_O_WORKDIR/xcreate_movie_AVS_DX_noattenuation < temp_input_movie
+# backup step 1 results (without attenuation)
+DIR_TEMP=$DIR_RESULT/step_1_no_attenuation
+rm -rf $DIR_TEMP
+mkdir -p $DIR_TEMP
+mv OUTPUT_FILES/* $DIR_TEMP/
+cp DATA/Par_file $DIR_TEMP/
+cp DATA/CMTSOLUTION $DIR_TEMP/
+cp DATA/STATIONS $DIR_TEMP/
+################################################################ step 2: ensemble forward wavefield
+cp DATA/Par_file_NOISE_2_noattenuation DATA/Par_file
+mpiexec $PBS_O_WORKDIR/xmeshfem3D_noattenuation
+mpiexec $PBS_O_WORKDIR/xspecfem3D_noattenuation
+echo "2" > temp_input_movie; echo "1" >> temp_input_movie; echo "$NSTEP" >> temp_input_movie
+$PBS_O_WORKDIR/xcreate_movie_AVS_DX_noattenuation < temp_input_movie
+# backup step 2 results (without attenuation)
+DIR_TEMP=$DIR_RESULT/step_2_no_attenuation
+rm -rf $DIR_TEMP
+mkdir -p $DIR_TEMP
+mv OUTPUT_FILES/* $DIR_TEMP/
+cp DATA/Par_file $DIR_TEMP/
+cp DATA/CMTSOLUTION $DIR_TEMP/
+cp DATA/STATIONS $DIR_TEMP/
+############################################################## step 3: ensemble kernels calculation
+cp DATA/Par_file_NOISE_3_noattenuation DATA/Par_file
+mpiexec $PBS_O_WORKDIR/xmeshfem3D_noattenuation
+mpiexec $PBS_O_WORKDIR/xspecfem3D_noattenuation
+echo "2" > temp_input_movie; echo "1" >> temp_input_movie; echo "$NSTEP" >> temp_input_movie
+$PBS_O_WORKDIR/xcreate_movie_AVS_DX_noattenuation < temp_input_movie
+
+## collect outputs (e.g., kernels, meshes) from local nodes to global disk
+d=`date`
+echo "Start collecting outputs from nodes, $d"
+pbsdsh $PBS_O_WORKDIR/NOISE_collect $LOCAL_PATH $PBS_O_WORKDIR/OUTPUT_FILES/
+d=`date`
+echo "Finish collecting outputs from nodes, $d"
+## combine_vol_data
+echo "0" > temp_input_movie
+for (( iproc=1;iproc<=$numnodes-1;iproc++ ))
+do
+ echo "$iproc" >> temp_input_movie
+done
+$PBS_O_WORKDIR/xcombine_vol_data temp_input_movie alpha_kernel $PBS_O_WORKDIR/OUTPUT_FILES $PBS_O_WORKDIR/OUTPUT_FILES $PBS_O_WORKDIR/OUTPUT_FILES 0 1
+rm temp_input_movie
+
+## backup step 3 results (without attenuation)
+DIR_TEMP=$DIR_RESULT/step_3_no_attenuation
+rm -rf $DIR_TEMP
+mkdir -p $DIR_TEMP
+mv OUTPUT_FILES/* $DIR_TEMP/
+cp DATA/Par_file $DIR_TEMP/
+cp DATA/CMTSOLUTION $DIR_TEMP/
+cp DATA/STATIONS $DIR_TEMP/
+
+## clean local nodes
+pbsdsh $PBS_O_WORKDIR/NOISE_clean_create $LOCAL_PATH
+d=`date`
+echo "Finish cleaning local disk on nodes, $d"
Property changes on: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/NOISE123.submit_atten
___________________________________________________________________
Name: svn:executable
+ *
Added: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/NOISE_clean_create
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/NOISE_clean_create (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/NOISE_clean_create 2011-04-01 15:57:37 UTC (rev 18161)
@@ -0,0 +1,10 @@
+#!/bin/bash -eu
+
+DIR_LOCAL=$1
+iproc=$PBS_VNODENUM
+
+if [ $(($iproc % 8)) = 0 ]; then
+ rm -rf $DIR_LOCAL
+ mkdir -p $DIR_LOCAL
+fi
+
Property changes on: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/NOISE_clean_create
___________________________________________________________________
Name: svn:executable
+ *
Added: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/NOISE_collect
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/NOISE_collect (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/NOISE_collect 2011-04-01 15:57:37 UTC (rev 18161)
@@ -0,0 +1,11 @@
+#!/bin/bash -eu
+
+DIR_LOCAL=$1
+DIR_GLOBAL=$2
+iproc=`printf "%06d" $PBS_VNODENUM`
+
+##### other files can be collected in the same way
+mv $DIR_LOCAL/*$iproc*reg1*kernel*.bin $DIR_GLOBAL/
+mv $DIR_LOCAL/*$iproc*reg1_*array_dims.txt $DIR_GLOBAL/
+mv $DIR_LOCAL/*$iproc*reg1_solver_data*.bin $DIR_GLOBAL/
+
Property changes on: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/NOISE_collect
___________________________________________________________________
Name: svn:executable
+ *
Added: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/Par_file_NOISE_1_attenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/Par_file_NOISE_1_attenuation (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/Par_file_NOISE_1_attenuation 2011-04-01 15:57:37 UTC (rev 18161)
@@ -0,0 +1,120 @@
+
+# forward or adjoint simulation
+SIMULATION_TYPE = 1
+NOISE_TOMOGRAPHY = 1 # flag of noise tomography, three steps (1,2,3). If earthquake tomography, set it to 0.
+SAVE_FORWARD = .true. # save last frame of forward simulation or not
+
+# number of chunks (1,2,3 or 6)
+NCHUNKS = 6
+
+# angular width of the first chunk (not used if full sphere with six chunks)
+ANGULAR_WIDTH_XI_IN_DEGREES = 90.d0 # angular size of a chunk
+ANGULAR_WIDTH_ETA_IN_DEGREES = 90.d0
+CENTER_LATITUDE_IN_DEGREES = 0.d0
+CENTER_LONGITUDE_IN_DEGREES = 15.d0
+GAMMA_ROTATION_AZIMUTH = 0.d0
+
+# number of elements at the surface along the two sides of the first chunk
+# (must be multiple of 16 and 8 * multiple of NPROC below)
+NEX_XI = 160
+NEX_ETA = 160
+
+# number of MPI processors along the two sides of the first chunk
+NPROC_XI = 10
+NPROC_ETA = 10
+
+# 1D models with real structure:
+# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
+#
+# 1D models with only one fictitious averaged crustal layer:
+# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135_onecrust
+#
+# fully 3D models:
+# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
+# s20rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
+# s29ea, s29ea,sea99_jp3d1994,sea99,jp3d1994,heterogen
+#
+# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
+# to take the 1D crustal model from the
+# associated reference model rather than the default 3D crustal model
+# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
+MODEL = 1D_isotropic_prem
+
+# parameters describing the Earth model
+OCEANS = .false.
+ELLIPTICITY = .false.
+TOPOGRAPHY = .false.
+GRAVITY = .false.
+ROTATION = .false.
+ATTENUATION = .true.
+
+# absorbing boundary conditions for a regional simulation
+ABSORBING_CONDITIONS = .false.
+
+# record length in minutes
+RECORD_LENGTH_IN_MINUTES = 1.0d0
+
+# save AVS or OpenDX movies
+#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
+#MOVIE_COARSE does not work with create_movie_AVS_DX
+MOVIE_SURFACE = .true.
+MOVIE_VOLUME = .false.
+MOVIE_COARSE = .false.
+NTSTEP_BETWEEN_FRAMES = 50
+HDUR_MOVIE = 0.d0
+
+# save movie in volume. Will save element if center of element is in prescribed volume
+# top/bottom: depth in KM, use MOVIE_TOP = -100 to make sure the surface is stored.
+# west/east: longitude, degrees East [-180/180] top/bottom: latitute, degrees North [-90/90]
+# start/stop: frames will be stored at MOVIE_START + i*NSTEP_BETWEEN_FRAMES, where i=(0,1,2..) and iNSTEP_BETWEEN_FRAMES <= MOVIE_STOP
+# movie_volume_type: 1=strain, 2=time integral of strain, 3=\mu*time integral of strain
+# type 4 saves the trace and deviatoric stress in the whole volume, 5=displacement, 6=velocity
+MOVIE_VOLUME_TYPE = 2
+MOVIE_TOP_KM = -100.0
+MOVIE_BOTTOM_KM = 1000.0
+MOVIE_WEST_DEG = -90.0
+MOVIE_EAST_DEG = 90.0
+MOVIE_NORTH_DEG = 90.0
+MOVIE_SOUTH_DEG = -90.0
+MOVIE_START = 0
+MOVIE_STOP = 40000
+
+# save mesh files to check the mesh
+SAVE_MESH_FILES = .false.
+
+# restart files (number of runs can be 1, 2 or 3, choose 1 for no restart files)
+NUMBER_OF_RUNS = 1
+NUMBER_OF_THIS_RUN = 1
+
+# path to store the local database files on each node
+LOCAL_PATH = /scratch/yangl/DATABASES_MPI
+#/scratch/yangl/DATABASES_MPI
+
+# interval at which we output time step info and max of norm of displacement
+NTSTEP_BETWEEN_OUTPUT_INFO = 10000
+
+# interval in time steps for temporary writing of seismograms
+NTSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000
+NTSTEP_BETWEEN_READ_ADJSRC = 1000
+
+# output format for the seismograms (one can use either or all of the three formats)
+OUTPUT_SEISMOS_ASCII_TEXT = .true.
+OUTPUT_SEISMOS_SAC_ALPHANUM = .false.
+OUTPUT_SEISMOS_SAC_BINARY = .false.
+
+# rotate seismograms to Radial-Transverse-Z or use default North-East-Z reference frame
+ROTATE_SEISMOGRAMS_RT = .false.
+
+# decide if master process writes all the seismograms or if all processes do it in parallel
+WRITE_SEISMOGRAMS_BY_MASTER = .true.
+
+# save all seismograms in one large combined file instead of one file per seismogram
+# to avoid overloading shared non-local file systems such as GPFS for instance
+SAVE_ALL_SEISMOS_IN_ONE_FILE = .false.
+USE_BINARY_FOR_LARGE_FILE = .false.
+
+# flag to impose receivers at the surface or allow them to be buried
+RECEIVERS_CAN_BE_BURIED = .true.
+
+# print source time function
+PRINT_SOURCE_TIME_FUNCTION = .true.
Added: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/Par_file_NOISE_1_noattenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/Par_file_NOISE_1_noattenuation (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/Par_file_NOISE_1_noattenuation 2011-04-01 15:57:37 UTC (rev 18161)
@@ -0,0 +1,120 @@
+
+# forward or adjoint simulation
+SIMULATION_TYPE = 1
+NOISE_TOMOGRAPHY = 1 # flag of noise tomography, three steps (1,2,3). If earthquake tomography, set it to 0.
+SAVE_FORWARD = .true. # save last frame of forward simulation or not
+
+# number of chunks (1,2,3 or 6)
+NCHUNKS = 6
+
+# angular width of the first chunk (not used if full sphere with six chunks)
+ANGULAR_WIDTH_XI_IN_DEGREES = 90.d0 # angular size of a chunk
+ANGULAR_WIDTH_ETA_IN_DEGREES = 90.d0
+CENTER_LATITUDE_IN_DEGREES = 0.d0
+CENTER_LONGITUDE_IN_DEGREES = 15.d0
+GAMMA_ROTATION_AZIMUTH = 0.d0
+
+# number of elements at the surface along the two sides of the first chunk
+# (must be multiple of 16 and 8 * multiple of NPROC below)
+NEX_XI = 160
+NEX_ETA = 160
+
+# number of MPI processors along the two sides of the first chunk
+NPROC_XI = 10
+NPROC_ETA = 10
+
+# 1D models with real structure:
+# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
+#
+# 1D models with only one fictitious averaged crustal layer:
+# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135_onecrust
+#
+# fully 3D models:
+# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
+# s20rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
+# s29ea, s29ea,sea99_jp3d1994,sea99,jp3d1994,heterogen
+#
+# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
+# to take the 1D crustal model from the
+# associated reference model rather than the default 3D crustal model
+# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
+MODEL = 1D_isotropic_prem
+
+# parameters describing the Earth model
+OCEANS = .false.
+ELLIPTICITY = .false.
+TOPOGRAPHY = .false.
+GRAVITY = .false.
+ROTATION = .false.
+ATTENUATION = .false.
+
+# absorbing boundary conditions for a regional simulation
+ABSORBING_CONDITIONS = .false.
+
+# record length in minutes
+RECORD_LENGTH_IN_MINUTES = 1.0d0
+
+# save AVS or OpenDX movies
+#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
+#MOVIE_COARSE does not work with create_movie_AVS_DX
+MOVIE_SURFACE = .true.
+MOVIE_VOLUME = .false.
+MOVIE_COARSE = .false.
+NTSTEP_BETWEEN_FRAMES = 50
+HDUR_MOVIE = 0.d0
+
+# save movie in volume. Will save element if center of element is in prescribed volume
+# top/bottom: depth in KM, use MOVIE_TOP = -100 to make sure the surface is stored.
+# west/east: longitude, degrees East [-180/180] top/bottom: latitute, degrees North [-90/90]
+# start/stop: frames will be stored at MOVIE_START + i*NSTEP_BETWEEN_FRAMES, where i=(0,1,2..) and iNSTEP_BETWEEN_FRAMES <= MOVIE_STOP
+# movie_volume_type: 1=strain, 2=time integral of strain, 3=\mu*time integral of strain
+# type 4 saves the trace and deviatoric stress in the whole volume, 5=displacement, 6=velocity
+MOVIE_VOLUME_TYPE = 2
+MOVIE_TOP_KM = -100.0
+MOVIE_BOTTOM_KM = 1000.0
+MOVIE_WEST_DEG = -90.0
+MOVIE_EAST_DEG = 90.0
+MOVIE_NORTH_DEG = 90.0
+MOVIE_SOUTH_DEG = -90.0
+MOVIE_START = 0
+MOVIE_STOP = 40000
+
+# save mesh files to check the mesh
+SAVE_MESH_FILES = .false.
+
+# restart files (number of runs can be 1, 2 or 3, choose 1 for no restart files)
+NUMBER_OF_RUNS = 1
+NUMBER_OF_THIS_RUN = 1
+
+# path to store the local database files on each node
+LOCAL_PATH = /scratch/yangl/DATABASES_MPI
+#/scratch/yangl/DATABASES_MPI
+
+# interval at which we output time step info and max of norm of displacement
+NTSTEP_BETWEEN_OUTPUT_INFO = 10000
+
+# interval in time steps for temporary writing of seismograms
+NTSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000
+NTSTEP_BETWEEN_READ_ADJSRC = 1000
+
+# output format for the seismograms (one can use either or all of the three formats)
+OUTPUT_SEISMOS_ASCII_TEXT = .true.
+OUTPUT_SEISMOS_SAC_ALPHANUM = .false.
+OUTPUT_SEISMOS_SAC_BINARY = .false.
+
+# rotate seismograms to Radial-Transverse-Z or use default North-East-Z reference frame
+ROTATE_SEISMOGRAMS_RT = .false.
+
+# decide if master process writes all the seismograms or if all processes do it in parallel
+WRITE_SEISMOGRAMS_BY_MASTER = .true.
+
+# save all seismograms in one large combined file instead of one file per seismogram
+# to avoid overloading shared non-local file systems such as GPFS for instance
+SAVE_ALL_SEISMOS_IN_ONE_FILE = .false.
+USE_BINARY_FOR_LARGE_FILE = .false.
+
+# flag to impose receivers at the surface or allow them to be buried
+RECEIVERS_CAN_BE_BURIED = .true.
+
+# print source time function
+PRINT_SOURCE_TIME_FUNCTION = .true.
Added: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/Par_file_NOISE_2_attenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/Par_file_NOISE_2_attenuation (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/Par_file_NOISE_2_attenuation 2011-04-01 15:57:37 UTC (rev 18161)
@@ -0,0 +1,120 @@
+
+# forward or adjoint simulation
+SIMULATION_TYPE = 1
+NOISE_TOMOGRAPHY = 2 # flag of noise tomography, three steps (1,2,3). If earthquake tomography, set it to 0.
+SAVE_FORWARD = .true. # save last frame of forward simulation or not
+
+# number of chunks (1,2,3 or 6)
+NCHUNKS = 6
+
+# angular width of the first chunk (not used if full sphere with six chunks)
+ANGULAR_WIDTH_XI_IN_DEGREES = 90.d0 # angular size of a chunk
+ANGULAR_WIDTH_ETA_IN_DEGREES = 90.d0
+CENTER_LATITUDE_IN_DEGREES = 0.d0
+CENTER_LONGITUDE_IN_DEGREES = 15.d0
+GAMMA_ROTATION_AZIMUTH = 0.d0
+
+# number of elements at the surface along the two sides of the first chunk
+# (must be multiple of 16 and 8 * multiple of NPROC below)
+NEX_XI = 160
+NEX_ETA = 160
+
+# number of MPI processors along the two sides of the first chunk
+NPROC_XI = 10
+NPROC_ETA = 10
+
+# 1D models with real structure:
+# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
+#
+# 1D models with only one fictitious averaged crustal layer:
+# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135_onecrust
+#
+# fully 3D models:
+# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
+# s20rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
+# s29ea, s29ea,sea99_jp3d1994,sea99,jp3d1994,heterogen
+#
+# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
+# to take the 1D crustal model from the
+# associated reference model rather than the default 3D crustal model
+# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
+MODEL = 1D_isotropic_prem
+
+# parameters describing the Earth model
+OCEANS = .false.
+ELLIPTICITY = .false.
+TOPOGRAPHY = .false.
+GRAVITY = .false.
+ROTATION = .false.
+ATTENUATION = .true.
+
+# absorbing boundary conditions for a regional simulation
+ABSORBING_CONDITIONS = .false.
+
+# record length in minutes
+RECORD_LENGTH_IN_MINUTES = 1.0d0
+
+# save AVS or OpenDX movies
+#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
+#MOVIE_COARSE does not work with create_movie_AVS_DX
+MOVIE_SURFACE = .true.
+MOVIE_VOLUME = .false.
+MOVIE_COARSE = .false.
+NTSTEP_BETWEEN_FRAMES = 50
+HDUR_MOVIE = 0.d0
+
+# save movie in volume. Will save element if center of element is in prescribed volume
+# top/bottom: depth in KM, use MOVIE_TOP = -100 to make sure the surface is stored.
+# west/east: longitude, degrees East [-180/180] top/bottom: latitute, degrees North [-90/90]
+# start/stop: frames will be stored at MOVIE_START + i*NSTEP_BETWEEN_FRAMES, where i=(0,1,2..) and iNSTEP_BETWEEN_FRAMES <= MOVIE_STOP
+# movie_volume_type: 1=strain, 2=time integral of strain, 3=\mu*time integral of strain
+# type 4 saves the trace and deviatoric stress in the whole volume, 5=displacement, 6=velocity
+MOVIE_VOLUME_TYPE = 2
+MOVIE_TOP_KM = -100.0
+MOVIE_BOTTOM_KM = 1000.0
+MOVIE_WEST_DEG = -90.0
+MOVIE_EAST_DEG = 90.0
+MOVIE_NORTH_DEG = 90.0
+MOVIE_SOUTH_DEG = -90.0
+MOVIE_START = 0
+MOVIE_STOP = 40000
+
+# save mesh files to check the mesh
+SAVE_MESH_FILES = .false.
+
+# restart files (number of runs can be 1, 2 or 3, choose 1 for no restart files)
+NUMBER_OF_RUNS = 1
+NUMBER_OF_THIS_RUN = 1
+
+# path to store the local database files on each node
+LOCAL_PATH = /scratch/yangl/DATABASES_MPI
+#/scratch/yangl/DATABASES_MPI
+
+# interval at which we output time step info and max of norm of displacement
+NTSTEP_BETWEEN_OUTPUT_INFO = 10000
+
+# interval in time steps for temporary writing of seismograms
+NTSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000
+NTSTEP_BETWEEN_READ_ADJSRC = 1000
+
+# output format for the seismograms (one can use either or all of the three formats)
+OUTPUT_SEISMOS_ASCII_TEXT = .true.
+OUTPUT_SEISMOS_SAC_ALPHANUM = .false.
+OUTPUT_SEISMOS_SAC_BINARY = .false.
+
+# rotate seismograms to Radial-Transverse-Z or use default North-East-Z reference frame
+ROTATE_SEISMOGRAMS_RT = .false.
+
+# decide if master process writes all the seismograms or if all processes do it in parallel
+WRITE_SEISMOGRAMS_BY_MASTER = .true.
+
+# save all seismograms in one large combined file instead of one file per seismogram
+# to avoid overloading shared non-local file systems such as GPFS for instance
+SAVE_ALL_SEISMOS_IN_ONE_FILE = .false.
+USE_BINARY_FOR_LARGE_FILE = .false.
+
+# flag to impose receivers at the surface or allow them to be buried
+RECEIVERS_CAN_BE_BURIED = .true.
+
+# print source time function
+PRINT_SOURCE_TIME_FUNCTION = .true.
Added: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/Par_file_NOISE_2_noattenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/Par_file_NOISE_2_noattenuation (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/Par_file_NOISE_2_noattenuation 2011-04-01 15:57:37 UTC (rev 18161)
@@ -0,0 +1,120 @@
+
+# forward or adjoint simulation
+SIMULATION_TYPE = 1
+NOISE_TOMOGRAPHY = 2 # flag of noise tomography, three steps (1,2,3). If earthquake tomography, set it to 0.
+SAVE_FORWARD = .true. # save last frame of forward simulation or not
+
+# number of chunks (1,2,3 or 6)
+NCHUNKS = 6
+
+# angular width of the first chunk (not used if full sphere with six chunks)
+ANGULAR_WIDTH_XI_IN_DEGREES = 90.d0 # angular size of a chunk
+ANGULAR_WIDTH_ETA_IN_DEGREES = 90.d0
+CENTER_LATITUDE_IN_DEGREES = 0.d0
+CENTER_LONGITUDE_IN_DEGREES = 15.d0
+GAMMA_ROTATION_AZIMUTH = 0.d0
+
+# number of elements at the surface along the two sides of the first chunk
+# (must be multiple of 16 and 8 * multiple of NPROC below)
+NEX_XI = 160
+NEX_ETA = 160
+
+# number of MPI processors along the two sides of the first chunk
+NPROC_XI = 10
+NPROC_ETA = 10
+
+# 1D models with real structure:
+# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
+#
+# 1D models with only one fictitious averaged crustal layer:
+# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135_onecrust
+#
+# fully 3D models:
+# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
+# s20rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
+# s29ea, s29ea,sea99_jp3d1994,sea99,jp3d1994,heterogen
+#
+# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
+# to take the 1D crustal model from the
+# associated reference model rather than the default 3D crustal model
+# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
+MODEL = 1D_isotropic_prem
+
+# parameters describing the Earth model
+OCEANS = .false.
+ELLIPTICITY = .false.
+TOPOGRAPHY = .false.
+GRAVITY = .false.
+ROTATION = .false.
+ATTENUATION = .false.
+
+# absorbing boundary conditions for a regional simulation
+ABSORBING_CONDITIONS = .false.
+
+# record length in minutes
+RECORD_LENGTH_IN_MINUTES = 1.0d0
+
+# save AVS or OpenDX movies
+#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
+#MOVIE_COARSE does not work with create_movie_AVS_DX
+MOVIE_SURFACE = .true.
+MOVIE_VOLUME = .false.
+MOVIE_COARSE = .false.
+NTSTEP_BETWEEN_FRAMES = 50
+HDUR_MOVIE = 0.d0
+
+# save movie in volume. Will save element if center of element is in prescribed volume
+# top/bottom: depth in KM, use MOVIE_TOP = -100 to make sure the surface is stored.
+# west/east: longitude, degrees East [-180/180] top/bottom: latitute, degrees North [-90/90]
+# start/stop: frames will be stored at MOVIE_START + i*NSTEP_BETWEEN_FRAMES, where i=(0,1,2..) and iNSTEP_BETWEEN_FRAMES <= MOVIE_STOP
+# movie_volume_type: 1=strain, 2=time integral of strain, 3=\mu*time integral of strain
+# type 4 saves the trace and deviatoric stress in the whole volume, 5=displacement, 6=velocity
+MOVIE_VOLUME_TYPE = 2
+MOVIE_TOP_KM = -100.0
+MOVIE_BOTTOM_KM = 1000.0
+MOVIE_WEST_DEG = -90.0
+MOVIE_EAST_DEG = 90.0
+MOVIE_NORTH_DEG = 90.0
+MOVIE_SOUTH_DEG = -90.0
+MOVIE_START = 0
+MOVIE_STOP = 40000
+
+# save mesh files to check the mesh
+SAVE_MESH_FILES = .false.
+
+# restart files (number of runs can be 1, 2 or 3, choose 1 for no restart files)
+NUMBER_OF_RUNS = 1
+NUMBER_OF_THIS_RUN = 1
+
+# path to store the local database files on each node
+LOCAL_PATH = /scratch/yangl/DATABASES_MPI
+#/scratch/yangl/DATABASES_MPI
+
+# interval at which we output time step info and max of norm of displacement
+NTSTEP_BETWEEN_OUTPUT_INFO = 10000
+
+# interval in time steps for temporary writing of seismograms
+NTSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000
+NTSTEP_BETWEEN_READ_ADJSRC = 1000
+
+# output format for the seismograms (one can use either or all of the three formats)
+OUTPUT_SEISMOS_ASCII_TEXT = .true.
+OUTPUT_SEISMOS_SAC_ALPHANUM = .false.
+OUTPUT_SEISMOS_SAC_BINARY = .false.
+
+# rotate seismograms to Radial-Transverse-Z or use default North-East-Z reference frame
+ROTATE_SEISMOGRAMS_RT = .false.
+
+# decide if master process writes all the seismograms or if all processes do it in parallel
+WRITE_SEISMOGRAMS_BY_MASTER = .true.
+
+# save all seismograms in one large combined file instead of one file per seismogram
+# to avoid overloading shared non-local file systems such as GPFS for instance
+SAVE_ALL_SEISMOS_IN_ONE_FILE = .false.
+USE_BINARY_FOR_LARGE_FILE = .false.
+
+# flag to impose receivers at the surface or allow them to be buried
+RECEIVERS_CAN_BE_BURIED = .true.
+
+# print source time function
+PRINT_SOURCE_TIME_FUNCTION = .true.
Added: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/Par_file_NOISE_3_noattenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/Par_file_NOISE_3_noattenuation (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/Par_file_NOISE_3_noattenuation 2011-04-01 15:57:37 UTC (rev 18161)
@@ -0,0 +1,120 @@
+
+# forward or adjoint simulation
+SIMULATION_TYPE = 3
+NOISE_TOMOGRAPHY = 3 # flag of noise tomography, three steps (1,2,3). If earthquake tomography, set it to 0.
+SAVE_FORWARD = .false. # save last frame of forward simulation or not
+
+# number of chunks (1,2,3 or 6)
+NCHUNKS = 6
+
+# angular width of the first chunk (not used if full sphere with six chunks)
+ANGULAR_WIDTH_XI_IN_DEGREES = 90.d0 # angular size of a chunk
+ANGULAR_WIDTH_ETA_IN_DEGREES = 90.d0
+CENTER_LATITUDE_IN_DEGREES = 0.d0
+CENTER_LONGITUDE_IN_DEGREES = 15.d0
+GAMMA_ROTATION_AZIMUTH = 0.d0
+
+# number of elements at the surface along the two sides of the first chunk
+# (must be multiple of 16 and 8 * multiple of NPROC below)
+NEX_XI = 160
+NEX_ETA = 160
+
+# number of MPI processors along the two sides of the first chunk
+NPROC_XI = 10
+NPROC_ETA = 10
+
+# 1D models with real structure:
+# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
+#
+# 1D models with only one fictitious averaged crustal layer:
+# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135_onecrust
+#
+# fully 3D models:
+# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
+# s20rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
+# s29ea, s29ea,sea99_jp3d1994,sea99,jp3d1994,heterogen
+#
+# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
+# to take the 1D crustal model from the
+# associated reference model rather than the default 3D crustal model
+# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
+MODEL = 1D_isotropic_prem
+
+# parameters describing the Earth model
+OCEANS = .false.
+ELLIPTICITY = .false.
+TOPOGRAPHY = .false.
+GRAVITY = .false.
+ROTATION = .false.
+ATTENUATION = .false.
+
+# absorbing boundary conditions for a regional simulation
+ABSORBING_CONDITIONS = .false.
+
+# record length in minutes
+RECORD_LENGTH_IN_MINUTES = 1.0d0
+
+# save AVS or OpenDX movies
+#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
+#MOVIE_COARSE does not work with create_movie_AVS_DX
+MOVIE_SURFACE = .true.
+MOVIE_VOLUME = .false.
+MOVIE_COARSE = .false.
+NTSTEP_BETWEEN_FRAMES = 50
+HDUR_MOVIE = 0.d0
+
+# save movie in volume. Will save element if center of element is in prescribed volume
+# top/bottom: depth in KM, use MOVIE_TOP = -100 to make sure the surface is stored.
+# west/east: longitude, degrees East [-180/180] top/bottom: latitute, degrees North [-90/90]
+# start/stop: frames will be stored at MOVIE_START + i*NSTEP_BETWEEN_FRAMES, where i=(0,1,2..) and iNSTEP_BETWEEN_FRAMES <= MOVIE_STOP
+# movie_volume_type: 1=strain, 2=time integral of strain, 3=\mu*time integral of strain
+# type 4 saves the trace and deviatoric stress in the whole volume, 5=displacement, 6=velocity
+MOVIE_VOLUME_TYPE = 2
+MOVIE_TOP_KM = -100.0
+MOVIE_BOTTOM_KM = 1000.0
+MOVIE_WEST_DEG = -90.0
+MOVIE_EAST_DEG = 90.0
+MOVIE_NORTH_DEG = 90.0
+MOVIE_SOUTH_DEG = -90.0
+MOVIE_START = 0
+MOVIE_STOP = 40000
+
+# save mesh files to check the mesh
+SAVE_MESH_FILES = .false.
+
+# restart files (number of runs can be 1, 2 or 3, choose 1 for no restart files)
+NUMBER_OF_RUNS = 1
+NUMBER_OF_THIS_RUN = 1
+
+# path to store the local database files on each node
+LOCAL_PATH = /scratch/yangl/DATABASES_MPI
+#/scratch/yangl/DATABASES_MPI
+
+# interval at which we output time step info and max of norm of displacement
+NTSTEP_BETWEEN_OUTPUT_INFO = 10000
+
+# interval in time steps for temporary writing of seismograms
+NTSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000
+NTSTEP_BETWEEN_READ_ADJSRC = 1000
+
+# output format for the seismograms (one can use either or all of the three formats)
+OUTPUT_SEISMOS_ASCII_TEXT = .true.
+OUTPUT_SEISMOS_SAC_ALPHANUM = .false.
+OUTPUT_SEISMOS_SAC_BINARY = .false.
+
+# rotate seismograms to Radial-Transverse-Z or use default North-East-Z reference frame
+ROTATE_SEISMOGRAMS_RT = .false.
+
+# decide if master process writes all the seismograms or if all processes do it in parallel
+WRITE_SEISMOGRAMS_BY_MASTER = .true.
+
+# save all seismograms in one large combined file instead of one file per seismogram
+# to avoid overloading shared non-local file systems such as GPFS for instance
+SAVE_ALL_SEISMOS_IN_ONE_FILE = .false.
+USE_BINARY_FOR_LARGE_FILE = .false.
+
+# flag to impose receivers at the surface or allow them to be buried
+RECEIVERS_CAN_BE_BURIED = .true.
+
+# print source time function
+PRINT_SOURCE_TIME_FUNCTION = .true.
Added: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/STATIONS_NOISE
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/STATIONS_NOISE (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/STATIONS_NOISE 2011-04-01 15:57:37 UTC (rev 18161)
@@ -0,0 +1,13 @@
+A1 II 0.0 0.0 0.0 0.0
+A2 II 0.0 5.0 0.0 0.0
+A3 II 0.0 10.0 0.0 0.0
+A4 II 0.0 15.0 0.0 0.0
+A5 II 0.0 20.0 0.0 0.0
+A6 II 0.0 25.0 0.0 0.0
+A7 II 0.0 30.0 0.0 0.0
+A8 II 0.0 35.0 0.0 0.0
+A9 II 0.0 40.0 0.0 0.0
+A10 II 0.0 45.0 0.0 0.0
+A11 II 0.0 50.0 0.0 0.0
+A12 II 0.0 55.0 0.0 0.0
+A13 II 0.0 60.0 0.0 0.0
Added: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/STATIONS_NOISE_ADJOINT
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/STATIONS_NOISE_ADJOINT (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/STATIONS_NOISE_ADJOINT 2011-04-01 15:57:37 UTC (rev 18161)
@@ -0,0 +1 @@
+A7 II 0.0 30.0 0.0 0.0
Added: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/S_squared
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/S_squared (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/S_squared 2011-04-01 15:57:37 UTC (rev 18161)
@@ -0,0 +1,799 @@
+ -7.5810000e+01 2.0083475e-19
+ -7.5620000e+01 1.9220480e-19
+ -7.5430000e+01 1.7529785e-19
+ -7.5240000e+01 1.5080677e-19
+ -7.5050000e+01 1.1973878e-19
+ -7.4860000e+01 8.3378339e-20
+ -7.4670000e+01 4.3239679e-20
+ -7.4480000e+01 1.0106876e-21
+ -7.4290000e+01 -4.1510205e-20
+ -7.4100000e+01 -8.2482300e-20
+ -7.3910000e+01 -1.2009421e-19
+ -7.3720000e+01 -1.5263679e-19
+ -7.3530000e+01 -1.7857550e-19
+ -7.3340000e+01 -1.9661979e-19
+ -7.3150000e+01 -2.0578729e-19
+ -7.2960000e+01 -2.0545992e-19
+ -7.2770000e+01 -1.9542883e-19
+ -7.2580000e+01 -1.7592525e-19
+ -7.2390000e+01 -1.4763374e-19
+ -7.2200000e+01 -1.1168541e-19
+ -7.2010000e+01 -6.9628508e-20
+ -7.1820000e+01 -2.3375690e-20
+ -7.1630000e+01 2.4872197e-20
+ -7.1440000e+01 7.2723698e-20
+ -7.1250000e+01 1.1771219e-19
+ -7.1060000e+01 1.5742382e-19
+ -7.0870000e+01 1.8963105e-19
+ -7.0680000e+01 2.1242398e-19
+ -7.0490000e+01 2.2433077e-19
+ -7.0300000e+01 2.2441895e-19
+ -7.0110000e+01 2.1236989e-19
+ -6.9920000e+01 1.8852037e-19
+ -6.9730000e+01 1.5386707e-19
+ -6.9540000e+01 1.1003221e-19
+ -6.9350000e+01 5.9191273e-20
+ -6.9160000e+01 3.9661551e-21
+ -6.8970000e+01 -5.2710589e-20
+ -6.8780000e+01 -1.0775270e-19
+ -6.8590000e+01 -1.5808949e-19
+ -6.8400000e+01 -2.0083774e-19
+ -6.8210000e+01 -2.3346641e-19
+ -6.8020000e+01 -2.5394490e-19
+ -6.7830000e+01 -2.6086654e-19
+ -6.7640000e+01 -2.5354037e-19
+ -6.7450000e+01 -2.3204589e-19
+ -6.7260000e+01 -1.9724741e-19
+ -6.7070000e+01 -1.5076633e-19
+ -6.6880000e+01 -9.4911984e-20
+ -6.6690000e+01 -3.2573180e-20
+ -6.6500000e+01 3.2925061e-20
+ -6.6310000e+01 9.7993561e-20
+ -6.6120000e+01 1.5897190e-19
+ -6.5930000e+01 2.1233279e-19
+ -6.5740000e+01 2.5488724e-19
+ -6.5550000e+01 2.8397889e-19
+ -6.5360000e+01 2.9765587e-19
+ -6.5170000e+01 2.9480847e-19
+ -6.4980000e+01 2.7526252e-19
+ -6.4790000e+01 2.3982016e-19
+ -6.4600000e+01 1.9024237e-19
+ -6.4410000e+01 1.2917147e-19
+ -6.4220000e+01 5.9995411e-20
+ -6.4030000e+01 -1.3339828e-20
+ -6.3840000e+01 -8.6559668e-20
+ -6.3650000e+01 -1.5531831e-19
+ -6.3460000e+01 -2.1546775e-19
+ -6.3270000e+01 -2.6332049e-19
+ -6.3080000e+01 -2.9588787e-19
+ -6.2890000e+01 -3.1107739e-19
+ -6.2700000e+01 -3.0783468e-19
+ -6.2510000e+01 -2.8622009e-19
+ -6.2320000e+01 -2.4741406e-19
+ -6.2130000e+01 -1.9365082e-19
+ -6.1940000e+01 -1.2808502e-19
+ -6.1750000e+01 -5.4600362e-20
+ -6.1560000e+01 2.2426535e-20
+ -6.1370000e+01 9.8391447e-20
+ -6.1180000e+01 1.6873968e-19
+ -6.0990000e+01 2.2923767e-19
+ -6.0800000e+01 2.7622279e-19
+ -6.0610000e+01 3.0681694e-19
+ -6.0420000e+01 3.1909258e-19
+ -6.0230000e+01 3.1218306e-19
+ -6.0040000e+01 2.8633230e-19
+ -5.9850000e+01 2.4288267e-19
+ -5.9660000e+01 1.8420267e-19
+ -5.9470000e+01 1.1355904e-19
+ -5.9280000e+01 3.4940085e-20
+ -5.9090000e+01 -4.7160837e-20
+ -5.8900000e+01 -1.2799106e-19
+ -5.8710000e+01 -2.0280178e-19
+ -5.8520000e+01 -2.6711570e-19
+ -5.8330000e+01 -3.1698554e-19
+ -5.8140000e+01 -3.4922950e-19
+ -5.7950000e+01 -3.6163018e-19
+ -5.7760000e+01 -3.5308443e-19
+ -5.7570000e+01 -3.2369325e-19
+ -5.7380000e+01 -2.7478380e-19
+ -5.7190000e+01 -2.0885893e-19
+ -5.7000000e+01 -1.2947405e-19
+ -5.6810000e+01 -4.1045823e-20
+ -5.6620000e+01 5.1397826e-20
+ -5.6430000e+01 1.4251234e-19
+ -5.6240000e+01 2.2695775e-19
+ -5.6050000e+01 2.9972579e-19
+ -5.5860000e+01 3.5645090e-19
+ -5.5670000e+01 3.9368369e-19
+ -5.5480000e+01 4.0910773e-19
+ -5.5290000e+01 4.0168438e-19
+ -5.5100000e+01 3.7171582e-19
+ -5.4910000e+01 3.2082240e-19
+ -5.4720000e+01 2.5183700e-19
+ -5.4530000e+01 1.6862509e-19
+ -5.4340000e+01 7.5844464e-20
+ -5.4150000e+01 -2.1337642e-20
+ -5.3960000e+01 -1.1754995e-19
+ -5.3770000e+01 -2.0752479e-19
+ -5.3580000e+01 -2.8638873e-19
+ -5.3390000e+01 -3.4991992e-19
+ -5.3200000e+01 -3.9475923e-19
+ -5.3010000e+01 -4.1856685e-19
+ -5.2820000e+01 -4.2011991e-19
+ -5.2630000e+01 -3.9935087e-19
+ -5.2440000e+01 -3.5732939e-19
+ -5.2250000e+01 -2.9619265e-19
+ -5.2060000e+01 -2.1903051e-19
+ -5.1870000e+01 -1.2973246e-19
+ -5.1680000e+01 -3.2803460e-20
+ -5.1490000e+01 6.6844511e-20
+ -5.1300000e+01 1.6411831e-19
+ -5.1110000e+01 2.5397989e-19
+ -5.0920000e+01 3.3168548e-19
+ -5.0730000e+01 3.9301904e-19
+ -5.0540000e+01 4.3451051e-19
+ -5.0350000e+01 4.5362837e-19
+ -5.0160000e+01 4.4893518e-19
+ -4.9970000e+01 4.2019525e-19
+ -4.9780000e+01 3.6842482e-19
+ -4.9590000e+01 2.9587708e-19
+ -4.9400000e+01 2.0595855e-19
+ -4.9210000e+01 1.0307694e-19
+ -4.9020000e+01 -7.5736812e-21
+ -4.8830000e+01 -1.2028009e-19
+ -4.8640000e+01 -2.2911445e-19
+ -4.8450000e+01 -3.2825766e-19
+ -4.8260000e+01 -4.1232107e-19
+ -4.8070000e+01 -4.7664596e-19
+ -4.7880000e+01 -5.1756148e-19
+ -4.7690000e+01 -5.3258546e-19
+ -4.7500000e+01 -5.2055661e-19
+ -4.7310000e+01 -4.8169230e-19
+ -4.7120000e+01 -4.1757152e-19
+ -4.6930000e+01 -3.3104772e-19
+ -4.6740000e+01 -2.2610035e-19
+ -4.6550000e+01 -1.0763689e-19
+ -4.6360000e+01 1.8741548e-20
+ -4.6170000e+01 1.4699816e-19
+ -4.5980000e+01 2.7092971e-19
+ -4.5790000e+01 3.8443735e-19
+ -4.5600000e+01 4.8179060e-19
+ -4.5410000e+01 5.5787848e-19
+ -4.5220000e+01 6.0844227e-19
+ -4.5030000e+01 6.3028490e-19
+ -4.4840000e+01 6.2145114e-19
+ -4.4650000e+01 5.8137141e-19
+ -4.4460000e+01 5.1096036e-19
+ -4.4270000e+01 4.1266036e-19
+ -4.4080000e+01 2.9041926e-19
+ -4.3890000e+01 1.4959319e-19
+ -4.3700000e+01 -3.2324241e-21
+ -4.3510000e+01 -1.6050562e-19
+ -4.3320000e+01 -3.1405071e-19
+ -4.3130000e+01 -4.5547972e-19
+ -4.2940000e+01 -5.7665400e-19
+ -4.2750000e+01 -6.7017014e-19
+ -4.2560000e+01 -7.2983910e-19
+ -4.2370000e+01 -7.5112406e-19
+ -4.2180000e+01 -7.3150248e-19
+ -4.1990000e+01 -6.7072072e-19
+ -4.1800000e+01 -5.7091573e-19
+ -4.1610000e+01 -4.3658770e-19
+ -4.1420000e+01 -2.7441787e-19
+ -4.1230000e+01 -9.2938263e-20
+ -4.1040000e+01 9.7928901e-20
+ -4.0850000e+01 2.8743343e-19
+ -4.0660000e+01 4.6462495e-19
+ -4.0470000e+01 6.1900506e-19
+ -4.0280000e+01 7.4116124e-19
+ -4.0090000e+01 8.2334194e-19
+ -3.9900000e+01 8.5993725e-19
+ -3.9710000e+01 8.4783677e-19
+ -3.9520000e+01 7.8664520e-19
+ -3.9330000e+01 6.7874560e-19
+ -3.9140000e+01 5.2920900e-19
+ -3.8950000e+01 3.4555773e-19
+ -3.8760000e+01 1.3739630e-19
+ -3.8570000e+01 -8.4071055e-20
+ -3.8380000e+01 -3.0661268e-19
+ -3.8190000e+01 -5.1760598e-19
+ -3.8000000e+01 -7.0470977e-19
+ -3.7810000e+01 -8.5654316e-19
+ -3.7620000e+01 -9.6333999e-19
+ -3.7430000e+01 -1.0175471e-18
+ -3.7240000e+01 -1.0143341e-18
+ -3.7050000e+01 -9.5198336e-19
+ -3.6860000e+01 -8.3212971e-19
+ -3.6670000e+01 -6.5982566e-19
+ -3.6480000e+01 -4.4341220e-19
+ -3.6290000e+01 -1.9418638e-19
+ -3.6100000e+01 7.4131800e-20
+ -3.5910000e+01 3.4611692e-19
+ -3.5720000e+01 6.0550846e-19
+ -3.5530000e+01 8.3618806e-19
+ -3.5340000e+01 1.0232050e-18
+ -3.5150000e+01 1.1537822e-18
+ -3.4960000e+01 1.2182298e-18
+ -3.4770000e+01 1.2106966e-18
+ -3.4580000e+01 1.1296951e-18
+ -3.4390000e+01 9.7835102e-19
+ -3.4200000e+01 7.6434700e-19
+ -3.4010000e+01 4.9955265e-19
+ -3.3820000e+01 1.9935691e-19
+ -3.3630000e+01 -1.1825671e-19
+ -3.3440000e+01 -4.3383120e-19
+ -3.3250000e+01 -7.2767627e-19
+ -3.3060000e+01 -9.8114099e-19
+ -3.2870000e+01 -1.1778295e-18
+ -3.2680000e+01 -1.3046756e-18
+ -3.2490000e+01 -1.3528028e-18
+ -3.2300000e+01 -1.3181151e-18
+ -3.2110000e+01 -1.2015820e-18
+ -3.1920000e+01 -1.0092053e-18
+ -3.1730000e+01 -7.5167615e-19
+ -3.1540000e+01 -4.4375100e-19
+ -3.1350000e+01 -1.0339019e-19
+ -3.1160000e+01 2.4928452e-19
+ -3.0970000e+01 5.9315068e-19
+ -3.0780000e+01 9.0729832e-19
+ -3.0590000e+01 1.1722309e-18
+ -3.0400000e+01 1.3709932e-18
+ -3.0210000e+01 1.4901705e-18
+ -3.0020000e+01 1.5207105e-18
+ -2.9830000e+01 1.4585276e-18
+ -2.9640000e+01 1.3048553e-18
+ -2.9450000e+01 1.0663233e-18
+ -2.9260000e+01 7.5474484e-19
+ -2.9070000e+01 3.8661012e-19
+ -2.8880000e+01 -1.7701900e-20
+ -2.8690000e+01 -4.3496996e-19
+ -2.8500000e+01 -8.4039233e-19
+ -2.8310000e+01 -1.2090055e-18
+ -2.8120000e+01 -1.5171856e-18
+ -2.7930000e+01 -1.7441435e-18
+ -2.7740000e+01 -1.8733199e-18
+ -2.7550000e+01 -1.8935795e-18
+ -2.7360000e+01 -1.8001133e-18
+ -2.7170000e+01 -1.5949702e-18
+ -2.6980000e+01 -1.2871610e-18
+ -2.6790000e+01 -8.9230633e-19
+ -2.6600000e+01 -4.3183467e-19
+ -2.6410000e+01 6.8230074e-20
+ -2.6220000e+01 5.7881722e-19
+ -2.6030000e+01 1.0695922e-18
+ -2.5840000e+01 1.5108511e-18
+ -2.5650000e+01 1.8753932e-18
+ -2.5460000e+01 2.1402407e-18
+ -2.5270000e+01 2.2880846e-18
+ -2.5080000e+01 2.3083610e-18
+ -2.4890000e+01 2.1978907e-18
+ -2.4700000e+01 1.9610503e-18
+ -2.4510000e+01 1.6094796e-18
+ -2.4320000e+01 1.1613619e-18
+ -2.4130000e+01 6.4034236e-19
+ -2.3940000e+01 7.4170036e-20
+ -2.3750000e+01 -5.0684259e-19
+ -2.3560000e+01 -1.0714449e-18
+ -2.3370000e+01 -1.5890677e-18
+ -2.3180000e+01 -2.0313591e-18
+ -2.2990000e+01 -2.3735726e-18
+ -2.2800000e+01 -2.5957604e-18
+ -2.2610000e+01 -2.6837397e-18
+ -2.2420000e+01 -2.6298080e-18
+ -2.2230000e+01 -2.4331948e-18
+ -2.2040000e+01 -2.1002378e-18
+ -2.1850000e+01 -1.6442742e-18
+ -2.1660000e+01 -1.0852384e-18
+ -2.1470000e+01 -4.4896117e-19
+ -2.1280000e+01 2.3382952e-19
+ -2.1090000e+01 9.2879539e-19
+ -2.0900000e+01 1.5995338e-18
+ -2.0710000e+01 2.2093363e-18
+ -2.0520000e+01 2.7230887e-18
+ -2.0330000e+01 3.1092217e-18
+ -2.0140000e+01 3.3416006e-18
+ -1.9950000e+01 3.4012322e-18
+ -1.9760000e+01 3.2776657e-18
+ -1.9570000e+01 2.9699741e-18
+ -1.9380000e+01 2.4872234e-18
+ -1.9190000e+01 1.8483709e-18
+ -1.9000000e+01 1.0815724e-18
+ -1.8810000e+01 2.2292353e-19
+ -1.8620000e+01 -6.8529585e-19
+ -1.8430000e+01 -1.5967693e-18
+ -1.8240000e+01 -2.4635475e-18
+ -1.8050000e+01 -3.2385574e-18
+ -1.7860000e+01 -3.8780665e-18
+ -1.7670000e+01 -4.3439600e-18
+ -1.7480000e+01 -4.6057158e-18
+ -1.7290000e+01 -4.6419909e-18
+ -1.7100000e+01 -4.4417680e-18
+ -1.6910000e+01 -4.0050424e-18
+ -1.6720000e+01 -3.3430568e-18
+ -1.6530000e+01 -2.4781055e-18
+ -1.6340000e+01 -1.4429418e-18
+ -1.6150000e+01 -2.7981391e-19
+ -1.5960000e+01 9.6084977e-19
+ -1.5770000e+01 2.2220983e-18
+ -1.5580000e+01 3.4424597e-18
+ -1.5390000e+01 4.5582867e-18
+ -1.5200000e+01 5.5064387e-18
+ -1.5010000e+01 6.2272659e-18
+ -1.4820000e+01 6.6678244e-18
+ -1.4630000e+01 6.7852074e-18
+ -1.4440000e+01 6.5498212e-18
+ -1.4250000e+01 5.9483931e-18
+ -1.4060000e+01 4.9864585e-18
+ -1.3870000e+01 3.6900596e-18
+ -1.3680000e+01 2.1064010e-18
+ -1.3490000e+01 3.0324576e-19
+ -1.3300000e+01 -1.6330909e-18
+ -1.3110000e+01 -3.6011971e-18
+ -1.2920000e+01 -5.4893960e-18
+ -1.2730000e+01 -7.1817592e-18
+ -1.2540000e+01 -8.5649837e-18
+ -1.2350000e+01 -9.5357145e-18
+ -1.2160000e+01 -1.0007822e-17
+ -1.1970000e+01 -9.9191087e-18
+ -1.1780000e+01 -9.2369373e-18
+ -1.1590000e+01 -7.9623161e-18
+ -1.1400000e+01 -6.1320895e-18
+ -1.1210000e+01 -3.8189941e-18
+ -1.1020000e+01 -1.1294951e-18
+ -1.0830000e+01 1.8005368e-18
+ -1.0640000e+01 4.8121066e-18
+ -1.0450000e+01 7.7310128e-18
+ -1.0260000e+01 1.0376990e-17
+ -1.0070000e+01 1.2573612e-17
+ -9.8800000e+00 1.4158452e-17
+ -9.6900000e+00 1.4993027e-17
+ -9.5000000e+00 1.4972037e-17
+ -9.3100000e+00 1.4031478e-17
+ -9.1200000e+00 1.2155218e-17
+ -8.9300000e+00 9.3796708e-18
+ -8.7400000e+00 5.7962211e-18
+ -8.5500000e+00 1.5510879e-18
+ -8.3600000e+00 -3.1576515e-18
+ -8.1700000e+00 -8.0861012e-18
+ -7.9800000e+00 -1.2953785e-17
+ -7.7900000e+00 -1.7456081e-17
+ -7.6000000e+00 -2.1279676e-17
+ -7.4100000e+00 -2.4120614e-17
+ -7.2200000e+00 -2.5704247e-17
+ -7.0300000e+00 -2.5806098e-17
+ -6.8400000e+00 -2.4272354e-17
+ -6.6500000e+00 -2.1038471e-17
+ -6.4600000e+00 -1.6144202e-17
+ -6.2700000e+00 -9.7433531e-18
+ -6.0800000e+00 -2.1066809e-18
+ -5.8900000e+00 6.3833442e-18
+ -5.7000000e+00 1.5246693e-17
+ -5.5100000e+00 2.3928993e-17
+ -5.3200000e+00 3.1833180e-17
+ -5.1300000e+00 3.8358117e-17
+ -4.9400000e+00 4.2941933e-17
+ -4.7500000e+00 4.5107075e-17
+ -4.5600000e+00 4.4503585e-17
+ -4.3700000e+00 4.0946839e-17
+ -4.1800000e+00 3.4446037e-17
+ -3.9900000e+00 2.5220178e-17
+ -3.8000000e+00 1.3698996e-17
+ -3.6100000e+00 5.0743105e-19
+ -3.4200000e+00 -1.3566495e-17
+ -3.2300000e+00 -2.7619073e-17
+ -3.0400000e+00 -4.0690808e-17
+ -2.8500000e+00 -5.1834784e-17
+ -2.6600000e+00 -6.0188387e-17
+ -2.4700000e+00 -6.5042697e-17
+ -2.2800000e+00 -6.5903659e-17
+ -2.0900000e+00 -6.2539553e-17
+ -1.9000000e+00 -5.5010154e-17
+ -1.7100000e+00 -4.3674315e-17
+ -1.5200000e+00 -2.9174442e-17
+ -1.3300000e+00 -1.2398213e-17
+ -1.1400000e+00 5.5801177e-18
+ -9.5000000e-01 2.3574673e-17
+ -7.6000000e-01 4.0373502e-17
+ -5.7000000e-01 5.4828016e-17
+ -3.8000000e-01 6.5939107e-17
+ -1.9000000e-01 7.2932921e-17
+ 0.0000000e+00 7.5320009e-17
+ 1.9000000e-01 7.2932921e-17
+ 3.8000000e-01 6.5939107e-17
+ 5.7000000e-01 5.4828016e-17
+ 7.6000000e-01 4.0373502e-17
+ 9.5000000e-01 2.3574673e-17
+ 1.1400000e+00 5.5801177e-18
+ 1.3300000e+00 -1.2398213e-17
+ 1.5200000e+00 -2.9174442e-17
+ 1.7100000e+00 -4.3674315e-17
+ 1.9000000e+00 -5.5010154e-17
+ 2.0900000e+00 -6.2539553e-17
+ 2.2800000e+00 -6.5903659e-17
+ 2.4700000e+00 -6.5042697e-17
+ 2.6600000e+00 -6.0188387e-17
+ 2.8500000e+00 -5.1834784e-17
+ 3.0400000e+00 -4.0690808e-17
+ 3.2300000e+00 -2.7619073e-17
+ 3.4200000e+00 -1.3566495e-17
+ 3.6100000e+00 5.0743105e-19
+ 3.8000000e+00 1.3698996e-17
+ 3.9900000e+00 2.5220178e-17
+ 4.1800000e+00 3.4446037e-17
+ 4.3700000e+00 4.0946839e-17
+ 4.5600000e+00 4.4503585e-17
+ 4.7500000e+00 4.5107075e-17
+ 4.9400000e+00 4.2941933e-17
+ 5.1300000e+00 3.8358117e-17
+ 5.3200000e+00 3.1833180e-17
+ 5.5100000e+00 2.3928993e-17
+ 5.7000000e+00 1.5246693e-17
+ 5.8900000e+00 6.3833442e-18
+ 6.0800000e+00 -2.1066809e-18
+ 6.2700000e+00 -9.7433531e-18
+ 6.4600000e+00 -1.6144202e-17
+ 6.6500000e+00 -2.1038471e-17
+ 6.8400000e+00 -2.4272354e-17
+ 7.0300000e+00 -2.5806098e-17
+ 7.2200000e+00 -2.5704247e-17
+ 7.4100000e+00 -2.4120614e-17
+ 7.6000000e+00 -2.1279676e-17
+ 7.7900000e+00 -1.7456081e-17
+ 7.9800000e+00 -1.2953785e-17
+ 8.1700000e+00 -8.0861012e-18
+ 8.3600000e+00 -3.1576515e-18
+ 8.5500000e+00 1.5510879e-18
+ 8.7400000e+00 5.7962211e-18
+ 8.9300000e+00 9.3796708e-18
+ 9.1200000e+00 1.2155218e-17
+ 9.3100000e+00 1.4031478e-17
+ 9.5000000e+00 1.4972037e-17
+ 9.6900000e+00 1.4993027e-17
+ 9.8800000e+00 1.4158452e-17
+ 1.0070000e+01 1.2573612e-17
+ 1.0260000e+01 1.0376990e-17
+ 1.0450000e+01 7.7310128e-18
+ 1.0640000e+01 4.8121066e-18
+ 1.0830000e+01 1.8005368e-18
+ 1.1020000e+01 -1.1294951e-18
+ 1.1210000e+01 -3.8189941e-18
+ 1.1400000e+01 -6.1320895e-18
+ 1.1590000e+01 -7.9623161e-18
+ 1.1780000e+01 -9.2369373e-18
+ 1.1970000e+01 -9.9191087e-18
+ 1.2160000e+01 -1.0007822e-17
+ 1.2350000e+01 -9.5357145e-18
+ 1.2540000e+01 -8.5649837e-18
+ 1.2730000e+01 -7.1817592e-18
+ 1.2920000e+01 -5.4893960e-18
+ 1.3110000e+01 -3.6011971e-18
+ 1.3300000e+01 -1.6330909e-18
+ 1.3490000e+01 3.0324576e-19
+ 1.3680000e+01 2.1064010e-18
+ 1.3870000e+01 3.6900596e-18
+ 1.4060000e+01 4.9864585e-18
+ 1.4250000e+01 5.9483931e-18
+ 1.4440000e+01 6.5498212e-18
+ 1.4630000e+01 6.7852074e-18
+ 1.4820000e+01 6.6678244e-18
+ 1.5010000e+01 6.2272659e-18
+ 1.5200000e+01 5.5064387e-18
+ 1.5390000e+01 4.5582867e-18
+ 1.5580000e+01 3.4424597e-18
+ 1.5770000e+01 2.2220983e-18
+ 1.5960000e+01 9.6084977e-19
+ 1.6150000e+01 -2.7981391e-19
+ 1.6340000e+01 -1.4429418e-18
+ 1.6530000e+01 -2.4781055e-18
+ 1.6720000e+01 -3.3430568e-18
+ 1.6910000e+01 -4.0050424e-18
+ 1.7100000e+01 -4.4417680e-18
+ 1.7290000e+01 -4.6419909e-18
+ 1.7480000e+01 -4.6057158e-18
+ 1.7670000e+01 -4.3439600e-18
+ 1.7860000e+01 -3.8780665e-18
+ 1.8050000e+01 -3.2385574e-18
+ 1.8240000e+01 -2.4635475e-18
+ 1.8430000e+01 -1.5967693e-18
+ 1.8620000e+01 -6.8529585e-19
+ 1.8810000e+01 2.2292353e-19
+ 1.9000000e+01 1.0815724e-18
+ 1.9190000e+01 1.8483709e-18
+ 1.9380000e+01 2.4872234e-18
+ 1.9570000e+01 2.9699741e-18
+ 1.9760000e+01 3.2776657e-18
+ 1.9950000e+01 3.4012322e-18
+ 2.0140000e+01 3.3416006e-18
+ 2.0330000e+01 3.1092217e-18
+ 2.0520000e+01 2.7230887e-18
+ 2.0710000e+01 2.2093363e-18
+ 2.0900000e+01 1.5995338e-18
+ 2.1090000e+01 9.2879539e-19
+ 2.1280000e+01 2.3382952e-19
+ 2.1470000e+01 -4.4896117e-19
+ 2.1660000e+01 -1.0852384e-18
+ 2.1850000e+01 -1.6442742e-18
+ 2.2040000e+01 -2.1002378e-18
+ 2.2230000e+01 -2.4331948e-18
+ 2.2420000e+01 -2.6298080e-18
+ 2.2610000e+01 -2.6837397e-18
+ 2.2800000e+01 -2.5957604e-18
+ 2.2990000e+01 -2.3735726e-18
+ 2.3180000e+01 -2.0313591e-18
+ 2.3370000e+01 -1.5890677e-18
+ 2.3560000e+01 -1.0714449e-18
+ 2.3750000e+01 -5.0684259e-19
+ 2.3940000e+01 7.4170036e-20
+ 2.4130000e+01 6.4034236e-19
+ 2.4320000e+01 1.1613619e-18
+ 2.4510000e+01 1.6094796e-18
+ 2.4700000e+01 1.9610503e-18
+ 2.4890000e+01 2.1978907e-18
+ 2.5080000e+01 2.3083610e-18
+ 2.5270000e+01 2.2880846e-18
+ 2.5460000e+01 2.1402407e-18
+ 2.5650000e+01 1.8753932e-18
+ 2.5840000e+01 1.5108511e-18
+ 2.6030000e+01 1.0695922e-18
+ 2.6220000e+01 5.7881722e-19
+ 2.6410000e+01 6.8230074e-20
+ 2.6600000e+01 -4.3183467e-19
+ 2.6790000e+01 -8.9230633e-19
+ 2.6980000e+01 -1.2871610e-18
+ 2.7170000e+01 -1.5949702e-18
+ 2.7360000e+01 -1.8001133e-18
+ 2.7550000e+01 -1.8935795e-18
+ 2.7740000e+01 -1.8733199e-18
+ 2.7930000e+01 -1.7441435e-18
+ 2.8120000e+01 -1.5171856e-18
+ 2.8310000e+01 -1.2090055e-18
+ 2.8500000e+01 -8.4039233e-19
+ 2.8690000e+01 -4.3496996e-19
+ 2.8880000e+01 -1.7701900e-20
+ 2.9070000e+01 3.8661012e-19
+ 2.9260000e+01 7.5474484e-19
+ 2.9450000e+01 1.0663233e-18
+ 2.9640000e+01 1.3048553e-18
+ 2.9830000e+01 1.4585276e-18
+ 3.0020000e+01 1.5207105e-18
+ 3.0210000e+01 1.4901705e-18
+ 3.0400000e+01 1.3709932e-18
+ 3.0590000e+01 1.1722309e-18
+ 3.0780000e+01 9.0729832e-19
+ 3.0970000e+01 5.9315068e-19
+ 3.1160000e+01 2.4928452e-19
+ 3.1350000e+01 -1.0339019e-19
+ 3.1540000e+01 -4.4375100e-19
+ 3.1730000e+01 -7.5167615e-19
+ 3.1920000e+01 -1.0092053e-18
+ 3.2110000e+01 -1.2015820e-18
+ 3.2300000e+01 -1.3181151e-18
+ 3.2490000e+01 -1.3528028e-18
+ 3.2680000e+01 -1.3046756e-18
+ 3.2870000e+01 -1.1778295e-18
+ 3.3060000e+01 -9.8114099e-19
+ 3.3250000e+01 -7.2767627e-19
+ 3.3440000e+01 -4.3383120e-19
+ 3.3630000e+01 -1.1825671e-19
+ 3.3820000e+01 1.9935691e-19
+ 3.4010000e+01 4.9955265e-19
+ 3.4200000e+01 7.6434700e-19
+ 3.4390000e+01 9.7835102e-19
+ 3.4580000e+01 1.1296951e-18
+ 3.4770000e+01 1.2106966e-18
+ 3.4960000e+01 1.2182298e-18
+ 3.5150000e+01 1.1537822e-18
+ 3.5340000e+01 1.0232050e-18
+ 3.5530000e+01 8.3618806e-19
+ 3.5720000e+01 6.0550846e-19
+ 3.5910000e+01 3.4611692e-19
+ 3.6100000e+01 7.4131800e-20
+ 3.6290000e+01 -1.9418638e-19
+ 3.6480000e+01 -4.4341220e-19
+ 3.6670000e+01 -6.5982566e-19
+ 3.6860000e+01 -8.3212971e-19
+ 3.7050000e+01 -9.5198336e-19
+ 3.7240000e+01 -1.0143341e-18
+ 3.7430000e+01 -1.0175471e-18
+ 3.7620000e+01 -9.6333999e-19
+ 3.7810000e+01 -8.5654316e-19
+ 3.8000000e+01 -7.0470977e-19
+ 3.8190000e+01 -5.1760598e-19
+ 3.8380000e+01 -3.0661268e-19
+ 3.8570000e+01 -8.4071055e-20
+ 3.8760000e+01 1.3739630e-19
+ 3.8950000e+01 3.4555773e-19
+ 3.9140000e+01 5.2920900e-19
+ 3.9330000e+01 6.7874560e-19
+ 3.9520000e+01 7.8664520e-19
+ 3.9710000e+01 8.4783677e-19
+ 3.9900000e+01 8.5993725e-19
+ 4.0090000e+01 8.2334194e-19
+ 4.0280000e+01 7.4116124e-19
+ 4.0470000e+01 6.1900506e-19
+ 4.0660000e+01 4.6462495e-19
+ 4.0850000e+01 2.8743343e-19
+ 4.1040000e+01 9.7928901e-20
+ 4.1230000e+01 -9.2938263e-20
+ 4.1420000e+01 -2.7441787e-19
+ 4.1610000e+01 -4.3658770e-19
+ 4.1800000e+01 -5.7091573e-19
+ 4.1990000e+01 -6.7072072e-19
+ 4.2180000e+01 -7.3150248e-19
+ 4.2370000e+01 -7.5112406e-19
+ 4.2560000e+01 -7.2983910e-19
+ 4.2750000e+01 -6.7017014e-19
+ 4.2940000e+01 -5.7665400e-19
+ 4.3130000e+01 -4.5547972e-19
+ 4.3320000e+01 -3.1405071e-19
+ 4.3510000e+01 -1.6050562e-19
+ 4.3700000e+01 -3.2324241e-21
+ 4.3890000e+01 1.4959319e-19
+ 4.4080000e+01 2.9041926e-19
+ 4.4270000e+01 4.1266036e-19
+ 4.4460000e+01 5.1096036e-19
+ 4.4650000e+01 5.8137141e-19
+ 4.4840000e+01 6.2145114e-19
+ 4.5030000e+01 6.3028490e-19
+ 4.5220000e+01 6.0844227e-19
+ 4.5410000e+01 5.5787848e-19
+ 4.5600000e+01 4.8179060e-19
+ 4.5790000e+01 3.8443735e-19
+ 4.5980000e+01 2.7092971e-19
+ 4.6170000e+01 1.4699816e-19
+ 4.6360000e+01 1.8741548e-20
+ 4.6550000e+01 -1.0763689e-19
+ 4.6740000e+01 -2.2610035e-19
+ 4.6930000e+01 -3.3104772e-19
+ 4.7120000e+01 -4.1757152e-19
+ 4.7310000e+01 -4.8169230e-19
+ 4.7500000e+01 -5.2055661e-19
+ 4.7690000e+01 -5.3258546e-19
+ 4.7880000e+01 -5.1756148e-19
+ 4.8070000e+01 -4.7664596e-19
+ 4.8260000e+01 -4.1232107e-19
+ 4.8450000e+01 -3.2825766e-19
+ 4.8640000e+01 -2.2911445e-19
+ 4.8830000e+01 -1.2028009e-19
+ 4.9020000e+01 -7.5736812e-21
+ 4.9210000e+01 1.0307694e-19
+ 4.9400000e+01 2.0595855e-19
+ 4.9590000e+01 2.9587708e-19
+ 4.9780000e+01 3.6842482e-19
+ 4.9970000e+01 4.2019525e-19
+ 5.0160000e+01 4.4893518e-19
+ 5.0350000e+01 4.5362837e-19
+ 5.0540000e+01 4.3451051e-19
+ 5.0730000e+01 3.9301904e-19
+ 5.0920000e+01 3.3168548e-19
+ 5.1110000e+01 2.5397989e-19
+ 5.1300000e+01 1.6411831e-19
+ 5.1490000e+01 6.6844511e-20
+ 5.1680000e+01 -3.2803460e-20
+ 5.1870000e+01 -1.2973246e-19
+ 5.2060000e+01 -2.1903051e-19
+ 5.2250000e+01 -2.9619265e-19
+ 5.2440000e+01 -3.5732939e-19
+ 5.2630000e+01 -3.9935087e-19
+ 5.2820000e+01 -4.2011991e-19
+ 5.3010000e+01 -4.1856685e-19
+ 5.3200000e+01 -3.9475923e-19
+ 5.3390000e+01 -3.4991992e-19
+ 5.3580000e+01 -2.8638873e-19
+ 5.3770000e+01 -2.0752479e-19
+ 5.3960000e+01 -1.1754995e-19
+ 5.4150000e+01 -2.1337642e-20
+ 5.4340000e+01 7.5844464e-20
+ 5.4530000e+01 1.6862509e-19
+ 5.4720000e+01 2.5183700e-19
+ 5.4910000e+01 3.2082240e-19
+ 5.5100000e+01 3.7171582e-19
+ 5.5290000e+01 4.0168438e-19
+ 5.5480000e+01 4.0910773e-19
+ 5.5670000e+01 3.9368369e-19
+ 5.5860000e+01 3.5645090e-19
+ 5.6050000e+01 2.9972579e-19
+ 5.6240000e+01 2.2695775e-19
+ 5.6430000e+01 1.4251234e-19
+ 5.6620000e+01 5.1397826e-20
+ 5.6810000e+01 -4.1045823e-20
+ 5.7000000e+01 -1.2947405e-19
+ 5.7190000e+01 -2.0885893e-19
+ 5.7380000e+01 -2.7478380e-19
+ 5.7570000e+01 -3.2369325e-19
+ 5.7760000e+01 -3.5308443e-19
+ 5.7950000e+01 -3.6163018e-19
+ 5.8140000e+01 -3.4922950e-19
+ 5.8330000e+01 -3.1698554e-19
+ 5.8520000e+01 -2.6711570e-19
+ 5.8710000e+01 -2.0280178e-19
+ 5.8900000e+01 -1.2799106e-19
+ 5.9090000e+01 -4.7160837e-20
+ 5.9280000e+01 3.4940085e-20
+ 5.9470000e+01 1.1355904e-19
+ 5.9660000e+01 1.8420267e-19
+ 5.9850000e+01 2.4288267e-19
+ 6.0040000e+01 2.8633230e-19
+ 6.0230000e+01 3.1218306e-19
+ 6.0420000e+01 3.1909258e-19
+ 6.0610000e+01 3.0681694e-19
+ 6.0800000e+01 2.7622279e-19
+ 6.0990000e+01 2.2923767e-19
+ 6.1180000e+01 1.6873968e-19
+ 6.1370000e+01 9.8391447e-20
+ 6.1560000e+01 2.2426535e-20
+ 6.1750000e+01 -5.4600362e-20
+ 6.1940000e+01 -1.2808502e-19
+ 6.2130000e+01 -1.9365082e-19
+ 6.2320000e+01 -2.4741406e-19
+ 6.2510000e+01 -2.8622009e-19
+ 6.2700000e+01 -3.0783468e-19
+ 6.2890000e+01 -3.1107739e-19
+ 6.3080000e+01 -2.9588787e-19
+ 6.3270000e+01 -2.6332049e-19
+ 6.3460000e+01 -2.1546775e-19
+ 6.3650000e+01 -1.5531831e-19
+ 6.3840000e+01 -8.6559668e-20
+ 6.4030000e+01 -1.3339828e-20
+ 6.4220000e+01 5.9995411e-20
+ 6.4410000e+01 1.2917147e-19
+ 6.4600000e+01 1.9024237e-19
+ 6.4790000e+01 2.3982016e-19
+ 6.4980000e+01 2.7526252e-19
+ 6.5170000e+01 2.9480847e-19
+ 6.5360000e+01 2.9765587e-19
+ 6.5550000e+01 2.8397889e-19
+ 6.5740000e+01 2.5488724e-19
+ 6.5930000e+01 2.1233279e-19
+ 6.6120000e+01 1.5897190e-19
+ 6.6310000e+01 9.7993561e-20
+ 6.6500000e+01 3.2925061e-20
+ 6.6690000e+01 -3.2573180e-20
+ 6.6880000e+01 -9.4911984e-20
+ 6.7070000e+01 -1.5076633e-19
+ 6.7260000e+01 -1.9724741e-19
+ 6.7450000e+01 -2.3204589e-19
+ 6.7640000e+01 -2.5354037e-19
+ 6.7830000e+01 -2.6086654e-19
+ 6.8020000e+01 -2.5394490e-19
+ 6.8210000e+01 -2.3346641e-19
+ 6.8400000e+01 -2.0083774e-19
+ 6.8590000e+01 -1.5808949e-19
+ 6.8780000e+01 -1.0775270e-19
+ 6.8970000e+01 -5.2710589e-20
+ 6.9160000e+01 3.9661551e-21
+ 6.9350000e+01 5.9191273e-20
+ 6.9540000e+01 1.1003221e-19
+ 6.9730000e+01 1.5386707e-19
+ 6.9920000e+01 1.8852037e-19
+ 7.0110000e+01 2.1236989e-19
+ 7.0300000e+01 2.2441895e-19
+ 7.0490000e+01 2.2433077e-19
+ 7.0680000e+01 2.1242398e-19
+ 7.0870000e+01 1.8963105e-19
+ 7.1060000e+01 1.5742382e-19
+ 7.1250000e+01 1.1771219e-19
+ 7.1440000e+01 7.2723698e-20
+ 7.1630000e+01 2.4872197e-20
+ 7.1820000e+01 -2.3375690e-20
+ 7.2010000e+01 -6.9628508e-20
+ 7.2200000e+01 -1.1168541e-19
+ 7.2390000e+01 -1.4763374e-19
+ 7.2580000e+01 -1.7592525e-19
+ 7.2770000e+01 -1.9542883e-19
+ 7.2960000e+01 -2.0545992e-19
+ 7.3150000e+01 -2.0578729e-19
+ 7.3340000e+01 -1.9661979e-19
+ 7.3530000e+01 -1.7857550e-19
+ 7.3720000e+01 -1.5263679e-19
+ 7.3910000e+01 -1.2009421e-19
+ 7.4100000e+01 -8.2482300e-20
+ 7.4290000e+01 -4.1510205e-20
+ 7.4480000e+01 1.0106876e-21
+ 7.4670000e+01 4.3239679e-20
+ 7.4860000e+01 8.3378339e-20
+ 7.5050000e+01 1.1973878e-19
+ 7.5240000e+01 1.5080677e-19
+ 7.5430000e+01 1.7529785e-19
+ 7.5620000e+01 1.9220480e-19
+ 7.5810000e+01 2.0083475e-19
Added: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/adj_traveltime_filter.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/adj_traveltime_filter.f90 (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/adj_traveltime_filter.f90 2011-04-01 15:57:37 UTC (rev 18161)
@@ -0,0 +1,381 @@
+!!!!! this subroutine is used for preparing adjoint sources in this example ONLY
+!!!!! generally, you should make your own measuremnts and calculate corresponding adjoint sources
+!!!!! FLEXWIN could be of help, or you can modify this code a little bit
+
+program adj_traveltime
+
+implicit none
+
+integer, parameter :: nstep = 31599
+double precision, parameter :: dt = 0.19d0
+double precision, parameter :: C_crit = 0.5d0
+
+double precision, parameter :: length_time_window = 2000.0d0
+integer, parameter :: taper_type = 1
+integer, parameter :: i_skip_window =1
+double precision, parameter :: factor = 5d-2
+
+
+integer, parameter :: nrec = 1
+double precision, parameter :: pi = 3.141592653589793d0
+
+double precision :: data_origin(nstep,nrec),syn_origin(nstep,nrec),adj(nstep,nrec),adj_density(nstep,nrec)
+double precision :: data_picked(nstep,nrec),syn_picked(nstep,nrec),data_filtered(nstep,nrec),syn_filtered(nstep,nrec)
+double precision :: data_temp(nstep),syn_temp(nstep),data_direct(nstep),syn_direct(nstep)
+double precision :: trace_data_max(nrec), trace_syn_max(nrec), t(nstep),data_max,syn_max
+integer :: flag(nrec), irec, itime, length_window, i_start_window, i_end_window
+character(len=150) :: station_name,file_data,file_syn,file_data_direct,file_syn_direct
+character(len=150) :: file_adj,file_adj_density, file_misfit, file_adj_BHX,file_adj_BHZ
+double precision :: c(nstep), data_trace(nstep), syn_trace(nstep)
+integer :: I(nstep)
+double precision :: Norm_data_crit,Norm_syn_crit,AMP,Norm_data,Norm_syn,Norm,Norm_adj_temp
+double precision :: c_current,c_final(100), w(nstep), misfit_traveltime, traveltime_delay
+integer :: I_current, n_current, index_current, l_current, l, lag
+integer :: I_final(100), index_final(100)
+
+double precision, allocatable :: taper(:),corre(:), adj_temp(:),data_trace_temp(:),syn_trace_temp(:)
+
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+integer,parameter :: filter_flag = 1
+integer :: ifreq, nfreq
+real :: F1,F2,D(8),G,DELT
+!real freq_low(2),freq_high(2)
+!data freq_low / 1.0d-4 , 1.0d-4/
+!data freq_high / 5.0d-1 , 5.0d-2/
+real freq_low(1),freq_high(1)
+data freq_low / 6.67d-3 /
+data freq_high / 3.33d-2 /
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+
+nfreq=size(freq_low)
+
+file_misfit = './OUTPUT_FILES/misfit_traveltime_delay'
+open(unit=1111,file=trim(file_misfit),status='unknown')
+
+DELT=dt*1.0d3
+adj=0.0d0
+misfit_traveltime = 0.0d0
+!!!! loading data and synthetics !!!!
+do irec = 1,nrec
+ file_data = './SEM/A7.II.MXZ.sem.ascii'
+ open(unit=1001,file=trim(file_data),status='old',action='read')
+ do itime = 1,nstep
+ !read(1001,*) t(itime),data_origin(itime,irec) ! original
+ read(1001,*) t(itime),data_origin(nstep-itime+1,irec) ! reversed
+ end do
+ close(1001)
+
+ data_temp(1)=0.0
+ data_temp(nstep)=0.0
+ do itime = 2,nstep-1
+ data_temp(itime)=( data_origin(itime+1,irec) - data_origin(itime-1,irec) )/(2*dt)
+ end do
+ data_origin(:,irec)=data_temp
+ trace_data_max(irec)=maxval(abs(data_temp))
+end do
+do itime = 1,nstep
+ t(itime) = t(itime)-t(1)
+end do
+data_max=maxval(trace_data_max)
+
+!!!! taper !!!!
+length_window = floor(length_time_window/dt)+1
+allocate(taper(length_window)); taper(:)=1.0d0
+allocate(adj_temp(length_window)); adj_temp(:)=0.0d0
+allocate(corre(2*length_window-1)); corre(:)=0.0d0
+allocate(data_trace_temp(length_window)); data_trace_temp(:)=0.0d0
+allocate(syn_trace_temp(length_window)); syn_trace_temp(:)=0.0d0
+if (taper_type .eq. 1) then ! cosine taper, otherwise using a constant (1.0) instead
+ do l=1,length_window
+ taper(l) = (1.0-cos(pi*2.0*(l-1)/(length_window-1)))/2.0
+ end do
+end if
+
+!!!! computing adj sources !!!!
+ if (filter_flag .eq. 0) then
+ nfreq=1
+ end if
+do irec = 1,nrec
+ do ifreq = 1,nfreq
+ data_filtered(:,irec)=data_origin(:,irec)
+ syn_filtered(:,irec)= syn_origin(:,irec)
+ !!!!!!!!!!!!!!!!!! to be implemented !!!!!!!!!!!!!!!!!
+ if (filter_flag .eq. 1) then
+ ! THIS SECTION CALCULATES THE FILTER AND MUST BE CALLED BEFORE
+ ! FILTER IS CALLED
+ F1=freq_low(ifreq)
+ F2=freq_high(ifreq)
+ call BNDPAS(F1,F2,DELT,D,G,nstep)
+ ! F1 = LOW FREQUENCY CUTOFF (6 DB DOWN)
+ ! F2 = HIGH FREQUENCY CUTOFF (6 DB DOWN)
+ ! DELT = SAMPLE INTERVAL IN MILLISECONDS
+ ! D = WILL CONTAIN 8 Z DOMAIN COEFICIENTS OF RECURSIVE FILTER
+ ! G = WILL CONTAIN THE GAIN OF THE FILTER,
+ call FILTER(data_filtered(:,irec),nstep,D,G,2)
+ call FILTER(syn_filtered(:,irec),nstep,D,G,2)
+ ! X = DATA VECTOR OF LENGTH N CONTAINING DATA TO BE FILTERED
+ ! D = FILTER COEFFICIENTS CALCULATED BY BNDPAS
+ ! G = FILTER GAIN
+ ! IG = 1 one pass
+ ! ig = 2 two passes
+ end if
+
+
+
+
+
+ c(:)=0.0
+ I(:)=(length_window-1)
+ data_trace=data_filtered(:,irec)
+ syn_trace= syn_filtered(:,irec)
+! if ( ( trace_data_max(irec) > (factor*data_max) ) .and. (trace_syn_max(irec) > (factor*syn_max)) ) then
+ if (1 .eq. 1) then
+ !!!! cross-correlation !!!!
+ do i_start_window = (nstep+1)/2, (nstep+1)/2 !the positive branch
+! do i_start_window = (nstep+1)/2 - length_window, (nstep+1)/2 - length_window !the negative branch
+ i_end_window = i_start_window + length_window - 1
+
+ data_trace_temp=data_trace(i_start_window:i_end_window)
+
+ Norm_data_crit=sqrt(DOT_PRODUCT(data_trace,data_trace))
+ AMP = Norm_data_crit
+
+ Norm_data=sqrt(DOT_PRODUCT( data_trace_temp,data_trace_temp ))
+ Norm = Norm_data*Norm_data
+
+ w(:)=0.0
+ w(i_start_window : i_end_window)=taper
+ end do
+ end if
+ data_trace=data_trace*w
+ data_filtered(:,irec)=data_trace
+ !!!! normal adjoint sources !!!!
+ adj_temp=data_filtered(i_start_window : i_end_window,irec)
+ Norm_adj_temp = - DOT_PRODUCT(adj_temp,adj_temp)*dt ! minus sign comes from integration by part !
+ !!!! choose which one to use !!!!
+ adj_temp=adj_temp/Norm_adj_temp ! normal adjoint sources
+ ! adj_temp=adj_temp/Norm_adj_temp ! ray density map DeltaT=+1
+ ! adj_temp=adj_temp/Norm_adj_temp * ((-1)**irec) ! ray density map DeltaT=+-1
+
+ adj(i_start_window : i_end_window ,irec )=adj_temp
+ traveltime_delay=dt
+ misfit_traveltime = misfit_traveltime + traveltime_delay * traveltime_delay / 2.0
+ write(1111,*) traveltime_delay
+
+ end do !do ifreq=1,nfreq
+
+ !!!! output !!!!
+ file_adj_BHZ = './SEM/A7.II.MXZ.adj'
+ open(unit=1002,file=trim(file_adj_BHZ),status='unknown')
+ do itime = 1,nstep
+ write(1002,*) t(itime), adj(nstep-itime+1,irec)
+! write(1002,*) t(itime), adj(itime,irec)
+ end do
+ close(1002)
+ file_adj_BHZ = './SEM/A7.II.MXZ.adj_reversed'
+ open(unit=1002,file=trim(file_adj_BHZ),status='unknown')
+ do itime = 1,nstep
+! write(1002,*) t(itime), adj(nstep-itime+1,irec)
+ write(1002,*) t(itime), adj(itime,irec)
+ end do
+ close(1002)
+
+
+ file_adj_BHZ = './SEM/A7.II.MXE.adj'
+ open(unit=1002,file=trim(file_adj_BHZ),status='unknown')
+ do itime = 1,nstep
+ write(1002,*) t(itime), 0.0
+ end do
+ close(1002)
+ file_adj_BHZ = './SEM/A7.II.MXN.adj'
+ open(unit=1002,file=trim(file_adj_BHZ),status='unknown')
+ do itime = 1,nstep
+ write(1002,*) t(itime), 0.0
+ end do
+ close(1002)
+
+enddo
+
+
+close(1111)
+
+file_misfit = './OUTPUT_FILES/misfit_traveltime'
+open(unit=1001,file=trim(file_misfit),status='unknown')
+write(1001,*) misfit_traveltime, traveltime_delay
+close(1001)
+
+end program adj_traveltime
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+SUBROUTINE BNDPAS(F1,F2,DELT,D,G,N)
+! RECURSIVE BUTTERWORTH BAND PASS FILTER (KANASEWICH, TIME SERIES
+! ANALYSIS IN GEOPHYSICS, UNIVERSITY OF ALBERTA PRESS, 1975; SHANKS,
+! JOHN L, RECURSION FILTERS FOR DIGITAL PROCESSING, GEOPHYSICS, V32,
+! FILTER. THE FILTER WILL HAVE 8 POLES IN THE S PLANE AND IS
+! APPLIED IN FORWARD AND REVERSE DIRECTIONS SO AS TO HAVE ZERO
+! PHASE SHIFT. THE GAIN AT THE TWO FREQUENCIES SPECIFIED AS
+! CUTOFF FREQUENCIES WILL BE -6DB AND THE ROLLOFF WILL BE ABOUT
+! THE FILTER TO PREVENT ALIASING PROBLEMS.
+ COMPLEX P(4),S(8),Z1,Z2
+ real D(8),XC(3),XD(3),XE(3)
+ double precision :: X(N)
+ DATA ISW/0/,TWOPI/6.2831853/
+! THIS SECTION CALCULATES THE FILTER AND MUST BE CALLED BEFORE
+! FILTER IS CALLED
+
+! F1 = LOW FREQUENCY CUTOFF (6 DB DOWN)
+! F2 = HIGH FREQUENCY CUTOFF (6 DB DOWN)
+! DELT = SAMPLE INTERVAL IN MILLISECONDS
+! D = WILL CONTAIN 8 Z DOMAIN COEFICIENTS OF RECURSIVE FILTER
+! G = WILL CONTAIN THE GAIN OF THE FILTER,
+
+ DT=DELT/1000.0
+ TDT=2.0/DT
+ FDT=4.0/DT
+ ISW=1
+ P(1)=CMPLX(-.3826834,.9238795)
+ P(2)=CMPLX(-.3826834,-.9238795)
+ P(3)=CMPLX(-.9238795,.3826834)
+ P(4)=CMPLX(-.9238795,-.3826834)
+ W1=TWOPI*F1
+ W2=TWOPI*F2
+ W1=TDT*TAN(W1/TDT)
+ W2=TDT*TAN(W2/TDT)
+ HWID=(W2-W1)/2.0
+ WW=W1*W2
+ DO 19 I=1,4
+ Z1=P(I)*HWID
+ Z2=Z1*Z1-WW
+ Z2=CSQRT(Z2)
+ S(I)=Z1+Z2
+ 19 S(I+4)=Z1-Z2
+ G=.5/HWID
+ G=G*G
+ G=G*G
+ DO 29 I=1,7,2
+ B=-2.0*REAL(S(I))
+ Z1=S(I)*S(I+1)
+ C=REAL(Z1)
+ A=TDT+B+C/TDT
+ G=G*A
+ D(I)=(C*DT-FDT)/A
+ 29 D(I+1)=(A-2.0*B)/A
+ G=G*G
+ 5 FORMAT ('-FILTER GAIN IS ', 9E12.6)
+ RETURN
+
+ ENTRY FILTER(X,N,D,G,IG)
+
+! X = DATA VECTOR OF LENGTH N CONTAINING DATA TO BE FILTERED
+! D = FILTER COEFFICIENTS CALCULATED BY BNDPAS
+! G = FILTER GAIN
+! IG = 1 one pass
+! ig = 2 two passes
+
+ IF (ISW.EQ.1) GO TO 31
+ WRITE (6,6)
+ 6 FORMAT ('1BNDPAS MUST BE CALLED BEFORE FILTER')
+ return
+
+! APPLY FILTER IN FORWARD DIRECTION
+
+ 31 XM2=X(1)
+ XM1=X(2)
+ XM=X(3)
+ XC(1)=XM2
+ XC(2)=XM1-D(1)*XC(1)
+ XC(3)=XM-XM2-D(1)*XC(2)-D(2)*XC(1)
+ XD(1)=XC(1)
+ XD(2)=XC(2)-D(3)*XD(1)
+ XD(3)=XC(3)-XC(1)-D(3)*XD(2)-D(4)*XD(1)
+ XE(1)=XD(1)
+ XE(2)=XD(2)-D(5)*XE(1)
+ XE(3)=XD(3)-XD(1)-D(5)*XE(2)-D(6)*XE(1)
+ X(1)=XE(1)
+ X(2)=XE(2)-D(7)*X(1)
+ X(3)=XE(3)-XE(1)-D(7)*X(2)-D(8)*X(1)
+ DO 39 I=4,N
+ XM2=XM1
+ XM1=XM
+ XM=X(I)
+ K=I-((I-1)/3)*3
+ GO TO (34,35,36),K
+ 34 M=1
+ M1=3
+ M2=2
+ GO TO 37
+ 35 M=2
+ M1=1
+ M2=3
+ GO TO 37
+ 36 M=3
+ M1=2
+ M2=1
+ 37 XC(M)=XM-XM2-D(1)*XC(M1)-D(2)*XC(M2)
+ XD(M)=XC(M)-XC(M2)-D(3)*XD(M1)-D(4)*XD(M2)
+ XE(M)=XD(M)-XD(M2)-D(5)*XE(M1)-D(6)*XE(M2)
+ 39 X(I)=XE(M)-XE(M2)-D(7)*X(I-1)-D(8)*X(I-2)
+!
+!
+ if(ig.eq.1) goto 3333
+ XM2=X(N)
+ XM1=X(N-1)
+ XM=X(N-2)
+ XC(1)=XM2
+ XC(2)=XM1-D(1)*XC(1)
+ XC(3)=XM-XM2-D(1)*XC(2)-D(2)*XC(1)
+ XD(1)=XC(1)
+ XD(2)=XC(2)-D(3)*XD(1)
+ XD(3)=XC(3)-XC(1)-D(3)*XD(2)-D(4)*XD(1)
+ XE(1)=XD(1)
+ XE(2)=XD(2)-D(5)*XE(1)
+ XE(3)=XD(3)-XD(1)-D(5)*XE(2)-D(6)*XE(1)
+ X(N)=XE(1)
+ X(N-1)=XE(2)-D(7)*X(1)
+ X(N-2)=XE(3)-XE(1)-D(7)*X(2)-D(8)*X(1)
+ DO 49 I=4,N
+ XM2=XM1
+ XM1=XM
+ J=N-I+1
+ XM=X(J)
+ K=I-((I-1)/3)*3
+ GO TO (44,45,46),K
+ 44 M=1
+ M1=3
+ M2=2
+ GO TO 47
+ 45 M=2
+ M1=1
+ M2=3
+ GO TO 47
+ 46 M=3
+ M1=2
+ M2=1
+ 47 XC(M)=XM-XM2-D(1)*XC(M1)-D(2)*XC(M2)
+ XD(M)=XC(M)-XC(M2)-D(3)*XD(M1)-D(4)*XD(M2)
+ XE(M)=XD(M)-XD(M2)-D(5)*XE(M1)-D(6)*XE(M2)
+ 49 X(J)=XE(M)-XE(M2)-D(7)*X(J+1)-D(8)*X(J+2)
+ 3333 continue
+ if (ig.eq.1) then
+ gg=sqrt(g) ! if only pass once, modify gain
+ else
+ gg=g
+ endif
+ DO 59 I=1,N
+ 59 X(I)=X(I)/gg
+ RETURN
+END
+
Added: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/irec_master_noise
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/irec_master_noise (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/irec_master_noise 2011-04-01 15:57:37 UTC (rev 18161)
@@ -0,0 +1 @@
+1
Added: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/nu_master
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/nu_master (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_global/nu_master 2011-04-01 15:57:37 UTC (rev 18161)
@@ -0,0 +1,3 @@
+0
+0
+1
Added: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/CMTSOLUTION_NOISE
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/CMTSOLUTION_NOISE (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/CMTSOLUTION_NOISE 2011-04-01 15:57:37 UTC (rev 18161)
@@ -0,0 +1,13 @@
+PDE 1994 6 9 0 33 16.40 -13.8300 -67.5600 637.0 6.9 6.8 NORTHERN BOLIVIA
+event name: 060994A
+time shift: 29.0000
+half duration: 30.0000
+latitude: 0.0000
+longitude: 0.0000
+depth: 0.0000
+Mrr: 0.00e+27
+Mtt: 0.00e+27
+Mpp: 0.00e+27
+Mrt: 0.00e+28
+Mrp: 0.00e+26
+Mtp: 0.00e+27
Added: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/NOISE123.submit_atten
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/NOISE123.submit_atten (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/NOISE123.submit_atten 2011-04-01 15:57:37 UTC (rev 18161)
@@ -0,0 +1,136 @@
+#!/bin/bash
+#PBS -N NOISE_test_regional
+#PBS -q tromp
+#PBS -o zzz_test_regional/job_info/output
+#PBS -e zzz_test_regional/job_info/error
+#PBS -l nodes=50:ppn=8,walltime=100:00:00
+
+DIR_RESULT="zzz_test_regional"
+NSTEP=799
+
+cd $PBS_O_WORKDIR
+######################################## pre-simulation ###########################################
+## gather simulation infomation
+NPROC_XI=`grep NPROC_XI DATA/Par_file | cut -d = -f 2 `
+NPROC_ETA=`grep NPROC_ETA DATA/Par_file | cut -d = -f 2`
+NCHUNKS=`grep NCHUNKS DATA/Par_file | cut -d = -f 2 `
+numnodes=$(( $NCHUNKS * $NPROC_XI * $NPROC_ETA ))
+cat $PBS_NODEFILE > $DIR_RESULT/job_info/nodes
+echo $PBS_JOBID > $DIR_RESULT/job_info/id
+
+## check number of procs
+NPROC_1=`wc -l < $DIR_RESULT/job_info/nodes` # number of procs assigned to the job
+NPROC_2=$numnodes # number of procs needed in SPECFEM
+if [ $NPROC_1 -ne $NPROC_2 ]; then
+ echo "inconsistent number of procs, check your job submission script & DATA/Par_file" 1>&2
+ echo "NPROC from job submission: $NPROC_1" 1>&2
+ echo "NPROC from DATA/Par_file: $NPROC_2" 1>&2
+ exit
+fi
+
+## clean local nodes
+LOCAL_PATH=`grep LOCAL_PATH DATA/Par_file | cut -d = -f 2 | sed 's/ //g'`
+d=`date`
+echo "Start cleaning local disk on nodes, $d"
+pbsdsh $PBS_O_WORKDIR/NOISE_clean_create $LOCAL_PATH
+d=`date`
+echo "Finish cleaning local disk on nodes, $d"
+
+## clean OUTPUT_FILES
+rm -rf OUTPUT_FILES
+mkdir -p OUTPUT_FILES
+######################################## simulation ###############################################
+################################################################### step 1: ensemble forward source
+cp DATA/Par_file_NOISE_1_attenuation DATA/Par_file
+mpiexec $PBS_O_WORKDIR/xmeshfem3D_attenuation
+mpiexec $PBS_O_WORKDIR/xspecfem3D_attenuation
+echo "2" > temp_input_movie; echo "1" >> temp_input_movie; echo "$NSTEP" >> temp_input_movie
+$PBS_O_WORKDIR/xcreate_movie_AVS_DX_attenuation < temp_input_movie
+# backup step 1 results (with attenuation)
+DIR_TEMP=$DIR_RESULT/step_1
+rm -rf $DIR_TEMP
+mkdir -p $DIR_TEMP
+mv OUTPUT_FILES/* $DIR_TEMP/
+cp DATA/Par_file $DIR_TEMP/
+cp DATA/CMTSOLUTION $DIR_TEMP/
+cp DATA/STATIONS $DIR_TEMP/
+################################################################ step 2: ensemble forward wavefield
+cp DATA/Par_file_NOISE_2_attenuation DATA/Par_file
+mpiexec $PBS_O_WORKDIR/xmeshfem3D_attenuation
+mpiexec $PBS_O_WORKDIR/xspecfem3D_attenuation
+echo "2" > temp_input_movie; echo "1" >> temp_input_movie; echo "$NSTEP" >> temp_input_movie
+$PBS_O_WORKDIR/xcreate_movie_AVS_DX_attenuation < temp_input_movie
+# backup step 2 results (with attenuation)
+DIR_TEMP=$DIR_RESULT/step_2
+rm -rf $DIR_TEMP
+mkdir -p $DIR_TEMP
+mv OUTPUT_FILES/* $DIR_TEMP/
+cp DATA/Par_file $DIR_TEMP/
+cp DATA/CMTSOLUTION $DIR_TEMP/
+cp DATA/STATIONS $DIR_TEMP/
+############################################################## prepare adjoint sources ############
+cp $DIR_RESULT/step_2/A7.II.MXZ.sem.ascii ./SEM/
+./NOISE_adj
+################################################################### step 1: ensemble forward source
+cp DATA/Par_file_NOISE_1_noattenuation DATA/Par_file
+mpiexec $PBS_O_WORKDIR/xmeshfem3D_noattenuation
+mpiexec $PBS_O_WORKDIR/xspecfem3D_noattenuation
+echo "2" > temp_input_movie; echo "1" >> temp_input_movie; echo "$NSTEP" >> temp_input_movie
+$PBS_O_WORKDIR/xcreate_movie_AVS_DX_noattenuation < temp_input_movie
+# backup step 1 results (without attenuation)
+DIR_TEMP=$DIR_RESULT/step_1_no_attenuation
+rm -rf $DIR_TEMP
+mkdir -p $DIR_TEMP
+mv OUTPUT_FILES/* $DIR_TEMP/
+cp DATA/Par_file $DIR_TEMP/
+cp DATA/CMTSOLUTION $DIR_TEMP/
+cp DATA/STATIONS $DIR_TEMP/
+################################################################ step 2: ensemble forward wavefield
+cp DATA/Par_file_NOISE_2_noattenuation DATA/Par_file
+mpiexec $PBS_O_WORKDIR/xmeshfem3D_noattenuation
+mpiexec $PBS_O_WORKDIR/xspecfem3D_noattenuation
+echo "2" > temp_input_movie; echo "1" >> temp_input_movie; echo "$NSTEP" >> temp_input_movie
+$PBS_O_WORKDIR/xcreate_movie_AVS_DX_noattenuation < temp_input_movie
+# backup step 2 results (without attenuation)
+DIR_TEMP=$DIR_RESULT/step_2_no_attenuation
+rm -rf $DIR_TEMP
+mkdir -p $DIR_TEMP
+mv OUTPUT_FILES/* $DIR_TEMP/
+cp DATA/Par_file $DIR_TEMP/
+cp DATA/CMTSOLUTION $DIR_TEMP/
+cp DATA/STATIONS $DIR_TEMP/
+############################################################## step 3: ensemble kernels calculation
+cp DATA/Par_file_NOISE_3_noattenuation DATA/Par_file
+mpiexec $PBS_O_WORKDIR/xmeshfem3D_noattenuation
+mpiexec $PBS_O_WORKDIR/xspecfem3D_noattenuation
+echo "2" > temp_input_movie; echo "1" >> temp_input_movie; echo "$NSTEP" >> temp_input_movie
+$PBS_O_WORKDIR/xcreate_movie_AVS_DX_noattenuation < temp_input_movie
+
+## collect outputs (e.g., kernels, meshes) from local nodes to global disk
+d=`date`
+echo "Start collecting outputs from nodes, $d"
+pbsdsh $PBS_O_WORKDIR/NOISE_collect $LOCAL_PATH $PBS_O_WORKDIR/OUTPUT_FILES/
+d=`date`
+echo "Finish collecting outputs from nodes, $d"
+## combine_vol_data
+echo "0" > temp_input_movie
+for (( iproc=1;iproc<=$numnodes-1;iproc++ ))
+do
+ echo "$iproc" >> temp_input_movie
+done
+$PBS_O_WORKDIR/xcombine_vol_data temp_input_movie alpha_kernel $PBS_O_WORKDIR/OUTPUT_FILES $PBS_O_WORKDIR/OUTPUT_FILES $PBS_O_WORKDIR/OUTPUT_FILES 0 1
+rm temp_input_movie
+
+## backup step 3 results (without attenuation)
+DIR_TEMP=$DIR_RESULT/step_3_no_attenuation
+rm -rf $DIR_TEMP
+mkdir -p $DIR_TEMP
+mv OUTPUT_FILES/* $DIR_TEMP/
+cp DATA/Par_file $DIR_TEMP/
+cp DATA/CMTSOLUTION $DIR_TEMP/
+cp DATA/STATIONS $DIR_TEMP/
+
+## clean local nodes
+pbsdsh $PBS_O_WORKDIR/NOISE_clean_create $LOCAL_PATH
+d=`date`
+echo "Finish cleaning local disk on nodes, $d"
Property changes on: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/NOISE123.submit_atten
___________________________________________________________________
Name: svn:executable
+ *
Added: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/NOISE_clean_create
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/NOISE_clean_create (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/NOISE_clean_create 2011-04-01 15:57:37 UTC (rev 18161)
@@ -0,0 +1,11 @@
+#!/bin/bash -eu
+
+DIR_LOCAL=$1
+iproc=$PBS_VNODENUM
+
+##### other files can be collected in the same way
+if [ $(($iproc % 8)) = 0 ]; then
+ rm -rf $DIR_LOCAL
+ mkdir -p $DIR_LOCAL
+fi
+
Property changes on: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/NOISE_clean_create
___________________________________________________________________
Name: svn:executable
+ *
Added: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/NOISE_collect
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/NOISE_collect (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/NOISE_collect 2011-04-01 15:57:37 UTC (rev 18161)
@@ -0,0 +1,11 @@
+#!/bin/bash -eu
+
+DIR_LOCAL=$1
+DIR_GLOBAL=$2
+iproc=`printf "%06d" $PBS_VNODENUM`
+
+##### other files can be collected in the same way
+mv $DIR_LOCAL/*$iproc*reg1*kernel*.bin $DIR_GLOBAL/
+mv $DIR_LOCAL/*$iproc*reg1_*array_dims.txt $DIR_GLOBAL/
+mv $DIR_LOCAL/*$iproc*reg1_solver_data*.bin $DIR_GLOBAL/
+
Property changes on: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/NOISE_collect
___________________________________________________________________
Name: svn:executable
+ *
Added: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_1_attenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_1_attenuation (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_1_attenuation 2011-04-01 15:57:37 UTC (rev 18161)
@@ -0,0 +1,120 @@
+
+# forward or adjoint simulation
+SIMULATION_TYPE = 1
+NOISE_TOMOGRAPHY = 1 # flag of noise tomography, three steps (1,2,3). If earthquake tomography, set it to 0.
+SAVE_FORWARD = .true. # save last frame of forward simulation or not
+
+# number of chunks (1,2,3 or 6)
+NCHUNKS = 1
+
+# angular width of the first chunk (not used if full sphere with six chunks)
+ANGULAR_WIDTH_XI_IN_DEGREES = 90.d0 # angular size of a chunk
+ANGULAR_WIDTH_ETA_IN_DEGREES = 90.d0
+CENTER_LATITUDE_IN_DEGREES = 0.d0
+CENTER_LONGITUDE_IN_DEGREES = 15.d0
+GAMMA_ROTATION_AZIMUTH = 0.d0
+
+# number of elements at the surface along the two sides of the first chunk
+# (must be multiple of 16 and 8 * multiple of NPROC below)
+NEX_XI = 160
+NEX_ETA = 160
+
+# number of MPI processors along the two sides of the first chunk
+NPROC_XI = 20
+NPROC_ETA = 20
+
+# 1D models with real structure:
+# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
+#
+# 1D models with only one fictitious averaged crustal layer:
+# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135_onecrust
+#
+# fully 3D models:
+# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
+# s20rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
+# s29ea, s29ea,sea99_jp3d1994,sea99,jp3d1994,heterogen
+#
+# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
+# to take the 1D crustal model from the
+# associated reference model rather than the default 3D crustal model
+# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
+MODEL = 1D_isotropic_prem
+
+# parameters describing the Earth model
+OCEANS = .false.
+ELLIPTICITY = .false.
+TOPOGRAPHY = .false.
+GRAVITY = .false.
+ROTATION = .false.
+ATTENUATION = .true.
+
+# absorbing boundary conditions for a regional simulation
+ABSORBING_CONDITIONS = .true.
+
+# record length in minutes
+RECORD_LENGTH_IN_MINUTES = 1.0d0
+
+# save AVS or OpenDX movies
+#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
+#MOVIE_COARSE does not work with create_movie_AVS_DX
+MOVIE_SURFACE = .true.
+MOVIE_VOLUME = .false.
+MOVIE_COARSE = .false.
+NTSTEP_BETWEEN_FRAMES = 50
+HDUR_MOVIE = 0.d0
+
+# save movie in volume. Will save element if center of element is in prescribed volume
+# top/bottom: depth in KM, use MOVIE_TOP = -100 to make sure the surface is stored.
+# west/east: longitude, degrees East [-180/180] top/bottom: latitute, degrees North [-90/90]
+# start/stop: frames will be stored at MOVIE_START + i*NSTEP_BETWEEN_FRAMES, where i=(0,1,2..) and iNSTEP_BETWEEN_FRAMES <= MOVIE_STOP
+# movie_volume_type: 1=strain, 2=time integral of strain, 3=\mu*time integral of strain
+# type 4 saves the trace and deviatoric stress in the whole volume, 5=displacement, 6=velocity
+MOVIE_VOLUME_TYPE = 2
+MOVIE_TOP_KM = -100.0
+MOVIE_BOTTOM_KM = 1000.0
+MOVIE_WEST_DEG = -90.0
+MOVIE_EAST_DEG = 90.0
+MOVIE_NORTH_DEG = 90.0
+MOVIE_SOUTH_DEG = -90.0
+MOVIE_START = 0
+MOVIE_STOP = 40000
+
+# save mesh files to check the mesh
+SAVE_MESH_FILES = .false.
+
+# restart files (number of runs can be 1, 2 or 3, choose 1 for no restart files)
+NUMBER_OF_RUNS = 1
+NUMBER_OF_THIS_RUN = 1
+
+# path to store the local database files on each node
+LOCAL_PATH = /scratch/yangl/DATABASES_MPI
+#/scratch/yangl/DATABASES_MPI
+
+# interval at which we output time step info and max of norm of displacement
+NTSTEP_BETWEEN_OUTPUT_INFO = 10000
+
+# interval in time steps for temporary writing of seismograms
+NTSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000
+NTSTEP_BETWEEN_READ_ADJSRC = 1000
+
+# output format for the seismograms (one can use either or all of the three formats)
+OUTPUT_SEISMOS_ASCII_TEXT = .true.
+OUTPUT_SEISMOS_SAC_ALPHANUM = .false.
+OUTPUT_SEISMOS_SAC_BINARY = .false.
+
+# rotate seismograms to Radial-Transverse-Z or use default North-East-Z reference frame
+ROTATE_SEISMOGRAMS_RT = .false.
+
+# decide if master process writes all the seismograms or if all processes do it in parallel
+WRITE_SEISMOGRAMS_BY_MASTER = .true.
+
+# save all seismograms in one large combined file instead of one file per seismogram
+# to avoid overloading shared non-local file systems such as GPFS for instance
+SAVE_ALL_SEISMOS_IN_ONE_FILE = .false.
+USE_BINARY_FOR_LARGE_FILE = .false.
+
+# flag to impose receivers at the surface or allow them to be buried
+RECEIVERS_CAN_BE_BURIED = .true.
+
+# print source time function
+PRINT_SOURCE_TIME_FUNCTION = .true.
Added: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_1_noattenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_1_noattenuation (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_1_noattenuation 2011-04-01 15:57:37 UTC (rev 18161)
@@ -0,0 +1,120 @@
+
+# forward or adjoint simulation
+SIMULATION_TYPE = 1
+NOISE_TOMOGRAPHY = 1 # flag of noise tomography, three steps (1,2,3). If earthquake tomography, set it to 0.
+SAVE_FORWARD = .true. # save last frame of forward simulation or not
+
+# number of chunks (1,2,3 or 6)
+NCHUNKS = 1
+
+# angular width of the first chunk (not used if full sphere with six chunks)
+ANGULAR_WIDTH_XI_IN_DEGREES = 90.d0 # angular size of a chunk
+ANGULAR_WIDTH_ETA_IN_DEGREES = 90.d0
+CENTER_LATITUDE_IN_DEGREES = 0.d0
+CENTER_LONGITUDE_IN_DEGREES = 15.d0
+GAMMA_ROTATION_AZIMUTH = 0.d0
+
+# number of elements at the surface along the two sides of the first chunk
+# (must be multiple of 16 and 8 * multiple of NPROC below)
+NEX_XI = 160
+NEX_ETA = 160
+
+# number of MPI processors along the two sides of the first chunk
+NPROC_XI = 20
+NPROC_ETA = 20
+
+# 1D models with real structure:
+# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
+#
+# 1D models with only one fictitious averaged crustal layer:
+# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135_onecrust
+#
+# fully 3D models:
+# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
+# s20rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
+# s29ea, s29ea,sea99_jp3d1994,sea99,jp3d1994,heterogen
+#
+# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
+# to take the 1D crustal model from the
+# associated reference model rather than the default 3D crustal model
+# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
+MODEL = 1D_isotropic_prem
+
+# parameters describing the Earth model
+OCEANS = .false.
+ELLIPTICITY = .false.
+TOPOGRAPHY = .false.
+GRAVITY = .false.
+ROTATION = .false.
+ATTENUATION = .false.
+
+# absorbing boundary conditions for a regional simulation
+ABSORBING_CONDITIONS = .true.
+
+# record length in minutes
+RECORD_LENGTH_IN_MINUTES = 1.0d0
+
+# save AVS or OpenDX movies
+#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
+#MOVIE_COARSE does not work with create_movie_AVS_DX
+MOVIE_SURFACE = .true.
+MOVIE_VOLUME = .false.
+MOVIE_COARSE = .false.
+NTSTEP_BETWEEN_FRAMES = 50
+HDUR_MOVIE = 0.d0
+
+# save movie in volume. Will save element if center of element is in prescribed volume
+# top/bottom: depth in KM, use MOVIE_TOP = -100 to make sure the surface is stored.
+# west/east: longitude, degrees East [-180/180] top/bottom: latitute, degrees North [-90/90]
+# start/stop: frames will be stored at MOVIE_START + i*NSTEP_BETWEEN_FRAMES, where i=(0,1,2..) and iNSTEP_BETWEEN_FRAMES <= MOVIE_STOP
+# movie_volume_type: 1=strain, 2=time integral of strain, 3=\mu*time integral of strain
+# type 4 saves the trace and deviatoric stress in the whole volume, 5=displacement, 6=velocity
+MOVIE_VOLUME_TYPE = 2
+MOVIE_TOP_KM = -100.0
+MOVIE_BOTTOM_KM = 1000.0
+MOVIE_WEST_DEG = -90.0
+MOVIE_EAST_DEG = 90.0
+MOVIE_NORTH_DEG = 90.0
+MOVIE_SOUTH_DEG = -90.0
+MOVIE_START = 0
+MOVIE_STOP = 40000
+
+# save mesh files to check the mesh
+SAVE_MESH_FILES = .false.
+
+# restart files (number of runs can be 1, 2 or 3, choose 1 for no restart files)
+NUMBER_OF_RUNS = 1
+NUMBER_OF_THIS_RUN = 1
+
+# path to store the local database files on each node
+LOCAL_PATH = /scratch/yangl/DATABASES_MPI
+#/scratch/yangl/DATABASES_MPI
+
+# interval at which we output time step info and max of norm of displacement
+NTSTEP_BETWEEN_OUTPUT_INFO = 10000
+
+# interval in time steps for temporary writing of seismograms
+NTSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000
+NTSTEP_BETWEEN_READ_ADJSRC = 1000
+
+# output format for the seismograms (one can use either or all of the three formats)
+OUTPUT_SEISMOS_ASCII_TEXT = .true.
+OUTPUT_SEISMOS_SAC_ALPHANUM = .false.
+OUTPUT_SEISMOS_SAC_BINARY = .false.
+
+# rotate seismograms to Radial-Transverse-Z or use default North-East-Z reference frame
+ROTATE_SEISMOGRAMS_RT = .false.
+
+# decide if master process writes all the seismograms or if all processes do it in parallel
+WRITE_SEISMOGRAMS_BY_MASTER = .true.
+
+# save all seismograms in one large combined file instead of one file per seismogram
+# to avoid overloading shared non-local file systems such as GPFS for instance
+SAVE_ALL_SEISMOS_IN_ONE_FILE = .false.
+USE_BINARY_FOR_LARGE_FILE = .false.
+
+# flag to impose receivers at the surface or allow them to be buried
+RECEIVERS_CAN_BE_BURIED = .true.
+
+# print source time function
+PRINT_SOURCE_TIME_FUNCTION = .true.
Added: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_2_attenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_2_attenuation (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_2_attenuation 2011-04-01 15:57:37 UTC (rev 18161)
@@ -0,0 +1,120 @@
+
+# forward or adjoint simulation
+SIMULATION_TYPE = 1
+NOISE_TOMOGRAPHY = 2 # flag of noise tomography, three steps (1,2,3). If earthquake tomography, set it to 0.
+SAVE_FORWARD = .true. # save last frame of forward simulation or not
+
+# number of chunks (1,2,3 or 6)
+NCHUNKS = 1
+
+# angular width of the first chunk (not used if full sphere with six chunks)
+ANGULAR_WIDTH_XI_IN_DEGREES = 90.d0 # angular size of a chunk
+ANGULAR_WIDTH_ETA_IN_DEGREES = 90.d0
+CENTER_LATITUDE_IN_DEGREES = 0.d0
+CENTER_LONGITUDE_IN_DEGREES = 15.d0
+GAMMA_ROTATION_AZIMUTH = 0.d0
+
+# number of elements at the surface along the two sides of the first chunk
+# (must be multiple of 16 and 8 * multiple of NPROC below)
+NEX_XI = 160
+NEX_ETA = 160
+
+# number of MPI processors along the two sides of the first chunk
+NPROC_XI = 20
+NPROC_ETA = 20
+
+# 1D models with real structure:
+# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
+#
+# 1D models with only one fictitious averaged crustal layer:
+# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135_onecrust
+#
+# fully 3D models:
+# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
+# s20rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
+# s29ea, s29ea,sea99_jp3d1994,sea99,jp3d1994,heterogen
+#
+# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
+# to take the 1D crustal model from the
+# associated reference model rather than the default 3D crustal model
+# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
+MODEL = 1D_isotropic_prem
+
+# parameters describing the Earth model
+OCEANS = .false.
+ELLIPTICITY = .false.
+TOPOGRAPHY = .false.
+GRAVITY = .false.
+ROTATION = .false.
+ATTENUATION = .true.
+
+# absorbing boundary conditions for a regional simulation
+ABSORBING_CONDITIONS = .true.
+
+# record length in minutes
+RECORD_LENGTH_IN_MINUTES = 1.0d0
+
+# save AVS or OpenDX movies
+#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
+#MOVIE_COARSE does not work with create_movie_AVS_DX
+MOVIE_SURFACE = .true.
+MOVIE_VOLUME = .false.
+MOVIE_COARSE = .false.
+NTSTEP_BETWEEN_FRAMES = 50
+HDUR_MOVIE = 0.d0
+
+# save movie in volume. Will save element if center of element is in prescribed volume
+# top/bottom: depth in KM, use MOVIE_TOP = -100 to make sure the surface is stored.
+# west/east: longitude, degrees East [-180/180] top/bottom: latitute, degrees North [-90/90]
+# start/stop: frames will be stored at MOVIE_START + i*NSTEP_BETWEEN_FRAMES, where i=(0,1,2..) and iNSTEP_BETWEEN_FRAMES <= MOVIE_STOP
+# movie_volume_type: 1=strain, 2=time integral of strain, 3=\mu*time integral of strain
+# type 4 saves the trace and deviatoric stress in the whole volume, 5=displacement, 6=velocity
+MOVIE_VOLUME_TYPE = 2
+MOVIE_TOP_KM = -100.0
+MOVIE_BOTTOM_KM = 1000.0
+MOVIE_WEST_DEG = -90.0
+MOVIE_EAST_DEG = 90.0
+MOVIE_NORTH_DEG = 90.0
+MOVIE_SOUTH_DEG = -90.0
+MOVIE_START = 0
+MOVIE_STOP = 40000
+
+# save mesh files to check the mesh
+SAVE_MESH_FILES = .false.
+
+# restart files (number of runs can be 1, 2 or 3, choose 1 for no restart files)
+NUMBER_OF_RUNS = 1
+NUMBER_OF_THIS_RUN = 1
+
+# path to store the local database files on each node
+LOCAL_PATH = /scratch/yangl/DATABASES_MPI
+#/scratch/yangl/DATABASES_MPI
+
+# interval at which we output time step info and max of norm of displacement
+NTSTEP_BETWEEN_OUTPUT_INFO = 10000
+
+# interval in time steps for temporary writing of seismograms
+NTSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000
+NTSTEP_BETWEEN_READ_ADJSRC = 1000
+
+# output format for the seismograms (one can use either or all of the three formats)
+OUTPUT_SEISMOS_ASCII_TEXT = .true.
+OUTPUT_SEISMOS_SAC_ALPHANUM = .false.
+OUTPUT_SEISMOS_SAC_BINARY = .false.
+
+# rotate seismograms to Radial-Transverse-Z or use default North-East-Z reference frame
+ROTATE_SEISMOGRAMS_RT = .false.
+
+# decide if master process writes all the seismograms or if all processes do it in parallel
+WRITE_SEISMOGRAMS_BY_MASTER = .true.
+
+# save all seismograms in one large combined file instead of one file per seismogram
+# to avoid overloading shared non-local file systems such as GPFS for instance
+SAVE_ALL_SEISMOS_IN_ONE_FILE = .false.
+USE_BINARY_FOR_LARGE_FILE = .false.
+
+# flag to impose receivers at the surface or allow them to be buried
+RECEIVERS_CAN_BE_BURIED = .true.
+
+# print source time function
+PRINT_SOURCE_TIME_FUNCTION = .true.
Added: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_2_noattenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_2_noattenuation (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_2_noattenuation 2011-04-01 15:57:37 UTC (rev 18161)
@@ -0,0 +1,120 @@
+
+# forward or adjoint simulation
+SIMULATION_TYPE = 1
+NOISE_TOMOGRAPHY = 2 # flag of noise tomography, three steps (1,2,3). If earthquake tomography, set it to 0.
+SAVE_FORWARD = .true. # save last frame of forward simulation or not
+
+# number of chunks (1,2,3 or 6)
+NCHUNKS = 1
+
+# angular width of the first chunk (not used if full sphere with six chunks)
+ANGULAR_WIDTH_XI_IN_DEGREES = 90.d0 # angular size of a chunk
+ANGULAR_WIDTH_ETA_IN_DEGREES = 90.d0
+CENTER_LATITUDE_IN_DEGREES = 0.d0
+CENTER_LONGITUDE_IN_DEGREES = 15.d0
+GAMMA_ROTATION_AZIMUTH = 0.d0
+
+# number of elements at the surface along the two sides of the first chunk
+# (must be multiple of 16 and 8 * multiple of NPROC below)
+NEX_XI = 160
+NEX_ETA = 160
+
+# number of MPI processors along the two sides of the first chunk
+NPROC_XI = 20
+NPROC_ETA = 20
+
+# 1D models with real structure:
+# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
+#
+# 1D models with only one fictitious averaged crustal layer:
+# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135_onecrust
+#
+# fully 3D models:
+# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
+# s20rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
+# s29ea, s29ea,sea99_jp3d1994,sea99,jp3d1994,heterogen
+#
+# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
+# to take the 1D crustal model from the
+# associated reference model rather than the default 3D crustal model
+# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
+MODEL = 1D_isotropic_prem
+
+# parameters describing the Earth model
+OCEANS = .false.
+ELLIPTICITY = .false.
+TOPOGRAPHY = .false.
+GRAVITY = .false.
+ROTATION = .false.
+ATTENUATION = .false.
+
+# absorbing boundary conditions for a regional simulation
+ABSORBING_CONDITIONS = .true.
+
+# record length in minutes
+RECORD_LENGTH_IN_MINUTES = 1.0d0
+
+# save AVS or OpenDX movies
+#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
+#MOVIE_COARSE does not work with create_movie_AVS_DX
+MOVIE_SURFACE = .true.
+MOVIE_VOLUME = .false.
+MOVIE_COARSE = .false.
+NTSTEP_BETWEEN_FRAMES = 50
+HDUR_MOVIE = 0.d0
+
+# save movie in volume. Will save element if center of element is in prescribed volume
+# top/bottom: depth in KM, use MOVIE_TOP = -100 to make sure the surface is stored.
+# west/east: longitude, degrees East [-180/180] top/bottom: latitute, degrees North [-90/90]
+# start/stop: frames will be stored at MOVIE_START + i*NSTEP_BETWEEN_FRAMES, where i=(0,1,2..) and iNSTEP_BETWEEN_FRAMES <= MOVIE_STOP
+# movie_volume_type: 1=strain, 2=time integral of strain, 3=\mu*time integral of strain
+# type 4 saves the trace and deviatoric stress in the whole volume, 5=displacement, 6=velocity
+MOVIE_VOLUME_TYPE = 2
+MOVIE_TOP_KM = -100.0
+MOVIE_BOTTOM_KM = 1000.0
+MOVIE_WEST_DEG = -90.0
+MOVIE_EAST_DEG = 90.0
+MOVIE_NORTH_DEG = 90.0
+MOVIE_SOUTH_DEG = -90.0
+MOVIE_START = 0
+MOVIE_STOP = 40000
+
+# save mesh files to check the mesh
+SAVE_MESH_FILES = .false.
+
+# restart files (number of runs can be 1, 2 or 3, choose 1 for no restart files)
+NUMBER_OF_RUNS = 1
+NUMBER_OF_THIS_RUN = 1
+
+# path to store the local database files on each node
+LOCAL_PATH = /scratch/yangl/DATABASES_MPI
+#/scratch/yangl/DATABASES_MPI
+
+# interval at which we output time step info and max of norm of displacement
+NTSTEP_BETWEEN_OUTPUT_INFO = 10000
+
+# interval in time steps for temporary writing of seismograms
+NTSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000
+NTSTEP_BETWEEN_READ_ADJSRC = 1000
+
+# output format for the seismograms (one can use either or all of the three formats)
+OUTPUT_SEISMOS_ASCII_TEXT = .true.
+OUTPUT_SEISMOS_SAC_ALPHANUM = .false.
+OUTPUT_SEISMOS_SAC_BINARY = .false.
+
+# rotate seismograms to Radial-Transverse-Z or use default North-East-Z reference frame
+ROTATE_SEISMOGRAMS_RT = .false.
+
+# decide if master process writes all the seismograms or if all processes do it in parallel
+WRITE_SEISMOGRAMS_BY_MASTER = .true.
+
+# save all seismograms in one large combined file instead of one file per seismogram
+# to avoid overloading shared non-local file systems such as GPFS for instance
+SAVE_ALL_SEISMOS_IN_ONE_FILE = .false.
+USE_BINARY_FOR_LARGE_FILE = .false.
+
+# flag to impose receivers at the surface or allow them to be buried
+RECEIVERS_CAN_BE_BURIED = .true.
+
+# print source time function
+PRINT_SOURCE_TIME_FUNCTION = .true.
Added: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_3_noattenuation
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_3_noattenuation (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/Par_file_NOISE_3_noattenuation 2011-04-01 15:57:37 UTC (rev 18161)
@@ -0,0 +1,120 @@
+
+# forward or adjoint simulation
+SIMULATION_TYPE = 3
+NOISE_TOMOGRAPHY = 3 # flag of noise tomography, three steps (1,2,3). If earthquake tomography, set it to 0.
+SAVE_FORWARD = .false. # save last frame of forward simulation or not
+
+# number of chunks (1,2,3 or 6)
+NCHUNKS = 1
+
+# angular width of the first chunk (not used if full sphere with six chunks)
+ANGULAR_WIDTH_XI_IN_DEGREES = 90.d0 # angular size of a chunk
+ANGULAR_WIDTH_ETA_IN_DEGREES = 90.d0
+CENTER_LATITUDE_IN_DEGREES = 0.d0
+CENTER_LONGITUDE_IN_DEGREES = 15.d0
+GAMMA_ROTATION_AZIMUTH = 0.d0
+
+# number of elements at the surface along the two sides of the first chunk
+# (must be multiple of 16 and 8 * multiple of NPROC below)
+NEX_XI = 160
+NEX_ETA = 160
+
+# number of MPI processors along the two sides of the first chunk
+NPROC_XI = 20
+NPROC_ETA = 20
+
+# 1D models with real structure:
+# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
+#
+# 1D models with only one fictitious averaged crustal layer:
+# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135_onecrust
+#
+# fully 3D models:
+# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
+# s20rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
+# s29ea, s29ea,sea99_jp3d1994,sea99,jp3d1994,heterogen
+#
+# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
+# to take the 1D crustal model from the
+# associated reference model rather than the default 3D crustal model
+# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
+MODEL = 1D_isotropic_prem
+
+# parameters describing the Earth model
+OCEANS = .false.
+ELLIPTICITY = .false.
+TOPOGRAPHY = .false.
+GRAVITY = .false.
+ROTATION = .false.
+ATTENUATION = .false.
+
+# absorbing boundary conditions for a regional simulation
+ABSORBING_CONDITIONS = .true.
+
+# record length in minutes
+RECORD_LENGTH_IN_MINUTES = 1.0d0
+
+# save AVS or OpenDX movies
+#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
+#MOVIE_COARSE does not work with create_movie_AVS_DX
+MOVIE_SURFACE = .true.
+MOVIE_VOLUME = .false.
+MOVIE_COARSE = .false.
+NTSTEP_BETWEEN_FRAMES = 50
+HDUR_MOVIE = 0.d0
+
+# save movie in volume. Will save element if center of element is in prescribed volume
+# top/bottom: depth in KM, use MOVIE_TOP = -100 to make sure the surface is stored.
+# west/east: longitude, degrees East [-180/180] top/bottom: latitute, degrees North [-90/90]
+# start/stop: frames will be stored at MOVIE_START + i*NSTEP_BETWEEN_FRAMES, where i=(0,1,2..) and iNSTEP_BETWEEN_FRAMES <= MOVIE_STOP
+# movie_volume_type: 1=strain, 2=time integral of strain, 3=\mu*time integral of strain
+# type 4 saves the trace and deviatoric stress in the whole volume, 5=displacement, 6=velocity
+MOVIE_VOLUME_TYPE = 2
+MOVIE_TOP_KM = -100.0
+MOVIE_BOTTOM_KM = 1000.0
+MOVIE_WEST_DEG = -90.0
+MOVIE_EAST_DEG = 90.0
+MOVIE_NORTH_DEG = 90.0
+MOVIE_SOUTH_DEG = -90.0
+MOVIE_START = 0
+MOVIE_STOP = 40000
+
+# save mesh files to check the mesh
+SAVE_MESH_FILES = .false.
+
+# restart files (number of runs can be 1, 2 or 3, choose 1 for no restart files)
+NUMBER_OF_RUNS = 1
+NUMBER_OF_THIS_RUN = 1
+
+# path to store the local database files on each node
+LOCAL_PATH = /scratch/yangl/DATABASES_MPI
+#/scratch/yangl/DATABASES_MPI
+
+# interval at which we output time step info and max of norm of displacement
+NTSTEP_BETWEEN_OUTPUT_INFO = 10000
+
+# interval in time steps for temporary writing of seismograms
+NTSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000
+NTSTEP_BETWEEN_READ_ADJSRC = 1000
+
+# output format for the seismograms (one can use either or all of the three formats)
+OUTPUT_SEISMOS_ASCII_TEXT = .true.
+OUTPUT_SEISMOS_SAC_ALPHANUM = .false.
+OUTPUT_SEISMOS_SAC_BINARY = .false.
+
+# rotate seismograms to Radial-Transverse-Z or use default North-East-Z reference frame
+ROTATE_SEISMOGRAMS_RT = .false.
+
+# decide if master process writes all the seismograms or if all processes do it in parallel
+WRITE_SEISMOGRAMS_BY_MASTER = .true.
+
+# save all seismograms in one large combined file instead of one file per seismogram
+# to avoid overloading shared non-local file systems such as GPFS for instance
+SAVE_ALL_SEISMOS_IN_ONE_FILE = .false.
+USE_BINARY_FOR_LARGE_FILE = .false.
+
+# flag to impose receivers at the surface or allow them to be buried
+RECEIVERS_CAN_BE_BURIED = .true.
+
+# print source time function
+PRINT_SOURCE_TIME_FUNCTION = .true.
Added: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/STATIONS_NOISE
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/STATIONS_NOISE (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/STATIONS_NOISE 2011-04-01 15:57:37 UTC (rev 18161)
@@ -0,0 +1,13 @@
+A1 II 0.0 0.0 0.0 0.0
+A2 II 0.0 5.0 0.0 0.0
+A3 II 0.0 10.0 0.0 0.0
+A4 II 0.0 15.0 0.0 0.0
+A5 II 0.0 20.0 0.0 0.0
+A6 II 0.0 25.0 0.0 0.0
+A7 II 0.0 30.0 0.0 0.0
+A8 II 0.0 35.0 0.0 0.0
+A9 II 0.0 40.0 0.0 0.0
+A10 II 0.0 45.0 0.0 0.0
+A11 II 0.0 50.0 0.0 0.0
+A12 II 0.0 55.0 0.0 0.0
+A13 II 0.0 60.0 0.0 0.0
Added: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/STATIONS_NOISE_ADJOINT
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/STATIONS_NOISE_ADJOINT (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/STATIONS_NOISE_ADJOINT 2011-04-01 15:57:37 UTC (rev 18161)
@@ -0,0 +1 @@
+A7 II 0.0 30.0 0.0 0.0
Added: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/S_squared
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/S_squared (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/S_squared 2011-04-01 15:57:37 UTC (rev 18161)
@@ -0,0 +1,799 @@
+ -7.5810000e+01 2.0083475e-19
+ -7.5620000e+01 1.9220480e-19
+ -7.5430000e+01 1.7529785e-19
+ -7.5240000e+01 1.5080677e-19
+ -7.5050000e+01 1.1973878e-19
+ -7.4860000e+01 8.3378339e-20
+ -7.4670000e+01 4.3239679e-20
+ -7.4480000e+01 1.0106876e-21
+ -7.4290000e+01 -4.1510205e-20
+ -7.4100000e+01 -8.2482300e-20
+ -7.3910000e+01 -1.2009421e-19
+ -7.3720000e+01 -1.5263679e-19
+ -7.3530000e+01 -1.7857550e-19
+ -7.3340000e+01 -1.9661979e-19
+ -7.3150000e+01 -2.0578729e-19
+ -7.2960000e+01 -2.0545992e-19
+ -7.2770000e+01 -1.9542883e-19
+ -7.2580000e+01 -1.7592525e-19
+ -7.2390000e+01 -1.4763374e-19
+ -7.2200000e+01 -1.1168541e-19
+ -7.2010000e+01 -6.9628508e-20
+ -7.1820000e+01 -2.3375690e-20
+ -7.1630000e+01 2.4872197e-20
+ -7.1440000e+01 7.2723698e-20
+ -7.1250000e+01 1.1771219e-19
+ -7.1060000e+01 1.5742382e-19
+ -7.0870000e+01 1.8963105e-19
+ -7.0680000e+01 2.1242398e-19
+ -7.0490000e+01 2.2433077e-19
+ -7.0300000e+01 2.2441895e-19
+ -7.0110000e+01 2.1236989e-19
+ -6.9920000e+01 1.8852037e-19
+ -6.9730000e+01 1.5386707e-19
+ -6.9540000e+01 1.1003221e-19
+ -6.9350000e+01 5.9191273e-20
+ -6.9160000e+01 3.9661551e-21
+ -6.8970000e+01 -5.2710589e-20
+ -6.8780000e+01 -1.0775270e-19
+ -6.8590000e+01 -1.5808949e-19
+ -6.8400000e+01 -2.0083774e-19
+ -6.8210000e+01 -2.3346641e-19
+ -6.8020000e+01 -2.5394490e-19
+ -6.7830000e+01 -2.6086654e-19
+ -6.7640000e+01 -2.5354037e-19
+ -6.7450000e+01 -2.3204589e-19
+ -6.7260000e+01 -1.9724741e-19
+ -6.7070000e+01 -1.5076633e-19
+ -6.6880000e+01 -9.4911984e-20
+ -6.6690000e+01 -3.2573180e-20
+ -6.6500000e+01 3.2925061e-20
+ -6.6310000e+01 9.7993561e-20
+ -6.6120000e+01 1.5897190e-19
+ -6.5930000e+01 2.1233279e-19
+ -6.5740000e+01 2.5488724e-19
+ -6.5550000e+01 2.8397889e-19
+ -6.5360000e+01 2.9765587e-19
+ -6.5170000e+01 2.9480847e-19
+ -6.4980000e+01 2.7526252e-19
+ -6.4790000e+01 2.3982016e-19
+ -6.4600000e+01 1.9024237e-19
+ -6.4410000e+01 1.2917147e-19
+ -6.4220000e+01 5.9995411e-20
+ -6.4030000e+01 -1.3339828e-20
+ -6.3840000e+01 -8.6559668e-20
+ -6.3650000e+01 -1.5531831e-19
+ -6.3460000e+01 -2.1546775e-19
+ -6.3270000e+01 -2.6332049e-19
+ -6.3080000e+01 -2.9588787e-19
+ -6.2890000e+01 -3.1107739e-19
+ -6.2700000e+01 -3.0783468e-19
+ -6.2510000e+01 -2.8622009e-19
+ -6.2320000e+01 -2.4741406e-19
+ -6.2130000e+01 -1.9365082e-19
+ -6.1940000e+01 -1.2808502e-19
+ -6.1750000e+01 -5.4600362e-20
+ -6.1560000e+01 2.2426535e-20
+ -6.1370000e+01 9.8391447e-20
+ -6.1180000e+01 1.6873968e-19
+ -6.0990000e+01 2.2923767e-19
+ -6.0800000e+01 2.7622279e-19
+ -6.0610000e+01 3.0681694e-19
+ -6.0420000e+01 3.1909258e-19
+ -6.0230000e+01 3.1218306e-19
+ -6.0040000e+01 2.8633230e-19
+ -5.9850000e+01 2.4288267e-19
+ -5.9660000e+01 1.8420267e-19
+ -5.9470000e+01 1.1355904e-19
+ -5.9280000e+01 3.4940085e-20
+ -5.9090000e+01 -4.7160837e-20
+ -5.8900000e+01 -1.2799106e-19
+ -5.8710000e+01 -2.0280178e-19
+ -5.8520000e+01 -2.6711570e-19
+ -5.8330000e+01 -3.1698554e-19
+ -5.8140000e+01 -3.4922950e-19
+ -5.7950000e+01 -3.6163018e-19
+ -5.7760000e+01 -3.5308443e-19
+ -5.7570000e+01 -3.2369325e-19
+ -5.7380000e+01 -2.7478380e-19
+ -5.7190000e+01 -2.0885893e-19
+ -5.7000000e+01 -1.2947405e-19
+ -5.6810000e+01 -4.1045823e-20
+ -5.6620000e+01 5.1397826e-20
+ -5.6430000e+01 1.4251234e-19
+ -5.6240000e+01 2.2695775e-19
+ -5.6050000e+01 2.9972579e-19
+ -5.5860000e+01 3.5645090e-19
+ -5.5670000e+01 3.9368369e-19
+ -5.5480000e+01 4.0910773e-19
+ -5.5290000e+01 4.0168438e-19
+ -5.5100000e+01 3.7171582e-19
+ -5.4910000e+01 3.2082240e-19
+ -5.4720000e+01 2.5183700e-19
+ -5.4530000e+01 1.6862509e-19
+ -5.4340000e+01 7.5844464e-20
+ -5.4150000e+01 -2.1337642e-20
+ -5.3960000e+01 -1.1754995e-19
+ -5.3770000e+01 -2.0752479e-19
+ -5.3580000e+01 -2.8638873e-19
+ -5.3390000e+01 -3.4991992e-19
+ -5.3200000e+01 -3.9475923e-19
+ -5.3010000e+01 -4.1856685e-19
+ -5.2820000e+01 -4.2011991e-19
+ -5.2630000e+01 -3.9935087e-19
+ -5.2440000e+01 -3.5732939e-19
+ -5.2250000e+01 -2.9619265e-19
+ -5.2060000e+01 -2.1903051e-19
+ -5.1870000e+01 -1.2973246e-19
+ -5.1680000e+01 -3.2803460e-20
+ -5.1490000e+01 6.6844511e-20
+ -5.1300000e+01 1.6411831e-19
+ -5.1110000e+01 2.5397989e-19
+ -5.0920000e+01 3.3168548e-19
+ -5.0730000e+01 3.9301904e-19
+ -5.0540000e+01 4.3451051e-19
+ -5.0350000e+01 4.5362837e-19
+ -5.0160000e+01 4.4893518e-19
+ -4.9970000e+01 4.2019525e-19
+ -4.9780000e+01 3.6842482e-19
+ -4.9590000e+01 2.9587708e-19
+ -4.9400000e+01 2.0595855e-19
+ -4.9210000e+01 1.0307694e-19
+ -4.9020000e+01 -7.5736812e-21
+ -4.8830000e+01 -1.2028009e-19
+ -4.8640000e+01 -2.2911445e-19
+ -4.8450000e+01 -3.2825766e-19
+ -4.8260000e+01 -4.1232107e-19
+ -4.8070000e+01 -4.7664596e-19
+ -4.7880000e+01 -5.1756148e-19
+ -4.7690000e+01 -5.3258546e-19
+ -4.7500000e+01 -5.2055661e-19
+ -4.7310000e+01 -4.8169230e-19
+ -4.7120000e+01 -4.1757152e-19
+ -4.6930000e+01 -3.3104772e-19
+ -4.6740000e+01 -2.2610035e-19
+ -4.6550000e+01 -1.0763689e-19
+ -4.6360000e+01 1.8741548e-20
+ -4.6170000e+01 1.4699816e-19
+ -4.5980000e+01 2.7092971e-19
+ -4.5790000e+01 3.8443735e-19
+ -4.5600000e+01 4.8179060e-19
+ -4.5410000e+01 5.5787848e-19
+ -4.5220000e+01 6.0844227e-19
+ -4.5030000e+01 6.3028490e-19
+ -4.4840000e+01 6.2145114e-19
+ -4.4650000e+01 5.8137141e-19
+ -4.4460000e+01 5.1096036e-19
+ -4.4270000e+01 4.1266036e-19
+ -4.4080000e+01 2.9041926e-19
+ -4.3890000e+01 1.4959319e-19
+ -4.3700000e+01 -3.2324241e-21
+ -4.3510000e+01 -1.6050562e-19
+ -4.3320000e+01 -3.1405071e-19
+ -4.3130000e+01 -4.5547972e-19
+ -4.2940000e+01 -5.7665400e-19
+ -4.2750000e+01 -6.7017014e-19
+ -4.2560000e+01 -7.2983910e-19
+ -4.2370000e+01 -7.5112406e-19
+ -4.2180000e+01 -7.3150248e-19
+ -4.1990000e+01 -6.7072072e-19
+ -4.1800000e+01 -5.7091573e-19
+ -4.1610000e+01 -4.3658770e-19
+ -4.1420000e+01 -2.7441787e-19
+ -4.1230000e+01 -9.2938263e-20
+ -4.1040000e+01 9.7928901e-20
+ -4.0850000e+01 2.8743343e-19
+ -4.0660000e+01 4.6462495e-19
+ -4.0470000e+01 6.1900506e-19
+ -4.0280000e+01 7.4116124e-19
+ -4.0090000e+01 8.2334194e-19
+ -3.9900000e+01 8.5993725e-19
+ -3.9710000e+01 8.4783677e-19
+ -3.9520000e+01 7.8664520e-19
+ -3.9330000e+01 6.7874560e-19
+ -3.9140000e+01 5.2920900e-19
+ -3.8950000e+01 3.4555773e-19
+ -3.8760000e+01 1.3739630e-19
+ -3.8570000e+01 -8.4071055e-20
+ -3.8380000e+01 -3.0661268e-19
+ -3.8190000e+01 -5.1760598e-19
+ -3.8000000e+01 -7.0470977e-19
+ -3.7810000e+01 -8.5654316e-19
+ -3.7620000e+01 -9.6333999e-19
+ -3.7430000e+01 -1.0175471e-18
+ -3.7240000e+01 -1.0143341e-18
+ -3.7050000e+01 -9.5198336e-19
+ -3.6860000e+01 -8.3212971e-19
+ -3.6670000e+01 -6.5982566e-19
+ -3.6480000e+01 -4.4341220e-19
+ -3.6290000e+01 -1.9418638e-19
+ -3.6100000e+01 7.4131800e-20
+ -3.5910000e+01 3.4611692e-19
+ -3.5720000e+01 6.0550846e-19
+ -3.5530000e+01 8.3618806e-19
+ -3.5340000e+01 1.0232050e-18
+ -3.5150000e+01 1.1537822e-18
+ -3.4960000e+01 1.2182298e-18
+ -3.4770000e+01 1.2106966e-18
+ -3.4580000e+01 1.1296951e-18
+ -3.4390000e+01 9.7835102e-19
+ -3.4200000e+01 7.6434700e-19
+ -3.4010000e+01 4.9955265e-19
+ -3.3820000e+01 1.9935691e-19
+ -3.3630000e+01 -1.1825671e-19
+ -3.3440000e+01 -4.3383120e-19
+ -3.3250000e+01 -7.2767627e-19
+ -3.3060000e+01 -9.8114099e-19
+ -3.2870000e+01 -1.1778295e-18
+ -3.2680000e+01 -1.3046756e-18
+ -3.2490000e+01 -1.3528028e-18
+ -3.2300000e+01 -1.3181151e-18
+ -3.2110000e+01 -1.2015820e-18
+ -3.1920000e+01 -1.0092053e-18
+ -3.1730000e+01 -7.5167615e-19
+ -3.1540000e+01 -4.4375100e-19
+ -3.1350000e+01 -1.0339019e-19
+ -3.1160000e+01 2.4928452e-19
+ -3.0970000e+01 5.9315068e-19
+ -3.0780000e+01 9.0729832e-19
+ -3.0590000e+01 1.1722309e-18
+ -3.0400000e+01 1.3709932e-18
+ -3.0210000e+01 1.4901705e-18
+ -3.0020000e+01 1.5207105e-18
+ -2.9830000e+01 1.4585276e-18
+ -2.9640000e+01 1.3048553e-18
+ -2.9450000e+01 1.0663233e-18
+ -2.9260000e+01 7.5474484e-19
+ -2.9070000e+01 3.8661012e-19
+ -2.8880000e+01 -1.7701900e-20
+ -2.8690000e+01 -4.3496996e-19
+ -2.8500000e+01 -8.4039233e-19
+ -2.8310000e+01 -1.2090055e-18
+ -2.8120000e+01 -1.5171856e-18
+ -2.7930000e+01 -1.7441435e-18
+ -2.7740000e+01 -1.8733199e-18
+ -2.7550000e+01 -1.8935795e-18
+ -2.7360000e+01 -1.8001133e-18
+ -2.7170000e+01 -1.5949702e-18
+ -2.6980000e+01 -1.2871610e-18
+ -2.6790000e+01 -8.9230633e-19
+ -2.6600000e+01 -4.3183467e-19
+ -2.6410000e+01 6.8230074e-20
+ -2.6220000e+01 5.7881722e-19
+ -2.6030000e+01 1.0695922e-18
+ -2.5840000e+01 1.5108511e-18
+ -2.5650000e+01 1.8753932e-18
+ -2.5460000e+01 2.1402407e-18
+ -2.5270000e+01 2.2880846e-18
+ -2.5080000e+01 2.3083610e-18
+ -2.4890000e+01 2.1978907e-18
+ -2.4700000e+01 1.9610503e-18
+ -2.4510000e+01 1.6094796e-18
+ -2.4320000e+01 1.1613619e-18
+ -2.4130000e+01 6.4034236e-19
+ -2.3940000e+01 7.4170036e-20
+ -2.3750000e+01 -5.0684259e-19
+ -2.3560000e+01 -1.0714449e-18
+ -2.3370000e+01 -1.5890677e-18
+ -2.3180000e+01 -2.0313591e-18
+ -2.2990000e+01 -2.3735726e-18
+ -2.2800000e+01 -2.5957604e-18
+ -2.2610000e+01 -2.6837397e-18
+ -2.2420000e+01 -2.6298080e-18
+ -2.2230000e+01 -2.4331948e-18
+ -2.2040000e+01 -2.1002378e-18
+ -2.1850000e+01 -1.6442742e-18
+ -2.1660000e+01 -1.0852384e-18
+ -2.1470000e+01 -4.4896117e-19
+ -2.1280000e+01 2.3382952e-19
+ -2.1090000e+01 9.2879539e-19
+ -2.0900000e+01 1.5995338e-18
+ -2.0710000e+01 2.2093363e-18
+ -2.0520000e+01 2.7230887e-18
+ -2.0330000e+01 3.1092217e-18
+ -2.0140000e+01 3.3416006e-18
+ -1.9950000e+01 3.4012322e-18
+ -1.9760000e+01 3.2776657e-18
+ -1.9570000e+01 2.9699741e-18
+ -1.9380000e+01 2.4872234e-18
+ -1.9190000e+01 1.8483709e-18
+ -1.9000000e+01 1.0815724e-18
+ -1.8810000e+01 2.2292353e-19
+ -1.8620000e+01 -6.8529585e-19
+ -1.8430000e+01 -1.5967693e-18
+ -1.8240000e+01 -2.4635475e-18
+ -1.8050000e+01 -3.2385574e-18
+ -1.7860000e+01 -3.8780665e-18
+ -1.7670000e+01 -4.3439600e-18
+ -1.7480000e+01 -4.6057158e-18
+ -1.7290000e+01 -4.6419909e-18
+ -1.7100000e+01 -4.4417680e-18
+ -1.6910000e+01 -4.0050424e-18
+ -1.6720000e+01 -3.3430568e-18
+ -1.6530000e+01 -2.4781055e-18
+ -1.6340000e+01 -1.4429418e-18
+ -1.6150000e+01 -2.7981391e-19
+ -1.5960000e+01 9.6084977e-19
+ -1.5770000e+01 2.2220983e-18
+ -1.5580000e+01 3.4424597e-18
+ -1.5390000e+01 4.5582867e-18
+ -1.5200000e+01 5.5064387e-18
+ -1.5010000e+01 6.2272659e-18
+ -1.4820000e+01 6.6678244e-18
+ -1.4630000e+01 6.7852074e-18
+ -1.4440000e+01 6.5498212e-18
+ -1.4250000e+01 5.9483931e-18
+ -1.4060000e+01 4.9864585e-18
+ -1.3870000e+01 3.6900596e-18
+ -1.3680000e+01 2.1064010e-18
+ -1.3490000e+01 3.0324576e-19
+ -1.3300000e+01 -1.6330909e-18
+ -1.3110000e+01 -3.6011971e-18
+ -1.2920000e+01 -5.4893960e-18
+ -1.2730000e+01 -7.1817592e-18
+ -1.2540000e+01 -8.5649837e-18
+ -1.2350000e+01 -9.5357145e-18
+ -1.2160000e+01 -1.0007822e-17
+ -1.1970000e+01 -9.9191087e-18
+ -1.1780000e+01 -9.2369373e-18
+ -1.1590000e+01 -7.9623161e-18
+ -1.1400000e+01 -6.1320895e-18
+ -1.1210000e+01 -3.8189941e-18
+ -1.1020000e+01 -1.1294951e-18
+ -1.0830000e+01 1.8005368e-18
+ -1.0640000e+01 4.8121066e-18
+ -1.0450000e+01 7.7310128e-18
+ -1.0260000e+01 1.0376990e-17
+ -1.0070000e+01 1.2573612e-17
+ -9.8800000e+00 1.4158452e-17
+ -9.6900000e+00 1.4993027e-17
+ -9.5000000e+00 1.4972037e-17
+ -9.3100000e+00 1.4031478e-17
+ -9.1200000e+00 1.2155218e-17
+ -8.9300000e+00 9.3796708e-18
+ -8.7400000e+00 5.7962211e-18
+ -8.5500000e+00 1.5510879e-18
+ -8.3600000e+00 -3.1576515e-18
+ -8.1700000e+00 -8.0861012e-18
+ -7.9800000e+00 -1.2953785e-17
+ -7.7900000e+00 -1.7456081e-17
+ -7.6000000e+00 -2.1279676e-17
+ -7.4100000e+00 -2.4120614e-17
+ -7.2200000e+00 -2.5704247e-17
+ -7.0300000e+00 -2.5806098e-17
+ -6.8400000e+00 -2.4272354e-17
+ -6.6500000e+00 -2.1038471e-17
+ -6.4600000e+00 -1.6144202e-17
+ -6.2700000e+00 -9.7433531e-18
+ -6.0800000e+00 -2.1066809e-18
+ -5.8900000e+00 6.3833442e-18
+ -5.7000000e+00 1.5246693e-17
+ -5.5100000e+00 2.3928993e-17
+ -5.3200000e+00 3.1833180e-17
+ -5.1300000e+00 3.8358117e-17
+ -4.9400000e+00 4.2941933e-17
+ -4.7500000e+00 4.5107075e-17
+ -4.5600000e+00 4.4503585e-17
+ -4.3700000e+00 4.0946839e-17
+ -4.1800000e+00 3.4446037e-17
+ -3.9900000e+00 2.5220178e-17
+ -3.8000000e+00 1.3698996e-17
+ -3.6100000e+00 5.0743105e-19
+ -3.4200000e+00 -1.3566495e-17
+ -3.2300000e+00 -2.7619073e-17
+ -3.0400000e+00 -4.0690808e-17
+ -2.8500000e+00 -5.1834784e-17
+ -2.6600000e+00 -6.0188387e-17
+ -2.4700000e+00 -6.5042697e-17
+ -2.2800000e+00 -6.5903659e-17
+ -2.0900000e+00 -6.2539553e-17
+ -1.9000000e+00 -5.5010154e-17
+ -1.7100000e+00 -4.3674315e-17
+ -1.5200000e+00 -2.9174442e-17
+ -1.3300000e+00 -1.2398213e-17
+ -1.1400000e+00 5.5801177e-18
+ -9.5000000e-01 2.3574673e-17
+ -7.6000000e-01 4.0373502e-17
+ -5.7000000e-01 5.4828016e-17
+ -3.8000000e-01 6.5939107e-17
+ -1.9000000e-01 7.2932921e-17
+ 0.0000000e+00 7.5320009e-17
+ 1.9000000e-01 7.2932921e-17
+ 3.8000000e-01 6.5939107e-17
+ 5.7000000e-01 5.4828016e-17
+ 7.6000000e-01 4.0373502e-17
+ 9.5000000e-01 2.3574673e-17
+ 1.1400000e+00 5.5801177e-18
+ 1.3300000e+00 -1.2398213e-17
+ 1.5200000e+00 -2.9174442e-17
+ 1.7100000e+00 -4.3674315e-17
+ 1.9000000e+00 -5.5010154e-17
+ 2.0900000e+00 -6.2539553e-17
+ 2.2800000e+00 -6.5903659e-17
+ 2.4700000e+00 -6.5042697e-17
+ 2.6600000e+00 -6.0188387e-17
+ 2.8500000e+00 -5.1834784e-17
+ 3.0400000e+00 -4.0690808e-17
+ 3.2300000e+00 -2.7619073e-17
+ 3.4200000e+00 -1.3566495e-17
+ 3.6100000e+00 5.0743105e-19
+ 3.8000000e+00 1.3698996e-17
+ 3.9900000e+00 2.5220178e-17
+ 4.1800000e+00 3.4446037e-17
+ 4.3700000e+00 4.0946839e-17
+ 4.5600000e+00 4.4503585e-17
+ 4.7500000e+00 4.5107075e-17
+ 4.9400000e+00 4.2941933e-17
+ 5.1300000e+00 3.8358117e-17
+ 5.3200000e+00 3.1833180e-17
+ 5.5100000e+00 2.3928993e-17
+ 5.7000000e+00 1.5246693e-17
+ 5.8900000e+00 6.3833442e-18
+ 6.0800000e+00 -2.1066809e-18
+ 6.2700000e+00 -9.7433531e-18
+ 6.4600000e+00 -1.6144202e-17
+ 6.6500000e+00 -2.1038471e-17
+ 6.8400000e+00 -2.4272354e-17
+ 7.0300000e+00 -2.5806098e-17
+ 7.2200000e+00 -2.5704247e-17
+ 7.4100000e+00 -2.4120614e-17
+ 7.6000000e+00 -2.1279676e-17
+ 7.7900000e+00 -1.7456081e-17
+ 7.9800000e+00 -1.2953785e-17
+ 8.1700000e+00 -8.0861012e-18
+ 8.3600000e+00 -3.1576515e-18
+ 8.5500000e+00 1.5510879e-18
+ 8.7400000e+00 5.7962211e-18
+ 8.9300000e+00 9.3796708e-18
+ 9.1200000e+00 1.2155218e-17
+ 9.3100000e+00 1.4031478e-17
+ 9.5000000e+00 1.4972037e-17
+ 9.6900000e+00 1.4993027e-17
+ 9.8800000e+00 1.4158452e-17
+ 1.0070000e+01 1.2573612e-17
+ 1.0260000e+01 1.0376990e-17
+ 1.0450000e+01 7.7310128e-18
+ 1.0640000e+01 4.8121066e-18
+ 1.0830000e+01 1.8005368e-18
+ 1.1020000e+01 -1.1294951e-18
+ 1.1210000e+01 -3.8189941e-18
+ 1.1400000e+01 -6.1320895e-18
+ 1.1590000e+01 -7.9623161e-18
+ 1.1780000e+01 -9.2369373e-18
+ 1.1970000e+01 -9.9191087e-18
+ 1.2160000e+01 -1.0007822e-17
+ 1.2350000e+01 -9.5357145e-18
+ 1.2540000e+01 -8.5649837e-18
+ 1.2730000e+01 -7.1817592e-18
+ 1.2920000e+01 -5.4893960e-18
+ 1.3110000e+01 -3.6011971e-18
+ 1.3300000e+01 -1.6330909e-18
+ 1.3490000e+01 3.0324576e-19
+ 1.3680000e+01 2.1064010e-18
+ 1.3870000e+01 3.6900596e-18
+ 1.4060000e+01 4.9864585e-18
+ 1.4250000e+01 5.9483931e-18
+ 1.4440000e+01 6.5498212e-18
+ 1.4630000e+01 6.7852074e-18
+ 1.4820000e+01 6.6678244e-18
+ 1.5010000e+01 6.2272659e-18
+ 1.5200000e+01 5.5064387e-18
+ 1.5390000e+01 4.5582867e-18
+ 1.5580000e+01 3.4424597e-18
+ 1.5770000e+01 2.2220983e-18
+ 1.5960000e+01 9.6084977e-19
+ 1.6150000e+01 -2.7981391e-19
+ 1.6340000e+01 -1.4429418e-18
+ 1.6530000e+01 -2.4781055e-18
+ 1.6720000e+01 -3.3430568e-18
+ 1.6910000e+01 -4.0050424e-18
+ 1.7100000e+01 -4.4417680e-18
+ 1.7290000e+01 -4.6419909e-18
+ 1.7480000e+01 -4.6057158e-18
+ 1.7670000e+01 -4.3439600e-18
+ 1.7860000e+01 -3.8780665e-18
+ 1.8050000e+01 -3.2385574e-18
+ 1.8240000e+01 -2.4635475e-18
+ 1.8430000e+01 -1.5967693e-18
+ 1.8620000e+01 -6.8529585e-19
+ 1.8810000e+01 2.2292353e-19
+ 1.9000000e+01 1.0815724e-18
+ 1.9190000e+01 1.8483709e-18
+ 1.9380000e+01 2.4872234e-18
+ 1.9570000e+01 2.9699741e-18
+ 1.9760000e+01 3.2776657e-18
+ 1.9950000e+01 3.4012322e-18
+ 2.0140000e+01 3.3416006e-18
+ 2.0330000e+01 3.1092217e-18
+ 2.0520000e+01 2.7230887e-18
+ 2.0710000e+01 2.2093363e-18
+ 2.0900000e+01 1.5995338e-18
+ 2.1090000e+01 9.2879539e-19
+ 2.1280000e+01 2.3382952e-19
+ 2.1470000e+01 -4.4896117e-19
+ 2.1660000e+01 -1.0852384e-18
+ 2.1850000e+01 -1.6442742e-18
+ 2.2040000e+01 -2.1002378e-18
+ 2.2230000e+01 -2.4331948e-18
+ 2.2420000e+01 -2.6298080e-18
+ 2.2610000e+01 -2.6837397e-18
+ 2.2800000e+01 -2.5957604e-18
+ 2.2990000e+01 -2.3735726e-18
+ 2.3180000e+01 -2.0313591e-18
+ 2.3370000e+01 -1.5890677e-18
+ 2.3560000e+01 -1.0714449e-18
+ 2.3750000e+01 -5.0684259e-19
+ 2.3940000e+01 7.4170036e-20
+ 2.4130000e+01 6.4034236e-19
+ 2.4320000e+01 1.1613619e-18
+ 2.4510000e+01 1.6094796e-18
+ 2.4700000e+01 1.9610503e-18
+ 2.4890000e+01 2.1978907e-18
+ 2.5080000e+01 2.3083610e-18
+ 2.5270000e+01 2.2880846e-18
+ 2.5460000e+01 2.1402407e-18
+ 2.5650000e+01 1.8753932e-18
+ 2.5840000e+01 1.5108511e-18
+ 2.6030000e+01 1.0695922e-18
+ 2.6220000e+01 5.7881722e-19
+ 2.6410000e+01 6.8230074e-20
+ 2.6600000e+01 -4.3183467e-19
+ 2.6790000e+01 -8.9230633e-19
+ 2.6980000e+01 -1.2871610e-18
+ 2.7170000e+01 -1.5949702e-18
+ 2.7360000e+01 -1.8001133e-18
+ 2.7550000e+01 -1.8935795e-18
+ 2.7740000e+01 -1.8733199e-18
+ 2.7930000e+01 -1.7441435e-18
+ 2.8120000e+01 -1.5171856e-18
+ 2.8310000e+01 -1.2090055e-18
+ 2.8500000e+01 -8.4039233e-19
+ 2.8690000e+01 -4.3496996e-19
+ 2.8880000e+01 -1.7701900e-20
+ 2.9070000e+01 3.8661012e-19
+ 2.9260000e+01 7.5474484e-19
+ 2.9450000e+01 1.0663233e-18
+ 2.9640000e+01 1.3048553e-18
+ 2.9830000e+01 1.4585276e-18
+ 3.0020000e+01 1.5207105e-18
+ 3.0210000e+01 1.4901705e-18
+ 3.0400000e+01 1.3709932e-18
+ 3.0590000e+01 1.1722309e-18
+ 3.0780000e+01 9.0729832e-19
+ 3.0970000e+01 5.9315068e-19
+ 3.1160000e+01 2.4928452e-19
+ 3.1350000e+01 -1.0339019e-19
+ 3.1540000e+01 -4.4375100e-19
+ 3.1730000e+01 -7.5167615e-19
+ 3.1920000e+01 -1.0092053e-18
+ 3.2110000e+01 -1.2015820e-18
+ 3.2300000e+01 -1.3181151e-18
+ 3.2490000e+01 -1.3528028e-18
+ 3.2680000e+01 -1.3046756e-18
+ 3.2870000e+01 -1.1778295e-18
+ 3.3060000e+01 -9.8114099e-19
+ 3.3250000e+01 -7.2767627e-19
+ 3.3440000e+01 -4.3383120e-19
+ 3.3630000e+01 -1.1825671e-19
+ 3.3820000e+01 1.9935691e-19
+ 3.4010000e+01 4.9955265e-19
+ 3.4200000e+01 7.6434700e-19
+ 3.4390000e+01 9.7835102e-19
+ 3.4580000e+01 1.1296951e-18
+ 3.4770000e+01 1.2106966e-18
+ 3.4960000e+01 1.2182298e-18
+ 3.5150000e+01 1.1537822e-18
+ 3.5340000e+01 1.0232050e-18
+ 3.5530000e+01 8.3618806e-19
+ 3.5720000e+01 6.0550846e-19
+ 3.5910000e+01 3.4611692e-19
+ 3.6100000e+01 7.4131800e-20
+ 3.6290000e+01 -1.9418638e-19
+ 3.6480000e+01 -4.4341220e-19
+ 3.6670000e+01 -6.5982566e-19
+ 3.6860000e+01 -8.3212971e-19
+ 3.7050000e+01 -9.5198336e-19
+ 3.7240000e+01 -1.0143341e-18
+ 3.7430000e+01 -1.0175471e-18
+ 3.7620000e+01 -9.6333999e-19
+ 3.7810000e+01 -8.5654316e-19
+ 3.8000000e+01 -7.0470977e-19
+ 3.8190000e+01 -5.1760598e-19
+ 3.8380000e+01 -3.0661268e-19
+ 3.8570000e+01 -8.4071055e-20
+ 3.8760000e+01 1.3739630e-19
+ 3.8950000e+01 3.4555773e-19
+ 3.9140000e+01 5.2920900e-19
+ 3.9330000e+01 6.7874560e-19
+ 3.9520000e+01 7.8664520e-19
+ 3.9710000e+01 8.4783677e-19
+ 3.9900000e+01 8.5993725e-19
+ 4.0090000e+01 8.2334194e-19
+ 4.0280000e+01 7.4116124e-19
+ 4.0470000e+01 6.1900506e-19
+ 4.0660000e+01 4.6462495e-19
+ 4.0850000e+01 2.8743343e-19
+ 4.1040000e+01 9.7928901e-20
+ 4.1230000e+01 -9.2938263e-20
+ 4.1420000e+01 -2.7441787e-19
+ 4.1610000e+01 -4.3658770e-19
+ 4.1800000e+01 -5.7091573e-19
+ 4.1990000e+01 -6.7072072e-19
+ 4.2180000e+01 -7.3150248e-19
+ 4.2370000e+01 -7.5112406e-19
+ 4.2560000e+01 -7.2983910e-19
+ 4.2750000e+01 -6.7017014e-19
+ 4.2940000e+01 -5.7665400e-19
+ 4.3130000e+01 -4.5547972e-19
+ 4.3320000e+01 -3.1405071e-19
+ 4.3510000e+01 -1.6050562e-19
+ 4.3700000e+01 -3.2324241e-21
+ 4.3890000e+01 1.4959319e-19
+ 4.4080000e+01 2.9041926e-19
+ 4.4270000e+01 4.1266036e-19
+ 4.4460000e+01 5.1096036e-19
+ 4.4650000e+01 5.8137141e-19
+ 4.4840000e+01 6.2145114e-19
+ 4.5030000e+01 6.3028490e-19
+ 4.5220000e+01 6.0844227e-19
+ 4.5410000e+01 5.5787848e-19
+ 4.5600000e+01 4.8179060e-19
+ 4.5790000e+01 3.8443735e-19
+ 4.5980000e+01 2.7092971e-19
+ 4.6170000e+01 1.4699816e-19
+ 4.6360000e+01 1.8741548e-20
+ 4.6550000e+01 -1.0763689e-19
+ 4.6740000e+01 -2.2610035e-19
+ 4.6930000e+01 -3.3104772e-19
+ 4.7120000e+01 -4.1757152e-19
+ 4.7310000e+01 -4.8169230e-19
+ 4.7500000e+01 -5.2055661e-19
+ 4.7690000e+01 -5.3258546e-19
+ 4.7880000e+01 -5.1756148e-19
+ 4.8070000e+01 -4.7664596e-19
+ 4.8260000e+01 -4.1232107e-19
+ 4.8450000e+01 -3.2825766e-19
+ 4.8640000e+01 -2.2911445e-19
+ 4.8830000e+01 -1.2028009e-19
+ 4.9020000e+01 -7.5736812e-21
+ 4.9210000e+01 1.0307694e-19
+ 4.9400000e+01 2.0595855e-19
+ 4.9590000e+01 2.9587708e-19
+ 4.9780000e+01 3.6842482e-19
+ 4.9970000e+01 4.2019525e-19
+ 5.0160000e+01 4.4893518e-19
+ 5.0350000e+01 4.5362837e-19
+ 5.0540000e+01 4.3451051e-19
+ 5.0730000e+01 3.9301904e-19
+ 5.0920000e+01 3.3168548e-19
+ 5.1110000e+01 2.5397989e-19
+ 5.1300000e+01 1.6411831e-19
+ 5.1490000e+01 6.6844511e-20
+ 5.1680000e+01 -3.2803460e-20
+ 5.1870000e+01 -1.2973246e-19
+ 5.2060000e+01 -2.1903051e-19
+ 5.2250000e+01 -2.9619265e-19
+ 5.2440000e+01 -3.5732939e-19
+ 5.2630000e+01 -3.9935087e-19
+ 5.2820000e+01 -4.2011991e-19
+ 5.3010000e+01 -4.1856685e-19
+ 5.3200000e+01 -3.9475923e-19
+ 5.3390000e+01 -3.4991992e-19
+ 5.3580000e+01 -2.8638873e-19
+ 5.3770000e+01 -2.0752479e-19
+ 5.3960000e+01 -1.1754995e-19
+ 5.4150000e+01 -2.1337642e-20
+ 5.4340000e+01 7.5844464e-20
+ 5.4530000e+01 1.6862509e-19
+ 5.4720000e+01 2.5183700e-19
+ 5.4910000e+01 3.2082240e-19
+ 5.5100000e+01 3.7171582e-19
+ 5.5290000e+01 4.0168438e-19
+ 5.5480000e+01 4.0910773e-19
+ 5.5670000e+01 3.9368369e-19
+ 5.5860000e+01 3.5645090e-19
+ 5.6050000e+01 2.9972579e-19
+ 5.6240000e+01 2.2695775e-19
+ 5.6430000e+01 1.4251234e-19
+ 5.6620000e+01 5.1397826e-20
+ 5.6810000e+01 -4.1045823e-20
+ 5.7000000e+01 -1.2947405e-19
+ 5.7190000e+01 -2.0885893e-19
+ 5.7380000e+01 -2.7478380e-19
+ 5.7570000e+01 -3.2369325e-19
+ 5.7760000e+01 -3.5308443e-19
+ 5.7950000e+01 -3.6163018e-19
+ 5.8140000e+01 -3.4922950e-19
+ 5.8330000e+01 -3.1698554e-19
+ 5.8520000e+01 -2.6711570e-19
+ 5.8710000e+01 -2.0280178e-19
+ 5.8900000e+01 -1.2799106e-19
+ 5.9090000e+01 -4.7160837e-20
+ 5.9280000e+01 3.4940085e-20
+ 5.9470000e+01 1.1355904e-19
+ 5.9660000e+01 1.8420267e-19
+ 5.9850000e+01 2.4288267e-19
+ 6.0040000e+01 2.8633230e-19
+ 6.0230000e+01 3.1218306e-19
+ 6.0420000e+01 3.1909258e-19
+ 6.0610000e+01 3.0681694e-19
+ 6.0800000e+01 2.7622279e-19
+ 6.0990000e+01 2.2923767e-19
+ 6.1180000e+01 1.6873968e-19
+ 6.1370000e+01 9.8391447e-20
+ 6.1560000e+01 2.2426535e-20
+ 6.1750000e+01 -5.4600362e-20
+ 6.1940000e+01 -1.2808502e-19
+ 6.2130000e+01 -1.9365082e-19
+ 6.2320000e+01 -2.4741406e-19
+ 6.2510000e+01 -2.8622009e-19
+ 6.2700000e+01 -3.0783468e-19
+ 6.2890000e+01 -3.1107739e-19
+ 6.3080000e+01 -2.9588787e-19
+ 6.3270000e+01 -2.6332049e-19
+ 6.3460000e+01 -2.1546775e-19
+ 6.3650000e+01 -1.5531831e-19
+ 6.3840000e+01 -8.6559668e-20
+ 6.4030000e+01 -1.3339828e-20
+ 6.4220000e+01 5.9995411e-20
+ 6.4410000e+01 1.2917147e-19
+ 6.4600000e+01 1.9024237e-19
+ 6.4790000e+01 2.3982016e-19
+ 6.4980000e+01 2.7526252e-19
+ 6.5170000e+01 2.9480847e-19
+ 6.5360000e+01 2.9765587e-19
+ 6.5550000e+01 2.8397889e-19
+ 6.5740000e+01 2.5488724e-19
+ 6.5930000e+01 2.1233279e-19
+ 6.6120000e+01 1.5897190e-19
+ 6.6310000e+01 9.7993561e-20
+ 6.6500000e+01 3.2925061e-20
+ 6.6690000e+01 -3.2573180e-20
+ 6.6880000e+01 -9.4911984e-20
+ 6.7070000e+01 -1.5076633e-19
+ 6.7260000e+01 -1.9724741e-19
+ 6.7450000e+01 -2.3204589e-19
+ 6.7640000e+01 -2.5354037e-19
+ 6.7830000e+01 -2.6086654e-19
+ 6.8020000e+01 -2.5394490e-19
+ 6.8210000e+01 -2.3346641e-19
+ 6.8400000e+01 -2.0083774e-19
+ 6.8590000e+01 -1.5808949e-19
+ 6.8780000e+01 -1.0775270e-19
+ 6.8970000e+01 -5.2710589e-20
+ 6.9160000e+01 3.9661551e-21
+ 6.9350000e+01 5.9191273e-20
+ 6.9540000e+01 1.1003221e-19
+ 6.9730000e+01 1.5386707e-19
+ 6.9920000e+01 1.8852037e-19
+ 7.0110000e+01 2.1236989e-19
+ 7.0300000e+01 2.2441895e-19
+ 7.0490000e+01 2.2433077e-19
+ 7.0680000e+01 2.1242398e-19
+ 7.0870000e+01 1.8963105e-19
+ 7.1060000e+01 1.5742382e-19
+ 7.1250000e+01 1.1771219e-19
+ 7.1440000e+01 7.2723698e-20
+ 7.1630000e+01 2.4872197e-20
+ 7.1820000e+01 -2.3375690e-20
+ 7.2010000e+01 -6.9628508e-20
+ 7.2200000e+01 -1.1168541e-19
+ 7.2390000e+01 -1.4763374e-19
+ 7.2580000e+01 -1.7592525e-19
+ 7.2770000e+01 -1.9542883e-19
+ 7.2960000e+01 -2.0545992e-19
+ 7.3150000e+01 -2.0578729e-19
+ 7.3340000e+01 -1.9661979e-19
+ 7.3530000e+01 -1.7857550e-19
+ 7.3720000e+01 -1.5263679e-19
+ 7.3910000e+01 -1.2009421e-19
+ 7.4100000e+01 -8.2482300e-20
+ 7.4290000e+01 -4.1510205e-20
+ 7.4480000e+01 1.0106876e-21
+ 7.4670000e+01 4.3239679e-20
+ 7.4860000e+01 8.3378339e-20
+ 7.5050000e+01 1.1973878e-19
+ 7.5240000e+01 1.5080677e-19
+ 7.5430000e+01 1.7529785e-19
+ 7.5620000e+01 1.9220480e-19
+ 7.5810000e+01 2.0083475e-19
Added: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/adj_traveltime_filter.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/adj_traveltime_filter.f90 (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/adj_traveltime_filter.f90 2011-04-01 15:57:37 UTC (rev 18161)
@@ -0,0 +1,381 @@
+!!!!! this subroutine is used for preparing adjoint sources in this example ONLY
+!!!!! generally, you should make your own measuremnts and calculate corresponding adjoint sources
+!!!!! FLEXWIN could be of help, or you can modify this code a little bit
+
+program adj_traveltime
+
+implicit none
+
+integer, parameter :: nstep = 799
+double precision, parameter :: dt = 0.19d0
+double precision, parameter :: C_crit = 0.5d0
+
+double precision, parameter :: length_time_window = 60.0d0
+integer, parameter :: taper_type = 1
+integer, parameter :: i_skip_window =1
+double precision, parameter :: factor = 5d-2
+
+
+integer, parameter :: nrec = 1
+double precision, parameter :: pi = 3.141592653589793d0
+
+double precision :: data_origin(nstep,nrec),syn_origin(nstep,nrec),adj(nstep,nrec),adj_density(nstep,nrec)
+double precision :: data_picked(nstep,nrec),syn_picked(nstep,nrec),data_filtered(nstep,nrec),syn_filtered(nstep,nrec)
+double precision :: data_temp(nstep),syn_temp(nstep),data_direct(nstep),syn_direct(nstep)
+double precision :: trace_data_max(nrec), trace_syn_max(nrec), t(nstep),data_max,syn_max
+integer :: flag(nrec), irec, itime, length_window, i_start_window, i_end_window
+character(len=150) :: station_name,file_data,file_syn,file_data_direct,file_syn_direct
+character(len=150) :: file_adj,file_adj_density, file_misfit, file_adj_BHX,file_adj_BHZ
+double precision :: c(nstep), data_trace(nstep), syn_trace(nstep)
+integer :: I(nstep)
+double precision :: Norm_data_crit,Norm_syn_crit,AMP,Norm_data,Norm_syn,Norm,Norm_adj_temp
+double precision :: c_current,c_final(100), w(nstep), misfit_traveltime, traveltime_delay
+integer :: I_current, n_current, index_current, l_current, l, lag
+integer :: I_final(100), index_final(100)
+
+double precision, allocatable :: taper(:),corre(:), adj_temp(:),data_trace_temp(:),syn_trace_temp(:)
+
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+integer,parameter :: filter_flag = 0
+integer :: ifreq, nfreq
+real :: F1,F2,D(8),G,DELT
+!real freq_low(2),freq_high(2)
+!data freq_low / 1.0d-4 , 1.0d-4/
+!data freq_high / 5.0d-1 , 5.0d-2/
+real freq_low(1),freq_high(1)
+data freq_low / 6.67d-3 /
+data freq_high / 3.33d-2 /
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+
+nfreq=size(freq_low)
+
+file_misfit = './OUTPUT_FILES/misfit_traveltime_delay'
+open(unit=1111,file=trim(file_misfit),status='unknown')
+
+DELT=dt*1.0d3
+adj=0.0d0
+misfit_traveltime = 0.0d0
+!!!! loading data and synthetics !!!!
+do irec = 1,nrec
+ file_data = './SEM/A7.II.MXZ.sem.ascii'
+ open(unit=1001,file=trim(file_data),status='old',action='read')
+ do itime = 1,nstep
+ !read(1001,*) t(itime),data_origin(itime,irec) ! original
+ read(1001,*) t(itime),data_origin(nstep-itime+1,irec) ! reversed
+ end do
+ close(1001)
+
+ data_temp(1)=0.0
+ data_temp(nstep)=0.0
+ do itime = 2,nstep-1
+ data_temp(itime)=( data_origin(itime+1,irec) - data_origin(itime-1,irec) )/(2*dt)
+ end do
+ data_origin(:,irec)=data_temp
+ trace_data_max(irec)=maxval(abs(data_temp))
+end do
+do itime = 1,nstep
+ t(itime) = t(itime)-t(1)
+end do
+data_max=maxval(trace_data_max)
+
+!!!! taper !!!!
+length_window = floor(length_time_window/dt)+1
+allocate(taper(length_window)); taper(:)=1.0d0
+allocate(adj_temp(length_window)); adj_temp(:)=0.0d0
+allocate(corre(2*length_window-1)); corre(:)=0.0d0
+allocate(data_trace_temp(length_window)); data_trace_temp(:)=0.0d0
+allocate(syn_trace_temp(length_window)); syn_trace_temp(:)=0.0d0
+if (taper_type .eq. 1) then ! cosine taper, otherwise using a constant (1.0) instead
+ do l=1,length_window
+ taper(l) = (1.0-cos(pi*2.0*(l-1)/(length_window-1)))/2.0
+ end do
+end if
+
+!!!! computing adj sources !!!!
+ if (filter_flag .eq. 0) then
+ nfreq=1
+ end if
+do irec = 1,nrec
+ do ifreq = 1,nfreq
+ data_filtered(:,irec)=data_origin(:,irec)
+ syn_filtered(:,irec)= syn_origin(:,irec)
+ !!!!!!!!!!!!!!!!!! to be implemented !!!!!!!!!!!!!!!!!
+ if (filter_flag .eq. 1) then
+ ! THIS SECTION CALCULATES THE FILTER AND MUST BE CALLED BEFORE
+ ! FILTER IS CALLED
+ F1=freq_low(ifreq)
+ F2=freq_high(ifreq)
+ call BNDPAS(F1,F2,DELT,D,G,nstep)
+ ! F1 = LOW FREQUENCY CUTOFF (6 DB DOWN)
+ ! F2 = HIGH FREQUENCY CUTOFF (6 DB DOWN)
+ ! DELT = SAMPLE INTERVAL IN MILLISECONDS
+ ! D = WILL CONTAIN 8 Z DOMAIN COEFICIENTS OF RECURSIVE FILTER
+ ! G = WILL CONTAIN THE GAIN OF THE FILTER,
+ call FILTER(data_filtered(:,irec),nstep,D,G,2)
+ call FILTER(syn_filtered(:,irec),nstep,D,G,2)
+ ! X = DATA VECTOR OF LENGTH N CONTAINING DATA TO BE FILTERED
+ ! D = FILTER COEFFICIENTS CALCULATED BY BNDPAS
+ ! G = FILTER GAIN
+ ! IG = 1 one pass
+ ! ig = 2 two passes
+ end if
+
+
+
+
+
+ c(:)=0.0
+ I(:)=(length_window-1)
+ data_trace=data_filtered(:,irec)
+ syn_trace= syn_filtered(:,irec)
+! if ( ( trace_data_max(irec) > (factor*data_max) ) .and. (trace_syn_max(irec) > (factor*syn_max)) ) then
+ if (1 .eq. 1) then
+ !!!! cross-correlation !!!!
+ do i_start_window = (nstep+1)/2, (nstep+1)/2 !the positive branch
+! do i_start_window = (nstep+1)/2 - length_window, (nstep+1)/2 - length_window !the negative branch
+ i_end_window = i_start_window + length_window - 1
+
+ data_trace_temp=data_trace(i_start_window:i_end_window)
+
+ Norm_data_crit=sqrt(DOT_PRODUCT(data_trace,data_trace))
+ AMP = Norm_data_crit
+
+ Norm_data=sqrt(DOT_PRODUCT( data_trace_temp,data_trace_temp ))
+ Norm = Norm_data*Norm_data
+
+ w(:)=0.0
+ w(i_start_window : i_end_window)=taper
+ end do
+ end if
+ data_trace=data_trace*w
+ data_filtered(:,irec)=data_trace
+ !!!! normal adjoint sources !!!!
+ adj_temp=data_filtered(i_start_window : i_end_window,irec)
+ Norm_adj_temp = - DOT_PRODUCT(adj_temp,adj_temp)*dt ! minus sign comes from integration by part !
+ !!!! choose which one to use !!!!
+ adj_temp=adj_temp/Norm_adj_temp ! normal adjoint sources
+ ! adj_temp=adj_temp/Norm_adj_temp ! ray density map DeltaT=+1
+ ! adj_temp=adj_temp/Norm_adj_temp * ((-1)**irec) ! ray density map DeltaT=+-1
+
+ adj(i_start_window : i_end_window ,irec )=adj_temp
+ traveltime_delay=dt
+ misfit_traveltime = misfit_traveltime + traveltime_delay * traveltime_delay / 2.0
+ write(1111,*) traveltime_delay
+
+ end do !do ifreq=1,nfreq
+
+ !!!! output !!!!
+ file_adj_BHZ = './SEM/A7.II.MXZ.adj'
+ open(unit=1002,file=trim(file_adj_BHZ),status='unknown')
+ do itime = 1,nstep
+ write(1002,*) t(itime), adj(nstep-itime+1,irec)
+! write(1002,*) t(itime), adj(itime,irec)
+ end do
+ close(1002)
+ file_adj_BHZ = './SEM/A7.II.MXZ.adj_reversed'
+ open(unit=1002,file=trim(file_adj_BHZ),status='unknown')
+ do itime = 1,nstep
+! write(1002,*) t(itime), adj(nstep-itime+1,irec)
+ write(1002,*) t(itime), adj(itime,irec)
+ end do
+ close(1002)
+
+
+ file_adj_BHZ = './SEM/A7.II.MXE.adj'
+ open(unit=1002,file=trim(file_adj_BHZ),status='unknown')
+ do itime = 1,nstep
+ write(1002,*) t(itime), 0.0
+ end do
+ close(1002)
+ file_adj_BHZ = './SEM/A7.II.MXN.adj'
+ open(unit=1002,file=trim(file_adj_BHZ),status='unknown')
+ do itime = 1,nstep
+ write(1002,*) t(itime), 0.0
+ end do
+ close(1002)
+
+enddo
+
+
+close(1111)
+
+file_misfit = './OUTPUT_FILES/misfit_traveltime'
+open(unit=1001,file=trim(file_misfit),status='unknown')
+write(1001,*) misfit_traveltime, traveltime_delay
+close(1001)
+
+end program adj_traveltime
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+SUBROUTINE BNDPAS(F1,F2,DELT,D,G,N)
+! RECURSIVE BUTTERWORTH BAND PASS FILTER (KANASEWICH, TIME SERIES
+! ANALYSIS IN GEOPHYSICS, UNIVERSITY OF ALBERTA PRESS, 1975; SHANKS,
+! JOHN L, RECURSION FILTERS FOR DIGITAL PROCESSING, GEOPHYSICS, V32,
+! FILTER. THE FILTER WILL HAVE 8 POLES IN THE S PLANE AND IS
+! APPLIED IN FORWARD AND REVERSE DIRECTIONS SO AS TO HAVE ZERO
+! PHASE SHIFT. THE GAIN AT THE TWO FREQUENCIES SPECIFIED AS
+! CUTOFF FREQUENCIES WILL BE -6DB AND THE ROLLOFF WILL BE ABOUT
+! THE FILTER TO PREVENT ALIASING PROBLEMS.
+ COMPLEX P(4),S(8),Z1,Z2
+ real D(8),XC(3),XD(3),XE(3)
+ double precision :: X(N)
+ DATA ISW/0/,TWOPI/6.2831853/
+! THIS SECTION CALCULATES THE FILTER AND MUST BE CALLED BEFORE
+! FILTER IS CALLED
+
+! F1 = LOW FREQUENCY CUTOFF (6 DB DOWN)
+! F2 = HIGH FREQUENCY CUTOFF (6 DB DOWN)
+! DELT = SAMPLE INTERVAL IN MILLISECONDS
+! D = WILL CONTAIN 8 Z DOMAIN COEFICIENTS OF RECURSIVE FILTER
+! G = WILL CONTAIN THE GAIN OF THE FILTER,
+
+ DT=DELT/1000.0
+ TDT=2.0/DT
+ FDT=4.0/DT
+ ISW=1
+ P(1)=CMPLX(-.3826834,.9238795)
+ P(2)=CMPLX(-.3826834,-.9238795)
+ P(3)=CMPLX(-.9238795,.3826834)
+ P(4)=CMPLX(-.9238795,-.3826834)
+ W1=TWOPI*F1
+ W2=TWOPI*F2
+ W1=TDT*TAN(W1/TDT)
+ W2=TDT*TAN(W2/TDT)
+ HWID=(W2-W1)/2.0
+ WW=W1*W2
+ DO 19 I=1,4
+ Z1=P(I)*HWID
+ Z2=Z1*Z1-WW
+ Z2=CSQRT(Z2)
+ S(I)=Z1+Z2
+ 19 S(I+4)=Z1-Z2
+ G=.5/HWID
+ G=G*G
+ G=G*G
+ DO 29 I=1,7,2
+ B=-2.0*REAL(S(I))
+ Z1=S(I)*S(I+1)
+ C=REAL(Z1)
+ A=TDT+B+C/TDT
+ G=G*A
+ D(I)=(C*DT-FDT)/A
+ 29 D(I+1)=(A-2.0*B)/A
+ G=G*G
+ 5 FORMAT ('-FILTER GAIN IS ', 9E12.6)
+ RETURN
+
+ ENTRY FILTER(X,N,D,G,IG)
+
+! X = DATA VECTOR OF LENGTH N CONTAINING DATA TO BE FILTERED
+! D = FILTER COEFFICIENTS CALCULATED BY BNDPAS
+! G = FILTER GAIN
+! IG = 1 one pass
+! ig = 2 two passes
+
+ IF (ISW.EQ.1) GO TO 31
+ WRITE (6,6)
+ 6 FORMAT ('1BNDPAS MUST BE CALLED BEFORE FILTER')
+ return
+
+! APPLY FILTER IN FORWARD DIRECTION
+
+ 31 XM2=X(1)
+ XM1=X(2)
+ XM=X(3)
+ XC(1)=XM2
+ XC(2)=XM1-D(1)*XC(1)
+ XC(3)=XM-XM2-D(1)*XC(2)-D(2)*XC(1)
+ XD(1)=XC(1)
+ XD(2)=XC(2)-D(3)*XD(1)
+ XD(3)=XC(3)-XC(1)-D(3)*XD(2)-D(4)*XD(1)
+ XE(1)=XD(1)
+ XE(2)=XD(2)-D(5)*XE(1)
+ XE(3)=XD(3)-XD(1)-D(5)*XE(2)-D(6)*XE(1)
+ X(1)=XE(1)
+ X(2)=XE(2)-D(7)*X(1)
+ X(3)=XE(3)-XE(1)-D(7)*X(2)-D(8)*X(1)
+ DO 39 I=4,N
+ XM2=XM1
+ XM1=XM
+ XM=X(I)
+ K=I-((I-1)/3)*3
+ GO TO (34,35,36),K
+ 34 M=1
+ M1=3
+ M2=2
+ GO TO 37
+ 35 M=2
+ M1=1
+ M2=3
+ GO TO 37
+ 36 M=3
+ M1=2
+ M2=1
+ 37 XC(M)=XM-XM2-D(1)*XC(M1)-D(2)*XC(M2)
+ XD(M)=XC(M)-XC(M2)-D(3)*XD(M1)-D(4)*XD(M2)
+ XE(M)=XD(M)-XD(M2)-D(5)*XE(M1)-D(6)*XE(M2)
+ 39 X(I)=XE(M)-XE(M2)-D(7)*X(I-1)-D(8)*X(I-2)
+!
+!
+ if(ig.eq.1) goto 3333
+ XM2=X(N)
+ XM1=X(N-1)
+ XM=X(N-2)
+ XC(1)=XM2
+ XC(2)=XM1-D(1)*XC(1)
+ XC(3)=XM-XM2-D(1)*XC(2)-D(2)*XC(1)
+ XD(1)=XC(1)
+ XD(2)=XC(2)-D(3)*XD(1)
+ XD(3)=XC(3)-XC(1)-D(3)*XD(2)-D(4)*XD(1)
+ XE(1)=XD(1)
+ XE(2)=XD(2)-D(5)*XE(1)
+ XE(3)=XD(3)-XD(1)-D(5)*XE(2)-D(6)*XE(1)
+ X(N)=XE(1)
+ X(N-1)=XE(2)-D(7)*X(1)
+ X(N-2)=XE(3)-XE(1)-D(7)*X(2)-D(8)*X(1)
+ DO 49 I=4,N
+ XM2=XM1
+ XM1=XM
+ J=N-I+1
+ XM=X(J)
+ K=I-((I-1)/3)*3
+ GO TO (44,45,46),K
+ 44 M=1
+ M1=3
+ M2=2
+ GO TO 47
+ 45 M=2
+ M1=1
+ M2=3
+ GO TO 47
+ 46 M=3
+ M1=2
+ M2=1
+ 47 XC(M)=XM-XM2-D(1)*XC(M1)-D(2)*XC(M2)
+ XD(M)=XC(M)-XC(M2)-D(3)*XD(M1)-D(4)*XD(M2)
+ XE(M)=XD(M)-XD(M2)-D(5)*XE(M1)-D(6)*XE(M2)
+ 49 X(J)=XE(M)-XE(M2)-D(7)*X(J+1)-D(8)*X(J+2)
+ 3333 continue
+ if (ig.eq.1) then
+ gg=sqrt(g) ! if only pass once, modify gain
+ else
+ gg=g
+ endif
+ DO 59 I=1,N
+ 59 X(I)=X(I)/gg
+ RETURN
+END
+
Added: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/irec_master_noise
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/irec_master_noise (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/irec_master_noise 2011-04-01 15:57:37 UTC (rev 18161)
@@ -0,0 +1 @@
+1
Added: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/nu_master
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/nu_master (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/test_regional/nu_master 2011-04-01 15:57:37 UTC (rev 18161)
@@ -0,0 +1,3 @@
+0
+0
+1
More information about the CIG-COMMITS
mailing list