Psmeca

From GMT Türkiye Wiki
(Sürümler arası farklar)
Jump to: navigation, search
 
(2 kullanıcı tarafından yapılan 8 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 mekanizması animasyonu ==
 
[[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, YMetin
# 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 çizimi ==
+
+
  
 
<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.
+
 
+
gmtset PAPER_MEDIA A4 MEASURE_UNIT inch
+
\rm harvard.ps
+
pscoast -Df -W  -R20/36/30.1/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>

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)

Personal tools