gfortran -o f90_fem_mesh4a f90_fem_mesh4a.f90
gfortran -o f90_fem_mesh4b f90_fem_mesh4b.f90
gfortran -o f90_fem_mesh4c f90_fem_mesh4c.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_mesh4a inp_test${i}Q.txt out_work_A.txt
    ./f90_fem_mesh4b out_work_A.txt out_work_B.txt
    ./f90_fem_mesh4c out_work_B.txt out_work_C.txt
    ./f90_gmt_mesh 3 out_work_A.txt _inp_gmt_A.txt
    ./f90_gmt_mesh 4 out_work_B.txt _inp_gmt_B.txt
    ./f90_gmt_mesh 4 out_work_C.txt _inp_gmt_C.txt
    cp out_work_C.txt out_test${i}Q.txt

    inp=_inp_gmt_A.txt
    fig1=fig_gmt_mesh_A${i}Q.eps
    gmt psxy $inp -R$range -JX$scale -W0.5,black -G255/255/191 -P -K > $fig1
    awk '10<=NR{print $0}' inp_test${i}Q.txt | gmt psxy -R -J -SC0.2 -G0 -N -O >> $fig1

    inp=_inp_gmt_B.txt
    fig2=fig_gmt_mesh_B${i}Q.eps
    gmt psxy $inp -R$range -JX$scale -W0.5,black -G255/255/191 -P -K > $fig2
    awk '10<=NR{print $0}' inp_test${i}Q.txt | gmt psxy -R -J -SC0.2 -G0 -N -O >> $fig2

    inp=_inp_gmt_C.txt
    fig3=fig_gmt_mesh_C${i}Q.eps
    gmt psxy $inp -R$range -JX$scale -W0.5,black -G255/255/191 -P -K > $fig3
    awk '10<=NR{print $0}' inp_test${i}Q.txt | gmt psxy -R -J -SC0.2 -G0 -N -O >> $fig3

    inp=_inp_gmt_C.txt
    fig4=fig_gmt_mesh_Cno${i}Q.eps
    gmt psxy $inp -R$range -JX$scale -W0.5,black -G255/255/191 -P -K > $fig4
    gmt pstext _inp_node.txt -R -J  -F+f+a+j -N -O >> $fig4

    inp=_inp_gmt_C.txt
    fig5=fig_gmt_mesh_Cne${i}Q.eps
    gmt psxy $inp -R$range -JX$scale -W0.5,black -G255/255/191 -P -K > $fig5
    gmt pstext _inp_nelt.txt -R -J -F+f+a+j -N -O >> $fig5
done
exit 0