[cig-commits] [commit] devel, master: change number of files defined in src/shared/write_c_binary.c (from 9 to 10), new index 9 corresponds to SURFACE_MOVIE used in NOISE_TOMOGRAPHY (d773215)

cig_noreply at geodynamics.org cig_noreply at geodynamics.org
Thu Nov 6 08:04:39 PST 2014 

Repository : https://github.com/geodynamics/specfem3d_globe

On branches: devel,master
Link       : https://github.com/geodynamics/specfem3d_globe/compare/bc58e579b3b0838a0968725a076f5904845437ca...be63f20cbb6f462104e949894dbe205d2398cd7f

>---------------------------------------------------------------

commit d7732153e2afd6d37e9f43d68627003c0dc4cc2f
Author: Yang Luo <yangl at princeton.edu>
Date:   Thu Mar 24 19:43:00 2011 +0000

    change number of files defined in src/shared/write_c_binary.c (from 9 to 10), new index 9 corresponds to SURFACE_MOVIE used in NOISE_TOMOGRAPHY


>---------------------------------------------------------------

d7732153e2afd6d37e9f43d68627003c0dc4cc2f
 .../global_long/adj_traveltime_filter.f90          | 39 +++++++++++-----------
 .../global_short/adj_traveltime_filter.f90         | 39 +++++++++++-----------
 noise_examples/pre-processing                      |  4 +--
 noise_examples/regional/adj_traveltime_filter.f90  | 39 +++++++++++-----------
 4 files changed, 62 insertions(+), 59 deletions(-)

diff --git a/noise_examples/global_long/adj_traveltime_filter.f90 b/noise_examples/global_long/adj_traveltime_filter.f90
index e6dd063..3fde181 100644
--- a/noise_examples/global_long/adj_traveltime_filter.f90
+++ b/noise_examples/global_long/adj_traveltime_filter.f90
@@ -1,3 +1,7 @@
+!!!!! 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
@@ -177,25 +181,22 @@ do irec = 1,nrec
       write(1002,*) t(itime), adj(itime,irec)
    end do
    close(1002)
-!   file_data         = './yang_data/'//trim(station_name)//'.AA.PRE.semp_picked'
-!   file_syn          = './yang_syn/'//trim(station_name)//'.AA.PRE.semp_picked'
-!   file_adj          = './yang_data/'//trim(station_name)//'.AA.PRE.semp_adj'
-!   file_adj_density  = './yang_syn/'//trim(station_name)//'.AA.PRE.semp_adjden'
-!   open(unit=1001,file=trim(file_data),status='unknown')
-!   open(unit=1002,file=trim(file_syn),status='unknown')
-!   open(unit=1003,file=trim(file_adj),status='unknown')
-!   open(unit=1004,file=trim(file_adj_density),status='unknown')
-!   do itime = 1,nstep
-!      write(1001,*) t(itime), data_origin(itime,irec)
-!      write(1002,*) t(itime), syn_origin(itime,irec)
-!      write(1003,*) t(itime), adj(itime,irec)
-!      write(1004,*) t(itime), adj_density(itime,irec)
-!   end do
-!   close(1001)
-!   close(1002)
-!   close(1003)
-!   close(1004)
-   enddo 
+
+
+   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)
diff --git a/noise_examples/global_short/adj_traveltime_filter.f90 b/noise_examples/global_short/adj_traveltime_filter.f90
index e6dd063..3fde181 100644
--- a/noise_examples/global_short/adj_traveltime_filter.f90
+++ b/noise_examples/global_short/adj_traveltime_filter.f90
@@ -1,3 +1,7 @@
+!!!!! 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
@@ -177,25 +181,22 @@ do irec = 1,nrec
       write(1002,*) t(itime), adj(itime,irec)
    end do
    close(1002)
