#!/bin/sh
### ====================================================================
### Synthetic seismograms in 2-D model UNCONFORMITY (Unix)
### ====================================================================
### Green function:
### ~~~~~~~~~~~~~~~
### 2-D two-point ray tracing ('crt.for' must not include 'scropc.for'):
echo "'u2d-crt.dat' /"|crt
### ....................................................................
### Converting output of CRT into Green function:
echo "' ' ' ' 's01.out' 's01i.out' 'green.out' /"|green
### --------------------------------------------------------------------
### Renaming already created PostScript plots:
### ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
for plot in `ls plot*.ps 2>/dev/null`
do
  file=`echo $plot|sed 's/plot/pl0t/'`
  mv $plot $file
done
### --------------------------------------------------------------------
### Time-domain seismograms:
### ~~~~~~~~~~~~~~~~~~~~~~~~
### Generating the source time function:
echo "'u2d-ss.dat' ' ' 'u2d-t0.gse' 'u2d-sst.out' 0 /"|ss
for plot in `ls plot*.ps`
do
  file=`echo $plot|sed 's/plot/sst0/'`
  mv $plot $file
done
### ....................................................................
### Calculating the synthetic seismograms (1st component):
echo "'green.out' 'source.dat' ' ' 'u2d-t0.gse' 'u2d-t1.gse' 1 /"|greenss
### Calculating the synthetic seismograms (3rd component):
echo "'green.out' 'source.dat' ' ' 'u2d-t0.gse' 'u2d-t3.gse' 3 /"|greenss
### --------------------------------------------------------------------
### Frequency-domain seismograms:
### ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
### Calculating the response functions (all components, no attenuation):
echo "1024 4. 200. 0.001 1.000 0. 100." >u2d-fr.tmp
echo "'green.out' 'source.dat' 'u2d-fr.tmp' ' ' 'u2d-rff.out' 0 /"|greenss
### Calculating the response functions (all components, noncausal att.):
echo "1024 4. 200. 0.001 1.000 1.   0." >u2d-fr.tmp
echo "'green.out' 'source.dat' 'u2d-fr.tmp' ' ' 'u2d-rfn.out' 0 /"|greenss
### Calculating the response functions (all components, attenuation):
echo "1024 4. 200. 0.001 1.000 1. 100." >u2d-fr.tmp
echo "'green.out' 'source.dat' 'u2d-fr.tmp' ' ' 'u2d-rfq.out' 0 /"|greenss
### ....................................................................
### Calculating the synthetic seismograms (1st component, no attenuation):
echo "'u2d-ss.dat' 'u2d-rff.out' 'u2d-f1.gse' 'u2d-ssf1.out' 1 /"|ss
for plot in `ls plot*.ps`
do
  file=`echo $plot|sed 's/plot/ssf1/'`
  mv $plot $file
done
### Calculating the synthetic seismograms (3rd component, no attenuation):
echo "'u2d-ss.dat' 'u2d-rff.out' 'u2d-f3.gse' 'u2d-ssf3.out' 3 /"|ss
for plot in `ls plot*.ps`
do
  file=`echo $plot|sed 's/plot/ssf3/'`
  mv $plot $file
done
### Calculating the synthetic seismograms (1st component, noncausal att.):
echo "'u2d-ss.dat' 'u2d-rfn.out' 'u2d-n1.gse' 'u2d-ssn1.out' 1 /"|ss
for plot in `ls plot*.ps`
do
  file=`echo $plot|sed 's/plot/ssn1/'`
  mv $plot $file
done
### Calculating the synthetic seismograms (3rd component, noncausal att.):
echo "'u2d-ss.dat' 'u2d-rfn.out' 'u2d-n3.gse' 'u2d-ssn3.out' 3 /"|ss
for plot in `ls plot*.ps`
do
  file=`echo $plot|sed 's/plot/ssn3/'`
  mv $plot $file
done
### Calculating the synthetic seismograms (1st component, attenuation):
echo "'u2d-ss.dat' 'u2d-rfq.out' 'u2d-q1.gse' 'u2d-ssq1.out' 1 /"|ss
for plot in `ls plot*.ps`
do
  file=`echo $plot|sed 's/plot/ssq1/'`
  mv $plot $file
done
### Calculating the synthetic seismograms (3rd component, attenuation):
echo "'u2d-ss.dat' 'u2d-rfq.out' 'u2d-q3.gse' 'u2d-ssq3.out' 3 /"|ss
for plot in `ls plot*.ps`
do
  file=`echo $plot|sed 's/plot/ssq3/'`
  mv $plot $file
