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

MAJ des fichiers p4sat.cmd et ckfitres.pic pour l'analyse/ajustement des...

MAJ des fichiers p4sat.cmd et ckfitres.pic pour l'analyse/ajustement des traces satellites et sources de LD14jan19, Reza 22/02/2019
parent c182b816
#########################################################################
### Script piapp pour verifier la sortie du programme trkacfit.cc
### Script piapp pour verifier la sortie du programme trkacxfit.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
setaxesatt 'font=helvetica,bold,16 fixedfontsize minorticks'
delobjs *
if ( $# < 1 ) then
echo 'ckfitres.pic / missing arguments '
echo ' Usage: exec ckfitres.pic PPFName_fit6cx '
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'
defscript ckacx ' Check trkacxfit.cc result : Auto-Correlation and Cross-Cor fit for individual baselines ckacx BaseFileName NumSeq NTrk'
if ( $# < 3 ) then
echo 'ckacx / missing arguments '
echo ' Usage: ckacx BaseFileName NumSeq NbTrk (or NbSource) '
echo ' BaseFileName without .ppf, will open ac_NumSeq_BaseFileName.ppf and cx_NumSeq_BaseFileName.ppf'
return
endif
setaxesatt 'font=helvetica,bold,16 fixedfontsize minorticks'
set flnm $1
set seq $2
set NTrk $2
JJMAX = $NTrk+1
set flnmA ac_${seq}_${flnm}.ppf
set flnmB ac_${seq}_${flnm}.ppf
echo ' BaseFileName=' $flnm ' NumSeq=' $seq ' NTrk=' $NTrk ' --> JJMAX=' $JJMAX
echo '---ckacx--- 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 '---ckacx--- 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
end
echo '--- Opening ' $flnmB
openppf $flnmB
for j 1:$JJMAX
endscript
defscript ck6cx ' Check trkacxfit.cc result : Simultaneous 6 Cross-Cor fit : ck6cx PPFFileName NumSeq NTrk'
if ( $# < 3 ) then
echo 'ckacx / missing arguments '
echo ' Usage ckacx BaseFileName NumSeq NTrk (or NSource) '
echo ' BaseFileName without .ppf, will open ac_NumSeq_BaseFileName.ppf and cx_NumSeq_BaseFileName.ppf'
return
endif
setaxesatt 'font=helvetica,bold,16 fixedfontsize minorticks'
defscript ck6 ' ck6 numseq (1...n) numtrak (1...n) '
set trk $2
set seq $1
echo '--- ck6 script: NumSeq=' $seq ' NumTrack=' $trk
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'
for i 1:7
disp cx_${i}_${trk}_${seq} 'cdreal '
disp simcx_${i}_${trk}_${seq} 'cdreal same blue'
disp cx_${i}_${trk}_${seq} 'cdimag red same'
disp simcx_${i}_${trk}_${seq} 'cdimag same orange'
end
end
endscript
set flnm6cx $1
echo '----ckfitres.pic---- Opening file:' $flnm6cx
openppf $flnm6cx
echo '======= Use scripts ckacx ck6 to plot obs.visi / expected after fit :'
echo '---- Cx6Fit : ck6 numseq (1...n) numtrak (1...n)'
echo '---- AC,Cx Fits: ckacx BaseFileName NumSeq NbTrk (or NbSource) '
######################################################################
####### Exemples de commandes pour analyses des satellites ds PAON4
####### Analyse PAON4 ,
######## 1ere version , R. Ansari, Novembre 2018
######## 1ere version , R. Ansari, Novembre 2018
### MAJ Decembre 2018 , Fev 2018
### Voir les commandes pour l'analyse la plus recente en fin de fichier
######################################################################
#####################
......@@ -125,3 +127,79 @@ ${PRGBASE}/Objs/tfm2dt ChezOlivier/Outputs/cor_B12_6dec18.ppf 1275,1282 dt_B12.p
${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
##################################################################################################
##### 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
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