gfortran -o f90_fem_domain f90_fem_domain.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

./f90_fem_domain 4 out_test3Q.txt inp_dom.txt out_dom_test3Q.txt
inp1=_mat_01.txt
inp2=_mat_02.txt
inp3=_mat_03.txt
fig=fig_gmt_mesh_col3Q.eps
gmt psxy $inp1 -R$range -JX$scale -W0.5,black -G255/255/191 -P -K > $fig
gmt psxy $inp2 -R -J -W0.5,black -G255/191/191 -O -K >> $fig
gmt psxy $inp3 -R -J -W0.5,black -G191/255/223 -O >> $fig

./f90_fem_domain 3 out_test3T.txt inp_dom.txt out_dom_test3T.txt
inp1=_mat_01.txt
inp2=_mat_02.txt
inp3=_mat_03.txt
fig=fig_gmt_mesh_col3T.eps
gmt psxy $inp1 -R$range -JX$scale -W0.5,black -G255/255/191 -P -K > $fig
gmt psxy $inp2 -R -J -W0.5,black -G255/191/191 -O -K >> $fig
gmt psxy $inp3 -R -J -W0.5,black -G191/255/223 -O >> $fig