gfortran -o f90_gmt_plnt f90_gmt_plnt.f90

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

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

# define of mesh colour (1:MATEL)
cat << EOT > _col_mesh.txt
 0,0,0
 255,192,203
 255,250,205
EOT

./f90_gmt_plnt out_4nod_ps.csv
cp _fig_gmt_mesh.eps fig_gmt_mesh_ps.eps
cp _fig_gmt_ps1.eps fig_gmt_ps1_ps.eps
cp _fig_gmt_ps2.eps fig_gmt_ps2_ps.eps

rm _*