gfortran -o f90_FEM_DOMAIN.exe f90_FEM_DOMAIN.f90

rem **************************************************************
rem program name           : f90_FEM_DOMAIN
rem node number per element: 4 
rem element and node file  : out_test3Q.txt
rem domain definition file : inp_dom.txt
rem output file            : out_dom_test3Q.txtren 
rem **************************************************************
f90_FEM_DOMAIN 4 out_test3Q.txt inp_dom.txt out_dom_test3Q.txt
set range=-0/30/0/10
set paper=12/4
set bxtext=a10f5:"x-direction"
set bytext=a10f5:"y-direction"
gmtset ANOT_FONT_SIZE 10
gmtset ANOT_OFFSET 0.2c
gmtset LABEL_FONT_SIZE 10
gmtset LABEL_OFFSET 0.1c
set inp1=_mat_01.txt
set inp2=_mat_02.txt
set inp3=_mat_03.txt
set fig=fig_gmt_mesh_col3Q.eps
psxy %inp1% -R%range% -JX%paper% -W0.5/black -G255/255/191 -m -P -K > %fig%
psxy %inp2% -R -J -W0.5/black -G255/191/191 -m -O -K >> %fig%
psxy %inp3% -R -J -W0.5/black -G191/255/223 -m -O >> %fig%


f90_FEM_DOMAIN 3 out_test3T.txt inp_dom.txt out_dom_test3T.txt
set range=-0/30/0/10
set paper=12/4
set bxtext=a10f5:"x-direction"
set bytext=a10f5:"y-direction"
gmtset ANOT_FONT_SIZE 10
gmtset ANOT_OFFSET 0.2c
gmtset LABEL_FONT_SIZE 10
gmtset LABEL_OFFSET 0.1c
set inp1=_mat_01.txt
set inp2=_mat_02.txt
set inp3=_mat_03.txt
set fig=fig_gmt_mesh_col3T.eps
psxy %inp1% -R%range% -JX%paper% -W0.5/black -G255/255/191 -m -P -K > %fig%
psxy %inp2% -R -J -W0.5/black -G255/191/191 -m -O -K >> %fig%
psxy %inp3% -R -J -W0.5/black -G191/255/223 -m -O >> %fig%