Psmeca

From GMT Türkiye Wiki
(Sürümler arası farklar)
Jump to: navigation, search
10. satır: 10. satır:
  
  
−
 
+
−
 
+
== Harvard CMT çizimi ==
−
 
+
Bunun için -Sd opsiyonu kullanılır
−
== Odak mekanizması animasyonu ==
+
−
[[Resim:focal_animate.gif]]
+
  
 
<bash>
 
<bash>
−
#!/bin/csh
+
psmeca << end -JM6 -R30/32/34/36 -Sd.3 -W2  -E255/255/255 -G0/0/255 >> harvard.ps
−
 
+
−
gmtset PAPER_MEDIA A6
+
−
 
+
−
# rake için bir while loop yapalım -180 +180 arası
+
−
set n = -180
+
−
while ($n <= 180)
+
−
 
+
−
# 0-360 derece cinsinden rake değerini al ve buna göre fay tipini belirle
+
−
# Aki ve Richard konvensiyonu
+
−
@ m = $n + 180
+
−
 
+
−
 
+
−
# Rake değerine göre fayın tipini belirleyip bunu fay isimli değişkene atıyalım.
+
−
# Ancak bunun tek satır olması gerekiyor, dikkat edin.
+
−
set fay = `echo $m | awk '{ if ($1<=10 && $1>=0) fay="right-lateral" ; else if ($1<80 && $1>10)
+
−
fay="oblique right-lateral normal"; else if ($1<=100 && $1>=80)fay="normal"; else if ($1<170 && $1>100)
+
−
fay="oblique left-lateral normal";else if ($1<=190 && $1>=170) fay="left-lateral" ; else if ($1<=260 &&
+
−
$1>190) fay="left-lateral reverse" ; else if ($1<=280 && $1>=260) fay="reverse" ; else if ($1<350 &&
+
−
$1>280) fay="oblique right-lateral reverse" ; else if ($1<=360 && $1>=350)  fay="right-lateral" ;
+
−
print fay}'`
+
−
 
+
−
# 90, 80, 60, 40, ve 20 derece eğimli düzlemler için mekanizma çizimi.
+
−
psmeca << son  -Sa0.4 -C -W -R39./40.1/38./38.9 -JM3  -G255/0/0 -P -K  > ! $m.ps
+
−
39.45 38.1 9.30505 45 90 $n 6
+
−
39.45 38.25 9.30505 45 80 $n 6
+
−
39.45 38.4 9.30505 45 60 $n 6
+
−
39.45 38.55 9.30505 45 40 $n 6
+
−
39.45 38.70 9.30505 45 20 $n 6
+
−
son
+
−
 
+
−
 
+
−
# rake değerini ve fay eğimlerini yanlara yazdır.
+
−
pstext << son -O -K -R -JM >> $m.ps
+
−
39.6 38.10 14 0 9 5 $n 
+
−
39.6 38.25 14 0 9 5 $n 
+
−
39.6 38.40 14 0 9 5 $n 
+
−
39.6 38.55 14 0 9 5 $n
+
−
39.6 38.70 14 0 9 5 $n
+
−
39.25 38.10 14 0 9 5 90 
+
−
39.25 38.25 14 0 9 5 80 
+
−
39.25 38.40 14 0 9 5 60 
+
−
39.25 38.55 14 0 9 5 40
+
−
39.25 38.70 14 0 9 5 20
+
−
son
+
−
 
+
−
# Fayın tipini üstte yazdır.
+
−
pstext << son -O -K -R -JM -G0/0/255>> $m.ps
+
−
39.55 38.87 12 0 9 6 $fay
+
−
son
+
−
 
+
−
# Eğim ve Rake başlığını yazdır
+
−
pstext << son -O -R -JM  -G255/0/0>> $m.ps
+
−
39.6  38.80 14 0 9 5 Rake
+
−
39.25  38.80 14 0 9 5 Dip
+
−
son
+
−
 
+
−
 
+
−
# rake değerini bir artır
+
−
@ n = $n + 1
+
−
end
+
−
 
+
−
# Imagemagic ile plot edilen şekilleri gif animasyonuna çevir
+
−
convert -delay 1 -dispose 2 -loop `ll [0-9]*ps | \
+
−
awk '{print $9}' | sort -g | awk '{printf "%s", $1" "}'` focal_animate.gif
+
−
 
+
−
# temizlik
+
−
\rm *ps
+
−
 
+
−
</bash>
+
−
 
