p4sat.cmd 6.39 KB
Newer Older
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
######################################################################
#######  Exemples de commandes pour analyses des satellites ds PAON4 
#######  Analyse PAON4 , 
######## 1ere version , R. Ansari, Novembre 2018 
######################################################################

#####################
## A/ On utilise le programme predictsatsgp4 pour generer des fichiers de traces (position en fonction du temps) 
##    pour les satellites 
#####################

###### B01 , 20 Oct 2018 
# Satellites vers 15:21 UT de B01   Galileo 41859  
~/Work/GIT_BAORadio/AnaPAON4/Satellites/Objs/predictsatsgp4 -T "2018/10/20 15:21:18" -H 0.,81.1 -K 41859 -S 10 -s 1,0,0,15 -p 0 TLE_20181123/galileo.txt 

# Satellites vers 23:40 UT de B01   Galileo 41550 
~/Work/GIT_BAORadio/AnaPAON4/Satellites/Objs/predictsatsgp4 -T "2018/10/20 23:39:59.0" -H 0.,81.1 -K 41550 -S 10 -s 1,0,0,15 -p 0 TLE_20181123/galileo.txt

###### B02 , 21 Oct 2018 
# Satellites vers 22:40 UT de B02   Galileo  40890
~/Work/GIT_BAORadio/AnaPAON4/Satellites/Objs/predictsatsgp4 -T "2018/10/21 22:02:0.0" -H 0.,82.1 -K 40890 -S 10 -s 1,0,0,15 -p 0 TLE_20181123/galileo.txt

###### B04 , 26 Oct 2018 
# Satellites vers 27/10/2018 00:18:00 UT de B04   
~/Work/GIT_BAORadio/AnaPAON4/Satellites/Objs/predictsatsgp4 -T "2018/10/27 00:18:0.0" -H 0.,78.1 -K 40889 -s 1,0,0,15 -p 0 galileo.txt


#####################
## B/ On convertit les fichiers de trace au format texte en DataTable PPF ou FITS (bin-table) avec trk2dt 
#####################

#### Converting to PPF DataTable track files
~/Work/GIT_BAORadio/AnaPAON4/Satellites/Objs/trk2dt trk_41859_20181020.txt trk_B01_41859.ppf
~/Work/GIT_BAORadio/AnaPAON4/Satellites/Objs/trk2dt trk_41550_20181020.txt trk_B01_41550.ppf
~/Work/GIT_BAORadio/AnaPAON4/Satellites/Objs/trk2dt trk_40890_20181021.txt trk_B02_40890.ppf
~/Work/GIT_BAORadio/AnaPAON4/Satellites/Objs/trk2dt trk_40889_20181027.txt trk_B04_40889.ppf


#####################
## C/ On calcule les visibilites PAON4 a partir des fichiers de traces format PPF 
##    a l'aide du programme JSkyMap trk2vis 
#####################
#### Computing visibility DataTable from track files
~/Work/Jiao/JSkyMap/Objs/trk2vis -paon4 -in trk_B01_41859.ppf -freq 1285 -sdec 9 -out trkvis_B01_41859.ppf
~/Work/Jiao/JSkyMap/Objs/trk2vis -paon4 -in trk_B01_41550.ppf -freq 1285 -sdec 9 -out trkvis_B01_41550.ppf
~/Work/Jiao/JSkyMap/Objs/trk2vis -paon4 -in trk_B02_40890.ppf -freq 1285 -sdec 9 -out trkvis_B02_40890.ppf
~/Work/Jiao/JSkyMap/Objs/trk2vis -paon4 -in trk_B04_40889.ppf -freq 1285 -sdec 11.9 -out trkvis_B04_40889.ppf
48 49 50


########################################
51 52
### Analyse de A21_14oct18 (et B04_26oct18  B12_6dec18 ) avec fit des parametres de geometrie (pointage) et 
### fit des phases (Dec. 2018 
53

54
set PRGBASE = /sps/baoradio/Reza/AnaP4/AnaPAON4/
55
###  A21_14oct18
56 57 58
# 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/*
59

60 61
${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
62

63 64 65 66 67 68 69 70
# 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 
71 72 73
~/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

74 75
##############################################
#### il est possible aussi  d'extraire le DataTable depuis cartes TFM ds piapp
76 77 78 79 80 81 82 83 84 85
delobjs * 
openppf tfm_A21_14oct18.ppf
FREQ_MIN = 1275 
FREQ_MAX = 1282 
TIM_START = 1320
TIM_END = 1480
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
86
##############################################
87

88
###  A19_12oct18  (z=-5.58 deg , CygA) 
89

90 91 92
${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* 
93

94 95 96
${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
97

98 99 100 101
## 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 
102

103 104
#  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
105

106 107
#  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 
108

109 110 111 112 113 114 115 116 117
###  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*
118

119 120
${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
121

122
${PRGBASE}/Objs/tfm2dt ChezOlivier/Outputs/cor_B12_6dec18.ppf 1275,1282 dt_B12.ppf 
123

124 125
## 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
126 127