The Perl Toolchain Summit needs more sponsors. If your company depends on Perl, please support this very important event.

SYNOPSIS

 FILE/PERL PROGRAM NAME: Synseis.pl

 AUTHOR:  Juan M. Lorenzo (sm), gllore@lsu.edu

 DATE:10/02/2017 05:33:04 PM

 DESCRIPTION:

 Version: 2.0

USE

NOTES for old V 1.0 This script creates the synthetic seismogram in an ascii file usage Synseis.sh Site number

 ***** SETTING SITE WATER DEPTHS******* NOTE 1 ***************************
water_depth=1123. # water depth at site 904 delay = 370 ms
water_depth=1; #1123.  #wqter depth incurrent zrhov model


 ********SETTING PATHS ******* NOTE 2 ***************************
 
 zrhov is the only file that Synseis reads on input 
 zrhov contains depth versus density and velosity in three columns
 values versus depth DO NOT HAVE  to be at regular intervals but
 commonly are.

----

# Setting paths path=./ zrhov_filename=./zrhov.904 # input filename zrhov_filename=./zrhov # input filename output_source=./source.out reflec_coef_time=./rc_t reflec_coef_depth=./rc_z reg_density_file=./zrho.reg reg_velocity_file=./zv.reg # # *********** OPTIONS ************* NOTE 3 ************************* The program has various options, e.g. 1) variable sampling rate 2) synthetic Ricker source 3A) real MCS source ready to roll! 3B) sources that need to be resampled at a finer sampling rate

 For each option several parameters must be turned off and others turned
 on.  At present this program works using an MCS source wavelet with
 a 2ms sampling interval taken from EW Line 1027 at CDP 1377 1.106-1.166ms
 (1027.source) found in /projects/projects5/Geol4068/Synseis_class/sources
 Be careful with the units as SUnix uses microseconds and you may
 like to think in terms of milliseconds or just seconds!

 ---------------- uncomment before next line if using resampled source
  OTHER OPTIONAL SOURCES with DIFFERENT (!!!) sampling intervals

 1ms SI; first source used, fromCDP 1400, line2, SCS
 input_source_filename=$path/sources/CDP1400.source
 time_sampling_interval=0.001  in seconds

 1ms SI;  CDP 1210 1086-1173ms 
 input_source_filename=$path/sources/line2.source 
 time_sampling_interval=0.001  in seconds

 ******************************************************
 WE USE THE FOLLOWING SOURCE IN THIS EXAMPLE
 *******************************************************
 2ms SI CDP 1377 1.106-1.166ms
 input_source_filename=./1027.source 
 time_sampling_interval=0.002  in seconds

depth_sampling_interval=0.4 in meters

 Remember that you can't use an MCS source AND a Ricker wavelet
 and a seismic source wavelet you would like to change within Synseis
 SIMLUTANEOUSLY.  Therefore, if you want to use on eof these options
 other options MUST be set off.  I advise you to not use the following
 options for the time being without help.

 If you want to use the following options
 ----------------- uncomment before next line for source for sampling
  source_resampling_interval=0.004
##  tmin_resampled_source=0.004 
  xstart4resampling_source=0.0
 -X0$xstart4resampling_source \
-S$input_source_filename

 --------------- uncomment below this line if using Ricker --------
 Ricker_endtime=0.15     s
 Ricker_file=$path2/Geol4068/users/$1/modeling/output/ricker.out
 Ricker_frequency=40.    Hz
 -AF$Ricker_frequency \
 -AE$Ricker_endtime \
 -Ao$Ricker_file 
 -----------------------------------------------------------------  
 ************* Running the program *********** NOTE   ***************
 -V allows output of all values during run time : not recommended
 unless you desire to inspect the gizzards of the beast

 -X0$xstart4resampling_source \
Synseis \
-S$input_source_filename \
-CZ$reflec_coef_depth \
-CT$reflec_coef_time \
-I$time_sampling_interval \
-IZ$depth_sampling_interval \
-LD$reg_density_file \
-LV$reg_velocity_file \
-Ro$output_source \
-Z$zrhov_filename \
-W$water_depth
 -V
 

NOTEs V 2.0

Examples

SEISMIC UNIX NOTES

CHANGES and their DATES

   upgrade ss plot to use suxwigb

Instantiate classes

       Create a new version of the package
       Personalize to give it a new name if you wish

Declare

  local variables 

default values for geomtery of boxes

Get configuration information input hash and an array reference

set the different parameters includes sugain variables

Set a2b

Set up

        data_in parameter values

Set a2b

Set a2b

Set a2b

Set a2b

Set sed

Set Synseis

DEFINE FLOW(S)

Reading and plotting data and synthetic seismogram

RUN FLOW(S)

flow 1 to create ss and accompanying files flow 2 create zrhoreg.bin

plot zrho_reg.bin

 xgraph 

Set

 xgraph plot zvreg.bin

Set

 xgraph plot reflection_coef_depth.bin

Set

 xgraph plot reflection_coef_time.bin

Plot

synthetic seismogram

LOG FLOW(S) TO SCREEN AND FILE