Procedure for making heliocentric phased, trailed spectra.
(example lists in HMXB/data/LPH040/RV_measurements)

1. Correct spectra to heliocentric using IRAF dopcor:
(
PACKAGE = kpnoslit
   TASK = dopcor

input   =       @speclist_fits  List of input spectra
output  =                       List of output spectra
redshift=              -vhelio  Redshift or velocity (Km/s)
(isveloc=                  yes) Is the redshift parameter a velocity?
(add    =                   no) Add to previous dispersion correction?
(dispers=                  yes) Apply dispersion correction?
(flux   =                   no) Apply flux correction?
(factor =                   3.) Flux correction factor (power of 1+z)
(apertur=                     ) List of apertures to correct
(verbose=                   no) Print corrections performed?
(mode   =                   ql)
)

2. Use (in IRAF) wspectext to convert fits spectra to ascii. 
(Input list=speclist_fits, output list=speclist_txt, 
converts *.ms.fits to *.ms.txt)

3. Use hselect ($I,HJD) to create hjdlist

4. For use later, edit hjdlist to make script addhjd to add HJD to gridded
spectra (convert *.ms.grd to *.ms.spec):
(
awk '{print "2454840.64208157 "$1" "$2}' con000083.ms.grd > con000083.ms.spec
awk '{print "2454840.73074613 "$1" "$2}' con000092.ms.grd > con000092.ms.spec
awk '{print "2454844.64130534 "$1" "$2}' con000102.ms.grd > con000102.ms.spec
awk '{print "2454844.79540599 "$1" "$2}' con000108.ms.grd > con000108.ms.spec
awk '{print "2454848.72885679 "$1" "$2}' con000114.ms.grd > con000114.ms.spec
awk '{print "2454848.88033089 "$1" "$2}' con000120.ms.grd > con000120.ms.spec
awk '{print "2454848.8968116 "$1" "$2}' con000122.ms.grd > con000122.ms.spec
awk '{print "2454848.91221596 "$1" "$2}' con000124.ms.grd > con000124.ms.spec
awk '{print "2454849.61680566 "$1" "$2}' con000143.ms.grd > con000143.ms.spec
awk '{print "2454849.63422338 "$1" "$2}' con000146.ms.grd > con000146.ms.spec
awk '{print "2454849.65409467 "$1" "$2}' con000149.ms.grd > con000149.ms.spec
awk '{print "2454850.60650588 "$1" "$2}' con000179.ms.grd > con000179.ms.spec
awk '{print "2454850.65164184 "$1" "$2}' con000185.ms.grd > con000185.ms.spec
awk '{print "2454850.69871042 "$1" "$2}' con000191.ms.grd > con000191.ms.spec
awk '{print "2454850.74756141 "$1" "$2}' con000197.ms.grd > con000197.ms.spec
awk '{print "2454850.77692329 "$1" "$2}' con000199.ms.grd > con000199.ms.spec
awk '{print "2454850.80627356 "$1" "$2}' con000201.ms.grd > con000201.ms.spec
awk '{print "2454850.8356703 "$1" "$2}' con000203.ms.grd > con000203.ms.spec
awk '{print "2454850.86522911 "$1" "$2}' con000205.ms.grd > con000205.ms.spec
awk '{print "2454850.8947764 "$1" "$2}' con000207.ms.grd > con000207.ms.spec
awk '{print "2454850.92431229 "$1" "$2}' con000209.ms.grd > con000209.ms.spec
awk '{print "2454442.82862695 "$1" "$2}' con015056.ms.grd > con015056.ms.spec
awk '{print "2454442.99308504 "$1" "$2}' con015072.ms.grd > con015072.ms.spec
awk '{print "2454446.99981414 "$1" "$2}' con015177.ms.grd > con015177.ms.spec
awk '{print "2454336.94159354 "$1" "$2}' con09128.ms.grd > con09128.ms.spec
awk '{print "2454342.97704712 "$1" "$2}' con09329.ms.grd > con09329.ms.spec
)

[Steps 5 to 9 can all be done with the script txt2sort:
 ./mkgrid
./addhjd
./speclist2
./tog
./phsrt
]

