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