mirror of
https://github.com/saitohirga/WSJT-X.git
synced 2024-11-05 08:51:19 -05:00
32 lines
625 B
FortranFixed
32 lines
625 B
FortranFixed
|
subroutine limit(x,jz)
|
||
|
|
|
||
|
|
real x(jz)
|
||
|
|
logical noping
|
||
|
|
common/limcom/ nslim2
|
||
|
|
|
||
|
|
noping=.false.
|
||
|
|
xlim=1.e30
|
||
|
|
if(nslim2.eq.1) xlim=3.0
|
||
|
|
if(nslim2.ge.2) xlim=1.0
|
||
|
|
if(nslim2.ge.3) noping=.true.
|
||
|
|
|
||
|
|
sq=0.
|
||
|
|
do i=1,jz
|
||
|
|
sq=sq+x(i)*x(i)
|
||
|
|
enddo
|
||
|
|
rms=sqrt(sq/jz)
|
||
|
|
rms0=14.5
|
||
|
|
x1=xlim*rms0
|
||
|
|
fac=1.0/xlim
|
||
|
|
if(fac.lt.1.0) fac=1.0
|
||
|
|
if(noping .and. rms.gt.20.0) fac=0.01 !Crude attempt at ping excision
|
||
|
|
|
||
|
|
do i=1,jz
|
||
|
|
if(x(i).lt.-x1) x(i)=-x1
|
||
|
|
if(x(i).gt.x1) x(i)=x1
|
||
|
|
x(i)=fac*x(i)
|
||
|
|
enddo
|
||
|
|
|
||
|
|
return
|
||
|
|
end
|