Commit dd96558c authored by Reza  ANSARI's avatar Reza ANSARI
Browse files

MAJ commandes ajustement signaux satellites et ajout script verification...

MAJ commandes ajustement signaux satellites et ajout script verification resultat fit par trkacfit.cc , Reza 27/12/2018
parent 7bb2abbb
#########################################################################
### Script piapp pour verifier la sortie du programme trkacfit.cc
### Fit de parametres de D_dish, pointage des antennes, phases
### sur auto-cor et cross-cor
### Usage : exec ckfitres BaseFileName NTrk
### BaseFileName --> fichiers BaseFileName.ppf et BaseFileName_cx.ppf
#########################################################################
if ( $# < 2 ) then
echo 'ckfitres / missing arguments '
echo ' Usage ckfitres BaseFileName NTrk (or NSat) '
echo ' BaseFileName without .ppf, will open BaseFileName.ppf and BaseFileName_cx.ppf'
return
endif
set flnm $1
set flnmA ${flnm}.ppf
set flnmB ${flnm}_cx.ppf
set NTrk $2
JJMAX = $NTrk+1
echo ' BaseFileName=' $flnm ' NTrk=' $NTrk ' --> JJMAX=' $JJMAX
delobjs *
setaxesatt 'font=helvetica,bold,16 fixedfontsize minorticks'
echo '--- Opening ' $flnmA
openppf $flnmA
for j 1:$JJMAX
newwin 2 2 600 500
for i 1:5
disp ac_${i}_${j} 'black'
disp simac_${i}_${j} 'blue same'
end
end
echo '--- Opening ' $flnmB
openppf $flnmB
for j 1:$JJMAX
newwin 3 2 800 500
for i 1:7
disp cx_${i}_${j} 'cdreal black'
disp simcx_${i}_${j} 'cdreal same blue'
disp cx_${i}_${j} 'cdimag same red'
disp simcx_${i}_${j} 'cdimag same orange'
end
end
......@@ -48,20 +48,31 @@
########################################
### Analyse de A21_14oct18 B04_26oct18 B12_6dec18
### Analyse de A21_14oct18 (et B04_26oct18 B12_6dec18 ) avec fit des parametres de geometrie (pointage) et
### fit des phases (Dec. 2018
set PRGBASE = /sps/baoradio/Reza/AnaP4/AnaPAON4/
### A21_14oct18
~/Work/GIT_BAORadio/AnaPAON4/Satellites/Objs/predictsatsgp4 -T "2018/10/14 23:00:00" -H 0.,78.58 -K 43057 TLE_20181123/*
~/Work/GIT_BAORadio/AnaPAON4/Satellites/Objs/predictsatsgp4 -T "2018/10/15 1:33:00" -H 0.,78.58 -K 43055 TLE_20181123/*
# 1/ Calcul et conversion de fichiers de traces satellites
${PRGBASE}/Satellites/Objs/predictsatsgp4 -T "2018/10/14 23:00:00" -H 0.,78.58 -K 43057 TLE_20181123/*
${PRGBASE}/Satellites/Objs/predictsatsgp4 -T "2018/10/15 1:33:00" -H 0.,78.58 -K 43055 TLE_20181123/*
~/Work/GIT_BAORadio/AnaPAON4/Satellites/Objs/trk2dt trk_43057_20181014.txt trk_A21_43057.ppf 2018/10/14
~/Work/GIT_BAORadio/AnaPAON4/Satellites/Objs/trk2dt trk_43055_20181015.txt trk_A21_43055.ppf 2018/10/14
${PRGBASE}/Satellites/Objs/trk2dt trk_43057_20181014.txt trk_A21_43057.ppf 2018/10/14
${PRGBASE}/Satellites/Objs/trk2dt trk_43055_20181015.txt trk_A21_43055.ppf 2018/10/14
# 2/ Extraction de DataTable Visi[Auto,Cross]=f(time) a partir de fichiers TFM
${PRGBASE}/Objs/tfm2dt ChezOlivier/Outputs/cor_A21_14oct18.ppf 1275,1282 dt_A21.ppf
# 3/ Ajustement geometrie/phases sur les auto / cross
${PRGBASE}/Objs/trkacfit -D 4.5 -sdec 10 -ckf tstA21.ppf -out A21_fac.txt dt_A21,trk_A21_43057,1358,1410,1278.5 dt_A21,trk_A21_43055,1512,1555,1278.5
# Etape utile pour faire des traces ds piapp
~/Work/GIT_BAORadio/AnaPAON4/Objs
~/Work/Jiao/JSkyMap/Objs/trk2vis -paon4 -in trk_A21_43057.ppf -freq 1278.5 -sdec 11.42 -out trkvis_A21_43057.ppf
~/Work/Jiao/JSkyMap/Objs/trk2vis -paon4 -in trk_A21_43055.ppf -freq 1278.5 -sdec 11.42 -out trkvis_A21_43055.ppf
#### ds piapp
##############################################
#### il est possible aussi d'extraire le DataTable depuis cartes TFM ds piapp
delobjs *
openppf tfm_A21_14oct18.ppf
FREQ_MIN = 1275
......@@ -72,45 +83,45 @@ c++args $FREQ_MIN $FREQ_MAX $TIM_START $TIM_END 1,1,1,1
c++execfrf anasat.cc
saveppf dtvis dt_A21.ppf
savefits dtvis dt_A21.fits
##############################################
### B04_26oct18
~/Work/GIT_BAORadio/AnaPAON4/Satellites/Objs/predictsatsgp4 -T "2018/10/27 00:18:00" -H 0.,78.08 -K 40889 TLE_20181123/*
~/Work/GIT_BAORadio/AnaPAON4/Satellites/Objs/trk2dt trk_40889_20181027.txt trk_B04_40889.ppf 2018/10/26
### A19_12oct18 (z=-5.58 deg , CygA)
#### ds piapp
delobjs *
openppf tfm_B04_26oct18.ppf
FREQ_MIN = 1275
FREQ_MAX = 1282
TIM_START = 1400
TIM_END = 1540
c++args $FREQ_MIN $FREQ_MAX $TIM_START $TIM_END 1,1,1,1
c++execfrf anasat.cc
saveppf dtvis dt_B04.ppf
savefits dtvis dt_B04.fits
${PRGBASE}/Satellites/Objs/predictsatsgp4 -T "2018/10/12 18: 5:00" -H 180.,84.42 -K 43002 TLE_20181123/bei*
${PRGBASE}/Satellites/Objs/predictsatsgp4 -T "2018/10/12 20:46:00" -H 180.,84.42 -K 40889 TLE_20181123/gal*
${PRGBASE}/Satellites/Objs/predictsatsgp4 -T "2018/10/13 3:24:00" -H 180.,84.42 -K 43647 TLE_20181123/bei*
${PRGBASE}/Satellites/Objs/trk2dt trk_43002_20181012.txt trk_A19_43002.ppf 2018/10/12
${PRGBASE}/Satellites/Objs/trk2dt trk_40889_20181012.txt trk_A19_40889.ppf 2018/10/12
${PRGBASE}/Satellites/Objs/trk2dt trk_43647_20181013.txt trk_A19_43647.ppf 2018/10/12
### B12_6dec18
~/Work/GIT_BAORadio/AnaPAON4/Satellites/Objs/predictsatsgp4 -T "2018/12/06 19:56:00" -H 0.,79.08 -K 43055 TLE_20181206/*
~/Work/GIT_BAORadio/AnaPAON4/Satellites/Objs/predictsatsgp4 -T "2018/12/07 4:18:00" -H 0.,79.08 -K 37847 TLE_20181206/*
## Bande de frequence de Beidou
${PRGBASE}/Objs/tfm2dt ChezOlivier/Outputs/cor_A19_12oct18.ppf 1268,1275 dt_A19_B.ppf
## Bande de frequence de Galileo
${PRGBASE}/Objs/tfm2dt ChezOlivier/Outputs/cor_A19_12oct18.ppf 1275,1282 dt_A19_G.ppf
~/Work/GIT_BAORadio/AnaPAON4/Satellites/Objs/trk2dt trk_43055_20181206.txt trk_B12_43055.ppf 2018/12/06
~/Work/GIT_BAORadio/AnaPAON4/Satellites/Objs/trk2dt trk_37847_20181207.txt trk_B12_37847.ppf 2018/12/06
# Fit avec satellites Beidou ne marche pas bien !
# ${PRGBASE}/Objs/trkacfit -D 4.5 -sdec -5 -ckf tstA19.ppf -out A19_fac.txt dt_A19_B,trk_A19_43002,1065,1105,1272 dt_A19_G,trk_A19_40889,1230,1260,1278.5 dt_A19_B,trk_A19_43647,1630,1660,1272
# On ne fit qu'avec le seul satellite Galileo ...
${PRGBASE}/Objs/trkacfit -D 4.5 -sdec -5 -ckf tstA19.ppf -out A19_fac.txt dt_A19_G,trk_A19_40889,1230,1260,1278.5
#### ds piapp
delobjs *
openppf tfm_B12_6dec18.ppf
FREQ_MIN = 1275
FREQ_MAX = 1282
TIM_START = 1175
TIM_END = 1225
c++args $FREQ_MIN $FREQ_MAX $TIM_START $TIM_END 1,1,1,1
c++execfrf anasat.cc
saveppf dtvis dt_B12.ppf
savefits dtvis dt_B12.fits
### B04_26oct18
${PRGBASE}/Satellites/Objs/predictsatsgp4 -T "2018/10/27 00:18:00" -H 0.,78.08 -K 40889 TLE_20181123/*
${PRGBASE}/Satellites/Objs/trk2dt trk_40889_20181027.txt trk_B04_40889.ppf 2018/10/26
### B12_6dec18 (z=+11.92 deg)
${PRGBASE}/Satellites/Objs/predictsatsgp4 -T "2018/12/06 19:56:00" -H 0.,79.08 -K 43055 TLE_20181206/gal*
${PRGBASE}/Satellites/Objs/predictsatsgp4 -T "2018/12/07 4:18:00" -H 0.,79.08 -K 37847 TLE_20181206/gal*
${PRGBASE}/Satellites/Objs/trk2dt trk_43055_20181206.txt trk_B12_43055.ppf 2018/12/06
${PRGBASE}/Satellites/Objs/trk2dt trk_37847_20181207.txt trk_B12_37847.ppf 2018/12/06
${PRGBASE}/Objs/tfm2dt ChezOlivier/Outputs/cor_B12_6dec18.ppf 1275,1282 dt_B12.ppf
## Ajustement
${PRGBASE}/Objs/trkacfit -D 4.5 -sdec 10 -ckf tstB12.ppf -out B12_fac.txt dt_B12,trk_B12_43055,1178,1215,1278.5 dt_B12,trk_B12_37847,1665,1695,1278.5
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment