[cig-commits] [commit] devel, master: two very short tests for NOISE_TOMOGRAPHY, just to make sure scripts work before submitting the real jobs (6b400df)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Thu Nov 6 08:04:48 PST 2014
Repository : https://github.com/geodynamics/specfem3d_globe
On branches: devel,master
Link : https://github.com/geodynamics/specfem3d_globe/compare/bc58e579b3b0838a0968725a076f5904845437ca...be63f20cbb6f462104e949894dbe205d2398cd7f
>---------------------------------------------------------------
commit 6b400df4565ad3610e0e4c96a1273831ade50946
Author: Yang Luo <yangl at princeton.edu>
Date: Fri Apr 1 15:55:46 2011 +0000
two very short tests for NOISE_TOMOGRAPHY, just to make sure scripts work before submitting the real jobs
>---------------------------------------------------------------
6b400df4565ad3610e0e4c96a1273831ade50946
noise_examples/global_long/NOISE123.submit_atten | 20 ++++++++++++++++++++
noise_examples/global_long/adj_traveltime_filter.f90 | 4 ++--
noise_examples/global_short/NOISE123.submit_atten | 20 ++++++++++++++++++++
.../global_short/adj_traveltime_filter.f90 | 2 +-
noise_examples/pre-processing | 6 ++++++
noise_examples/regional/NOISE123.submit_atten | 20 ++++++++++++++++++++
noise_examples/regional/NOISE_clean_create | 1 -
noise_examples/regional/adj_traveltime_filter.f90 | 2 +-
8 files changed, 70 insertions(+), 5 deletions(-)
diff --git a/noise_examples/global_long/NOISE123.submit_atten b/noise_examples/global_long/NOISE123.submit_atten
index 1a51ecb..331b633 100755
--- a/noise_examples/global_long/NOISE123.submit_atten
+++ b/noise_examples/global_long/NOISE123.submit_atten
@@ -6,6 +6,7 @@
#PBS -l nodes=75:ppn=8,walltime=100:00:00
DIR_RESULT="zzz_global_long"
+NSTEP=31599
cd $PBS_O_WORKDIR
######################################## pre-simulation ###########################################
@@ -43,6 +44,8 @@ mkdir -p OUTPUT_FILES
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
@@ -55,6 +58,8 @@ cp DATA/STATIONS $DIR_TEMP/
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
@@ -64,11 +69,14 @@ 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
@@ -81,6 +89,8 @@ cp DATA/STATIONS $DIR_TEMP/
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
@@ -93,6 +103,8 @@ cp DATA/STATIONS $DIR_TEMP/
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`
@@ -100,6 +112,14 @@ 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
diff --git a/noise_examples/global_long/adj_traveltime_filter.f90 b/noise_examples/global_long/adj_traveltime_filter.f90
index 110fe24..a6031b1 100644
--- a/noise_examples/global_long/adj_traveltime_filter.f90
+++ b/noise_examples/global_long/adj_traveltime_filter.f90
@@ -6,7 +6,7 @@ program adj_traveltime
implicit none
-integer, parameter :: nstep = 31599
+integer, parameter :: nstep = 129599
double precision, parameter :: dt = 0.19d0
double precision, parameter :: C_crit = 0.5d0
@@ -57,7 +57,7 @@ adj=0.0d0
misfit_traveltime = 0.0d0
!!!! loading data and synthetics !!!!
do irec = 1,nrec
- file_data = './zzz_global_long/step_2/A7.II.MXZ.sem.ascii'
+ 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
diff --git a/noise_examples/global_short/NOISE123.submit_atten b/noise_examples/global_short/NOISE123.submit_atten
index d16d097..08ab721 100755
--- a/noise_examples/global_short/NOISE123.submit_atten
+++ b/noise_examples/global_short/NOISE123.submit_atten
@@ -6,6 +6,7 @@
#PBS -l nodes=75:ppn=8,walltime=100:00:00
DIR_RESULT="zzz_global_short"
+NSTEP=31599
cd $PBS_O_WORKDIR
######################################## pre-simulation ###########################################
@@ -43,6 +44,8 @@ mkdir -p OUTPUT_FILES
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
@@ -55,6 +58,8 @@ cp DATA/STATIONS $DIR_TEMP/
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
@@ -64,11 +69,14 @@ 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
@@ -81,6 +89,8 @@ cp DATA/STATIONS $DIR_TEMP/
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
@@ -93,6 +103,8 @@ cp DATA/STATIONS $DIR_TEMP/
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`
@@ -100,6 +112,14 @@ 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
diff --git a/noise_examples/global_short/adj_traveltime_filter.f90 b/noise_examples/global_short/adj_traveltime_filter.f90
index 079d72e..8935bc2 100644
--- a/noise_examples/global_short/adj_traveltime_filter.f90
+++ b/noise_examples/global_short/adj_traveltime_filter.f90
@@ -57,7 +57,7 @@ adj=0.0d0
misfit_traveltime = 0.0d0
!!!! loading data and synthetics !!!!
do irec = 1,nrec
- file_data = './zzz_global_short/step_2/A7.II.MXZ.sem.ascii'
+ 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
diff --git a/noise_examples/pre-processing b/noise_examples/pre-processing
index 27314de..9ca3088 100755
--- a/noise_examples/pre-processing
+++ b/noise_examples/pre-processing
@@ -32,14 +32,20 @@ cp DATA/Par_file_NOISE_1_attenuation DATA/Par_file
rm -rf OUTPUT_FILES/*
make xmeshfem3D
make xspecfem3D
+make xcreate_movie_AVS_DX
mv bin/xmeshfem3D xmeshfem3D_attenuation
mv bin/xspecfem3D xspecfem3D_attenuation
+mv bin/xcreate_movie_AVS_DX xcreate_movie_AVS_DX_attenuation
cp DATA/Par_file_NOISE_1_noattenuation DATA/Par_file
rm -rf OUTPUT_FILES/*
make xmeshfem3D
make xspecfem3D
+make xcreate_movie_AVS_DX
+make xcombine_vol_data
mv bin/xmeshfem3D xmeshfem3D_noattenuation
mv bin/xspecfem3D xspecfem3D_noattenuation
+mv bin/xcreate_movie_AVS_DX xcreate_movie_AVS_DX_noattenuation
+mv bin/xcombine_vol_data xcombine_vol_data
qsub NOISE123.submit_atten
diff --git a/noise_examples/regional/NOISE123.submit_atten b/noise_examples/regional/NOISE123.submit_atten
index 64b291b..3e00fee 100755
--- a/noise_examples/regional/NOISE123.submit_atten
+++ b/noise_examples/regional/NOISE123.submit_atten
@@ -6,6 +6,7 @@
#PBS -l nodes=50:ppn=8,walltime=100:00:00
DIR_RESULT="zzz_regional"
+NSTEP=31599
cd $PBS_O_WORKDIR
######################################## pre-simulation ###########################################
@@ -43,6 +44,8 @@ mkdir -p OUTPUT_FILES
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
@@ -55,6 +58,8 @@ cp DATA/STATIONS $DIR_TEMP/
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
@@ -64,11 +69,14 @@ 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
@@ -81,6 +89,8 @@ cp DATA/STATIONS $DIR_TEMP/
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
@@ -93,6 +103,8 @@ cp DATA/STATIONS $DIR_TEMP/
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`
@@ -100,6 +112,14 @@ 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
diff --git a/noise_examples/regional/NOISE_clean_create b/noise_examples/regional/NOISE_clean_create
index bb03b61..f3231e0 100755
--- a/noise_examples/regional/NOISE_clean_create
+++ b/noise_examples/regional/NOISE_clean_create
@@ -3,7 +3,6 @@
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
diff --git a/noise_examples/regional/adj_traveltime_filter.f90 b/noise_examples/regional/adj_traveltime_filter.f90
index ed68bd7..8935bc2 100644
--- a/noise_examples/regional/adj_traveltime_filter.f90
+++ b/noise_examples/regional/adj_traveltime_filter.f90
@@ -57,7 +57,7 @@ adj=0.0d0
misfit_traveltime = 0.0d0
!!!! loading data and synthetics !!!!
do irec = 1,nrec
- file_data = './zzz_regional/step_2/A7.II.MXZ.sem.ascii'
+ 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
More information about the CIG-COMMITS
mailing list