gfortran -o f90_GMTplnt.exe f90_GMTplnt.f90

makecpt -Cno_green -T-10/10/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 mesh colour (1:MATEL)
echo 255,250,205 > _col_mesh.txt
echo 255,192,203 >> _col_mesh.txt

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 1 1 > _inp_base.txt
echo 12.0 0.5 >> _inp_base.txt
echo 30.0 10.0 >> _inp_base.txt
echo  1.0 5.0 >> _inp_base.txt
echo "x-distance (mm)" >> _inp_base.txt
echo "y-distance (mm)" >> _inp_base.txt
echo  0 >> _inp_base.txt

f90_GMTplnt out4rand19h_m.csv
copy fig_gmt_mesh.eps fig_gmt_mesh_hole.eps
copy fig_gmt_ps1.eps fig_gmt_ps1_hole.eps
copy fig_gmt_ps2.eps fig_gmt_ps2_hole.eps

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