Commit 9065c96c authored by Reza  ANSARI's avatar Reza ANSARI
Browse files

MAJ des scripts suite analyse H,V de LD14jan19, Reza 30/07/2019

parent 85a0adb4
...@@ -202,6 +202,8 @@ ${PRGBASE}/Objs/tfm2dt ../ChezOlivier/LDOutputs/cor2_LD14jan19.ppf 1342,1350 dt ...@@ -202,6 +202,8 @@ ${PRGBASE}/Objs/tfm2dt ../ChezOlivier/LDOutputs/cor2_LD14jan19.ppf 1342,1350 dt
${PRGBASE}/Objs/tfm2dt ../ChezOlivier/LDOutputs/cor2_LD14jan19.ppf 1392,1400 dtsky1396_LD14jan19.ppf ${PRGBASE}/Objs/tfm2dt ../ChezOlivier/LDOutputs/cor2_LD14jan19.ppf 1392,1400 dtsky1396_LD14jan19.ppf
## Extraction du signal pour les sources du ciel entre 1410-1418 MHz -> FreqCentral=1416 MHz ## Extraction du signal pour les sources du ciel entre 1410-1418 MHz -> FreqCentral=1416 MHz
${PRGBASE}/Objs/tfm2dt ../ChezOlivier/LDOutputs/cor2_LD14jan19.ppf 1410,1418 dtsky1416_LD14jan19.ppf ${PRGBASE}/Objs/tfm2dt ../ChezOlivier/LDOutputs/cor2_LD14jan19.ppf 1410,1418 dtsky1416_LD14jan19.ppf
## Extraction du signal pour les sources du ciel entre 1435-1442 MHz -> FreqCentral=1439 MHz
${PRGBASE}/Objs/tfm2dt ../ChezOlivier/LDOutputs/cor2_LD14jan19.ppf 1435,1442 dtsky1439_LD14jan19.ppf
#### Fabrication des DataTables track a partir des fichiers texte fourni par JEC #### Fabrication des DataTables track a partir des fichiers texte fourni par JEC
foreach f ( CasA_survey_2019-1-14 CasA_survey_2019-1-15 CasA_survey_2019-1-16 CygA_survey_2019-1-17 CygA_survey_2019-1-18 ) foreach f ( CasA_survey_2019-1-14 CasA_survey_2019-1-15 CasA_survey_2019-1-16 CygA_survey_2019-1-17 CygA_survey_2019-1-18 )
...@@ -343,6 +345,8 @@ setenv DATAS /sps/baoradio/Olivier/PAON4/Work/LDOutputs/ ...@@ -343,6 +345,8 @@ setenv DATAS /sps/baoradio/Olivier/PAON4/Work/LDOutputs/
$TFM2DT -doVV $DATAS/cor2_LD14jan19.ppf 1275,1282 dt_LD14jan19_V.ppf $TFM2DT -doVV $DATAS/cor2_LD14jan19.ppf 1275,1282 dt_LD14jan19_V.ppf
## Extraction du signal pour les satellites entre 1270-1274 MHz (GPS) raie fine a 1272 MHz ## Extraction du signal pour les satellites entre 1270-1274 MHz (GPS) raie fine a 1272 MHz
$TFM2DT -doVV $DATAS/cor2_LD14jan19.ppf 1270,1274 dt1272_LD14jan19_V.ppf $TFM2DT -doVV $DATAS/cor2_LD14jan19.ppf 1270,1274 dt1272_LD14jan19_V.ppf
## Extraction du signal pour les satellites entre 1258-1266 MHz (GPS) en evitant la raie fine a 1272 MHz
$TFM2DT -doVV $DATAS/cor2_LD14jan19.ppf 1258,1266 dt1262_LD14jan19_V.ppf
## Extraction du signal pour les satellites entre 1318-1328 MHz (GPS et Galileo faible intensite mais au milieu de la bande de frequence) ## Extraction du signal pour les satellites entre 1318-1328 MHz (GPS et Galileo faible intensite mais au milieu de la bande de frequence)
$TFM2DT -doVV $DATAS/cor2_LD14jan19.ppf 1318,1328 dt1323_LD14jan19_V.ppf $TFM2DT -doVV $DATAS/cor2_LD14jan19.ppf 1318,1328 dt1323_LD14jan19_V.ppf
...@@ -352,6 +356,8 @@ $TFM2DT -doVV $DATAS/cor2_LD14jan19.ppf 1342,1350 dtsky1346_LD14jan19_V.ppf ...@@ -352,6 +356,8 @@ $TFM2DT -doVV $DATAS/cor2_LD14jan19.ppf 1342,1350 dtsky1346_LD14jan19_V.ppf
$TFM2DT -doVV $DATAS/cor2_LD14jan19.ppf 1392,1400 dtsky1396_LD14jan19_V.ppf $TFM2DT -doVV $DATAS/cor2_LD14jan19.ppf 1392,1400 dtsky1396_LD14jan19_V.ppf
## Extraction du signal pour les sources du ciel entre 1410-1418 MHz -> FreqCentral=1416 MHz ## Extraction du signal pour les sources du ciel entre 1410-1418 MHz -> FreqCentral=1416 MHz
$TFM2DT -doVV $DATAS/cor2_LD14jan19.ppf 1410,1418 dtsky1416_LD14jan19_V.ppf $TFM2DT -doVV $DATAS/cor2_LD14jan19.ppf 1410,1418 dtsky1416_LD14jan19_V.ppf
## Extraction du signal pour les sources du ciel entre 1435-1442 MHz -> FreqCentral=1439 MHz
$TFM2DT -doVV $DATAS/cor2_LD14jan19.ppf 1435,1442 dtsky1439_LD14jan19_V.ppf
#### Les datacards pour les traces vers cygA et casA #### Les datacards pour les traces vers cygA et casA
# On fabrique ldtrk_casA_V.d et ldtrk_cygA_V.d a partir de ldtrk_casA.d et ldtrk_cygA.d # On fabrique ldtrk_casA_V.d et ldtrk_cygA_V.d a partir de ldtrk_casA.d et ldtrk_cygA.d
...@@ -396,3 +402,6 @@ $TRKFIT -phifreq -D 4.5 -out cygA_V.txt -ckf cygA_V.ppf ldtrk_cygA_V.d > resfitc ...@@ -396,3 +402,6 @@ $TRKFIT -phifreq -D 4.5 -out cygA_V.txt -ckf cygA_V.ppf ldtrk_cygA_V.d > resfitc
$TRKFIT -phifreq -D 4.5 -out casA.txt -ckf casA.ppf ldtrk_casA.d > resfitcasA.log $TRKFIT -phifreq -D 4.5 -out casA.txt -ckf casA.ppf ldtrk_casA.d > resfitcasA.log
$TRKFIT -phifreq -D 4.5 -out cygA.txt -ckf cygA.ppf ldtrk_cygA.d > resfitcygA.log $TRKFIT -phifreq -D 4.5 -out cygA.txt -ckf cygA.ppf ldtrk_cygA.d > resfitcygA.log
set f = cygA
$TRKFIT -D 4.5 -out $f.txt -ckf $f.ppf ldtrk_$f.d
...@@ -191,3 +191,44 @@ setaxelabels 'Zenith Angle (deg)' 'Phase (deg)' 'font=helvetica,bolditalic,16 b ...@@ -191,3 +191,44 @@ setaxelabels 'Zenith Angle (deg)' 'Phase (deg)' 'font=helvetica,bolditalic,16 b
w2ps w2ps
######################################################################################
### Analyse ajustement phases en fonction de la frequence sur LD14jan19
######################################################################################
setaxesatt 'font=helvetica,bold,16 fixedfontsize minorticks'
newnt phase pol src az freq p2 p3 p4 p23 p24 p34
ntfrascii phase phase_HV.txt
print phase
set CLNCUT (fabs(freq-1323)>1.)
set PHI0_H ( 168. 65. 50. )
set APHI_H ( -95. 382. 439. )
set PHI0_V ( 95. 335 390. )
set APHI_V ( 0. -260. -180. )
newwin 3 2 1200 800
for i 1:3
if ( $i == 1 ) then
set TIT 'Polar H , fcircle->CasA, box->CygA'
set PHI0 ( $PHI0_H )
set APHI ( $APHI_H )
else
set TIT 'Polar V , fcircle->CasA, box->CygA'
set PHI0 ( $PHI0_V )
set APHI ( $APHI_V )
endif
for j 0:3
jj = $j+2
set VP p$jj
if ( $j == 2 ) then
set VP ((p4<50.)?p4+360.:p4)
endif
set phi0 $PHI0[j]
set aphi $APHI[j]
plot2d phase freq $VP (pol==$i)&&(src==1)&&$CLNCUT 'marker=fcircle,9 navyblue notit nsta'
plot2d phase freq $VP (pol==$i)&&(src==2) 'marker=box,11 same navyblue notit nsta'
func $phi0+(x-1250.)/250.*$aphi 1250. 1500. 10 'line=solid,2 same red'
setaxelabels 'Freq (MHz)' " Phase Phi_1$jj (deg)" 'font=helvetica,bolditalic,16 black'
settitle "$TIT" ' ' 'font=helvetica,bolditalic,16 black'
end
end
...@@ -79,3 +79,24 @@ end ...@@ -79,3 +79,24 @@ end
w2ps w2ps
######################################################################################
### Analyse Pointage H-V sur LD14jan19
######################################################################################
setaxesatt 'font=helvetica,bold,16 fixedfontsize minorticks'
newnt pointhv src az num D_H tet_H phi_H dtet_H dphi_H D_V tet_V phi_V dtet_V dphi_V
ntfrascii pointhv pointageHV_LD14jan19.txt
print pointhv
newwin 2 2 1200 800
for i 1:5
set PHI (dphi_H*M_PI/180.)
plot2dw pointhv dtet_H*cos($PHI) dtet_H*sin($PHI) az (num==$i) "marker=fcircle,15 colbr128 cmapdrw.font=helvetica,bold,14 nosizescale colscale notit nsta showcmap=right xylimits=-1.5,1.5,-1.5,1.5"
set PHI (dphi_V*M_PI/180.)
plot2dw pointhv dtet_V*cos($PHI) dtet_V*sin($PHI) az (num==$i) "marker=ftriangle,15 colbr128 nosizescale colscale nsta same"
settitle "Ant: $i Col->za H=circle V=triangle" ' ' 'font=helvetica,bold,16 black'
setaxelabels 'EW-Point-shift (deg)' 'NS-Point-Shift (deg) 'font=helvetica,bolditalic,16 black'
end
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