Commit 6f519221 authored by tilman.metz's avatar tilman.metz

translated input variables to english

parent 7d5596b0
......@@ -34,16 +34,16 @@ timestep_(in_seconds)_(DT) = 5.0e-5
#
#--------------------Source---------------------------------------
# Shape_of_source-signal:
(ricker=1;fumue=2;from_SIGNAL_FILE=3;SIN**3=4;deltapulse=5)_(QUELLART) = 4
point_source_(explosive=1;force_in_x=2;in_y=3;in_z=4;custom=5)_(QUELLTYP) = 4
# If QUELLTYP <5 the following two lines are ignored
(ricker=1;fumue=2;from_SIGNAL_FILE=3;SIN**3=4;deltapulse=5)_(SOURCE_SHAPE) = 4
point_source_(explosive=1;force_in_x=2;in_y=3;in_z=4;custom=5)_(SOURCE_TYPE) = 4
# If SOURCE_TYPE <5 the following two lines are ignored
force_angle_between_x_y_(in_degree)_(ALPHA) = 45.0
force_angle_between_x_z_(in_degree)_(BETA) = 45.0
# Plane wave (PW) excitation,if PLANE_WAVE_DEPTH>0, SRCREC is treated as 0
depth_of_PW_excitation_(no<=0)_(in_meter)_(PLANE_WAVE_DEPTH) = 0.0direct.com/
dip_of_PW_from_vertical_(in_degrees)_(PHI) =0.0
duration_of_source-signal_PW_(in_seconds)_(TS) = 0.0033
# External signal input instead of QUELLART
# External signal input instead of SOURCE_SHAPE
SIGNAL_FILE = ?
read_source_positions_from_SOURCE_FILE_(yes=1)_(SRCREC) = 1
SOURCE_FILE = ./sources/sources_toy.dat
......
......@@ -34,16 +34,16 @@ timestep_(in_seconds)_(DT) = 5.0e-5
#
#--------------------Source---------------------------------------
# Shape_of_source-signal:
(ricker=1;fumue=2;from_SIGNAL_FILE=3;SIN**3=4;deltapulse=5)_(QUELLART) = 4
point_source_(explosive=1;force_in_x=2;in_y=3;in_z=4;custom=5)_(QUELLTYP) = 4
# If QUELLTYP <5 the following two lines are ignored
(ricker=1;fumue=2;from_SIGNAL_FILE=3;SIN**3=4;deltapulse=5)_(SOURCE_SHAPE) = 4
point_source_(explosive=1;force_in_x=2;in_y=3;in_z=4;custom=5)_(SOURCE_TYPE) = 4
# If SOURCE_TYPE <5 the following two lines are ignored
force_angle_between_x_y_(in_degree)_(ALPHA) = 45.0
force_angle_between_x_z_(in_degree)_(BETA) = 45.0
# Plane wave (PW) excitation,if PLANE_WAVE_DEPTH>0, SRCREC is treated as 0
depth_of_PW_excitation_(no<=0)_(in_meter)_(PLANE_WAVE_DEPTH) = 0.0direct.com/
dip_of_PW_from_vertical_(in_degrees)_(PHI) =0.0
duration_of_source-signal_PW_(in_seconds)_(TS) = 0.0033
# External signal input instead of QUELLART
# External signal input instead of SOURCE_SHAPE
SIGNAL_FILE = ?
read_source_positions_from_SOURCE_FILE_(yes=1)_(SRCREC) = 1
SOURCE_FILE = ./sources/sources_toy.dat
......
......@@ -25,7 +25,7 @@
void exchange_par(void){
/* declaration of extern variables */
extern int NX, NY, NZ, QUELLART, QUELLTYP, SNAP, SNAP_FORMAT, SNAP_PLANE;
extern int NX, NY, NZ, SOURCE_SHAPE, SOURCE_TYPE, SNAP, SNAP_FORMAT, SNAP_PLANE;
extern int DRX, DRY, L, SRCREC, FDORDER;
extern float DX, DY, DZ, TIME, DT, *FL, TS, TAU, PLANE_WAVE_DEPTH, PHI;
extern float XREC1, XREC2, YREC1, YREC2, ZREC1, ZREC2;
......@@ -101,8 +101,8 @@ void exchange_par(void){
idum[5] = NX;
idum[6] = NY;
idum[7] = NZ;
idum[8] = QUELLART;
idum[9] = QUELLTYP;
idum[8] = SOURCE_SHAPE;
idum[9] = SOURCE_TYPE;
idum[10] = READMOD;
idum[11] = L;
idum[12] = FREE_SURF;
......@@ -229,8 +229,8 @@ void exchange_par(void){
NX = idum[5];
NY = idum[6];
NZ = idum[7];
QUELLART = idum[8];
QUELLTYP = idum[9];
SOURCE_SHAPE = idum[8];
SOURCE_TYPE = idum[9];
READMOD = idum[10];
L = idum[11];
FREE_SURF = idum[12];
......
......@@ -29,7 +29,7 @@ float TSNAP1, TSNAP2, TSNAPINC, *FL, TAU, REC_ARRAY_DEPTH, REC_ARRAY_DIST;
float XREC1, XREC2, YREC1, YREC2, ZREC1=0.0, ZREC2=0.0;
float REFREC[4]={0.0, 0.0, 0.0, 0.0}, DAMPING=8.0, VPPML, FPML;
int SEISMO, NDT, NDTSHIFT, NGEOPH, SEIS_FORMAT[6], FREE_SURF, READMOD, MOD_FORMAT[6], READREC, REC_ARRAY, LOG, FDORDER, FW=0, ABS_TYPE, BLOCK;
int NX, NY, NZ=1, NT, QUELLART, QUELLTYP, SNAP, SNAP_FORMAT, BOUNDARY, SRCREC, SNAP_PLANE;
int NX, NY, NZ=1, NT, SOURCE_SHAPE, SOURCE_TYPE, SNAP, SNAP_FORMAT, BOUNDARY, SRCREC, SNAP_PLANE;
float ALPHA, BETA;
int NXG, NYG, NZG, IDX, IDY, IDZ, L=1, NX1, NX2, NY1, NY2, NZ1, NZ2, DRX, DRY, RUN_MULTIPLE_SHOTS, FDCOEFF;
char SNAP_FILE[STRING_SIZE], SOURCE_FILE[STRING_SIZE], SIGNAL_FILE[STRING_SIZE], INV_FILE[STRING_SIZE];
......
......@@ -591,7 +591,7 @@ MPI_Barrier(MPI_COMM_WORLD);
MPI_Barrier(MPI_COMM_WORLD);
srcpos= fmatrix(1,7,1,nsrc);
sources(fpsrc,&nsrc,srcpos);
/*originally, QUELLTY=stype is defined in the source file, if not, QUELLTYP is taken from the input file */
/*originally, SOURCE_TYPE (stype) is defined in the source file, if not, SOURCE_TYPE is taken from the input file */
/*if (stype==NULL) printf("PE%d: Source type(s) undefined?! \n",MYID);*/
}
......@@ -712,7 +712,7 @@ CPML_coeff(K_x,alpha_prime_x,a_x,b_x,K_x_half,alpha_prime_x_half,a_x_half,b_x_ha
}
}
/*printf("source=%e,%e,%e,%e,%e,%e,%e\n",srcpos_loc[1][1],srcpos_loc[2][1],srcpos_loc[3][1],srcpos_loc[4][1],srcpos_loc[5][1],srcpos_loc[6][1],srcpos_loc[7][1]);*/
if(nsrc_loc>0){ wavelet(srcpos_loc,nsrc_loc,QUELLART,signals);
if(nsrc_loc>0){ wavelet(srcpos_loc,nsrc_loc,SOURCE_SHAPE,signals);
/*if (SEISMO){
saveseis(FP,signals,sectionvy,sectionvz,sectionp,sectioncurl,sectiondiv,recpos,recpos_loc,ntr,srcpos,ishot,ns,0,11);
}*/
......@@ -1044,7 +1044,7 @@ CPML_coeff(K_x,alpha_prime_x,a_x,b_x,K_x_half,alpha_prime_x_half,a_x_half,b_x_ha
}
if(nsrc_loc>0){
wavelet(srcpos_loc,nsrc_loc,QUELLART,signals);
wavelet(srcpos_loc,nsrc_loc,SOURCE_SHAPE,signals);
if(FILT==1) filt_seis(signals,nsrc_loc,NT,finv[nf-1]);}
zero_wavefield(NX,NY,NZ,vx,vy,vz,sxx,syy,szz,sxy,syz,sxz,rxx,ryy,rzz,rxy,ryz,rxz,psi_sxx_x,psi_sxy_x,psi_sxz_x,psi_sxy_y,psi_syy_y,psi_syz_y,psi_sxz_z,psi_syz_z,psi_szz_z,psi_vxx,psi_vyx,psi_vzx,psi_vxy,psi_vyy,psi_vzy,psi_vxz,psi_vyz,psi_vzz);
......
......@@ -27,7 +27,7 @@ int read_par(FILE *fp_in){
/* declaration of extern variables */
extern int NX, NY, NZ, QUELLART, QUELLTYP, SNAP, SNAP_FORMAT, SNAP_PLANE;
extern int NX, NY, NZ, SOURCE_SHAPE, SOURCE_TYPE, SNAP, SNAP_FORMAT, SNAP_PLANE;
extern int DRX, DRY, L, SRCREC, FDORDER, FW, FDCOEFF;
extern float DX, DY, DZ, TIME, DT, TS, *FL, TAU, PLANE_WAVE_DEPTH, PHI;
extern float XREC1, XREC2, YREC1, YREC2, ZREC1, ZREC2, ALPHA, BETA;
......@@ -107,11 +107,11 @@ int read_par(FILE *fp_in){
sscanf(cline,"%s =%f",s,&DT);
break;
case 14 :
nvarin=sscanf(cline,"%s =%i , %f , %f",s,&QUELLART,&FC,&AMP);
nvarin=sscanf(cline,"%s =%i , %f , %f",s,&SOURCE_SHAPE,&FC,&AMP);
switch(nvarin){
case 0:
case 1: QUELLART=0;
/* fprintf(stderr," Caution: QUELLART set to %d.\n",QUELLART); */
case 1: SOURCE_SHAPE=0;
/* fprintf(stderr," Caution: SOURCE_SHAPE set to %d.\n",SOURCE_SHAPE); */
case 2: FC=0.0;
/* fprintf(stderr," FC set to %d.\n",FC); */
case 3: AMP=1.0;
......@@ -121,7 +121,7 @@ int read_par(FILE *fp_in){
}
break;
case 15 :
sscanf(cline,"%s =%i",s,&QUELLTYP);
sscanf(cline,"%s =%i",s,&SOURCE_TYPE);
break;
case 16 :
sscanf(cline,"%s =%f",s,&ALPHA);
......
......@@ -27,7 +27,7 @@
void sources(FILE * fpsrc, int *nsrc, float **srcpos){
/* declaration of extern variables */
extern int MYID,SRC_MF,QUELLTYP;
extern int MYID,SRC_MF,SOURCE_TYPE;
extern FILE *FP;
extern float TS,REFSRC[3],SRCTSHIFT,FC,AMP;
......@@ -48,7 +48,7 @@ void sources(FILE * fpsrc, int *nsrc, float **srcpos){
case 3: tshift=0.0;
case 4: srcpos[5][l]=FC;
case 5: srcpos[6][l]=AMP;
case 6: srcpos[7][l]=QUELLTYP;
case 6: srcpos[7][l]=SOURCE_TYPE;
}
/*change stype to srcpos[7]*/
/*note that internally "y" is used for the vertical coordinate,
......@@ -107,7 +107,7 @@ void pwsources(int *nsrc, float **srcpos){ /* plane wave excitation */
/* declaration of extern variables */
extern float PLANE_WAVE_DEPTH, TS, DX, DY, PHI;
extern int MYID, NXG, NYG, SRCREC, FW,QUELLTYP;
extern int MYID, NXG, NYG, SRCREC, FW,SOURCE_TYPE;
extern FILE *FP;
float x, y, z, tan_phi;
......@@ -158,7 +158,7 @@ void pwsources(int *nsrc, float **srcpos){ /* plane wave excitation */
srcpos[4][isrc]=0.0;
srcpos[5][isrc]=1.0/TS;
srcpos[6][isrc]=1.0;
srcpos[7][isrc]=QUELLTYP;
srcpos[7][isrc]=SOURCE_TYPE;
}
/*double check if number of receivers match, counted within loop*/
......
......@@ -30,7 +30,7 @@ void wavelet(float ** srcpos_loc, int nsrc, int quelltype, float **signals){
/* extern variables */
extern int QUELLART, NT, MYID;
extern int SOURCE_SHAPE, NT, MYID;
extern float DT;
extern char SIGNAL_FILE[STRING_SIZE];
extern FILE *FP;
......@@ -94,6 +94,6 @@ void wavelet(float ** srcpos_loc, int nsrc, int quelltype, float **signals){
fprintf(FP," have been assigned with a source signal. \n");
if (QUELLART==3) free_vector(psource,1,NT);
if (SOURCE_SHAPE==3) free_vector(psource,1,NT);
}
......@@ -27,7 +27,7 @@
void writepar(FILE *fp, int ns){
/* declaration of extern variables */
extern int NX, NY, NZ, NT, QUELLART, QUELLTYP, FDORDER, RUN_MULTIPLE_SHOTS;
extern int NX, NY, NZ, NT, SOURCE_SHAPE, SOURCE_TYPE, FDORDER, RUN_MULTIPLE_SHOTS;
extern int SNAP, SNAP_FORMAT, REC_ARRAY, L, SNAP_PLANE,FW;
extern float DX, DY, DZ, TIME, DT, TS, *FL, TAU, PLANE_WAVE_DEPTH;
extern float XREC1, XREC2, YREC1, YREC2, ZREC1, ZREC2;
......@@ -99,7 +99,7 @@ void writepar(FILE *fp, int ns){
fprintf(fp," Wavelet of source:");
switch (QUELLART){
switch (SOURCE_SHAPE){
case 1 :
fprintf(fp," Ricker\n");
break;
......@@ -122,7 +122,7 @@ void writepar(FILE *fp, int ns){
}
fprintf(fp," Default type of source:");
switch (QUELLTYP){
switch (SOURCE_TYPE){
case 1 :
fprintf(fp," explosive point source (concentrated at a single gridpoint)\n");
break;
......@@ -141,8 +141,8 @@ void writepar(FILE *fp, int ns){
fprintf(fp," Angle between x and z direticon (BETA): %f\n", BETA);
break;
default :
fprintf(fp," WARNING: Default type of source ('%d') not available -> changed to explosive! ", QUELLTYP);
QUELLTYP=1;
fprintf(fp," WARNING: Default type of source ('%d') not available -> changed to explosive! ", SOURCE_TYPE);
SOURCE_TYPE=1;
}
fprintf(fp,"\n");
fprintf(fp,"\n Source file: %s \n", SOURCE_FILE);
......
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