2022-12-12 13:09:07 -05:00
|
|
|
subroutine q65c(itimer)
|
2022-12-09 19:19:59 -05:00
|
|
|
|
|
|
|
|
use timer_module, only: timer
|
2022-12-10 09:24:24 -05:00
|
|
|
use timer_impl, only: fini_timer !, limtrace
|
2022-12-09 19:19:59 -05:00
|
|
|
use, intrinsic :: iso_c_binding, only: C_NULL_CHAR
|
|
|
|
|
use FFTW3
|
|
|
|
|
use q65
|
|
|
|
|
use q65_decode
|
|
|
|
|
|
|
|
|
|
parameter (NFFT=32768)
|
|
|
|
|
include 'njunk.f90'
|
|
|
|
|
real*8 fcenter
|
|
|
|
|
integer nparams0(NJUNK+3),nparams(NJUNK+3)
|
2022-12-13 12:58:52 -05:00
|
|
|
logical first
|
2024-01-01 15:48:24 -05:00
|
|
|
logical*1 bAlso30
|
2022-12-09 19:19:59 -05:00
|
|
|
character*12 mycall,hiscall
|
|
|
|
|
character*6 mygrid,hisgrid
|
|
|
|
|
character*20 datetime
|
|
|
|
|
|
2023-12-07 11:43:25 -05:00
|
|
|
common/datcom2/dd(2,5760000),ss(400,NFFT),savg(NFFT),nparams0
|
2022-12-10 19:21:19 -05:00
|
|
|
|
|
|
|
|
!### REMEMBER that /npar/ is not updated until nparams=nparams0 is executed. ###
|
2023-05-13 13:43:04 -04:00
|
|
|
common/npar/fcenter,nutc,fselected,mousedf,mousefqso,nagain, &
|
2024-01-01 09:46:09 -05:00
|
|
|
ndepth,ndiskdat,ntx60,newdat,nn1,nn2,nfcal,nfshift, &
|
|
|
|
|
ntx30a,ntx30b,ntol,nxant,nCFOM,nfsample,nxpol,nmode, &
|
2024-01-01 13:26:26 -05:00
|
|
|
ndop00,nsave,nn3,nn4,nhsym,mycall,mygrid,hiscall,hisgrid, &
|
2023-12-14 13:37:28 -05:00
|
|
|
datetime,junk1,junk2,bAlso30
|
2022-12-09 19:19:59 -05:00
|
|
|
equivalence (nparams,fcenter)
|
|
|
|
|
data first/.true./
|
2022-12-10 19:21:19 -05:00
|
|
|
save first
|
2022-12-09 19:19:59 -05:00
|
|
|
|
2022-12-10 19:21:19 -05:00
|
|
|
nparams=nparams0 !Copy parameters into common/npar/
|
2023-12-08 12:06:20 -05:00
|
|
|
datetime(12:)='00 '
|
2023-12-01 08:46:40 -05:00
|
|
|
|
2022-12-10 09:24:24 -05:00
|
|
|
if(itimer.ne.0) then
|
|
|
|
|
call timer('decode0 ',101)
|
|
|
|
|
call fini_timer
|
|
|
|
|
return
|
|
|
|
|
endif
|
2022-12-10 19:21:19 -05:00
|
|
|
|
2022-12-09 19:19:59 -05:00
|
|
|
npatience=1
|
2023-01-16 18:51:57 -05:00
|
|
|
newdat=1 !Always on ??
|
2022-12-17 08:28:48 -05:00
|
|
|
|
2024-01-01 15:48:24 -05:00
|
|
|
call chkstat(dd,nhsym,dbdiff)
|
|
|
|
|
! if(nhsym.eq.390) write(*,3001) nutc,nhsym,dbdiff
|
|
|
|
|
!3001 format(i4.4,i6,f7.1)
|
|
|
|
|
|
|
|
|
|
if(dbdiff.gt.5.0) ntx30b=30
|
|
|
|
|
if(dbdiff.lt.-5.0) ntx30a=30
|
|
|
|
|
|
2024-01-01 13:56:39 -05:00
|
|
|
if(ntx30a.gt.5) then
|
|
|
|
|
dd(1:2,1:30*96000)=0.
|
|
|
|
|
ss(1:200,1:NFFT)=0.
|
2024-01-01 15:48:24 -05:00
|
|
|
do i=1,NFFT
|
|
|
|
|
savg(i)=sum(ss(201:400,i))
|
|
|
|
|
enddo
|
2024-01-01 13:56:39 -05:00
|
|
|
endif
|
|
|
|
|
if(ntx30b.gt.5) then
|
|
|
|
|
dd(1:2,30*96000+1:60*96000)=0.
|
|
|
|
|
ss(201:400,1:NFFT)=0.
|
2024-01-01 15:48:24 -05:00
|
|
|
do i=1,NFFT
|
|
|
|
|
savg(i)=sum(ss(1:200,i))
|
|
|
|
|
enddo
|
2024-01-01 13:26:26 -05:00
|
|
|
endif
|
2023-12-16 11:48:06 -05:00
|
|
|
|
2022-12-10 09:24:24 -05:00
|
|
|
call timer('decode0 ',0)
|
2022-12-12 15:23:51 -05:00
|
|
|
call decode0(dd,ss,savg)
|
2022-12-10 09:24:24 -05:00
|
|
|
call timer('decode0 ',1)
|
2022-12-09 19:19:59 -05:00
|
|
|
|
|
|
|
|
return
|
2022-12-12 13:09:07 -05:00
|
|
|
end subroutine q65c
|
