gfortran -o f90_fem_mesh3 f90_fem_mesh3.f90
gfortran -o f90_gmt_mesh f90_gmt_mesh.f90

range=0/30/0/10
scale=12/4

gmt set FONT_ANNOT_PRIMARY 10
gmt set MAP_ANNOT_OFFSET_PRIMARY 0.3c
gmt set FONT_LABEL 10p
gmt set MAP_LABEL_OFFSET 0.1c
gmt set MAP_TICK_LENGTH_PRIMARY -0.2c

for i in 1 2 3
do
    ./f90_fem_mesh3 inp_test${i}T.txt out_test${i}T.txt
    ./f90_gmt_mesh 3 out_test${i}T.txt _inp_gmt_${i}T.txt
    inp=_inp_gmt_${i}T.txt
    fig1=fig_gmt_mesh_${i}T.eps
    fig2=fig_gmt_mesh_${i}Tno.eps
    fig3=fig_gmt_mesh_${i}Tne.eps
    gmt psxy $inp -R$range -JX$scale -W0.5,black -G255/255/191 -P > $fig1
    gmt psxy $inp -R$range -JX$scale -W0.5,black -G255/255/191 -P -K > $fig2
    gmt pstext _inp_node.txt -R -J -F+f+a+j -N -O >> $fig2
    gmt psxy $inp -R$range -JX$scale -W0.5,black -G255/255/191 -P -K > $fig3
    gmt pstext _inp_nelt.txt -R -J -F+f+a+j -N -O >> $fig3
done
exit 0