Rabu, 22 Mei 2013

Cross Section of South Sumatra


#!/bin/csh
gmtset ANNOT_FONT_SIZE_PRIMARY 12p HEADER_FONT_SIZE 18p PLOT_DEGREE_FORMAT ddd:mm:ssF
#set region = 92/106/-7.5/6
set region = 96/103/-6/1
set size = M12c
set psfile = xxxx_all_3.ps
set grdfile=./indonesia.nc
set cptfile=color.cpt
makecpt -Csealand -Z > $cptfile
set AZ = 55
set pointboxlon_1 = 96.5
set pointboxlat_1 = -4.5
echo $pointboxlon_1 $pointboxlat_1 S1 >! S1.d
awk '{print $1+5.5, $2+4, $3"`" }' S1.d >! S11.d
cat S1.d S11.d >! garis_S11.d

#Draw bathymetry and basemap
grdgradient $grdfile -A30/270 -Gintens.grd -Nt0.30 -V
gmtset ANOT_FONT_SIZE 10
grdimage $grdfile -R$region -J$size -C$cptfile -Iintens.grd -B2g1/2g1wnSE -K -P -X0.7 -Y3.5 >! $psfile
pscoast -R$region -J$size -Dh -W1 -K -P -O >> $psfile

# Draw tectonical structures  ;   trech3 contains jasum & nusampu
awk -F, '{ print $1, $2}' trench.gmt    | psxy -R$region -J$size -Sf0.20i/0.04irt -W1.5/black -Gblack -M -N -O -K >> $psfile
awk -F, '{ print $1, $2}' all_transform.txt | psxy -R$region -J$size -W1.5 -M -N -O -K >> $psfile
awk -F, '{ print $1, $2}' all_ridge.txt     | psxy -R$region -J$size -W1.5 -M -N -O -K >> $psfile
awk -F, '{ print $1, $2}' pad_rel.txt    | psxy -J$size -R$region -Sc0.25c -W1.5/brown -O -K >> $psfile
awk -F, '{ print $1, $2}' pad_kat.txt    | psxy -J$size -R$region -Sc0.25c -W1.5/black -O -K >> $psfile

#Projection with the center
awk -F, '{ print $1, $2, "0"}' trench.gmt    | project -Q -C$pointboxlon_1/$pointboxlat_1 -A$AZ -Fxyzpqrs -L0/800 -W-200/200  >! xyzpqrs_trench1
awk -F, '{ print $1, $2, $3}' all_transform.txt | project -Q -C$pointboxlon_1/$pointboxlat_1 -A$AZ -Fxyzpqrs -L0/800 -W-200/200  >! xyzpqrs_transform1
awk -F, '{ print $1, $2, $3}' all_ridge.txt     | project -Q -C$pointboxlon_1/$pointboxlat_1 -A$AZ -Fxyzpqrs -L0/800 -W-200/200  >! xyzpqrs_ridge1
awk -F, '{ print $1, $2, $3}' pad_rel.txt  | project -Q -C$pointboxlon_1/$pointboxlat_1 -A$AZ -Fxyzpqrs -L0/800 -W-200/200  >! xyzpqrs.tabelrelokasi1
awk -F, '{ print $1, $2, $3}' pad_kat.txt  | project -Q -C$pointboxlon_1/$pointboxlat_1 -A$AZ -Fxyzpqrs -L0/800 -W-200/200  >! xyzpqrs.tabelkatalog1
awk '{print $1, $2}' xyzpqrs.tabelrelokasi1 | psxy -J$size -R$region -Sc0.25c -W1.5/darkblue -O -K >> $psfile
awk '{print $1, $2}' xyzpqrs.tabelkatalog1 | psxy -J$size -R$region -Sc0.25c -W1.5/red -O -K >> $psfile

#Draw lines
awk '{print $1, $2, 12, 0, 1, "LT", $3 }' S1.d | pstext -J$size -R$region -Gblack -Swhite -P -O -K >> $psfile
awk '{print $1, $2, 12, 0, 1, "LT", $3 }' S11.d | pstext -J$size -R$region -Gblack -Swhite -P -O -K >> $psfile
awk '{print $1, $2}' garis_S11.d | psxy -J$size -R$region -W1p,- -P -O -K >> $psfile

gmtset ANOT_FONT_SIZE 7
psbasemap -R$region -J$size -Tf97/3/0.4i/2 --HEADER_FONT_SIZE=6p -Lf97/1.6/25/200 --LABEL_FONT_SIZE=8p -O -K >> $psfile
pscoast  -R90/140/-12/8 -JM5.5c -Ggreen -Swhite -B20f10/10f5nwSE -Dl -W1 -P -O -K -X2.6 -Y5.1 >> $psfile
psxy -R -JM -W2 -P -O -K <<EOF>> $psfile
96    -6
96     1
103    1
103   -6
96    -6
EOF

#DRAW CROSSEC DIAGRAM
gmtset ANOT_FONT_SIZE 8
psbasemap -JX4.5/-2 -R0/850/-40/200 -B100/50SWne -P -O -K -X-2.5 -Y-7.8 >> $psfile
pstext -JX -R  -Gblack -P -O -K <<EOF>> $psfile
10   10   8 0 1 LT   S1
760  10   8 0 1 LT   S1'
200  -30   8 0 0.8 LT  Trench
650  -30   8 0 0.8 LT  Fault
EOF
awk '{print $4, $3}' xyzpqrs_trench1    | psxy -JX -R  -Gblack -St0.1 -P -O -K >> $psfile
awk '{print $4, $3}' xyzpqrs_transform1 | psxy -JX -R  -Gblack -St0.1 -P -O -K >> $psfile
awk '{print $4, $3}' xyzpqrs_ridge1     | psxy -JX -R  -Gblack -St0.1 -P -O -K >> $psfile
awk '{print $4, $3}' xyzpqrs.tabelrelokasi1  | psxy -JX -R  -Sc0.25c -W1.5/darkblue -P -O -K >> $psfile
awk '{print $4, $3}' xyzpqrs.tabelkatalog1  | psxy -JX -R  -Sc0.25c -W1.5/red -P -O -K >> $psfile

psscale -C$cptfile -Iintens.grd -D0.2c/1/1.2/0.1 -B2000 -O -X1.5 -Y7.3 >> $psfile

gs device $psfile
rm xyzpqrs* *.d

Tidak ada komentar:

Poskan Komentar