Tweak the averaging scheme: now straight avg up to N=4, then expenential with ntc=4.

This commit is contained in:
Joe Taylor 2021-02-10 11:14:46 -05:00
parent 50469f898e
commit 6395087e6f

View File

@ -234,11 +234,11 @@ subroutine q65_symspec(iwave,nmax,iz,jz,s1)
enddo
enddo
if(lnewdat) then
ntc=3
navg(iseq)=navg(iseq) + 1
ntc=min(navg(iseq),4) !Averaging time constant in sequences
u=1.0/ntc
s1a(:,:,iseq)=u*s1 + (1.0-u)*s1a(:,:,iseq)
navg(iseq)=navg(iseq) + 1
endif
endif
return
end subroutine q65_symspec