[cig-commits] r18160 - in seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples: . global_long global_short regional
yangl at geodynamics.org
yangl at geodynamics.org
Fri Apr 1 08:55:46 PDT 2011
Author: yangl
Date: 2011-04-01 08:55:46 -0700 (Fri, 01 Apr 2011)
New Revision: 18160
Modified:
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_long/NOISE123.submit_atten
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_long/adj_traveltime_filter.f90
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_short/NOISE123.submit_atten
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_short/adj_traveltime_filter.f90
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/pre-processing
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/regional/NOISE123.submit_atten
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/regional/NOISE_clean_create
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/regional/adj_traveltime_filter.f90
Log:
two very short tests for NOISE_TOMOGRAPHY, just to make sure scripts work before submitting the real jobs
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_long/NOISE123.submit_atten
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_long/NOISE123.submit_atten 2011-04-01 14:32:55 UTC (rev 18159)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_long/NOISE123.submit_atten 2011-04-01 15:55:46 UTC (rev 18160)
@@ -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 @@
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/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/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/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/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 @@
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
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_long/adj_traveltime_filter.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_long/adj_traveltime_filter.f90 2011-04-01 14:32:55 UTC (rev 18159)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_long/adj_traveltime_filter.f90 2011-04-01 15:55:46 UTC (rev 18160)
@@ -6,7 +6,7 @@
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 @@
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
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_short/NOISE123.submit_atten
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_short/NOISE123.submit_atten 2011-04-01 14:32:55 UTC (rev 18159)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_short/NOISE123.submit_atten 2011-04-01 15:55:46 UTC (rev 18160)
@@ -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 @@
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/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/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/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/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 @@
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
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_short/adj_traveltime_filter.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_short/adj_traveltime_filter.f90 2011-04-01 14:32:55 UTC (rev 18159)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/global_short/adj_traveltime_filter.f90 2011-04-01 15:55:46 UTC (rev 18160)
@@ -57,7 +57,7 @@
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
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/pre-processing
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/pre-processing 2011-04-01 14:32:55 UTC (rev 18159)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/pre-processing 2011-04-01 15:55:46 UTC (rev 18160)
@@ -32,14 +32,20 @@
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
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/regional/NOISE123.submit_atten
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/regional/NOISE123.submit_atten 2011-04-01 14:32:55 UTC (rev 18159)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/regional/NOISE123.submit_atten 2011-04-01 15:55:46 UTC (rev 18160)
@@ -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 @@
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/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/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/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/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 @@
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
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/regional/NOISE_clean_create
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/regional/NOISE_clean_create 2011-04-01 14:32:55 UTC (rev 18159)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/regional/NOISE_clean_create 2011-04-01 15:55:46 UTC (rev 18160)
@@ -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
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/regional/adj_traveltime_filter.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/regional/adj_traveltime_filter.f90 2011-04-01 14:32:55 UTC (rev 18159)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/noise_examples/regional/adj_traveltime_filter.f90 2011-04-01 15:55:46 UTC (rev 18160)
@@ -57,7 +57,7 @@
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