ajout script piapp trace de traces satellites ds le plan angulaire , Reza 15/03/2019

Scripts/pltsattrack.pic

+#### Figures avec traces satellites / sources
+## script piapp 
+## Voir fichier p4sat.cmd - Reza , 15 Mars 2019 
+
+setaxesatt 'font=helvetica,bold,16 fixedfontsize minorticks'
+set xylim ( -5 5 -5 5 )
+set zsrc ( +11.4 -6.65 -25.4 -35 )
+set y0 ( +11 -6 -25 -35 )
+set srcname ( CasA CygA Crab/M1 VirA/M87 )
+
+set cols ( red orange siennared )
+set ANGX (90.-elevation)*cos((90.-azimuth)*M_PI/180.)
+set ANGY (90.-elevation)*sin((90.-azimuth)*M_PI/180.)
+
+defscript plttrk
+  ymin = $xylim[2]+$y0[J]
+  ymax = $xylim[3]+$y0[J]
+
+  drpanel $xylim[0] $xylim[1] $ymin $ymax 'stdaxes'
+  addcirc 0. $zsrc[J] 2.5 'line=solid,4 black'
+  addmarker 0. $zsrc[J] 'marker=fcircle,15 black'
+
+ for i 0:$#files 
+    del trk
+    openppf $files[i].ppf
+    rename $files[i] trk
+    plot2d trk $ANGX $ANGY 1 "same line=solid,2 $cols[i] cpts nsta notit"
+  end
+  del trk
+  openppf $SRC.ppf
+  rename $SRC trk
+  plot2d trk $ANGX $ANGY 1 "same line=dashed,2 navyblue cpts nsta notit"
+  setaxelabels 'EW-angle (degree)' 'NS-angle (degree)' 'font=helvetica,bolditalic,16 black'
+  settitle "$TRKTITLE" ' ' 'font=helvetica,bold,16 black'
+endscript
+
+J = 0
+set files ( trk_40889_20190114 trk_43566_20190115 trk_43055_20190115 )
+set SRC trk_CasA
+set TRKTITLE "LDJan2019  CasA , Galileo & CasA tracks "
+plttrk
+
+J = 1
+set files ( trk_43057_20190116 trk_41175_20190118 )
+set SRC trk_CygA
+set TRKTITLE "LDJan2019  CygA , Galileo & CygA tracks "
+plttrk
+
+J = 2
+set files ( trk_43055_20190119 trk_40534_20190119 )
+set SRC trk_M1 
+set TRKTITLE "LDJan2019  Crab/M1 , Galileo (red) & GPS (orange) M1/Crab tracks "
+plttrk
+
+J = 3
+set files ( trk_40890_20190122  trk_41859_20190122 )
+set SRC trk_M87 
+set TRKTITLE "LDJan2019  VirgoA/M87 , Galileo M87/VirgoA tracks "
+plttrk
+
+trk_40890_20190122
+@trk  dt_LD14jan19  12370,12410  1278.5  trk_41859_20190122