Join/Login

Business Software

Open Source Software

For Vendors

Blog

About

More

Articles

Create

SourceForge Podcast

Site Documentation

Subscribe to our Newsletter

Support Request

For Vendors Help Create Join Login

Business Software

Open Source Software

SourceForge Podcast

Resources

Articles

Case Studies

Blog

Menu

Help

Create

Join

Login

Home

Browse

WSJT

Echo

WSJT Echo

Weak signal ham radio communication

Brought to you by: bsomervi, chetfenn, dg2ycb, k1jt

Summary

Files

Reviews

Support

Mailing Lists

Git ▾

WSJT-X

WSPR

WSPR-X

WSJT

Echo

FMT

MAP65

MAP65IQ

WSJT-X Superbuild

More...

Old SVN

Menu ▾ ▴

Browse Commits

Fork

Merge Requests 0

Forks 1

Branches

master

[7a9056]: / f90 / avecho65.f90  Maximize  Restore  History

Download this file

123 lines (106 with data), 3.1 kB

  1
  2
  3
  4
  5
  6
  7
  8
  9
 10
 11
 12
 13
 14
 15
 16
 17
 18
 19
 20
 21
 22
 23
 24
 25
 26
 27
 28
 29
 30
 31
 32
 33
 34
 35
 36
 37
 38
 39
 40
 41
 42
 43
 44
 45
 46
 47
 48
 49
 50
 51
 52
 53
 54
 55
 56
 57
 58
 59
 60
 61
 62
 63
 64
 65
 66
 67
 68
 69
 70
 71
 72
 73
 74
 75
 76
 77
 78
 79
 80
 81
 82
 83
 84
 85
 86
 87
 88
 89
 90
 91
 92
 93
 94
 95
 96
 97
 98
 99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122

subroutine avecho65(cc,nutc,naz,nel,dop,nn,techo,fspread,fsample,i00,dphi,  &

     t0,f1a,dl,dc,pol,delta,rms1,rms2,snr,sigdb,dfreq,width,red,blue)

  parameter (NZ=520000,NZH=NZ/2,NTX=27*8192)

  parameter (NFFT=256*1024)

  complex cc(2,NZ)

  complex cx(0:NFFT-1),cy(0:NFFT-1)

  complex csx(-1000:1000),csy(-1000:1000)

  complex z

  real*4 sx(-1000:1000),sy(-1000:1000)

  real blue(2000),red(2000)

  integer ipkv(1)

  equivalence (ipk,ipkv)

  save sx,sy,nhdr

  abs2(z)=real(z)*real(z) + aimag(z)*aimag(z)

!  print*,'A',nn,dop,techo,fspread,fsample,i00,dphi

  if(nn.eq.0) then

     sx=0.

     sy=0.

     nhdr=1

  endif

  nn=nn+1

  if(nn.gt.100) nn=1

  cx(0:NZH-1)=cc(1,1:NZH)

  cy(0:NZH-1)=cc(2,1:NZH)

  call txtone(cx,tx,f1x,smx)

  call txtone(cy,ty,f1y,smy)

  t0=tx

  f1a=f1x

  if(smy.gt.smx) then

     t0=ty

     f1a=f1y

  endif

  write(*,3101) smy/smx,tx,ty,f1x,f1y,f1a

3101 format('B',f8.1,2f8.3,3f9.1)

  istart=nint((t0+techo)*fsample)

  iz2=min(NZ-1,istart+NTX-1)

  iz1=iz2-istart

  cx(0:iz1)=cc(1,istart:iz2)

  cx(iz1+1:)=0.

  cy(0:iz1)=cc(2,istart:iz2)

  cy(iz1+1:)=0.

  fdop=f1a+dop

  sq1=0.

  sq2=0.

  do i=0,NTX

     sq1=sq1 + abs2(cx(i))

     sq2=sq2 + abs2(cy(i))

  enddo

  rms1=sqrt(0.5*sq1/NTX)

  rms2=sqrt(0.5*sq2/NTX)

  call cspec(cx,fdop,csx)

  call cspec(cy,fdop,csy)

!  open(21,file='emecho.dat',status='unknown',access='stream',position='append')

!  write(21) nutc,naz,nel,dop,techo,fspread,csx,csy

!  close(21)

!### Do the following only when nn=1 ??  Or only on "first" call ??

  smax=0.

  do i=-200,200

     sx(i)=sx(i) + abs2(csx(i))

     sy(i)=sy(i) + abs2(csy(i))

     if(sx(i).gt.smax .or. sy(i).gt.smax) then

        smax=max(sx(i),sy(i))

        i0=i

     endif

  enddo

!###  i0=i00

  call polfit(csx,csy,nn,i0,dphi,dl,dc,pol,delta,red,blue)

  ave=min(sum(red(1101:1500))/400.,sum(red(1501:1900))/400.)

  sq1=dot_product(red(1101:1500)-ave,red(1101:1500)-ave)

  sq2=dot_product(red(1501:1900)-ave,red(1501:1900)-ave)

  rms=sqrt(min(sq1,sq2)/399.0)

  redmax=maxval(red(950:1050))

  snr=min(999.9,(redmax-ave)/rms)

  ipkv=maxloc(red(951:1050))+950

  df=fsample/NFFT

  dfreq=(ipk-1001)*df

  s=0.

  do i=ipk,ipk+300

     if(red(i).lt.1.0) exit

     s=s+(red(i)-ave)

  enddo

  do i=ipk-1,ipk-300,-1

     if(red(i).lt.1.0) exit

     s=s+(red(i)-ave)

  enddo

  bins=s/(red(ipk)-ave)

  width=df*bins

  sigdb=-99.0

  if(ave.gt.0.0) sigdb=10.0*log10(redmax/ave - 1.0) - 35.7

  open(20,file='emecho.txt',status='unknown',access='append')

  ndphi=nint(dphi)

  nf1a=nint(f1a)

  npol=nint(pol)

  ndelta=nint(delta)

  ndb=nint(sigdb)

  nsnr=nint(snr)

  if(nhdr.eq.1) write(20,1000) ndphi

1000 format(/'dPhi:',i5/'  UTC    N  Az  El  dB  S/N   DF    W  Pol   d    Lin  Circ Teme   Dop   Spread T0   F1a  Rx Ry'/  &

          95('-'))

  nhdr=0

  write(20,1010) nutc,nn,naz,nel,ndb,nsnr,dfreq,width,npol,ndelta,dl,dc,   &

       techo,dop,fspread,t0,0.001*f1a,nint(rms1),nint(rms2)

1010 format(i6.6,i4,4i4,2f6.1,i4,i5,2f6.2,f5.2,f8.1,f6.1,f5.2,f6.1,2i3)

  close(20)

  return

end subroutine avecho65

SourceForge

Create a Project

Open Source Software

Business Software

Top Downloaded Projects

Company

About

Team

SourceForge Headquarters
225 Broadway Suite 1600
San Diego, CA 92101
+1 (858) 422-6466

Resources

Support

Site Documentation

Site Status

SourceForge Reviews

© 2025 Slashdot Media. All Rights Reserved.

Terms Privacy Opt Out Advertise

×
Want the latest updates on software, tech news, and AI?

Get latest updates about software, tech news, and AI from SourceForge directly in your inbox once a month.

Submit