Psmeca

From GMT Türkiye Wiki
(Sürümler arası farklar)
Jump to: navigation, search
 
(2 kullanıcı tarafından yapılan 12 ara revizyon gösterilmiyor)
1. satır: 1. satır:
−
==Aki Convention==
+
Fay düzlemi (kuvvet çifti veya moment tensör) çözümlerinin harita üzerine çizdirilmesini sağlar. Çeşitli formatlarda veri girişine imkan verir. Aki ve Richard parametre formatı, Harvard-CMT formatı, moment tensör bileşenleri gibi verileri okuyabilmektedir.
  
−
* 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
 
−
<bash>
 
−
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
 
−
</bash>
 
  
 +
'''ÖNEMLİ NOT:''' GMT 4.4.0 versiyonunda bazı sürümlerdeki hata nedeniyle mekanizma küreleri yanlış çizilmektedir. Konu ile ilgili O. Tan'nın GMT-Türkiye  mail grubuna gönderdiği mesaj şöyledir:
 +
<pre>
 +
Merhaba,
  
 +
psmeca  ile ilgili problem, programin  son versiyonundaki kaynak kodunda yapilan bir hatadan kaynaklaniyor.
 +
Bu durumla karşılaşan arkadaşların, psmeca'nın bir önceki versiyonunu kullanmaları uygun olur. Aksi halde
 +
utilmeca.c  içinde hatalı yapılan radyan-derece dönüşümü  ( ilk satır 497 dir, benzer hatali satirlar vardir)
 +
str2 = d_atan2(sp2, cp2) * R2D;
 +
şeklinde düzeltilmeli ve tekrar derlenmeli  :)
  
 +
Konu ile iligili GMT listesindeki mailleriden biri şu adreste bulunabilir.
 +
https://listserv.hawaii.edu/cgi-bin/wa?A2=ind0902&L=GMT-HELP&D=0&P=77915
  
 +
iyi çalışmalar
 +
onur
 +
</pre>
  
−
== ODAK MEKANİZMA ANİMASYONU ==
 
−
[[Resim:focal_animate.gif]]
 
  
−
<bash>
+
== Aki ve Richard konvensiyonu ==
−
#!/bin/csh
+
  
−
gmtset PAPER_MEDIA A6
+
* Sa parametresi ile seçilir. Bunu ölçek ve metin boyu takip edebilir.  
−
 
+
* Gerekli veriler: Boylam, Enlem, Derinlik, Doğrultu, Eğim, Kayma Açısı, Büyüklük, X, Y,  Metin
−
# rake için bir while loop yapalım -180 +180 arası
+
X ve Y odak küresinin harita üzerindeki boylam ve enlem değerleridir.
−
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>
 +
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
 
</bash>
 
</bash>
  
−
 
+
== Harvard CMT çizimi ==
−
 
+
Bunun için -Sd opsiyonu kullanılır
−
'''== Harvard moment tensor plotu =='''
+
−
+
  
 
<bash>
 
<bash>
−
# Harvard sitesinden indirirken formatı GMT seçerseniz
+
psmeca << end -JM6 -R30/32/34/36 -Sd.3 -W2  -E255/255/255 -G0/0/255 > harvard.ps
−
# aşağıdaki formatta verilir tensör. Bu da -Sd seçeneği ile plot edilir.
+
−
 
+
−
 
+
−
\rm harvard.ps
+
−
pscoast -Df -W  -R20/36/30.1/45 -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>

09:52, 18 Mayıs 2009 itibarı ile sayfanın şu anki hâli

Fay düzlemi (kuvvet çifti veya moment tensör) çözümlerinin harita üzerine çizdirilmesini sağlar. Çeşitli formatlarda veri girişine imkan verir. Aki ve Richard parametre formatı, Harvard-CMT formatı, moment tensör bileşenleri gibi verileri okuyabilmektedir.


ÖNEMLİ NOT: GMT 4.4.0 versiyonunda bazı sürümlerdeki hata nedeniyle mekanizma küreleri yanlış çizilmektedir. Konu ile ilgili O. Tan'nın GMT-Türkiye mail grubuna gönderdiği mesaj şöyledir: 

Merhaba,

psmeca  ile ilgili problem, programin  son versiyonundaki kaynak kodunda yapilan bir hatadan kaynaklaniyor.
Bu durumla karşılaşan arkadaşların, psmeca'nın bir önceki versiyonunu kullanmaları uygun olur. Aksi halde 
 utilmeca.c  içinde hatalı yapılan radyan-derece dönüşümü   ( ilk satır 497 dir, benzer hatali satirlar vardir) 
	str2 = d_atan2(sp2, cp2) * R2D;
şeklinde düzeltilmeli ve tekrar derlenmeli  :)

Konu ile iligili GMT listesindeki mailleriden biri şu adreste bulunabilir.
https://listserv.hawaii.edu/cgi-bin/wa?A2=ind0902&L=GMT-HELP&D=0&P=77915

iyi çalışmalar
onur


[değiştir] Aki ve Richard konvensiyonu

  • Sa parametresi ile seçilir. Bunu ölçek ve metin boyu takip edebilir.
  • Gerekli veriler: Boylam, Enlem, Derinlik, Doğrultu, Eğim, Kayma Açısı, Büyüklük, X, Y, Metin

X ve Y odak küresinin harita üzerindeki boylam ve enlem değerleridir. 

 
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
 

[değiştir] 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=1869"
Personal tools