Elastik deformasyon çizimi
From GMT Türkiye Wiki
set atim = 3 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 while ($W <= 25)
set renk = `perl -e 'printf "%d %d %d ", int(rand(255)), int(rand(255)), int(rand(100))' | awk '{print $1"/"$2"/"$3}'`
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 echo $renk end
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