gfortran -o f90_GMTplnt.exe f90_GMTplnt.f90

makecpt -Cno_green -T-1/1/0.1 -I > 1_cpt.cpt
gawk "{a[NR]=$0}END{for(i=1;i<=NR-1;i++){print a[i]};gsub(/(128)/,\"255\",a[NR]);print a[NR]}" 1_cpt.cpt > 0_cpt.cpt

rem define of drawing data
rem integer::ind,ibc,ifc,isc
rem real(4)::baselength,grdspc
rem real(4)::bav,bfv
rem real(4)::conint,annint
rem character::strx*100,stry*100
rem integer::ncut,numcut(10)
echo 0 1 0 1 > _inp_base.txt
echo 12.0 50 >> _inp_base.txt
echo 5000.0 5000.0 >> _inp_base.txt
echo  0.5 0.5 >> _inp_base.txt
echo "x-distance (mm)" >> _inp_base.txt
echo "y-distance (mm)" >> _inp_base.txt
echo  1 1 >> _inp_base.txt

rem define of mesh colour (1:MATEL)
echo 0,0,0 > _col_mesh.txt
echo 255,192,203 >> _col_mesh.txt
echo 255,250,205 >> _col_mesh.txt

f90_GMTplnt out4nodPS0_L.csv
copy fig_gmt_mesh.eps fig_gmt_mesh_ps.eps
copy fig_gmt_ps1.eps fig_gmt_ps1_ps.eps
copy fig_gmt_ps2.eps fig_gmt_ps2_ps.eps

del fig_gmt_mesh.eps
del fig_gmt_ps1.eps
del fig_gmt_ps2.eps
del _*