Dear GMT Forum,
Thank you for the time reading the post ,my GMT is 5.4.1, I have the following script (it plots some geological faults):
region="-R-69.76/-62.81/-23.17/-13.0"
regionea="-R-68.8451/-67.4718/-16.889/-16.1095"
grdfile="MT_bedrock.grd"
proy="-JM5i"
proyea="-JM2i"
flags="-Xc -Yc -B2.5 -K -P -V"
flagsea="-X6.7c -Y15.85c -Bf -K -P -V"
pts="-Wthinnest -V"
color="GMT_relief.cpt"
archi="10_Fig2a_Bolfaults.ps"
psbasemap $region $proy $flags > $archi
grdgradient $grdfile -A180 -Ne0.7 -Gbol_MT.int
grdimage $grdfile $region -Ibol_MT.int $proy $flags -C$color -O -K -n+c>> $archi
pscoast $region $proy -B2.5 -B+t"Faults" -N1 -Df -W2/1 -Scornflowerblue -O -V -K >> $archi
awk '{if (NR > 1) print $2,$1}' Ftunari.txt | psxy -JM -R -Sf0.25i/8.5p+l+f -W1.3,red -O -K >> $archi
#another map inside the "big map"
#it is a small set of faults that are so near
grdimage $grdfile $regionea -Ibol_MT.int $proyea $flagsea -C$color -O -K -n+c>> $archi
pscoast $regionea $proyea -Bf -N1 -Df -W2/1 -X6.7c -Y15.85c -O -V -K >> $archi
#this is the attemp, but no draw
awk '{if (NR > 1) print $2,$1}' Fllojeta.txt | psxy -JM -R -Sf0.9/3.5p+l+f -W1.3,red -O -K -P >> $archi
psxy -R -J /dev/null -O >> $archi
gv $archi &
The map is plotted, no problem at all, however I would like to add a small map (a zoom of a region) and populate it with psxy data
grdimage $grdfile $regionea -Ibol_MT.int $proyea $flagsea -C$color -O -K -n+c>> $archi
pscoast $regionea $proyea -Bf -N1 -Df -W2/1 -X6.7c -Y15.85c -O -V -K >> $archi
#this is the attemp, but no draw
awk '{if (NR > 1) print $2,$1}' Fllojeta.txt | psxy -JM -R -Sf0.9/3.5p+l+f -W1.3,red -O -K -P >> $archi
Maybe the scale is not correct `-Sf0.9/3.5p+l+f` , does anyone make a plot like the one I present?, thank you for the time.
