gfortran -o f90_fem_ini4 f90_fem_ini4.f90
./f90_fem_ini4 0.0 150.0 0.0 inp_div_tunnel.csv out_div_tunnel.prn