do-spherical-vids.sh

#!/bin/bash
##
##  "DO"-script for a Radiative3D run:
##
##  Copy and edit this file to easily specify (and remember)
##  parameters used for a particular run.
##
##  If $1 is set, it is used as an extra identifier token in the
##  output directory name.
##
##  (Note: Do not edit line 3.  It is used by subsequent scripts to
##  identify this file as a "do script". (See wikify.sh.))
##
source scripts/do-fundamentals.sh
##
##  This DO-SCRIPT sets up a VIDEO run based on the Crust Pinch WCG model
##

SIMTARGET="video"             # Choice: 'waveform' or 'video'. Affects
                              # defaults not otherwise specified.
MODIDX=16  # Spherical Earth  # Model Index: Selects from custom coded models.
                              # 1: Lop Nor (base or moho depends on COMPARGS),
                              # 5: North Sea Crust Pinch model,
                              # 8: Crust Upthrust model

#SOURCETYP=$(EventParams "Selby" 0)
        # Lookup source parameters. Choices include 'expl', 'eq', 'Selby', etc.
        # Second parameter is usually isofrac (range [-1.0,1.0] or [-90,90]).
SOURCETYP=EXPL #SDR,22.5,90,0       # Or, specify source directly.
SOURCELOC=0,0,-10             # Event Location: Null Island: 0,0,elev
                              #                 North Pole:  10007,0,elev


FREQ=2.0                      # Phonon frequency to model
NUMPHONS=2K                   # Number of phonons to spray
RECTIME=12000                 # Recording duration of seismometers.
GATHER=40.0                   # Terminal gather radius, in kilometers.

SCAT1=0.8,0.01,0.5,0.2,50     # Scat Args (nu,eps,a,kappa) and Q in sedi's
SCAT2=0.8,0.01,0.5,0.3,1000   # Scat Args (nu,eps,a,kappa) and Q in crust
SCAT3=0.8,0.01,0.7,0.5,300    # Scat Args (nu,eps,a,kappa) and Q in mantle
                              # (Note: Q's specified are Q_s values.
                              # Q_p is computed from assumption that
                              # Q_kappa is infinite.)
COMPARGS=$SCAT1,$SCAT2,$SCAT3

MFPOVERRIDE=--overridemfp=50,40     # Uncomment to pin P,S Mean-Free-Paths.
NODEFLECT=--nodeflect               # Uncomment to make scattering
                                    # non-delfectionary.

ADDITIONAL="--ocsraw"                 # Additional params. (Such as --ocsraw,
                              # or --earthrad=xxxx) Radii: Earth: 6371 km,
                              # Mars: 3389, Mercury: 2440, Earth's
                              # moon: 1737

SEISORIG1=0,0,0               # Array locations: (Range, Azimuth, Z)
SEISORIG2=0,90,0              #
SEISDEST1=1000,0,0            #
SEISDEST2=1000,90,0           #
SEIS1=--seis-p2p=$SEISORIG1,$SEISDEST1,1.0,2.0,$GATHER,16
SEIS2=--seis-p2p=$SEISORIG2,$SEISDEST2,1.0,2.0,$GATHER,16
BINSIZE=60.0                  # Large time-bins, not useful for trace analy-
                              # sis, but makes seis files smaller. (We only
                              # want them to plot a map as the video
                              # backdrop.)

## RUN THE SIMULATION:
##
##  Set up output directory and run the sim:
##

PopDefaults $SIMTARGET        ## Defined in do-fundamentals.sh
CheckBuildStatus              ##  ''
CreateOutputDirectory $@      ##  ''
RunSimulation                 ##  ''

echo Begin Figure Generation: >> "$LOGFILE"
rwd=`pwd`       # Switch to output directory - the rest of our work
cd "$outdir"    # will be done there.
## (Do not edit/remove this comment block.)
#!/bin/bash
## ___FIG_GEN_START___
##

##                                      # Note: this part of the do-script will
## GENERATE VIDEOS:                     #   be carved off into 'do-figsonly.sh'
##                                      #   to facilitate re-running and/or
                                        #   customizing figure output after the
                                        #   initial run.