+
−
 
+
−
 
+
−
== Harvard moment tensör çizimi ==
+
−
 
+
−
[[Resim:harvardcmt.jpg]]
+
−
 
+
−
<bash>
+
−
# Harvard sitesinden indirirken formatı GMT seçerseniz
+
−
# aşağıdaki formatta verilir tensör. Bu da -Sd seçeneği ile plot edilir.
+
−
 
+
−
gmtset PAPER_MEDIA A4 MEASURE_UNIT inch
+
−
\rm harvard.ps
+
−
pscoast -Df -W  -R20/36/35/42 -JM8 -B5  -K > harvard.ps
+
−
 
+
−
psmeca << end -JM -R  -O -Sd.3 -W2  -E255/255/255 -G0/0/255 >> harvard.ps
+
 
31.48 35.79 59 0.67 -2.86 2.19 -0.26 -2.33 0.34 24 X Y 060177A
 
31.48 35.79 59 0.67 -2.86 2.19 -0.26 -2.33 0.34 24 X Y 060177A
−
27.94 35.18 85 0.60 -0.30 -0.30 0.63 1.15 0.98 24 X Y 112877A
 
−
25.45 34.19 42 0.59 -0.90 0.31 2.13 -0.21 -0.32 24 X Y 030778D
 
−
26.40 37.58 15 -4.61 5.45 -0.84 1.90 2.35 -0.85 24 X Y 061479A
 
−
26.77 38.46 15 -0.87 1.19 -0.32 0.32 0.36 -0.28 24 X Y 061679A
 
−
29.19 39.44 15 -0.80 1.16 -0.36 0.43 -0.12 -0.53 24 X Y 071879B
 
−
26.75 35.29 15 -2.02 0.55 1.47 -0.60 2.16 1.44 24 X Y 072379A
 
−
27.89 35.75 94 -3.89 6.47 -2.58 8.55 -0.27 1.85 23 X Y 082279C
 
−
30.16 35.99 15 -1.60 2.53 -0.93 7.11 0.44 -0.12 24 X Y 050280C
 
−
26.72 34.12 33 0.81 -1.80 0.99 1.34 -0.72 -0.20 24 X Y 052286B
 
−
29.16 34.65 15 -6.82 7.91 -1.08 2.63 3.71 2.10 23 X Y 100286A
 
−
28.71 37.31 15 -2.14 3.31 -1.17 -0.53 -1.51 -1.53 24 X Y 101186A
 
−
26.17 40.17 15 -0.10 0.51 -0.41 0.21 -0.18 -0.76 24 X Y 070603G
 
−
25.86 40.33 18 -0.01 0.31 -0.29 -0.63 0.09 -1.71 23 X Y 070903G
 
−
28.77 37.88 15 -0.81 0.98 -0.17 0.52 -0.34 -0.52 24 X Y 072303A
 
−
29.05 38.03 15 -0.94 1.86 -0.92 0.13 -0.42 -0.28 24 X Y 072603F
 
−
27.00 35.93 160 1.58 0.71 -2.29 -0.82 6.44 5.59 23 X Y 091303E
 
 
end
 
end
  
−
xv harvard.ps
 
  
 
</bash>
 
</bash>

Sayfanın 14:49, 2 Nisan 2009 tarihindeki hali

Aki ve Richard konvensiyonu

  • Sa parametresi ile seçilir. Bunu ölçek ve metin boyu takip edebilir.
  • Kolonlar: Boylam, Enlem, Derinlik, strike, eğim-dip, rake, Magnitüd, Topun çizdirileceği boylam, Topun çizdirileceği enlem, Metin 
 
psmeca -Sa1/10 -R -J -O -C -W3 -G0/0/0 -W2/0/0/0 -E255/255/255 << END>> x1.ps
35.05696  40.8025 5 43.869 89.9 179 5 34.9 40.88 Mw=5
END
 


Harvard CMT çizimi

Bunun için -Sd opsiyonu kullanılır 

 
psmeca << end -JM6 -R30/32/34/36 -Sd.3 -W2  -E255/255/255 -G0/0/255 >> harvard.ps
31.48 35.79 59 0.67 -2.86 2.19 -0.26 -2.33 0.34 24 X Y 060177A
end
 
 
 

--Ziyadin 08:44, 2 Nisan 2009 (CEST)

Retrieved from "https://atag.itu.edu.tr/w/index.php?title=Psmeca&oldid=1655"
Personal tools