p4sat.cmd 9.89 KB
Newer Older
1 2 3
######################################################################
#######  Exemples de commandes pour analyses des satellites ds PAON4 
#######  Analyse PAON4 , 
4 5 6
######## 1ere version , R. Ansari, Novembre 2018
###   MAJ Decembre 2018 , Fev 2018
###   Voir les commandes pour l'analyse la plus recente en fin de fichier 
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
######################################################################

#####################
## 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
50 51 52


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

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

62 63
${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
64

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

76 77
##############################################
#### il est possible aussi  d'extraire le DataTable depuis cartes TFM ds piapp
78 79 80 81 82 83 84 85 86 87
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
88
##############################################
89

90
###  A19_12oct18  (z=-5.58 deg , CygA) 
91

92 93 94
${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* 
95

96 97 98
${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
99

100 101 102 103
## 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 
104

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

108 109
#  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 
110

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

121 122
${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
123

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

126 127
## 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
128 129


130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205
##################################################################################################
##### Fevrier 2019 ################
### Analyse acquisition longue duree Janv 2019 , a partir du fichier LD_goodsat.txt d'Olivier
##################################################################################################

### AU CC-IN2P3 : se mettre ds /sps/baoradio/Reza/AnaP4/SatLD14jan19
source /pbs/throng/baoradio/Library/Sophya/env.csh  
cd /sps/baoradio/Reza/AnaP4/SatLD14jan19
set PRGBASE = /sps/baoradio/Reza/AnaP4/AnaPAON4 


${PRGBASE}/Objs/tfm2dt ChezOlivier/LDOutputs/cor2_LD14jan19.ppf  1275,1282 dt_LD14jan19.ppf

## Temps en minutes 1115,1160
${PRGBASE}/Satellites/Objs/predictsatsgp4  -T "2019/01/14 18:56:00" -H 0.,79.08 -K 40889 TLE_20190125/gal*
${PRGBASE}/Satellites/Objs/trk2dt trk_40889_20190114.txt trk_40889_20190114.ppf 2019/01/14

## Temps en minutes 1460,1500
${PRGBASE}/Satellites/Objs/predictsatsgp4  -T "2019/01/15 0:43:00" -H 0.,79.08 -K 43566 TLE_20190125/gal*
${PRGBASE}/Satellites/Objs/trk2dt trk_43566_20190115.txt trk_43566_20190115.ppf 2019/01/14

## Temps en minutes 2455,2495 
${PRGBASE}/Satellites/Objs/predictsatsgp4  -T "2019/01/15 17:16:00" -H 0.,79.08 -K 43055 TLE_20190125/gal*
${PRGBASE}/Satellites/Objs/trk2dt trk_43055_20190115.txt trk_43055_20190115.ppf 2019/01/14

## Temps en minutes 4130,4175
${PRGBASE}/Satellites/Objs/predictsatsgp4  -T "2019/01/16 21:15:00" -H 180.,83.35 -K 43057 TLE_20190125/gal*
${PRGBASE}/Satellites/Objs/trk2dt trk_43057_20190116.txt trk_43057_20190116.ppf 2019/01/14

## Temps en minutes 6120,6160
${PRGBASE}/Satellites/Objs/predictsatsgp4  -T "2019/01/18 6:20:00" -H 180.,83.35 -K 41175 TLE_20190125/gal*
${PRGBASE}/Satellites/Objs/trk2dt trk_41175_20190118.txt trk_41175_20190118.ppf 2019/01/14

## Temps en minutes 11520+720,+760 : 12240,12280
${PRGBASE}/Satellites/Objs/predictsatsgp4  -T "2019/01/22 12:20:00" -H 180.,55 -K 40890 TLE_20190125/gal*
${PRGBASE}/Satellites/Objs/trk2dt trk_40890_20190122.txt trk_40890_20190122.ppf 2019/01/14

## Temps en minutes 11520+850,+890 : 12370,12410
${PRGBASE}/Satellites/Objs/predictsatsgp4  -T "2019/01/22 14:35:00" -H 180.,55 -K 41859 TLE_20190125/gal*
${PRGBASE}/Satellites/Objs/trk2dt trk_41859_20190122.txt trk_41859_20190122.ppf 2019/01/14

#########################
####### Fichiers .d : Tracks-list  
#### ldtrk_casA.d
--- Observation Longue duree Janvier 2019 
Sources avec observation vers CasA  (z=+11.42)
@zenang 10 
@trk  dt_LD14jan19  1115,1160  1278.5  trk_40889_20190114
@trk  dt_LD14jan19  1460,1500  1278.5  trk_43566_20190115
@trk  dt_LD14jan19  2455,2495  1278.5  trk_43055_20190115

#### ldtrk_cygA.d
--- Observation Longue duree Janvier 2019 
Sources avec observation vers CygA  (z=-6.5)
@zenang -6.
@trk  dt_LD14jan19  4130,4175  1278.5  trk_43057_20190116
@trk  dt_LD14jan19  6120,6160  1278.5  trk_41175_20190118

#### ldtrk_virA.d
--- Observation Longue duree Janvier 2019 
Sources avec observation vers VirgoA  (z=-35)
@zenang -35 
@trk  dt_LD14jan19  12240,12280  1278.5  trk_40890_20190122
@trk  dt_LD14jan19  12370,12410  1278.5  trk_41859_20190122


##### Exemple d'execution de commandes
${PRGBASE}/Objs/trkacxfit -docx6s -docx6f -D 4.5 -out ldj.txt -ckf ldj.ppf ldtrk_casA.d ldtrk_cygA.d ldtrk_virA.d
 
#### Utiliser le script piapp : ckfitres.pic  et es scripts definis dedans pour verifier les resultats
piapp> exec ckfitres.pic cxb6_ldj.ppf
piapp> ck6 1 1
piapp> ck6 1 2
piapp> ckacx 1 1
piapp> ckacx 1 1