gfortran -o f90_wave_phase f90_wave_phase.f90

awk 'NR<=3{print $0}4<=NR{print $1}' inp_ran1.csv > inp_test1.csv
awk 'NR<=3{print $0}4<=NR{print $1}' acc_ran1.csv > inp_test2.csv

./f90_wave_phase inp_test1.csv out_fps1.csv out_pw1.csv
./f90_wave_phase inp_test2.csv out_fps2.csv out_pw2.csv

bash a_gmt_phase.sh