set inp=inp_test1Q.txt f90_FEM_MESH4A %inp% 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 out_work_A.txt 3 _inp_gmt_A.txt f90_GMT_MESH out_work_B.txt 4 _inp_gmt_B.txt f90_GMT_MESH out_work_C.txt 4 _inp_gmt_C.txt copy out_work_C.txt out_test1Q.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 inp=_inp_gmt_A.txt set fig=fig_gmt_mesh_A1Q.eps psxy %inp% -R%range% -JX%paper% -W0.5/black -G255/255/191 -m -P -K > %fig% gawk "10<=NR{print $0}" inp_test1Q.txt | psxy -R -J -SC0.2 -G0 -N -O >> %fig% set inp=_inp_gmt_B.txt set fig=fig_gmt_mesh_B1Q.eps psxy %inp% -R%range% -JX%paper% -W0.5/black -G255/255/191 -m -P -K > %fig% gawk "10<=NR{print $0}" inp_test1Q.txt | psxy -R -J -SC0.2 -G0 -N -O >> %fig% set inp=_inp_gmt_C.txt set fig=fig_gmt_mesh_C1Q.eps psxy %inp% -R%range% -JX%paper% -W0.5/black -G255/255/191 -m -P -K > %fig% gawk "10<=NR{print $0}" inp_test1Q.txt | psxy -R -J -SC0.2 -G0 -N -O >> %fig% set inp=_inp_gmt_C.txt set fig=fig_gmt_mesh_Cno1Q.eps psxy %inp% -R%range% -JX%paper% -W0.5/black -G255/255/191 -m -P -K > %fig% pstext _inp_node.txt -R -J -Gblue -N -O >> %fig% set inp=_inp_gmt_C.txt set fig=fig_gmt_mesh_Cne1Q.eps psxy %inp% -R%range% -JX%paper% -W0.5/black -G255/255/191 -m -P -K > %fig% pstext _inp_nelt.txt -R -J -Gred -N -O >> %fig%