done
### --------------------------------------------------------------------
### Plotting seismograms:
### ~~~~~~~~~~~~~~~~~~~~~
### Plotting the synthetic seismograms (1st component, time-domain):
echo "'u2d-sp-1.dat' 'u2d-t1.gse' 1 /"|sp
mv plot00.ps u2d-t1.ps
### Plotting the synthetic seismograms (3rd component, time-domain):
echo "'u2d-sp-1.dat' 'u2d-t3.gse' 3 /"|sp
mv plot00.ps u2d-t3.ps
### Plotting the synthetic seismograms (1st component, no attenuation):
echo "'u2d-sp-1.dat' 'u2d-f1.gse' 1 /"|sp
mv plot00.ps u2d-f1.ps
### Plotting the synthetic seismograms (3rd component, no attenuation):
echo "'u2d-sp-1.dat' 'u2d-f3.gse' 3 /"|sp
mv plot00.ps u2d-f3.ps
### Plotting the synthetic seismograms (1st component, noncausal att.):
echo "'u2d-sp-4.dat' 'u2d-n1.gse' 1 /"|sp
mv plot00.ps u2d-n1.ps
### Plotting the synthetic seismograms (3rd component, noncausal att.):
echo "'u2d-sp-4.dat' 'u2d-n3.gse' 3 /"|sp
mv plot00.ps u2d-n3.ps
### Plotting the synthetic seismograms (1st component, attenuation):
echo "'u2d-sp-4.dat' 'u2d-q1.gse' 1 /"|sp
mv plot00.ps u2d-q1.ps
### Plotting the synthetic seismograms (3rd component, attenuation):
echo "'u2d-sp-4.dat' 'u2d-q3.gse' 3 /"|sp
mv plot00.ps u2d-q3.ps
### --------------------------------------------------------------------
### Rerenaming old PostScript plots:
### ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
for plot in `ls pl0t*.ps 2>/dev/null`
do
  file=`echo $plot|sed 's/pl0t/plot/'`
  mv $plot $file
done
### --------------------------------------------------------------------
### In the output PostScript files 'u2d-*.ps' of the 'sp' program, the
### Setup section may be edited in the following way:
###  .48 setlinewidth 1 setlinecap 1 setlinejoin
### For the dvips TeX driver, line
### %%BoundingBox: (atend)
### should be replaced with the %%BoundingBox line from the %%Trailer of
### the PostScript file.
###
### Manual editing of the PostScript seismograms using the PE2 editor:
### pause "Press Crtl-C not to edit the output PostScript files"
### vi u2d-t1.ps
### vi u2d-t3.ps
### vi u2d-f1.ps
### vi u2d-f3.ps
### vi u2d-n1.ps
### vi u2d-n3.ps
### vi u2d-q1.ps
### vi u2d-q3.ps
###
### List of output PostScript figures:
### ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
###           Synthetic seismograms:
### u2d-t1.ps Time-domain, no attenuation, 1st component.
### u2d-t3.ps Time-domain, no attenuation, 3rd component.
### u2d-f1.ps Frequency-domain, no attenuation, 1st component.
### u2d-f3.ps Frequency-domain, no attenuation, 3rd component.
### u2d-n1.ps Frequency-domain, noncausal attenuation, 1st component.
### u2d-n3.ps Frequency-domain, noncausal attenuation, 3rd component.
### u2d-q1.ps Frequency-domain, attenuation, 1st component.
### u2d-q3.ps Frequency-domain, attenuation, 3rd component.
###
### Auxiliary PostScript figures:
### ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
### sst0??.ps Sketch of source time function (output of SS).
###           Sketches of synthetic seismograms (output of SS):
### sst1??.ps Time-domain, no attenuation, 1st component.
### sst3??.ps Time-domain, no attenuation, 3rd component.
### ssf1??.ps Frequency-domain, no attenuation, 1st component.
### ssf3??.ps Frequency-domain, no attenuation, 3rd component.
### ssn1??.ps Frequency-domain, noncausal attenuation, 1st component.
### ssn3??.ps Frequency-domain, noncausal attenuation, 3rd component.
### ssq1??.ps Frequency-domain, attenuation, 1st component.
### ssq3??.ps Frequency-domain, attenuation, 3rd component.
### ====================================================================