# Plotting travel times in 2-D model UNCONFORMITY # =============================================== # Input files required chk.pl: "data/u2d/" "u2d-mod.dat" chk.pl: "data/u2d/" "u2d-src.dat" chk.pl: "data/u2d/" "u2d-rec.dat" chk.pl: "data/u2d/" "u2d-crt.dat" chk.pl: "data/u2d/" "u2d-rpa.dat" chk.pl: "data/u2d/" "u2d-writ.dat" chk.pl: "data/u2d/" "u2d-ttc.dat" # Model UNCONFORMITY (Cormier and Mellen 1984) MODEL='u2d-mod.dat' # Source and receivers SRC='u2d-src.dat' REC='u2d-rec.dat' # Data filenames for CRT (for details see 'crt.for') DCRT='u2d-crt.dat' CODE='u2d-crt.dat' RPAR='u2d-rpa.dat' WRIT='u2d-writ.dat' # output files for plotting travel-time curves CRTLOG='u2d-log.out' # Kind of ray parameters (for details see 'init.for') INIPAR=2 # geographic-like spherical coordinates # Data common for travel-time and seismogram plotting (program SP) KODESP=3 SPCHRH=0.23 SPTMIN=0.300 SPTMAX=1.300 SPTLEN=16.00 SPTDIV=10 SPTSUB=1 SPXMIN=0.136 SPXMAX=-1.241 SPXLEN=13.77 SPXDIV=9 SPXSUB=3 # Linewidth for plotting (program SP) CALCOPS='0.48 setlinewidth' # Two-point ray tracing crt: # Generating input files CRTOUT for individual elementary waves echo.pl: "' ' 's01.out' 's01i.out' /" "> crt01.out" echo.pl: "' ' 's02.out' 's02i.out' /" "> crt02.out" echo.pl: "' ' 's03.out' 's03i.out' /" "> crt03.out" echo.pl: "' ' 's04.out' 's04i.out' /" "> crt04.out" echo.pl: "' ' 's05.out' 's05i.out' /" "> crt05.out" echo.pl: "' ' 's06.out' 's06i.out' /" "> crt06.out" echo.pl: "' ' 's07.out' 's07i.out' /" "> crt07.out" echo.pl: "' ' 's08.out' 's08i.out' /" "> crt08.out" echo.pl: "' ' 's09.out' 's09i.out' /" "> crt09.out" echo.pl: "' ' 's10.out' 's10i.out' /" "> crt10.out" echo.pl: "' ' 's11.out' 's11i.out' /" "> crt11.out" echo.pl: "' ' 's12.out' 's12i.out' /" "> crt12.out" echo.pl: "' ' 's13.out' 's13i.out' /" "> crt13.out" echo.pl: "' ' 's14.out' 's14i.out' /" "> crt14.out" echo.pl: "' ' 's15.out' 's15i.out' /" "> crt15.out" echo.pl: "' ' 's16.out' 's16i.out' /" "> crt16.out" echo.pl: "' ' 's17.out' 's17i.out' /" "> crt17.out" echo.pl: "' ' 's18.out' 's18i.out' /" "> crt18.out" echo.pl: "' ' 's19.out' 's19i.out' /" "> crt19.out" echo.pl: "' ' 's20.out' 's20i.out' /" "> crt20.out" echo.pl: "' ' 's21.out' 's21i.out' /" "> crt21.out" echo.pl: "' ' 's22.out' 's22i.out' /" "> crt22.out" # Generating files with travel times of individual elementary waves SRC=' ' NQ=4 KALL=0 KREC=1 CRTOUT='crt01.out' PTS='tt01.out' crtpts: CRTOUT='crt02.out' PTS='tt02.out' crtpts: CRTOUT='crt03.out' PTS='tt03.out' crtpts: CRTOUT='crt04.out' PTS='tt04.out' crtpts: CRTOUT='crt05.out' PTS='tt05.out' crtpts: CRTOUT='crt06.out' PTS='tt06.out' crtpts: CRTOUT='crt07.out' PTS='tt07.out' crtpts: CRTOUT='crt08.out' PTS='tt08.out' crtpts: CRTOUT='crt09.out' PTS='tt09.out' crtpts: CRTOUT='crt10.out' PTS='tt10.out' crtpts: CRTOUT='crt11.out' PTS='tt11.out' crtpts: CRTOUT='crt12.out' PTS='tt12.out' crtpts: CRTOUT='crt13.out' PTS='tt13.out' crtpts: CRTOUT='crt14.out' PTS='tt14.out' crtpts: CRTOUT='crt15.out' PTS='tt15.out' crtpts: CRTOUT='crt16.out' PTS='tt16.out' crtpts: CRTOUT='crt17.out' PTS='tt17.out' crtpts: CRTOUT='crt18.out' PTS='tt18.out' crtpts: CRTOUT='crt19.out' PTS='tt19.out' crtpts: CRTOUT='crt20.out' PTS='tt20.out' crtpts: CRTOUT='crt21.out' PTS='tt21.out' crtpts: CRTOUT='crt22.out' PTS='tt22.out' crtpts: # Data to control travel-time plotting (program SP) SP1='u2d-tt.ps' SP2=' ' SP3=' ' SS=' ' PTS=' ' FTT=' ' SPTTC='u2d-ttc.dat' KOLORTD=-1 SPSYMH=0.11 # Plotting travel times sp: # Comments: # ~~~~~~~~~ # 'u2d-t1.ps'... output PostScript file with travel-time curves. # Legenda: # black + unconverted P wave (both refracted and once reflected) # red + unconverted S wave (both refracted and once reflected) # green + once converted PS wave (both refracted and once reflected) # blue + once converted SP wave (both refracted and once reflected) # red x twice converted SPS wave (both refracted and once reflected) # blue x three time converted SPSP wave (once reflected)