Fix msk144sim - thanks Joe.

git-svn-id: svn+ssh://svn.code.sf.net/p/wsjt/wsjt/branches/wsjtx@6797 ab8295b8-cf94-4d9e-aec4-7959e3be5d79

lib/detectmsk144.f90

@@ -169,6 +169,8 @@ subroutine detectmsk144(cbig,n,pchk_file,lines,nmessages,nutc)
 
   do ip=1,ndet  !run through the candidates and try to sync/demod/decode
     imid=times(ip)*fs
+    if( imid .lt. NPTS/2 ) imid=NPTS/2
+    if( imid .gt. n-NPTS/2 ) imid=n-NPTS/2
     t0=times(ip)
     cdat=cbig(imid-NPTS/2+1:imid+NPTS/2)
     ferr=ferrs(ip)

lib/msk144sim.f90

@@ -3,7 +3,7 @@ program msk144sim
   use wavhdr
   parameter (NMAX=15*12000)
   real pings(0:NMAX-1)
-  real waveform(0:864-1)
+  real waveform(0:NMAX-1)
   character arg*8,msg*22,msgsent*22,fname*40
   character*512 pchk_file
   character*3 rpt(0:7)
@@ -12,8 +12,6 @@ program msk144sim
   type(hdr) h                          !Header for .wav file
   integer*2 iwave(0:NMAX-1)
   integer itone(144)                   !Message bits
-  integer b11(11)                      !Barker-11 code
-  data b11/1,1,1,0,0,0,1,0,0,1,0/
   data rpt /'26 ','27 ','28 ','R26','R27','R28','RRR','73 '/
 
   pchk_file='./peg-128-80-reg3.pchk'
@@ -47,7 +45,10 @@ program msk144sim
   dphi1=twopi*(freq+500)/12000.0
   phi=0.0
   indx=0
+  nreps=NMAX/(144*6)
+  do jrep=1,nreps
     do i=1,144
+      jj=mod(i,864)
       if( itone(i) .eq. 0 ) then
         dphi=dphi0
       else
@@ -59,6 +60,7 @@ program msk144sim
         phi=mod(phi+dphi,twopi)
       enddo 
     enddo
+  enddo 
 
   if(itype.lt.1 .or. itype.gt.7) then
      print*,'Illegal message'
@@ -73,12 +75,12 @@ program msk144sim
 1002 format('000000_',i6.6)
      open(10,file=fname(1:13)//'.wav',access='stream',status='unknown')
 
+     wave=0.0
+     iwave=0
      fac=sqrt(6000.0/2500.0)
-     j=-1
      do i=0,NMAX-1
-        j=mod(j+1,864)
         xx=gran()
-        wave(i)=pings(i)*waveform(j) + fac*xx
+        wave(i)=pings(i)*waveform(i) + fac*xx
         iwave(i)=30.0*wave(i)
      enddo