gfortran -o f90_gmt_seep f90_gmt_seep.f90

mv -f 0_cpt_tunnel.cpt 0_cpt.cpt

# define of mesh colour (1:MATEL)
cat << EOT > _col_mesh.txt
 255 255 204
   0 255 255
 255   0 255
   0 255 255
 255 255 204
EOT
# integer::ind,ibc,iph,ivc,isc
# real(4)::baselength,grdspc
# real(4)::bav,bfv
# real(4)::conint,annint
# character::strx*100,strz*100
cat << EOT > _inp_base.txt
0 0 0 2 1
12.0 0.5
50.0 10.0
10.0 50.0
"Distance (m)"
"Altitude (m)"
EOT
./f90_gmt_seep out_tunnel_0m.csv
mv -f _fig_gmt_conp.eps fig_gmt_conp_tunnel0.eps
mv -f _fig_gmt_conh.eps fig_gmt_conh_tunnel0.eps
mv -f _fig_gmt_vect.eps fig_gmt_vect_tunnel0.eps

./f90_gmt_seep out_tunnel_1m.csv
mv -f _fig_gmt_conp.eps fig_gmt_conp_tunnel1.eps
mv -f _fig_gmt_conh.eps fig_gmt_conh_tunnel1.eps
mv -f _fig_gmt_vect.eps fig_gmt_vect_tunnel1.eps

./f90_gmt_seep out_tunnel_2m.csv
mv -f _fig_gmt_conp.eps fig_gmt_conp_tunnel2.eps
mv -f _fig_gmt_conh.eps fig_gmt_conh_tunnel2.eps
mv -f _fig_gmt_vect.eps fig_gmt_vect_tunnel2.eps

./f90_gmt_seep out_tunnel_3m.csv
mv -f _fig_gmt_conp.eps fig_gmt_conp_tunnel3.eps
mv -f _fig_gmt_conh.eps fig_gmt_conh_tunnel3.eps
mv -f _fig_gmt_vect.eps fig_gmt_vect_tunnel3.eps

./f90_gmt_seep out_tunnel_4m.csv
mv -f _fig_gmt_conp.eps fig_gmt_conp_tunnel4.eps
mv -f _fig_gmt_conh.eps fig_gmt_conh_tunnel4.eps
mv -f _fig_gmt_vect.eps fig_gmt_vect_tunnel4.eps

mv -f _fig_gmt_mesh.eps fig_gmt_mesh_tunnel.eps

rm _*