| 
									
										
										
										
											2016-05-17 17:19:27 +00:00
										 |  |  | subroutine extract(s3,nadd,mode65,ntrials,naggressive,ndepth,nflip,     &
 | 
					
						
							|  |  |  |      mycall_12,hiscall_12,hisgrid,nexp_decode,ncount,nhist,decoded,     &
 | 
					
						
							|  |  |  |      ltext,nft,qual)
 | 
					
						
							| 
									
										
										
										
											2014-12-03 00:06:54 +00:00
										 |  |  | 
 | 
					
						
							|  |  |  | ! Input:
 | 
					
						
							|  |  |  | !   s3       64-point spectra for each of 63 data symbols
 | 
					
						
							|  |  |  | !   nadd     number of spectra summed into s3
 | 
					
						
							| 
									
										
										
										
											2015-02-05 17:43:43 +00:00
										 |  |  | !   nqd      0/1 to indicate decode attempt at QSO frequency
 | 
					
						
							| 
									
										
										
										
											2014-12-03 00:06:54 +00:00
										 |  |  | 
 | 
					
						
							|  |  |  | ! Output:
 | 
					
						
							| 
									
										
										
										
											2015-04-22 17:48:03 +00:00
										 |  |  | !   ncount   number of symbols requiring correction (-1 for no KV decode)
 | 
					
						
							| 
									
										
										
										
											2014-12-03 00:06:54 +00:00
										 |  |  | !   nhist    maximum number of identical symbol values
 | 
					
						
							|  |  |  | !   decoded  decoded message (if ncount >=0)
 | 
					
						
							|  |  |  | !   ltext    true if decoded message is free text
 | 
					
						
							| 
									
										
										
										
											2015-12-09 21:02:37 +00:00
										 |  |  | !   nft      0=no decode; 1=FT decode; 2=hinted decode
 | 
					
						
							| 
									
										
										
										
											2014-12-03 00:06:54 +00:00
										 |  |  | 
 | 
					
						
							|  |  |  |   use prog_args                       !shm_key, exe_dir, data_dir
 | 
					
						
							| 
									
										
										
										
											2015-04-22 17:48:03 +00:00
										 |  |  |   use packjt
 | 
					
						
							| 
									
										
										
										
											2016-03-10 14:25:22 +00:00
										 |  |  |   use jt65_mod
 | 
					
						
							| 
									
										
										
										
											2015-12-27 15:40:57 +00:00
										 |  |  |   use timer_module, only: timer
 | 
					
						
							| 
									
										
										
										
											2014-12-03 00:06:54 +00:00
										 |  |  | 
 | 
					
						
							|  |  |  |   real s3(64,63)
 | 
					
						
							|  |  |  |   character decoded*22
 | 
					
						
							| 
									
										
										
										
											2015-12-16 19:31:12 +00:00
										 |  |  |   character*12 mycall_12,hiscall_12
 | 
					
						
							|  |  |  |   character*6 mycall,hiscall,hisgrid
 | 
					
						
							| 
									
										
										
										
											2014-12-03 00:06:54 +00:00
										 |  |  |   integer dat4(12)
 | 
					
						
							|  |  |  |   integer mrsym(63),mr2sym(63),mrprob(63),mr2prob(63)
 | 
					
						
							| 
									
										
										
										
											2015-11-18 01:28:12 +00:00
										 |  |  |   integer correct(63),tmp(63)
 | 
					
						
							| 
									
										
										
										
											2015-12-15 21:24:22 +00:00
										 |  |  |   logical ltext
 | 
					
						
							| 
									
										
										
										
											2015-11-18 01:28:12 +00:00
										 |  |  |   common/chansyms65/correct
 | 
					
						
							| 
									
										
										
										
											2014-12-03 00:06:54 +00:00
										 |  |  |   save
 | 
					
						
							|  |  |  | 
 | 
					
						
							| 
									
										
										
										
											2016-06-10 14:18:10 +00:00
										 |  |  |   if(mode65.eq.-99) stop                   !Silence compiler warning
 | 
					
						
							| 
									
										
										
										
											2015-12-16 19:31:12 +00:00
										 |  |  |   mycall=mycall_12(1:6)
 | 
					
						
							|  |  |  |   hiscall=hiscall_12(1:6)
 | 
					
						
							| 
									
										
										
										
											2015-12-09 21:02:37 +00:00
										 |  |  |   qual=0.
 | 
					
						
							| 
									
										
										
										
											2015-11-18 01:28:12 +00:00
										 |  |  |   nbirdie=20
 | 
					
						
							|  |  |  |   npct=50
 | 
					
						
							|  |  |  |   afac1=1.1
 | 
					
						
							| 
									
										
										
										
											2015-12-09 21:02:37 +00:00
										 |  |  |   nft=0
 | 
					
						
							| 
									
										
										
										
											2014-12-03 00:06:54 +00:00
										 |  |  |   nfail=0
 | 
					
						
							|  |  |  |   decoded='                      '
 | 
					
						
							|  |  |  |   call pctile(s3,4032,npct,base)
 | 
					
						
							|  |  |  |   s3=s3/base
 | 
					
						
							| 
									
										
										
										
											2015-12-15 21:24:22 +00:00
										 |  |  |   s3a=s3                                            !###
 | 
					
						
							|  |  |  | 
 | 
					
						
							| 
									
										
										
										
											2014-12-03 00:06:54 +00:00
										 |  |  | ! Get most reliable and second-most-reliable symbol values, and their
 | 
					
						
							|  |  |  | ! probabilities
 | 
					
						
							|  |  |  | 1 call demod64a(s3,nadd,afac1,mrsym,mrprob,mr2sym,mr2prob,ntest,nlow)
 | 
					
						
							|  |  |  | 
 | 
					
						
							|  |  |  |   call chkhist(mrsym,nhist,ipk)       !Test for birdies and QRM
 | 
					
						
							|  |  |  |   if(nhist.ge.nbirdie) then
 | 
					
						
							|  |  |  |      nfail=nfail+1
 | 
					
						
							|  |  |  |      call pctile(s3,4032,npct,base)
 | 
					
						
							|  |  |  |      s3(ipk,1:63)=base
 | 
					
						
							|  |  |  |      if(nfail.gt.30) then
 | 
					
						
							|  |  |  |         decoded='                      '
 | 
					
						
							|  |  |  |         ncount=-1
 | 
					
						
							|  |  |  |         go to 900
 | 
					
						
							|  |  |  |      endif
 | 
					
						
							|  |  |  |      go to 1
 | 
					
						
							|  |  |  |   endif
 | 
					
						
							|  |  |  | 
 | 
					
						
							| 
									
										
										
										
											2015-12-09 21:02:37 +00:00
										 |  |  |   mrs=mrsym
 | 
					
						
							|  |  |  |   mrs2=mr2sym
 | 
					
						
							|  |  |  | 
 | 
					
						
							| 
									
										
										
										
											2015-04-22 17:48:03 +00:00
										 |  |  |   call graycode65(mrsym,63,-1)        !Remove gray code 
 | 
					
						
							|  |  |  |   call interleave63(mrsym,-1)         !Remove interleaving
 | 
					
						
							| 
									
										
										
										
											2014-12-03 00:06:54 +00:00
										 |  |  |   call interleave63(mrprob,-1)
 | 
					
						
							|  |  |  | 
 | 
					
						
							|  |  |  |   call graycode65(mr2sym,63,-1)      !Remove gray code and interleaving
 | 
					
						
							|  |  |  |   call interleave63(mr2sym,-1)       !from second-most-reliable symbols
 | 
					
						
							|  |  |  |   call interleave63(mr2prob,-1)
 | 
					
						
							| 
									
										
										
										
											2015-11-24 19:49:04 +00:00
										 |  |  |   ntry=0
 | 
					
						
							| 
									
										
										
										
											2014-12-03 00:06:54 +00:00
										 |  |  | 
 | 
					
						
							| 
									
										
										
										
											2015-12-07 13:44:00 +00:00
										 |  |  |   call timer('ftrsd   ',0)
 | 
					
						
							| 
									
										
										
										
											2015-12-15 21:24:22 +00:00
										 |  |  |   param=0
 | 
					
						
							|  |  |  |   call ftrsd2(mrsym,mrprob,mr2sym,mr2prob,ntrials,correct,param,ntry)
 | 
					
						
							| 
									
										
										
										
											2015-12-07 13:44:00 +00:00
										 |  |  |   call timer('ftrsd   ',1)
 | 
					
						
							| 
									
										
										
										
											2015-11-18 01:28:12 +00:00
										 |  |  |   ncandidates=param(0)
 | 
					
						
							|  |  |  |   nhard=param(1)
 | 
					
						
							|  |  |  |   nsoft=param(2)
 | 
					
						
							| 
									
										
										
										
											2015-12-06 04:37:14 +00:00
										 |  |  |   nerased=param(3)
 | 
					
						
							| 
									
										
										
										
											2015-12-31 01:30:31 +00:00
										 |  |  |   rtt=0.001*param(4)
 | 
					
						
							| 
									
										
										
										
											2015-12-09 21:02:37 +00:00
										 |  |  |   ntotal=param(5)
 | 
					
						
							| 
									
										
										
										
											2015-12-15 21:24:22 +00:00
										 |  |  |   qual=0.001*param(7)
 | 
					
						
							| 
									
										
										
										
											2015-12-31 01:30:31 +00:00
										 |  |  |   nd0=81
 | 
					
						
							|  |  |  |   r0=0.87
 | 
					
						
							|  |  |  |   if(naggressive.eq.10) then
 | 
					
						
							|  |  |  |      nd0=83
 | 
					
						
							|  |  |  |      r0=0.90
 | 
					
						
							|  |  |  |   endif
 | 
					
						
							|  |  |  |   if(ntotal.le.nd0 .and. rtt.le.r0) nft=1
 | 
					
						
							| 
									
										
										
										
											2015-12-09 21:02:37 +00:00
										 |  |  | 
 | 
					
						
							| 
									
										
										
										
											2016-03-21 16:03:11 +00:00
										 |  |  |   if(nft.eq.0 .and. iand(ndepth,32).eq.32) then
 | 
					
						
							| 
									
										
										
										
											2016-03-22 14:12:59 +00:00
										 |  |  |      qmin=2.0 - 0.1*naggressive
 | 
					
						
							| 
									
										
										
										
											2016-03-07 20:54:12 +00:00
										 |  |  |      call timer('hint65  ',0)
 | 
					
						
							| 
									
										
										
										
											2016-06-10 14:18:10 +00:00
										 |  |  |      call hint65(s3,mrs,mrs2,nadd,nflip,mycall,hiscall,hisgrid,qual,decoded)
 | 
					
						
							| 
									
										
										
										
											2015-12-15 21:24:22 +00:00
										 |  |  |      if(qual.ge.qmin) then
 | 
					
						
							| 
									
										
										
										
											2015-12-09 21:02:37 +00:00
										 |  |  |         nft=2
 | 
					
						
							| 
									
										
										
										
											2016-04-13 13:16:10 +00:00
										 |  |  |         ncount=0
 | 
					
						
							| 
									
										
										
										
											2015-12-09 21:02:37 +00:00
										 |  |  |      else
 | 
					
						
							| 
									
										
										
										
											2015-12-18 20:00:59 +00:00
										 |  |  |         decoded='                      '
 | 
					
						
							| 
									
										
										
										
											2015-12-09 21:02:37 +00:00
										 |  |  |         ntry=0
 | 
					
						
							| 
									
										
										
										
											2015-11-24 19:49:04 +00:00
										 |  |  |      endif
 | 
					
						
							| 
									
										
										
										
											2016-01-08 21:05:00 +00:00
										 |  |  |      call timer('hint65  ',1)
 | 
					
						
							| 
									
										
										
										
											2015-12-09 21:02:37 +00:00
										 |  |  |      go to 900
 | 
					
						
							| 
									
										
										
										
											2015-11-24 19:49:04 +00:00
										 |  |  |   endif
 | 
					
						
							|  |  |  | 
 | 
					
						
							| 
									
										
										
										
											2015-11-18 01:28:12 +00:00
										 |  |  |   ncount=-1
 | 
					
						
							| 
									
										
										
										
											2014-12-03 00:06:54 +00:00
										 |  |  |   decoded='                      '
 | 
					
						
							|  |  |  |   ltext=.false.
 | 
					
						
							| 
									
										
										
										
											2015-12-09 21:02:37 +00:00
										 |  |  |   if(nft.gt.0) then
 | 
					
						
							| 
									
										
										
										
											2015-12-15 21:24:22 +00:00
										 |  |  | ! Turn the corrected symbol array into channel symbols for subtraction;
 | 
					
						
							|  |  |  | ! pass it back to jt65a via common block "chansyms65".
 | 
					
						
							| 
									
										
										
										
											2015-11-24 19:49:04 +00:00
										 |  |  |      do i=1,12
 | 
					
						
							|  |  |  |         dat4(i)=correct(13-i)
 | 
					
						
							|  |  |  |      enddo
 | 
					
						
							|  |  |  |      do i=1,63
 | 
					
						
							| 
									
										
										
										
											2015-11-18 01:28:12 +00:00
										 |  |  |        tmp(i)=correct(64-i)
 | 
					
						
							|  |  |  |      enddo
 | 
					
						
							|  |  |  |      correct(1:63)=tmp(1:63)
 | 
					
						
							|  |  |  |      call interleave63(correct,63,1)
 | 
					
						
							|  |  |  |      call graycode65(correct,63,1)
 | 
					
						
							| 
									
										
										
										
											2014-12-03 00:06:54 +00:00
										 |  |  |      call unpackmsg(dat4,decoded)     !Unpack the user message
 | 
					
						
							| 
									
										
										
										
											2015-11-18 01:28:12 +00:00
										 |  |  |      ncount=0
 | 
					
						
							| 
									
										
										
										
											2014-12-03 00:06:54 +00:00
										 |  |  |      if(iand(dat4(10),8).ne.0) ltext=.true.
 | 
					
						
							|  |  |  |   endif
 | 
					
						
							| 
									
										
										
										
											2015-02-13 14:22:54 +00:00
										 |  |  | 900 continue
 | 
					
						
							| 
									
										
										
										
											2015-12-09 21:02:37 +00:00
										 |  |  |   if(nft.eq.1 .and. nhard.lt.0) decoded='                      '
 | 
					
						
							| 
									
										
										
										
											2015-12-15 21:24:22 +00:00
										 |  |  | 
 | 
					
						
							| 
									
										
										
										
											2015-02-13 14:22:54 +00:00
										 |  |  |   return
 | 
					
						
							| 
									
										
										
										
											2014-12-03 00:06:54 +00:00
										 |  |  | end subroutine extract
 | 
					
						
							| 
									
										
										
										
											2015-12-15 21:24:22 +00:00
										 |  |  | 
 | 
					
						
							|  |  |  | subroutine getpp(workdat,p)
 | 
					
						
							|  |  |  | 
 | 
					
						
							| 
									
										
										
										
											2016-03-10 14:25:22 +00:00
										 |  |  |   use jt65_mod
 | 
					
						
							| 
									
										
										
										
											2015-12-15 21:24:22 +00:00
										 |  |  |   integer workdat(63)
 | 
					
						
							|  |  |  |   integer a(63)
 | 
					
						
							|  |  |  | 
 | 
					
						
							|  |  |  |   a(1:63)=workdat(63:1:-1)
 | 
					
						
							|  |  |  |   call interleave63(a,1)
 | 
					
						
							|  |  |  |   call graycode(a,63,1,a)
 | 
					
						
							|  |  |  | 
 | 
					
						
							|  |  |  |   psum=0.
 | 
					
						
							|  |  |  |   do j=1,63
 | 
					
						
							|  |  |  |      i=a(j)+1
 | 
					
						
							|  |  |  |      x=s3a(i,j)
 | 
					
						
							|  |  |  |      s3a(i,j)=0.
 | 
					
						
							|  |  |  |      psum=psum + x
 | 
					
						
							|  |  |  |      s3a(i,j)=x
 | 
					
						
							|  |  |  |   enddo
 | 
					
						
							| 
									
										
										
										
											2015-12-31 01:30:31 +00:00
										 |  |  |   p=psum/63.0
 | 
					
						
							| 
									
										
										
										
											2015-12-15 21:24:22 +00:00
										 |  |  | 
 | 
					
						
							|  |  |  |   return
 | 
					
						
							|  |  |  | end subroutine getpp
 |