2012-11-21 12:42:53 -05:00
|
|
|
subroutine jt9a
|
|
|
|
|
|
|
|
|
|
! NB: this interface block is required by g95, but must be omitted
|
|
|
|
|
! for gfortran. (????)
|
|
|
|
|
|
|
|
|
|
interface
|
|
|
|
|
function address_jt9()
|
2013-03-23 11:35:32 -04:00
|
|
|
integer*1, pointer :: address_jt9
|
2012-11-21 12:42:53 -05:00
|
|
|
end function address_jt9
|
|
|
|
|
end interface
|
|
|
|
|
|
|
|
|
|
integer*1 attach_jt9,lock_jt9,unlock_jt9
|
|
|
|
|
integer size_jt9
|
2013-03-23 11:35:32 -04:00
|
|
|
integer*1, pointer :: p_jt9
|
2012-11-21 12:42:53 -05:00
|
|
|
character*80 cwd
|
|
|
|
|
logical fileExists
|
|
|
|
|
common/tracer/limtrace,lu
|
|
|
|
|
|
|
|
|
|
call getcwd(cwd)
|
2012-11-24 10:18:17 -05:00
|
|
|
open(12,file='timer.out',status='unknown')
|
2012-11-21 12:42:53 -05:00
|
|
|
|
|
|
|
|
limtrace=0
|
|
|
|
|
lu=12
|
|
|
|
|
i1=attach_jt9()
|
|
|
|
|
|
|
|
|
|
10 inquire(file=trim(cwd)//'/.lock',exist=fileExists)
|
|
|
|
|
if(fileExists) then
|
|
|
|
|
call sleep_msec(100)
|
|
|
|
|
go to 10
|
|
|
|
|
endif
|
|
|
|
|
|
|
|
|
|
inquire(file=trim(cwd)//'/.quit',exist=fileExists)
|
|
|
|
|
if(fileExists) then
|
|
|
|
|
! call ftnquit
|
|
|
|
|
i=detach_jt9()
|
|
|
|
|
go to 999
|
|
|
|
|
endif
|
|
|
|
|
|
|
|
|
|
nbytes=size_jt9()
|
|
|
|
|
if(nbytes.le.0) then
|
|
|
|
|
print*,'jt9a: Shared memory mem_jt9 does not exist.'
|
|
|
|
|
print*,"Must start 'jt9 -s' from within WSJT-X."
|
|
|
|
|
go to 999
|
|
|
|
|
endif
|
|
|
|
|
p_jt9=>address_jt9()
|
|
|
|
|
call jt9b(p_jt9,nbytes)
|
|
|
|
|
|
|
|
|
|
! write(*,1010)
|
|
|
|
|
!1010 format('<jt9aFinished>')
|
|
|
|
|
! flush(6)
|
|
|
|
|
|
|
|
|
|
100 inquire(file=trim(cwd)//'/.lock',exist=fileExists)
|
|
|
|
|
if(fileExists) go to 10
|
2013-03-26 10:00:51 -04:00
|
|
|
call sleep_msec(100)
|
2012-11-21 12:42:53 -05:00
|
|
|
go to 100
|
|
|
|
|
|
|
|
|
|
999 return
|
|
|
|
|
end subroutine jt9a
|
|
|
|
|
|
|
|
|
|
subroutine jt9b(jt9com,nbytes)
|
|
|
|
|
integer*1 jt9com(0:nbytes-1)
|
|
|
|
|
kss=0
|
|
|
|
|
ksavg=kss + 4*184*22000
|
|
|
|
|
kc0=ksavg + 4*22000
|
|
|
|
|
kid2=kc0 + 2*4*1800*1500
|
|
|
|
|
knutc=kid2 + 2*1800*12000
|
|
|
|
|
call jt9c(jt9com(kss),jt9com(ksavg),jt9com(kc0),jt9com(kid2),jt9com(knutc))
|
|
|
|
|
return
|
|
|
|
|
end subroutine jt9b
|
|
|
|
|
|
|
|
|
|
subroutine jt9c(ss,savg,c0,id2,nparams0)
|
|
|
|
|
parameter (NSMAX=22000)
|
|
|
|
|
integer*1 detach_jt9
|
|
|
|
|
real*4 ss(184*NSMAX),savg(NSMAX)
|
|
|
|
|
complex c0(1800*1500)
|
|
|
|
|
integer*2 id2(1800*12000)
|
|
|
|
|
|
|
|
|
|
integer nparams0(21),nparams(21)
|
|
|
|
|
character*20 datetime
|
|
|
|
|
common/npar/nutc,ndiskdat,ntrperiod,nfqso,newdat,npts8,nfa,nfb,ntol, &
|
|
|
|
|
kin,nzhsym,nsave,nagain,ndepth,nrxlog,nfsample,datetime
|
|
|
|
|
equivalence (nparams,nutc)
|
|
|
|
|
|
|
|
|
|
nparams=nparams0 !Copy parameters into common/npar/
|
|
|
|
|
npatience=1
|
|
|
|
|
|
2012-11-21 16:08:17 -05:00
|
|
|
call flush(6)
|
2013-04-16 12:25:14 -04:00
|
|
|
if(sum(nparams).ne.0) call decoder(ss,c0,0)
|
2012-11-21 12:42:53 -05:00
|
|
|
|
|
|
|
|
return
|
|
|
|
|
end subroutine jt9c
|
