### Calcul et trace de lobes synthetises 
### BAORadio --  R. Ansari  -   Avril 2008 
### ===> Fichiers pa.pic + pa.cc + pa.h
defscript go 
  delobjs *
  c++execfrf pa.cc pa.h 
endscript

go

setaxesatt 'font=times,bolditalic,18 fixedfontsize minorticks'

defscript tracelobe
  set gropt "connectpoints xylimits=-90,90,50,280000 notit nsta"
  if ( $# > 0 ) then
    gropt = "$gropt $1" 
  endif
  plot2d lobes teta*180/M_PI wmod2 m==0 "$gropt" 
  plot2d lobes teta*180/M_PI wmod2 m==50 'connectpoints same blue nsta'
  plot2d lobes teta*180/M_PI wmod2 m==150 'connectpoints same orange nsta'
  plot2d lobes teta*180/M_PI wmod2 m==200 'connectpoints same red nsta'
  setaxelabels ' Angle Theta (degree) ' 'Synth. lobes' 'font=times,bold,18' 
  settitle 'Reconstructed lobes (FFT-like synthesis)' ' '  'font=times,bold,18'
  set txt ( 'Theta=0 deg' 'Theta=+12 deg' 'Theta=+36 deg' 'Theta=48 deg' )
  set cols ( 'black' 'blue' 'orange' 'red' )
  textdrawer txt cols 'font=times,bolditalic,18 frame' 
endscript 


defscript tracelobe2
  plot2d lobes teta*180/M_PI wmod2 m==0 "black connectpoints notit nsta" 
  plot2d lobes teta*180/M_PI wmod2 m==1 "blue same connectpoints notit nsta" 
  plot2d lobes teta*180/M_PI wmod2 m==2 "magenta same connectpoints notit nsta" 
  plot2d lobes teta*180/M_PI wmod2 m==3 "red same connectpoints notit nsta" 
  setaxelabels ' Angle Theta (degree) ' 'Lobe' 'font=times,bold,18' 
  settitle 'Reconstructed lobes - 7 Lambda/2 dipole sum' ' '  'font=times,bold,18'
endscript 

defscript tracelobe3
  plot2d lobes wy wx m==0 'black connectpoints'
  plot2d lobes wy wx m==2 'red same connectpoints'
  plot2d lobes wy wx m==-2 'red same connectpoints'
  plot2d lobes wy wx m==1 'orange same connectpoints'
  plot2d lobes wy wx m==-1 'orange same connectpoints'
  plot2d lobes wy wx m==3 'blue same connectpoints'
  plot2d lobes wy wx m==-3 'blue same connectpoints'
  plot2d lobes wy wx m==4 'violet same connectpoints'
  plot2d lobes wy wx m==-4 'violet same connectpoints'
endscript