Commit 2519bd09 authored by laura.gassner's avatar laura.gassner

adaption of seismogram output options, output of adjoint sources only if VERBOSE=1

parent f6697382
......@@ -56,7 +56,7 @@
"k_max_PML" : "8.0",
"Receiver" : "comment",
"SEISMO" : "5",
"SEISMO" : "2",
"READREC" : "1",
"REC_FILE" : "./receiver/rec_toy_example_ac.dat",
"REFRECX, REFRECY" : "0.0 , 0.0",
......
......@@ -2092,14 +2092,19 @@ int main(int argc, char **argv){
itestshot+=TESTSHOT_INCR;
}
/* Write differences between measured and synthetic seismogramms to disk */
if (SEISMO){
/* Write differences between measured and synthetic seismogramms (adjoint sources) to disk */
if (SEISMO&&VERBOSE){
if(WAVETYPE==1 || WAVETYPE==3){
if ((QUELLTYPB==1)|| (QUELLTYPB==3)){
catseis(sectionvxdiff, fulldata_vx, recswitch, ntr_glob, MPI_COMM_NTR);
}
if ((QUELLTYPB==1)|| (QUELLTYPB==2)){
catseis(sectionvydiff, fulldata_vy, recswitch, ntr_glob, MPI_COMM_NTR);
if(ACOUSTIC)
}
if (QUELLTYPB==4){
catseis(sectionpdiff, fulldata_p, recswitch, ntr_glob, MPI_COMM_NTR);
}
}
if(WAVETYPE==2 || WAVETYPE==3){
catseis(sectionvzdiff, fulldata_vz, recswitch, ntr_glob, MPI_COMM_NTR);
}
......@@ -2111,11 +2116,16 @@ int main(int argc, char **argv){
/* Write measured filtered seismogramms to disk */
if (SEISMO && TIME_FILT && WRITE_FILTERED_DATA){
if(WAVETYPE==1 || WAVETYPE==3){
if ((QUELLTYPB==1)|| (QUELLTYPB==3)){
catseis(sectionvxdata, fulldata_vx, recswitch, ntr_glob, MPI_COMM_NTR);
}
if ((QUELLTYPB==1)|| (QUELLTYPB==2)){
catseis(sectionvydata, fulldata_vy, recswitch, ntr_glob, MPI_COMM_NTR);
if(ACOUSTIC)
}
if (QUELLTYPB==4){
catseis(sectionpdata, fulldata_p, recswitch, ntr_glob, MPI_COMM_NTR);
}
}
if(WAVETYPE==2 || WAVETYPE==3){
catseis(sectionvzdata, fulldata_vz, recswitch, ntr_glob, MPI_COMM_NTR);
}
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment