Compute a mean radial velocity given an input file listing observed wavelength, restframe wavelength, and feature name.
% redshift.sh --help
usage: redshift.sh dat.1 [-v] [-h]
arg1 - input file listing Wave_obs, Wave_rest, Feature_Name
Additional options:
-v = print verbose comments and run in debug mode
-h = just show usage message
Here is a simple example of a run:
% cat dat.1
4486.194 4340.470 H_gamma
5024.542 4861.330 H_beta_L
% redshift.sh dat.1
% cat redshift.out
Initial velocity values:
Wave_obs Wave_rest Feautre Vobs(km/s)
4486.194 4340.470 H_gamma 10072.0
5024.542 4861.330 H_beta_L 10072.1
Mean Velocity = 10072.012 -+ 0.041