****filter the real data*****************
./filter_qin
input the name of read data to be filtered in format .asc
CAN__VHZ.asc
input the name of soure wavelet (source.gnu)
source.gnu

the result is save in:
fCAN__VHZ.asc for display using xmgr
fdCAN__VHZ.asc for input of phase velocity computation

*****calculate the station reponse***********
station_response
Entrer dt et NBT:
10 900
input the name of read data in format .asc
CAN__VHZ.asc
input the name of synthetic data to be applied by station response
UZ__CAN

the result is saved in sCAN__UZ used for phase velocity calculation

******display real data with header************* 
 awk 'BEGIN{i=0} {if($1 == "delta:") dt=$2; if(NR>79) {i++; print (i-1)*dt, $1}}' CAN__VHZ.asc >xmCAN_VHZ.asc
 xmgr xmCAN_VHZ.asc

*****calculate the phase velocity*******
./vphase_qin
name of filtered data file in format .asc
fdCAN__VHZ.asc
name of synthetic seismogram file with station response
sCAN__UZ
minimum and maximum group velocity  and number  wavetrain
2.7 5.1 1
minimal and maximal frequency in Hertz
0.00333 0.0166 
name of phase velocity file and de groupe correspondant au synthetic
vph_vgrSprem_n0


the result is save in:
 vitphas (theoretical and pratical phase velocity)
 phase   (theoretical and pratical pahse)