5. grid ascii spectra to uniform grid using my interp1D (ccd #21)
(input *.ms.txt to output *.ms.grd, script mkgrid does a list of
spectra using interp1D:
interp1D con000083.ms.txt con000083.ms.grd 4300 4520 0.25
interp1D con000092.ms.txt con000092.ms.grd 4300 4520 0.25
interp1D con000102.ms.txt con000102.ms.grd 4300 4520 0.25
interp1D con000108.ms.txt con000108.ms.grd 4300 4520 0.25
interp1D con000114.ms.txt con000114.ms.grd 4300 4520 0.25
interp1D con000120.ms.txt con000120.ms.grd 4300 4520 0.25
interp1D con000122.ms.txt con000122.ms.grd 4300 4520 0.25
interp1D con000124.ms.txt con000124.ms.grd 4300 4520 0.25
interp1D con000143.ms.txt con000143.ms.grd 4300 4520 0.25
interp1D con000146.ms.txt con000146.ms.grd 4300 4520 0.25
interp1D con000149.ms.txt con000149.ms.grd 4300 4520 0.25
interp1D con000179.ms.txt con000179.ms.grd 4300 4520 0.25
interp1D con000185.ms.txt con000185.ms.grd 4300 4520 0.25
interp1D con000191.ms.txt con000191.ms.grd 4300 4520 0.25
interp1D con000197.ms.txt con000197.ms.grd 4300 4520 0.25
interp1D con000199.ms.txt con000199.ms.grd 4300 4520 0.25
interp1D con000201.ms.txt con000201.ms.grd 4300 4520 0.25
interp1D con000203.ms.txt con000203.ms.grd 4300 4520 0.25
interp1D con000205.ms.txt con000205.ms.grd 4300 4520 0.25
interp1D con000207.ms.txt con000207.ms.grd 4300 4520 0.25
interp1D con000209.ms.txt con000209.ms.grd 4300 4520 0.25
interp1D con015056.ms.txt con015056.ms.grd 4300 4520 0.25
interp1D con015072.ms.txt con015072.ms.grd 4300 4520 0.25
interp1D con015177.ms.txt con015177.ms.grd 4300 4520 0.25
interp1D con09128.ms.txt con09128.ms.grd 4300 4520 0.25
interp1D con09329.ms.txt con09329.ms.grd 4300 4520 0.25
) 

6. Run the script, made in step 4, addhjd to add HJD to the ascii spectra.
input = *.ms.txt, output = *.ms.spec

7. Use script speclist2 to sub HJD with phase:
(
awk '{print$2" "($1/0.67 - int($1/0.67))" "$3}' con000083.ms.spec > con000083.ms.phase
awk '{print$2" "($1/0.67 - int($1/0.67))" "$3}' con000092.ms.spec > con000092.ms.phase
awk '{print$2" "($1/0.67 - int($1/0.67))" "$3}' con000102.ms.spec > con000102.ms.phase
awk '{print$2" "($1/0.67 - int($1/0.67))" "$3}' con000108.ms.spec > con000108.ms.phase
awk '{print$2" "($1/0.67 - int($1/0.67))" "$3}' con000114.ms.spec > con000114.ms.phase
awk '{print$2" "($1/0.67 - int($1/0.67))" "$3}' con000120.ms.spec > con000120.ms.phase
awk '{print$2" "($1/0.67 - int($1/0.67))" "$3}' con000122.ms.spec > con000122.ms.phase
awk '{print$2" "($1/0.67 - int($1/0.67))" "$3}' con000124.ms.spec > con000124.ms.phase
awk '{print$2" "($1/0.67 - int($1/0.67))" "$3}' con000143.ms.spec > con000143.ms.phase
awk '{print$2" "($1/0.67 - int($1/0.67))" "$3}' con000146.ms.spec > con000146.ms.phase
awk '{print$2" "($1/0.67 - int($1/0.67))" "$3}' con000149.ms.spec > con000149.ms.phase
awk '{print$2" "($1/0.67 - int($1/0.67))" "$3}' con000179.ms.spec > con000179.ms.phase
awk '{print$2" "($1/0.67 - int($1/0.67))" "$3}' con000185.ms.spec > con000185.ms.phase
awk '{print$2" "($1/0.67 - int($1/0.67))" "$3}' con000191.ms.spec > con000191.ms.phase
awk '{print$2" "($1/0.67 - int($1/0.67))" "$3}' con000197.ms.spec > con000197.ms.phase
awk '{print$2" "($1/0.67 - int($1/0.67))" "$3}' con000199.ms.spec > con000199.ms.phase
awk '{print$2" "($1/0.67 - int($1/0.67))" "$3}' con000201.ms.spec > con000201.ms.phase
awk '{print$2" "($1/0.67 - int($1/0.67))" "$3}' con000203.ms.spec > con000203.ms.phase
awk '{print$2" "($1/0.67 - int($1/0.67))" "$3}' con000205.ms.spec > con000205.ms.phase
awk '{print$2" "($1/0.67 - int($1/0.67))" "$3}' con000207.ms.spec > con000207.ms.phase
awk '{print$2" "($1/0.67 - int($1/0.67))" "$3}' con000209.ms.spec > con000209.ms.phase
awk '{print$2" "($1/0.67 - int($1/0.67))" "$3}' con015056.ms.spec > con015056.ms.phase
awk '{print$2" "($1/0.67 - int($1/0.67))" "$3}' con015072.ms.spec > con015072.ms.phase
awk '{print$2" "($1/0.67 - int($1/0.67))" "$3}' con015177.ms.spec > con015177.ms.phase
awk '{print$2" "($1/0.67 - int($1/0.67))" "$3}' con09128.ms.spec > con09128.ms.phase
awk '{print$2" "($1/0.67 - int($1/0.67))" "$3}' con09329.ms.spec > con09329.ms.phase
)

8. cat phase spectra together to spectra3dtemp.txt 
with script tog:
(
cat con000083.ms.phase \
con000092.ms.phase \
con000102.ms.phase \
con000108.ms.phase \
con000114.ms.phase \
con000120.ms.phase \
con000122.ms.phase \
con000124.ms.phase \
con000143.ms.phase \
con000146.ms.phase \
con000149.ms.phase \
con000179.ms.phase \
con000185.ms.phase \
con000191.ms.phase \
con000197.ms.phase \
con000199.ms.phase \
con000201.ms.phase \
con000203.ms.phase \
con000205.ms.phase \
con000207.ms.phase \
con000209.ms.phase \
con015056.ms.phase \
con015072.ms.phase \
con015177.ms.phase \
con09128.ms.phase \
con09329.ms.phase > spectra3dtemp.txt
)

9. sort phased spectra with phsort:
(
awk '{print$2" "$1" "$3}' spectra3dtemp.txt > temp1
sort temp1 > temp2
awk '{print$2" "$1" "$3}' temp2 > spectra3d.txt
rm temp1
rm temp2
)

10. plot phased spectra with plotphase.gnu:
(
reset
set term x11
#set term pbm color
#set out "phase.pbm"
set xlabel "phase"
set ylabel "Heliocentric RV (km/s)"
set multiplot layout 3,3 title "Heliocentric radial velocities folded at 0.67 day period"
plot "LPH040_table_rv.txt" u ($3/0.67 - int($3/0.67)):14:15 w err title "Hgam blue"
plot "LPH040_table_rv.txt" u ($3/0.67 - int($3/0.67)):16:17 w err title "Hgam red"
plot "LPH040_table_rv.txt" u ($3/0.67 - int($3/0.67)):18:19 w err title "Hgam centre"
plot "LPH040_table_rv.txt" u ($3/0.67 - int($3/0.67)):($16-$18):(sqrt($17*$17+$19*$19)) w err title "Hgam red-blue"
plot "LPH040_table_rv.txt" u ($3/0.67 - int($3/0.67)):20:21 w err title "HeI 4388"
plot "LPH040_table_rv.txt" u ($3/0.67 - int($3/0.67)):22:23 w err title "HeI 4477"
plot "LPH040_table_rv.txt" u ($3/0.67 - int($3/0.67)):24:25 w err title "Mg 4481"
unset multiplot
)