[cig-commits] [commit] devel, master: fix EXAMPLES/noise_tomography to fit updates made in other subroutines. write SURFACE_MOVIE using the same c rountine for absorbing boundaries (presumably much faster and more effective, since there is only one file on each PROC for the whole simulation instead of one file for each time step) (cfee3ec)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Thu Nov 6 08:04:35 PST 2014
Repository : https://github.com/geodynamics/specfem3d_globe
On branches: devel,master
Link : https://github.com/geodynamics/specfem3d_globe/compare/bc58e579b3b0838a0968725a076f5904845437ca...be63f20cbb6f462104e949894dbe205d2398cd7f
>---------------------------------------------------------------
commit cfee3ec4fab14e22a97d5a9b68f354b7bbfe3462
Author: Yang Luo <yangl at princeton.edu>
Date: Tue Mar 22 20:44:44 2011 +0000
fix EXAMPLES/noise_tomography to fit updates made in other subroutines. write SURFACE_MOVIE using the same c rountine for absorbing boundaries (presumably much faster and more effective, since there is only one file on each PROC for the whole simulation instead of one file for each time step)
>---------------------------------------------------------------
cfee3ec4fab14e22a97d5a9b68f354b7bbfe3462
noise_examples/global_long/adj_traveltime_filter.f90 | 10 +++++-----
noise_examples/global_short/adj_traveltime_filter.f90 | 10 +++++-----
noise_examples/pre-processing | 14 +++++++-------
noise_examples/regional/adj_traveltime_filter.f90 | 6 +++---
4 files changed, 20 insertions(+), 20 deletions(-)
diff --git a/noise_examples/global_long/adj_traveltime_filter.f90 b/noise_examples/global_long/adj_traveltime_filter.f90
index 02e3ec6..e6dd063 100644
--- a/noise_examples/global_long/adj_traveltime_filter.f90
+++ b/noise_examples/global_long/adj_traveltime_filter.f90
@@ -54,11 +54,11 @@ adj=0.0d0
misfit_traveltime = 0.0d0
!!!! loading data and synthetics !!!!
do irec = 1,nrec
- file_data = './ZZZ_2/A7.II.LHZ.sem.ascii'
+ file_data = './ZZZ_2/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(nstep-itime+1,irec) ! original
- read(1001,*) t(itime),data_origin(nstep-itime+1,irec) ! reversed
+ !read(1001,*) t(itime),data_origin(itime,irec) ! original
+ read(1001,*) t(itime),data_origin(nstep-itime+1,irec) ! reversed
end do
close(1001)
@@ -163,14 +163,14 @@ do irec = 1,nrec
end do !do ifreq=1,nfreq
!!!! output !!!!
- file_adj_BHZ = './SEM/A7.II.LHZ.adj'
+ 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.LHZ.adj_reversed'
+ 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)
diff --git a/noise_examples/global_short/adj_traveltime_filter.f90 b/noise_examples/global_short/adj_traveltime_filter.f90
index 744b95d..e6dd063 100644
--- a/noise_examples/global_short/adj_traveltime_filter.f90
+++ b/noise_examples/global_short/adj_traveltime_filter.f90
@@ -54,11 +54,11 @@ adj=0.0d0
misfit_traveltime = 0.0d0
!!!! loading data and synthetics !!!!
do irec = 1,nrec
- file_data = './ZZZ_2/A7.II.LHZ.sem.ascii'
+ file_data = './ZZZ_2/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
+ !read(1001,*) t(itime),data_origin(itime,irec) ! original
+ read(1001,*) t(itime),data_origin(nstep-itime+1,irec) ! reversed
end do
close(1001)
@@ -163,14 +163,14 @@ do irec = 1,nrec
end do !do ifreq=1,nfreq
!!!! output !!!!
- file_adj_BHZ = './SEM/A7.II.LHZ.adj'
+ 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.LHZ.adj_reversed'
+ 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)
diff --git a/noise_examples/pre-processing b/noise_examples/pre-processing
index adfa886..86fd0e2 100755
--- a/noise_examples/pre-processing
+++ b/noise_examples/pre-processing
@@ -1,4 +1,4 @@
-#!/bin/bash
+#!/bin/bash -eu
example=$1
dir=$PWD/$example
@@ -9,10 +9,10 @@ echo ""; echo "";
cd ../../
-mkdir -p SEM zzz_job_info
+mkdir -p SEM zzz_job_info NOISE_TOMOGRAPHY
-cp $dir/A0.II.LHZ.adj ./SEM/
-cp $dir/cp_adj_sources ./SEM/
+#cp $dir/A0.II.LHZ.adj ./SEM/
+#cp $dir/cp_adj_sources ./SEM/
cp $dir/adj_traveltime_filter.f90 ./
cp $dir/NOISE123.submit_atten ./
@@ -28,9 +28,9 @@ cp $dir/Par_file_NOISE_* ./DATA/
cp $dir/CMTSOLUTION_NOISE ./DATA/
cp $dir/STATIONS_NOISE ./DATA/
-cd SEM
-./cp_adj_sources
-cd ..
+#cd SEM
+#./cp_adj_sources
+#cd ..
ifort -o NOISE_adj adj_traveltime_filter.f90
diff --git a/noise_examples/regional/adj_traveltime_filter.f90 b/noise_examples/regional/adj_traveltime_filter.f90
index b4b8c15..e6dd063 100644
--- a/noise_examples/regional/adj_traveltime_filter.f90
+++ b/noise_examples/regional/adj_traveltime_filter.f90
@@ -54,7 +54,7 @@ adj=0.0d0
misfit_traveltime = 0.0d0
!!!! loading data and synthetics !!!!
do irec = 1,nrec
- file_data = './ZZZ_2/A7.II.LHZ.sem.ascii'
+ file_data = './ZZZ_2/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
@@ -163,14 +163,14 @@ do irec = 1,nrec
end do !do ifreq=1,nfreq
!!!! output !!!!
- file_adj_BHZ = './SEM/A7.II.LHZ.adj'
+ 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.LHZ.adj_reversed'
+ 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)
More information about the CIG-COMMITS
mailing list