# Organize seisfiles
mkdir -p seisfiles  # Stash seis files in a subdirectory (aesthetics)
mv seis_???.octv seisfiles/ 2> /dev/null || : # (suppress error msg and code)

# Text values for Wiki script:
cat > wikify.incl <<EOF
  STA="NOR"
  STB="EAS"
  STC="xxx"
  MODELCODE="SPH"
EOF

# Determine RunID from name of do-script:   (for labeling figs)
SHFILE=RUNID
for file in *.sh; do
    if grep -qx '##  "DO"-script for a Radiative3D run:' "$file"
    then                # Catch ONLY file with R3D do marker
        SHFILE="$file"  # Also catches no more than one do- file
    fi                  #
done
RUNID=${SHFILE%.sh}
RUNID=${RUNID#do-}
echo "Found RunID: $RUNID"

# Some octave/gnuplot set-up and platform-independence stuff:
alias octave="octave --no-window-system"    # Prevents trouble on some systems.
export GNUTERM="dumb"                   # Ascii, but we don't disply so is OK.
shopt -s expand_aliases                 # Aliases in scripts... needs enabled.
reverser="tac"                                  # Linux has tac (backwords cat)
[ "`uname`" == "Darwin" ] && reverser="tail -r" # But OS X has 'tail -r',...
revtail() { $reverser | tail $@ | $reverser; }  # 'head' with neg line count
                                                # does not work on OS X, so
                                                # use revtail in its place.

##
## Pre-process:
##
echo "Processing output for visualization:"
./preprocess.sh
cat stdout.txt | \
    grep -A 10000 "#  R3D_GRID:" | \
    grep -B 10000 "#  END R3D_GRID" | \
    revtail -n +2 | \
    tail -n +9 | \
    sed 's/\*\*\*/000/g' > griddump.txt

##
## Elevation (Side-View) Video:
##
if false; then
echo "Making Scatter-Vid Elevation View:"
mkdir -p framecache
octave -qf <<EOF
  scattervid_model(0,
            struct(
              "destination","seisfiles/seis_031.octv", # 
              "gridfile",  "griddump.txt",             # A 3x1xN Cylinder grid
              "gridpair",  [1 2],                      # (1=LOP, 2=MAK, 3=WUS)
              "mparams",   "out_mparams.octv",         # Metadata
              "window",    [-0.05 1.4 10 -130],        # [LRTB], L,R fractional
              "seispat",   "seisfiles/seis_%03d.octv", # Seismometers files
              "seisrange", [],                         # ...which ones to plot
              "txtkey",    {{"LOP","MAK","WUS"}},      # Location names
              "loctext",   "Lop Nor",                  # Source verbose name
              "base",      [492.310, -263.650, 1.050], # Location of LOP
              "duration",  1000,                        # Time, (if<PhononTTL)
              "numframes", 100                         # Total frames in video
              )
            );
EOF
mv framecache/movie.mp4 scatvid-elev.mp4
if [ -d framecache-elev ]; then
    rm -rf framecache-elev-old
    mv framecache-elev framecache-elev-old
fi
mv framecache framecache-elev
fi


##
## Plan-View (from above) Video:
##
echo "Making Scatter-Vid Above View:"
mkdir -p framecache
octave -qf <<EOF
  scattervid_axial(0,
            struct(
              "seispat",   "seisfiles/seis_%03d.octv",
              "seisrange", [],
              "loctext",   "Null Island",
              "numframes", 600,
              "window", [-8600 8600 6400 -6400]
            )
          );
EOF
mv framecache/movie.mp4 scatvid-above.mp4
if [ -d framecache-above ]; then
    rm -rf framecache-above-old
    mv framecache-above framecache-above-old
fi
mv framecache framecache-above


## (Do not edit/remove this comment block.)
## ___FIG_GEN_END___
##
cd "$rwd"
echo End Figure Generation: `date +"%Y.%m.%d-%H:%M:%S"` >> "$LOGFILE"
echo '///'
echo "Figure generation complete. To re-run figures, run 'do-figsonly.sh' in"
echo "the output directory, which may be modified to customize figure output."
echo "Output has been placed in $outdir."
echo "Logfile contents:"
cat "$outdir"/logfile
## END
##