[cig-commits] r18178 - in seismo/3D/SPECFEM3D_GLOBE/trunk: EXAMPLES/global_s362ani/SEM EXAMPLES/regional_MiddleEast/SEM UTILS/adjoint_sources UTILS/adjoint_sources/amplitude UTILS/adjoint_sources/traveltime UTILS/lib src/specfem3D
danielpeter at geodynamics.org
danielpeter at geodynamics.org
Tue Apr 5 11:02:33 PDT 2011
Author: danielpeter
Date: 2011-04-05 11:02:33 -0700 (Tue, 05 Apr 2011)
New Revision: 18178
Added:
seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/lib/drw_ascfile.h
seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/lib/rw_fortran_wrapper.c
Removed:
seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/amplitude/drw_ascfile.h
seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/amplitude/rw_ascfile_c.c
seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/amplitude/rw_ascfile_f.f90
seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/amplitude/rw_fortran_wrapper.c
seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/traveltime/drw_ascfile.h
seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/traveltime/rw_ascfile_c.c
seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/traveltime/rw_ascfile_f.f90
seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/traveltime/rw_fortran_wrapper.c
Modified:
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/global_s362ani/SEM/create_adjoint_source.sh
seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/regional_MiddleEast/SEM/create_adjoint_source.sh
seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/README
seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/amplitude/Makefile
seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/amplitude/compile_cut
seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/traveltime/Makefile
seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/traveltime/compile_cut
seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/traveltime/create_adjsrc_traveltime.f90
seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/noise_tomography.f90
Log:
updates Makefiles for UTILS/adjoint_sources; adds a check in routine read_parameters_noise()
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/global_s362ani/SEM/create_adjoint_source.sh
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/global_s362ani/SEM/create_adjoint_source.sh 2011-04-05 16:25:05 UTC (rev 18177)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/global_s362ani/SEM/create_adjoint_source.sh 2011-04-05 18:02:33 UTC (rev 18178)
@@ -14,7 +14,7 @@
t2="755.0"
#window out single phase arrival on vertical component between t1 to t2 :
-~/SPECFEM3D_GLOBE/UTILS/cut_velocity/xcut_velocity $t1 $t2 3 ../REF_SEIS/$station.$network.MX*
+~/SPECFEM3D_GLOBE/UTILS/adjoint_sources/traveltime/xcreate_adjsrc_traveltime $t1 $t2 3 ../REF_SEIS/$station.$network.MX*
mv ../REF_SEIS/$station.$network.MX*adj ./
# rename adjoint source files:
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/regional_MiddleEast/SEM/create_adjoint_source.sh
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/regional_MiddleEast/SEM/create_adjoint_source.sh 2011-04-05 16:25:05 UTC (rev 18177)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/EXAMPLES/regional_MiddleEast/SEM/create_adjoint_source.sh 2011-04-05 18:02:33 UTC (rev 18178)
@@ -14,7 +14,7 @@
t2="225."
#window out single phase arrival on vertical component between t1 to t2 :
-~/SPECFEM3D_GLOBE/UTILS/cut_velocity/xcut_velocity $t1 $t2 3 ../REF_SEIS/$station.$network.MX*
+~/SPECFEM3D_GLOBE/UTILS/adjoint_sources/traveltime/xcreate_adjsrc_traveltime $t1 $t2 3 ../REF_SEIS/$station.$network.MX*
mv ../REF_SEIS/$station.$network.MX*adj ./
# rename adjoint source files:
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/README
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/README 2011-04-05 16:25:05 UTC (rev 18177)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/README 2011-04-05 18:02:33 UTC (rev 18178)
@@ -8,13 +8,13 @@
this program cuts certain portion of the seismograms and converts them into
the adjoint sources for generating banana-dougnut kernels.
- (see Tromp et al. (2005), eq. (45))
+ (see Tromp et al. (2005), eq. 45)
- amplitude/ subdirectory: amplitude adjoint sources
this program cuts a certain portion of displacement seismograms and
converts them into adjoint sources for generating classical amplitude
- kernels following Tromp et al. (2005) eq.67.
+ kernels following Tromp et al. (2005) eq. 67.
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/amplitude/Makefile
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/amplitude/Makefile 2011-04-05 16:25:05 UTC (rev 18177)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/amplitude/Makefile 2011-04-05 18:02:33 UTC (rev 18178)
@@ -21,8 +21,8 @@
create_adjsrc_amplitude.o: create_adjsrc_amplitude.f90
${F90} -Wall -c create_adjsrc_amplitude.f90
-rw_ascfile_c.o: ../lib/rw_ascfile_c.c
- ${CC} -c -o rw_ascfile_c.o ../lib/rw_ascfile_c.c
+rw_ascfile_c.o: ../../lib/rw_ascfile_c.c
+ ${CC} -c -o rw_ascfile_c.o ../../lib/rw_ascfile_c.c
clean:
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/amplitude/compile_cut
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/amplitude/compile_cut 2011-04-05 16:25:05 UTC (rev 18177)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/amplitude/compile_cut 2011-04-05 18:02:33 UTC (rev 18178)
@@ -10,7 +10,3 @@
echo
echo
-# obsolete...
-#gcc -c -o rw_ascfile_c.o rw_ascfile_c.c
-#gcc -c -o rw_fortran_wrapper.o rw_fortran_wrapper.c
-#ifort -O3 -o xcreate_adjsrc_amplitude create_adjsrc_amplitude.f90 rw_ascfile_f.f90 rw_ascfile_c.o rw_fortran_wrapper.o
Deleted: seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/amplitude/drw_ascfile.h
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/amplitude/drw_ascfile.h 2011-04-05 16:25:05 UTC (rev 18177)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/amplitude/drw_ascfile.h 2011-04-05 18:02:33 UTC (rev 18178)
@@ -1,29 +0,0 @@
-#ifndef _drw_ascfile_h
-#define _drw_ascfile_h
-
-void dread_ascfile(const char *ascfile,
- double *t0, double *dt, int *n,
- double *data);
-void dwrite_ascfile(const char *ascfile,
- double t0, double dt, int n,
- const double *data);
-
-void dread_ascfile_c(const char *ascfile,
- double *t0, double *dt, int *n,
- double *data);
-void dwrite_ascfile_c(const char *ascfile,
- const double *t0, const double *dt, const int *n,
- const double *data);
-
-void DREAD_ASCFILE_C(const char *ascfile,
- double *t0, double *dt, int *n,
- double *data);
-void DWRITE_ASCFILE_C(const char *ascfile,
- const double *t0, const double *dt, const int *n,
- const double *data);
-
-
-
-
-
-#endif
Deleted: seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/amplitude/rw_ascfile_c.c
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/amplitude/rw_ascfile_c.c 2011-04-05 16:25:05 UTC (rev 18177)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/amplitude/rw_ascfile_c.c 2011-04-05 18:02:33 UTC (rev 18178)
@@ -1,50 +0,0 @@
-#include <stdio.h>
-#include "drw_ascfile.h"
-
-void dread_ascfile(const char *ascfile,
- double *t0, double *dt, int *n,
- double *data)
-
-{
- FILE *fd;
- double junk,junk1;
- int i;
-
- if ((fd = fopen(ascfile,"r")) == NULL) {
- printf(" file %s cannot be opened\n",ascfile);
- exit(1);
- }
- i = 0;
- while ( fscanf(fd,"%lf %lf\n",&junk, data+i) != EOF ) {
- if (i == 0) junk1 = junk;
- if (i == 1) *dt = junk - junk1;
- i++;}
- *t0 = junk1;
- *n = i;
- if (fclose(fd) != 0) {
- printf(" file %s cannot be closed\n",ascfile);
- exit(1);}
-
-}
-
-void dwrite_ascfile(const char *ascfile,
- double t0, double dt, int n,
- const double *data)
-
-{
- FILE *fd;
- int i;
-
- if ((fd = fopen(ascfile,"w")) == NULL) {
- printf(" file %s cannot be opened to write\n",ascfile);
- exit(1);
- }
- i = 0;
- for (i=0; i<n; i++) {
- fprintf(fd,"%14.7g %18.7g\n", t0+i*dt, data[i]);
- }
- if (fclose(fd) != 0) {
- printf("file %s cannot be closed\n",ascfile);
- exit(1);}
-
-}
Deleted: seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/amplitude/rw_ascfile_f.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/amplitude/rw_ascfile_f.f90 2011-04-05 16:25:05 UTC (rev 18177)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/amplitude/rw_ascfile_f.f90 2011-04-05 18:02:33 UTC (rev 18178)
@@ -1,29 +0,0 @@
-
-
-! these functions can be used directly in fortran
-
-subroutine dread_ascfile_f(name,t0,dt,n,data)
-
- implicit none
-
- character(len=*) :: name
- real*8 :: t0, dt, data(*)
- integer :: n
-
- call dread_ascfile_c(trim(name)//char(0), t0, dt, n, data)
-
-end subroutine dread_ascfile_f
-
-
-subroutine dwrite_ascfile_f(name,t0,dt,n,data)
-
- implicit none
-
- character(len=*) :: name
- real*8 :: t0, dt, data(*)
- integer :: n
-
- call dwrite_ascfile_c(trim(name)//char(0), t0, dt, n, data)
-
-end subroutine dwrite_ascfile_f
-
Deleted: seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/amplitude/rw_fortran_wrapper.c
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/amplitude/rw_fortran_wrapper.c 2011-04-05 16:25:05 UTC (rev 18177)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/amplitude/rw_fortran_wrapper.c 2011-04-05 18:02:33 UTC (rev 18178)
@@ -1,33 +0,0 @@
-#include "drw_ascfile.h"
-
-void dread_ascfile_c(const char *name, double *t0, double *dt, int * n, double *data) {
-
- dread_ascfile(name,t0,dt,n,data);}
-
-void DREAD_ASCFILE_C(const char *name, double *t0, double *dt, int * n, double *data) {
-
- dread_ascfile(name,t0,dt,n,data);}
-
-void dread_ascfile_c_(const char *name, double *t0, double *dt, int * n, double *data) {
-
- dread_ascfile(name,t0,dt,n,data);}
-
-void DREAD_ASCFILE_C_(const char *name, double *t0, double *dt, int * n, double *data) {
-
- dread_ascfile(name,t0,dt,n,data);}
-
-void dwrite_ascfile_c(const char *name, const double *t0, const double *dt, const int *n, const double *data) {
-
- dwrite_ascfile(name,*t0,*dt,*n,data);}
-
-void DWRITE_ASCFILE_C(const char *name, const double *t0, const double *dt, const int *n, const double *data) {
-
- dwrite_ascfile(name,*t0,*dt,*n,data);}
-
-void dwrite_ascfile_c_(const char *name, const double *t0, const double *dt, const int *n, const double *data) {
-
- dwrite_ascfile(name,*t0,*dt,*n,data);}
-
-void DWRITE_ASCFILE_C_(const char *name, const double *t0, const double *dt, const int *n, const double *data) {
-
- dwrite_ascfile(name,*t0,*dt,*n,data);}
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/traveltime/Makefile
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/traveltime/Makefile 2011-04-05 16:25:05 UTC (rev 18177)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/traveltime/Makefile 2011-04-05 18:02:33 UTC (rev 18178)
@@ -21,8 +21,8 @@
create_adjsrc_traveltime.o: create_adjsrc_traveltime.f90
${F90} -Wall -c create_adjsrc_traveltime.f90
-rw_ascfile_c.o: ../lib/rw_ascfile_c.c
- ${CC} -c -o rw_ascfile_c.o ../lib/rw_ascfile_c.c
+rw_ascfile_c.o: ../../lib/rw_ascfile_c.c
+ ${CC} -c -o rw_ascfile_c.o ../../lib/rw_ascfile_c.c
clean:
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/traveltime/compile_cut
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/traveltime/compile_cut 2011-04-05 16:25:05 UTC (rev 18177)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/traveltime/compile_cut 2011-04-05 18:02:33 UTC (rev 18178)
@@ -10,7 +10,3 @@
echo
echo
-# obsolete...
-#gcc -c -o rw_ascfile_c.o rw_ascfile_c.c
-#gcc -c -o rw_fortran_wrapper.o rw_fortran_wrapper.c
-#ifort -O3 -o xcreate_adjsrc_traveltime create_adjsrc_traveltime.f90 rw_ascfile_f.f90 rw_ascfile_c.o rw_fortran_wrapper.o
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/traveltime/create_adjsrc_traveltime.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/traveltime/create_adjsrc_traveltime.f90 2011-04-05 16:25:05 UTC (rev 18177)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/traveltime/create_adjsrc_traveltime.f90 2011-04-05 18:02:33 UTC (rev 18178)
@@ -3,6 +3,7 @@
! this program cuts certain portion of the seismograms and converts them into
! the adjoint sources for generating banana-dougnut kernels.
! Qinya Liu, Caltech, May 2007
+! (renamed from cut_velocity.f90 (2011)
!
! call by: ./xcreate_adjsrc_traveltime t1 t2 ifile[0-5] E/N/Z-ascii-files [baz]
!
Deleted: seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/traveltime/drw_ascfile.h
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/traveltime/drw_ascfile.h 2011-04-05 16:25:05 UTC (rev 18177)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/traveltime/drw_ascfile.h 2011-04-05 18:02:33 UTC (rev 18178)
@@ -1,29 +0,0 @@
-#ifndef _drw_ascfile_h
-#define _drw_ascfile_h
-
-void dread_ascfile(const char *ascfile,
- double *t0, double *dt, int *n,
- double *data);
-void dwrite_ascfile(const char *ascfile,
- double t0, double dt, int n,
- const double *data);
-
-void dread_ascfile_c(const char *ascfile,
- double *t0, double *dt, int *n,
- double *data);
-void dwrite_ascfile_c(const char *ascfile,
- const double *t0, const double *dt, const int *n,
- const double *data);
-
-void DREAD_ASCFILE_C(const char *ascfile,
- double *t0, double *dt, int *n,
- double *data);
-void DWRITE_ASCFILE_C(const char *ascfile,
- const double *t0, const double *dt, const int *n,
- const double *data);
-
-
-
-
-
-#endif
Deleted: seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/traveltime/rw_ascfile_c.c
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/traveltime/rw_ascfile_c.c 2011-04-05 16:25:05 UTC (rev 18177)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/traveltime/rw_ascfile_c.c 2011-04-05 18:02:33 UTC (rev 18178)
@@ -1,50 +0,0 @@
-#include <stdio.h>
-#include "drw_ascfile.h"
-
-void dread_ascfile(const char *ascfile,
- double *t0, double *dt, int *n,
- double *data)
-
-{
- FILE *fd;
- double junk,junk1;
- int i;
-
- if ((fd = fopen(ascfile,"r")) == NULL) {
- printf(" file %s cannot be opened\n",ascfile);
- exit(1);
- }
- i = 0;
- while ( fscanf(fd,"%lf %lf\n",&junk, data+i) != EOF ) {
- if (i == 0) junk1 = junk;
- if (i == 1) *dt = junk - junk1;
- i++;}
- *t0 = junk1;
- *n = i;
- if (fclose(fd) != 0) {
- printf(" file %s cannot be closed\n",ascfile);
- exit(1);}
-
-}
-
-void dwrite_ascfile(const char *ascfile,
- double t0, double dt, int n,
- const double *data)
-
-{
- FILE *fd;
- int i;
-
- if ((fd = fopen(ascfile,"w")) == NULL) {
- printf(" file %s cannot be opened to write\n",ascfile);
- exit(1);
- }
- i = 0;
- for (i=0; i<n; i++) {
- fprintf(fd,"%14.7g %18.7g\n", t0+i*dt, data[i]);
- }
- if (fclose(fd) != 0) {
- printf("file %s cannot be closed\n",ascfile);
- exit(1);}
-
-}
Deleted: seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/traveltime/rw_ascfile_f.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/traveltime/rw_ascfile_f.f90 2011-04-05 16:25:05 UTC (rev 18177)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/traveltime/rw_ascfile_f.f90 2011-04-05 18:02:33 UTC (rev 18178)
@@ -1,29 +0,0 @@
-
-
-! these functions can be used directly in fortran
-
-subroutine dread_ascfile_f(name,t0,dt,n,data)
-
- implicit none
-
- character(len=*) :: name
- real*8 :: t0, dt, data(*)
- integer :: n
-
- call dread_ascfile_c(trim(name)//char(0), t0, dt, n, data)
-
-end subroutine dread_ascfile_f
-
-
-subroutine dwrite_ascfile_f(name,t0,dt,n,data)
-
- implicit none
-
- character(len=*) :: name
- real*8 :: t0, dt, data(*)
- integer :: n
-
- call dwrite_ascfile_c(trim(name)//char(0), t0, dt, n, data)
-
-end subroutine dwrite_ascfile_f
-
Deleted: seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/traveltime/rw_fortran_wrapper.c
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/traveltime/rw_fortran_wrapper.c 2011-04-05 16:25:05 UTC (rev 18177)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/traveltime/rw_fortran_wrapper.c 2011-04-05 18:02:33 UTC (rev 18178)
@@ -1,33 +0,0 @@
-#include "drw_ascfile.h"
-
-void dread_ascfile_c(const char *name, double *t0, double *dt, int * n, double *data) {
-
- dread_ascfile(name,t0,dt,n,data);}
-
-void DREAD_ASCFILE_C(const char *name, double *t0, double *dt, int * n, double *data) {
-
- dread_ascfile(name,t0,dt,n,data);}
-
-void dread_ascfile_c_(const char *name, double *t0, double *dt, int * n, double *data) {
-
- dread_ascfile(name,t0,dt,n,data);}
-
-void DREAD_ASCFILE_C_(const char *name, double *t0, double *dt, int * n, double *data) {
-
- dread_ascfile(name,t0,dt,n,data);}
-
-void dwrite_ascfile_c(const char *name, const double *t0, const double *dt, const int *n, const double *data) {
-
- dwrite_ascfile(name,*t0,*dt,*n,data);}
-
-void DWRITE_ASCFILE_C(const char *name, const double *t0, const double *dt, const int *n, const double *data) {
-
- dwrite_ascfile(name,*t0,*dt,*n,data);}
-
-void dwrite_ascfile_c_(const char *name, const double *t0, const double *dt, const int *n, const double *data) {
-
- dwrite_ascfile(name,*t0,*dt,*n,data);}
-
-void DWRITE_ASCFILE_C_(const char *name, const double *t0, const double *dt, const int *n, const double *data) {
-
- dwrite_ascfile(name,*t0,*dt,*n,data);}
Copied: seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/lib/drw_ascfile.h (from rev 18177, seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/traveltime/drw_ascfile.h)
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/lib/drw_ascfile.h (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/lib/drw_ascfile.h 2011-04-05 18:02:33 UTC (rev 18178)
@@ -0,0 +1,29 @@
+#ifndef _drw_ascfile_h
+#define _drw_ascfile_h
+
+void dread_ascfile(const char *ascfile,
+ double *t0, double *dt, int *n,
+ double *data);
+void dwrite_ascfile(const char *ascfile,
+ double t0, double dt, int n,
+ const double *data);
+
+void dread_ascfile_c(const char *ascfile,
+ double *t0, double *dt, int *n,
+ double *data);
+void dwrite_ascfile_c(const char *ascfile,
+ const double *t0, const double *dt, const int *n,
+ const double *data);
+
+void DREAD_ASCFILE_C(const char *ascfile,
+ double *t0, double *dt, int *n,
+ double *data);
+void DWRITE_ASCFILE_C(const char *ascfile,
+ const double *t0, const double *dt, const int *n,
+ const double *data);
+
+
+
+
+
+#endif
Copied: seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/lib/rw_fortran_wrapper.c (from rev 18177, seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/adjoint_sources/traveltime/rw_fortran_wrapper.c)
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/lib/rw_fortran_wrapper.c (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/lib/rw_fortran_wrapper.c 2011-04-05 18:02:33 UTC (rev 18178)
@@ -0,0 +1,33 @@
+#include "drw_ascfile.h"
+
+void dread_ascfile_c(const char *name, double *t0, double *dt, int * n, double *data) {
+
+ dread_ascfile(name,t0,dt,n,data);}
+
+void DREAD_ASCFILE_C(const char *name, double *t0, double *dt, int * n, double *data) {
+
+ dread_ascfile(name,t0,dt,n,data);}
+
+void dread_ascfile_c_(const char *name, double *t0, double *dt, int * n, double *data) {
+
+ dread_ascfile(name,t0,dt,n,data);}
+
+void DREAD_ASCFILE_C_(const char *name, double *t0, double *dt, int * n, double *data) {
+
+ dread_ascfile(name,t0,dt,n,data);}
+
+void dwrite_ascfile_c(const char *name, const double *t0, const double *dt, const int *n, const double *data) {
+
+ dwrite_ascfile(name,*t0,*dt,*n,data);}
+
+void DWRITE_ASCFILE_C(const char *name, const double *t0, const double *dt, const int *n, const double *data) {
+
+ dwrite_ascfile(name,*t0,*dt,*n,data);}
+
+void dwrite_ascfile_c_(const char *name, const double *t0, const double *dt, const int *n, const double *data) {
+
+ dwrite_ascfile(name,*t0,*dt,*n,data);}
+
+void DWRITE_ASCFILE_C_(const char *name, const double *t0, const double *dt, const int *n, const double *data) {
+
+ dwrite_ascfile(name,*t0,*dt,*n,data);}
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/noise_tomography.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/noise_tomography.f90 2011-04-05 16:25:05 UTC (rev 18177)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/noise_tomography.f90 2011-04-05 18:02:33 UTC (rev 18178)
@@ -129,6 +129,11 @@
close(IOUT_NOISE)
endif
+ ! checks master irec
+ if( irec_master_noise < 1 .or. irec_master_noise > nrec ) then
+ call exit_MPI(myrank,'error noise tomography: irec_master_noise is not in range of given number of receivers')
+ endif
+
! compute source arrays for "ensemble forward source", which is source of "ensemble forward wavefield"
if(myrank == islice_selected_rec(irec_master_noise) .OR. myrank == 0) then ! myrank == 0 is used for output only
call compute_arrays_source_noise(myrank, &
More information about the CIG-COMMITS
mailing list