Reorder the ibw loop in q65_loops.f90.

This commit is contained in:
Joe Taylor 2020-12-20 21:43:30 -05:00
parent ddf727f351
commit 545079e951

View File

@ -29,8 +29,8 @@ subroutine q65_loops(c00,npts2,nsps,mode,mode_q65,nsubmode,nFadingModel, &
idfmax=3 idfmax=3
idtmax=3 idtmax=3
ibwmin=0 ibwmin=1
ibwmax=2 ibwmax=3
maxdist=5 maxdist=5
if(iand(ndepth,3).ge.2) then if(iand(ndepth,3).ge.2) then
idfmax=5 idfmax=5
@ -66,13 +66,13 @@ subroutine q65_loops(c00,npts2,nsps,mode,mode_q65,nsubmode,nFadingModel, &
where(s3(1:LL*NN)>s3lim) s3(1:LL*NN)=s3lim where(s3(1:LL*NN)>s3lim) s3(1:LL*NN)=s3lim
endif endif
do ibw=ibwmin,ibwmax do ibw=ibwmin,ibwmax
nbw=ibw nbw=ibw/2
ndist=ndf**2 + ndt**2 + ((nbw-2))**2 if(mod(ibw,2).eq.0) nbw=-nbw
ndist=ndf**2 + ndt**2 + nbw**2
if(ndist.gt.maxdist) cycle if(ndist.gt.maxdist) cycle
b90=3.0**nbw !### Mult by ~baud/3.33 ??? ### xx=1.885*log(3.0*width)+nbw
if(nbw.eq.0) b90=baud b90=1.7**xx
if(b90.lt.0.5*width) cycle if(b90.gt.345.0) cycle
if(b90.gt.230.0) cycle
call timer('q65_intr',0) call timer('q65_intr',0)
b90ts = b90/baud b90ts = b90/baud
call q65_intrinsics_ff(s3,nsubmode,b90ts,nFadingModel,s3prob) call q65_intrinsics_ff(s3,nsubmode,b90ts,nFadingModel,s3prob)