gfortran -o f90_gmt_ini4_con f90_gmt_ini4_con.f90

cat << EOT > _cpt_ps.cpt
-200  204   0   0   -150  255   0 102
-150  255 153 102   -100  255 153 102
-100  255 204 102    -50  255 204 102
 -50  255 255 170      0  255 255 170
   0    0 255 255     25    0 255 255
  25    0 204 244     50    0 204 255
  50    0 153 255     75    0 153 255
  75    0   0 255    100    0   0 136
B 204   0   0
F   0   0 136
N 255 255 255
EOT

cat << EOT > _cpt_sf.cpt
0.5 255   0   0  1.0 255   0   0
1.0 255 255 170  1.2 255 255 170
1.2   0 255   0  1.5   0 255   0
1.5   0 255 255  2.0   0 128 255
B 128   0   0
F   0 128 255
N 255 255 255
EOT

# 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 1 1
 12.0 0.05
 5.0 5.0
 50 50
 "x-distance (m)"
 "y-distance (m)"
 0 0
EOT

# define of mesh colour (1:MATEL)
cat << EOT > _col_mesh.txt
 255,255,204
EOT

./f90_gmt_ini4_con out_div_tunnel.prn
cp _fig_gmt_mesh.eps fig_gmt_mesh.eps
cp _fig_gmt_ps1.eps fig_gmt_ps1.eps
cp _fig_gmt_ps2.eps fig_gmt_ps2.eps
cp _fig_gmt_sf.eps fig_gmt_sf.eps

#rm _*