-!   file_data         = './yang_data/'//trim(station_name)//'.AA.PRE.semp_picked'
-!   file_syn          = './yang_syn/'//trim(station_name)//'.AA.PRE.semp_picked'
-!   file_adj          = './yang_data/'//trim(station_name)//'.AA.PRE.semp_adj'
-!   file_adj_density  = './yang_syn/'//trim(station_name)//'.AA.PRE.semp_adjden'
-!   open(unit=1001,file=trim(file_data),status='unknown')
-!   open(unit=1002,file=trim(file_syn),status='unknown')
-!   open(unit=1003,file=trim(file_adj),status='unknown')
-!   open(unit=1004,file=trim(file_adj_density),status='unknown')
-!   do itime = 1,nstep
-!      write(1001,*) t(itime), data_origin(itime,irec)
-!      write(1002,*) t(itime), syn_origin(itime,irec)
-!      write(1003,*) t(itime), adj(itime,irec)
-!      write(1004,*) t(itime), adj_density(itime,irec)
-!   end do
-!   close(1001)
-!   close(1002)
-!   close(1003)
-!   close(1004)
-   enddo 
+
+
+   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)
diff --git a/noise_examples/pre-processing b/noise_examples/pre-processing
index efe55e8..0111d50 100755
--- a/noise_examples/pre-processing
+++ b/noise_examples/pre-processing
@@ -29,14 +29,14 @@ cp $dir/STATIONS_NOISE_ADJOINT           ./DATA/STATIONS_ADJOINT
 ifort -o NOISE_adj adj_traveltime_filter.f90
 
 cp DATA/Par_file_NOISE_1_attenuation DATA/Par_file
-rm OUTPUT_FILES/*
+rm -rf OUTPUT_FILES/*
 make xmeshfem3D
 make xspecfem3D
 mv bin/xmeshfem3D xmeshfem3D_attenuation
 mv bin/xspecfem3D xspecfem3D_attenuation
 
 cp DATA/Par_file_NOISE_1_noattenuation DATA/Par_file
-rm OUTPUT_FILES/*
+rm -rf OUTPUT_FILES/*
 make xmeshfem3D
 make xspecfem3D
 mv bin/xmeshfem3D xmeshfem3D_noattenuation
diff --git a/noise_examples/regional/adj_traveltime_filter.f90 b/noise_examples/regional/adj_traveltime_filter.f90
index e6dd063..3fde181 100644
--- a/noise_examples/regional/adj_traveltime_filter.f90
+++ b/noise_examples/regional/adj_traveltime_filter.f90
@@ -1,3 +1,7 @@
+!!!!! 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
@@ -177,25 +181,22 @@ do irec = 1,nrec
       write(1002,*) t(itime), adj(itime,irec)
    end do
    close(1002)
-!   file_data         = './yang_data/'//trim(station_name)//'.AA.PRE.semp_picked'
-!   file_syn          = './yang_syn/'//trim(station_name)//'.AA.PRE.semp_picked'
-!   file_adj          = './yang_data/'//trim(station_name)//'.AA.PRE.semp_adj'
-!   file_adj_density  = './yang_syn/'//trim(station_name)//'.AA.PRE.semp_adjden'
-!   open(unit=1001,file=trim(file_data),status='unknown')
-!   open(unit=1002,file=trim(file_syn),status='unknown')
-!   open(unit=1003,file=trim(file_adj),status='unknown')
-!   open(unit=1004,file=trim(file_adj_density),status='unknown')
-!   do itime = 1,nstep
-!      write(1001,*) t(itime), data_origin(itime,irec)
-!      write(1002,*) t(itime), syn_origin(itime,irec)
-!      write(1003,*) t(itime), adj(itime,irec)
-!      write(1004,*) t(itime), adj_density(itime,irec)
-!   end do
-!   close(1001)
-!   close(1002)
-!   close(1003)
-!   close(1004)
-   enddo 
+
+
+   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)

More information about the CIG-COMMITS mailing list