Tiden os, og den verden vi lever i!

Se stjernebilledet af himlen på et vilkårligt klokkeslet. HIMMELRUM

ing GMT-JOB ----- set wpath=/data10/shusen/tomochile set inpath=$wpath/3dtomo/FATOMO/COARSE/r1/VP # use colorscale velabs4.5-9.cpt set cpt=$wpath/gmt/velabs4.5-9.6.cpt set psfile=$wpath/gmt/3D-PSFILES/latsec1-3b.ftm.ps # file with final hypocenter locations set cnvfile=$inpath2/hypo.xyz # file containing grid nodes set nodefile=$inpath/latvert_nodes.xyz set title=" " # variables for earthquakes & which latsec set latsec=1 set lat=24.00 set lat1=-24.25 set lat2=-23.75 set Rstring=-71.2727/-69.2168/-3.0/79.0 set Jstring='x2.55/-0.025' set infile1=$inpath/latsec$latsec.xyz # column where to find Vp in the tomo2gmt outputfile set pos=6 psbasemap -J$Jstring -R$Rstring -Baf0.5:"longitude [degrees]":/a10f5:"depth [km]":WenS -K -V -P -X1 -Y2 >$psfile # create grid and plot with colorscale nawk 'substr($1,1,1)!="#" {print $3,$5,$pos}' pos=$pos $infile1 | \ xyz2grd -I$Istring -Gvpabs.grd -R -V # mask for simulps ( all nodes with DWS<10 ) nawk 'substr($1,1,1)!="#" {if ($11<1.0) print $3,$5,"9999"; \ if ($11>=1.0) print $3,$5,"0"}' $infile1 |\ xyz2grd -I$Istring -Gclip.grd -R -V grdmath vpabs.grd clip.grd ADD = plot.grd grdimage plot.grd -J$Jstring -R$Rstring -C$cpt -O -K -P -V >> $psfile # trace vp contours grdcontour vpabs.grd -J$Jstring -R -C0.5 -L5/8 -A0.5f10 -Wc1/200 -G3/20 -S10 -O -K -P -V >>$psfile # plot nodes with crosses nawk '{print $1,$2}' $nodefile |\ psxy -J$Jstring -R$Rstring -Sx0.01 -W25/25/25/0.5ta -O -K -P -V >>$psfile # plot quakes nawk '$2lat1 {print $1,$3}' lat1=$lat1 lat2=$lat2 $cnvfile |\ psxy -J$Jstring -R -Sc0.05 -G255 -L -W0.5 -P -O -K >>$psfile pstext -J$Jstring -R -O -P -K -V -N <>$psfile -71 78 10 0 1 1 Latitude Depth Section along $lat S END #plot colorscale psscale -C$cpt -D4.5/1.5/5/0.25h -L -B:."absolute velocity [km/s]": \ -O -P -X-2 -Y-2.3 >> $psfile ##-------------- printing gs $psfile # cleanup /bin/rm vpabs.grd clip.grd plot.grd sprd.grd echo "-- printing file, to abort --" pause lpr -Pphaser560-8d_213 $psfile exit