Use exponential message averaging (ntc=3) rather than simple sums in Q65.

This commit is contained in:
Joe Taylor
2021-02-08 11:10:02 -05:00
parent 3ccf949691
commit a6779ecda3
2 changed files with 4 additions and 1 deletions
+3 -1
View File
@@ -234,7 +234,9 @@ subroutine q65_symspec(iwave,nmax,iz,jz,s1)
enddo
enddo
if(lnewdat) then
s1a(:,:,iseq)=s1a(:,:,iseq) + s1
ntc=3
u=1.0/ntc
s1a(:,:,iseq)=u*s1 + (1.0-u)*s1a(:,:,iseq)
navg(iseq)=navg(iseq) + 1
endif