Add a comment to clarify units for 'drift'.

This commit is contained in:
Joe Taylor 2021-05-29 13:18:11 -04:00
parent 40f1af5af2
commit ccf5e0daa7

View File

@ -52,6 +52,7 @@ subroutine q65_loops(c00,npts2,nsps2,nsubmode,ndepth,jpk0, &
if(mod(idf,2).eq.0) ndf=-ndf
a=0.
a(1)=-(f0+0.5*baud*ndf)
! Variable 'drift' is frequency increase over full TxT. Therefore we want:
a(2)=-0.5*drift
call twkfreq(c00,c0,npts2,6000.0,a)
do idt=1,idtmax