further cleanup

src/IFOS2D.c

@@ -1990,6 +1990,11 @@ int main(int argc, char **argv){
                                                    }
                                                    h++;
                                                }
+//                                                /* for testing */
+//                                                sprintf(pf,"%s_pup.su.shot%i",SEIS_FILE,ishot);
+//                                                if(myid_ntr==0){
+//                                                    outseis_glob(FP,fopen(pf,"w"), 0, fulldata_p,recpos,recpos_loc,ntr_glob,srcpos,1,ns,SEIS_FORMAT,ishot,0);
+//                                                 }
                                            }
                                        }
                                 }
@@ -2073,8 +2078,6 @@ int main(int argc, char **argv){
                                         if ((TIME_FILT==1 )|| (TIME_FILT==2)){
                                             timedomain_filt(sectionread,F_LOW_PASS,ORDER,ntr_glob,ns,1);
                                         }
-                                        sprintf(pf,"%s_pup_data.su.shot%i",SEIS_FILE,ishot);
-                                            outseis_glob(FP,fopen(pf,"w"), 0, sectionread,recpos,recpos_loc,ntr_glob,srcpos,1,ns,SEIS_FORMAT,ishot,1);
                                         h=1;
                                         for(i=1;i<=ntr;i++){
                                             for(j=1;j<=ns;j++){

src/pup.c

@@ -194,6 +194,7 @@ void	pup(float **data_p, float **data_vy, FILE *fp, int ntr_glob, int  **recpos,
 		for (i=1; i<=ns; i++) {
 			data_p[k][i]=data_pup[i][k];
 		}
+		data_p[k][1]=0;		/*first sample need to be set to 0 to avoid artefacts after integrating from velocity to displacement*/
 	}
 
 // 	/* save pup*/