Elastik deformasyon çizimi
From GMT Türkiye Wiki
(Sürümler arası farklar)
| 1. satır: | 1. satır: | ||
''' Doğrultu atımlı faylar boyunca intersismik dönemde ve deprem anında meydana gelen yamulmanın hesaplanıp çizdirilmesi''' | ''' Doğrultu atımlı faylar boyunca intersismik dönemde ve deprem anında meydana gelen yamulmanın hesaplanıp çizdirilmesi''' | ||
| − | + | # Kayma hızı (intersismik cm; eşsismik metre) | |
set atim = 3 | set atim = 3 | ||
| + | |||
| + | # kitlenme derinliği | ||
set W = 5 | set W = 5 | ||
psbasemap -R-100/100/-2/2 -B20f10:"Faya dik mesafe (km)":Sn/.5:"Kayma hizi (cm/y)":We:."Deprem arasi donem": -JX22/6 -K> ! x.ps | psbasemap -R-100/100/-2/2 -B20f10:"Faya dik mesafe (km)":Sn/.5:"Kayma hizi (cm/y)":We:."Deprem arasi donem": -JX22/6 -K> ! x.ps | ||
| + | |||
| + | #Kitlenme derinliğini 25 km ye gelene kadar 5 km arttır | ||
while ($W <= 25) | while ($W <= 25) | ||
| + | # Her bir eğri için random renk üret | ||
set renk = `perl -e 'printf "%d %d %d ", int(rand(255)), int(rand(255)), int(rand(100))' | awk '{print $1"/"$2"/"$3}'` | set renk = `perl -e 'printf "%d %d %d ", int(rand(255)), int(rand(255)), int(rand(100))' | awk '{print $1"/"$2"/"$3}'` | ||
| + | # Plot et | ||
awk 'BEGIN {a=3.1416/180; for (i=-200;i<=200;i=i+1) print (i,('$atim'/3.141 *(atan2(i,'$W'))))}' |\ | awk 'BEGIN {a=3.1416/180; for (i=-200;i<=200;i=i+1) print (i,('$atim'/3.141 *(atan2(i,'$W'))))}' |\ | ||
psxy -R -W5/$renk -JX -M -O -K>> x.ps | psxy -R -W5/$renk -JX -M -O -K>> x.ps | ||
@ W = $W + 5 | @ W = $W + 5 | ||
| − | |||
end | end | ||
| + | #Deprem anı deformasyonu plot et | ||
awk 'BEGIN {a=3.1416/180; for (i=-200;i<200;i=i+1) if (i <0) print i,('$atim'/3.141 *(atan2(i,'$W')))+(('$atim')/2); | awk 'BEGIN {a=3.1416/180; for (i=-200;i<200;i=i+1) if (i <0) print i,('$atim'/3.141 *(atan2(i,'$W')))+(('$atim')/2); | ||
else print i,('$atim'/3.141 *(atan2(i,'$W')))-(('$atim')/2)}' | psxy -R-100/100/-2/2 -B50:"":Sn/.5:"Kayma (m)":We:."100 sene sonra deprem ani": | else print i,('$atim'/3.141 *(atan2(i,'$W')))-(('$atim')/2)}' | psxy -R-100/100/-2/2 -B50:"":Sn/.5:"Kayma (m)":We:."100 sene sonra deprem ani": | ||
-W10/'red' -JX22/6 -O -Y10>> x.ps | -W10/'red' -JX22/6 -O -Y10>> x.ps | ||
Sayfanın 23:43, 16 Haziran 2009 tarihindeki hali
Doğrultu atımlı faylar boyunca intersismik dönemde ve deprem anında meydana gelen yamulmanın hesaplanıp çizdirilmesi
# Kayma hızı (intersismik cm; eşsismik metre) set atim = 3
# kitlenme derinliği set W = 5 psbasemap -R-100/100/-2/2 -B20f10:"Faya dik mesafe (km)":Sn/.5:"Kayma hizi (cm/y)":We:."Deprem arasi donem": -JX22/6 -K> ! x.ps #Kitlenme derinliğini 25 km ye gelene kadar 5 km arttır while ($W <= 25)
# Her bir eğri için random renk üret
set renk = `perl -e 'printf "%d %d %d ", int(rand(255)), int(rand(255)), int(rand(100))' | awk '{print $1"/"$2"/"$3}'`
# Plot et
awk 'BEGIN {a=3.1416/180; for (i=-200;i<=200;i=i+1) print (i,('$atim'/3.141 *(atan2(i,'$W'))))}' |\
psxy -R -W5/$renk -JX -M -O -K>> x.ps
@ W = $W + 5
end
#Deprem anı deformasyonu plot et
awk 'BEGIN {a=3.1416/180; for (i=-200;i<200;i=i+1) if (i <0) print i,('$atim'/3.141 *(atan2(i,'$W')))+(('$atim')/2);
else print i,('$atim'/3.141 *(atan2(i,'$W')))-(('$atim')/2)}' | psxy -R-100/100/-2/2 -B50:"":Sn/.5:"Kayma (m)":We:."100 sene sonra deprem ani":
-W10/'red' -JX22/6 -O -Y10>> x.ps
