Skip to content
Snippets Groups Projects
Commit ce4d2167 authored by Lars Houpt's avatar Lars Houpt
Browse files

Suffix (.su) moved to the end of filename in I/O of seismic data.

parent 8f7849f3
No related branches found
No related tags found
No related merge requests found
......@@ -45,11 +45,11 @@ void readrec(int ishot, int **topo, st_acquisition *acq, GlobVar *gv, int rec_lo
switch (gv->READREC) { /* read receiver positions from file */
case 3: /* read receiver positions from SU file headers */
if (gv->RUN_MULTIPLE_SHOTS) {
sprintf(shotnum, ".shot%d", shot_id);
sprintf(shotnum, "shot%d", shot_id);
} else {
sprintf(shotnum, "%s", "");
}
sprintf(filename, "%s_%s.su%s", gv->SEIS_OBS_FILE, component[gv->ADJOINT_TYPE-1], shotnum);
sprintf(filename, "%s_%s_%s.su", gv->SEIS_OBS_FILE, component[gv->ADJOINT_TYPE-1], shotnum);
log_info("Reading receiver positions from file %s.\n", filename);
//acq->recpos = scan_suhead_rec(&filename, gv, topo, &ntr);
......
......@@ -51,12 +51,12 @@ void readseis(int shot_id, float **section, const st_acquisition *acq, int ntr,
}
if (gv->RUN_MULTIPLE_SHOTS) {
sprintf(shotnum, ".shot%d", shot_id);
sprintf(shotnum, "shot%d", shot_id);
} else {
sprintf(shotnum, "%s", "");
}
sprintf(data, "%s_%s.su%s", gv->SEIS_OBS_FILE, component[comp-1], shotnum);
sprintf(data, "%s_%s_%s.su", gv->SEIS_OBS_FILE, component[comp-1], shotnum);
if(gv->MYID_SHOT==0) log_info("Reading file %s.\n", data);
if (ntr > 0) {
......
......@@ -88,12 +88,12 @@ void saveseis(st_seismogram *section, st_acquisition *acq, int ntr_glob, int ish
}
if (gv->RUN_MULTIPLE_SHOTS) {
sprintf(vxf, "%s_vx%s.%s.shot%d", seisfile, file_iter, file_ext, shot_id);
sprintf(vyf, "%s_vy%s.%s.shot%d", seisfile, file_iter, file_ext, shot_id);
sprintf(vzf, "%s_vz%s.%s.shot%d", seisfile, file_iter, file_ext, shot_id);
sprintf(curlf, "%s_rot%s.%s.shot%d", seisfile, file_iter, file_ext, shot_id);
sprintf(divf, "%s_div%s.%s.shot%d", seisfile, file_iter, file_ext, shot_id);
sprintf(pf, "%s_p%s.%s.shot%d", seisfile, file_iter, file_ext, shot_id);
sprintf(vxf, "%s_vx%s_shot%d.%s", seisfile, file_iter, shot_id, file_ext);
sprintf(vyf, "%s_vy%s_shot%d.%s", seisfile, file_iter, shot_id, file_ext);
sprintf(vzf, "%s_vz%s_shot%d.%s", seisfile, file_iter, shot_id, file_ext);
sprintf(curlf, "%s_rot%s_shot%d.%s", seisfile, file_iter, shot_id, file_ext);
sprintf(divf, "%s_div%s_shot%d.%s", seisfile, file_iter, shot_id, file_ext);
sprintf(pf, "%s_p%s_shot%d.%s", seisfile, file_iter, shot_id, file_ext);
} else {
sprintf(vxf, "%s_vx%s.%s", seisfile, file_iter, file_ext);
sprintf(vyf, "%s_vy%s.%s", seisfile, file_iter, file_ext);
......
......@@ -66,7 +66,7 @@ void savesig(float **signals, st_acquisition *acq, int nsrc_loc, int ishot, int
switch (fswitch) {
case 0:
if (gv->RUN_MULTIPLE_SHOTS) {
sprintf(sigf, "%s_srcsig%s.%s.shot%d", seisfile, file_iter, file_ext, shot_id);
sprintf(sigf, "%s_srcsig%s_shot%d.%s", seisfile, file_iter, shot_id, file_ext);
nsig = 1;
} else {
sprintf(sigf, "%s_srcsig%s.%s", seisfile, file_iter, file_ext);
......@@ -75,7 +75,7 @@ void savesig(float **signals, st_acquisition *acq, int nsrc_loc, int ishot, int
break;
case 1:
if (gv->RUN_MULTIPLE_SHOTS) {
sprintf(sigf, "%s_srcsig_stfi%s.%s.shot%d", seisfile, file_iter, file_ext, shot_id);
sprintf(sigf, "%s_srcsig_stfi%s_shot%d.%s", seisfile, file_iter, shot_id, file_ext);
nsig = 1;
} else {
sprintf(sigf, "%s_srcsig_stfi%s.%s", seisfile, file_iter, file_ext);
......@@ -84,7 +84,7 @@ void savesig(float **signals, st_acquisition *acq, int nsrc_loc, int ishot, int
break;
case 2:
if (gv->RUN_MULTIPLE_SHOTS) {
sprintf(sigf, "%s_srcsig_stfi_stepl%s.%s.shot%d", seisfile, file_iter, file_ext, shot_id);
sprintf(sigf, "%s_srcsig_stfi_stepl%s_shot%d.%s", seisfile, file_iter, shot_id, file_ext);
nsig = 1;
} else {
sprintf(sigf, "%s_srcsig_stfi_stepl%s.%s", seisfile, file_iter, file_ext);
......
......@@ -72,7 +72,7 @@ void time_window(float **sectiondata1, float **sectiondata2, int **recpos_loc, i
}
}
//log_info("Read mute from pick file %s, starting at line %d.\n", pickfile, gv->TW_START_READING_AT_LINE);
log_info("Read mute from pick file %s, starting at line %d.\n", pickfile, gv->TW_START_READING_AT_LINE);
for(i=1;i<=gv->TW_START_READING_AT_LINE-1;i++){
fscanf(fptime,"%s %s\n", dummy1, dummy2);
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment