ASTRO> time 12:00:00.000 02-004-2022
ASTRO> obs Pico
I-OBSERVATORY,  Selected PICOVELETA observatory
I-OBSERVATORY,  PICO status is set to year 2021
   >> Current Outer baseband frequency range (as of 2021)
   >> Folding frequency = 16.00 GHz
ASTRO> let frequency 86.1
ASTRO> planet

 MARS         RA 21:33:03.2682 Dec -15:50:06.161 Az     46.59298 El     22.72529 Sun.D.   52.3 Vel.  -10.443
 DE= 1.80 DS= 1.42 Maj=  5.21 Min=  5.18 PA=356.38 TB=214.09 S(  86.1)=   23.95
 Frequency   86.1: Beam   27.9  Tmb     5.03  Flux     23.66  Size   28.0

URANUS       RA 02:42:31.3682 Dec 15:22:47.939 Az    -62.53318 El     53.46409 Sun.D.   30.2 Vel.   15.024
 DE=20.57 DS=19.71 Maj=  3.43 Min=  3.35 PA=265.32 TB=139.10 S(  86.1)=    6.55
 Frequency   86.1: Beam   27.9  Tmb     1.39  Flux      6.52  Size   27.9


+++++ WITH URANUS +++++

E090H Pol. H

Uranus_PolH -> Ta/gc = 1.00 +/- 0.07 K

S/Ta*_E090H = (6.52+/-10% Jy) / (1.00+/-0.07 K) = 6.52+/-0.79 Jy/K  -> G_E090H = 1./(S/Ta*) = 0.153 +/- 0.019 K/Jy

E090 Pol. V

Uranus_PolV -> Ta/gc = 1.01 +/- 0.06 K

S/Ta*_E090V = (6.52+/-10% Jy) / (1.01+/-0.06 K)) = 6.45+/-0.75 Jy/K -> G_E090V = 1/(S/Ta*) = 0.155 +/- 0.018 K/Jy 


///// WITH MARS ////

E090 Pol. H

Mars_PolH -> Ta/gc = 3.42 +/- 0.06 K

S/Ta*_E090H = (23.66+-/10% Jy) / (3.42+/-0.06 K) = 6.92+/-0.70 Jy/K  -> G_E090H = 1/(S/Ta*) = 0.145 +/- 0.015 K/Jy 

E090 Pol. V

Mars_PolV -> Ta/gc = 3.50 +/- 0.07 K

S/Ta*_E090H = (23.66+/-10% Jy) / (3.50+/-0.07 K) = 6.76+/-0.69 Jy/K  -> G_E090H = 1/(S/Ta*) = 0.148 +/- 0.015 K/Jy


======= AVERAGE URANUS & MARS =======

E090 Pol H.:

S/Ta*_E090H = 6.72 +/- 0.53 Jy/K  -> G_H = 0.149 +/- 0.012 K/ Jy

In [13]: np.mean([SM/Ta_H_Mars, SU/Ta_H])
Out[13]: 6.7190643274853805+/-0.5307402477633616

In [14]: 1/np.mean([SM/Ta_H_Mars, SU/Ta_H])
Out[14]: 0.14883024648377663+/-0.011756131217610174

E090 Pol. V:

In [26]: np.mean([STa_MV, STa_UV])
Out[26]: 6.607722772277228+/-0.5096687119092774

In [27]: 1/np.mean([STa_MV, STa_UV])
Out[27]: 0.15133806826695437+/-0.011673050001441681








Full log:

(base) torne@PMac ~$ ipython
Python 3.9.7 (default, Sep 16 2021, 08:50:36) 
Type 'copyright', 'credits' or 'license' for more information
IPython 7.29.0 -- An enhanced Interactive Python. Type '?' for help.

In [1]: from uncertainties import ufloat

In [2]: SU = ufloat(6.52, 6.52*0.1)

In [3]: SU
Out[3]: 6.52+/-0.652

In [4]: Ta_H = ufloat(1.00, 0.07)

In [5]: 

In [5]: SU/Ta_H
Out[5]: 6.52+/-0.7958674261458375

In [6]: 1./(SU/Ta_H)
Out[6]: 0.15337423312883436+/-0.018721711067076232

In [7]: SM = ufloat(23.66, 23.66*0.1)

In [8]: 

In [8]: Ta_H_Mars = ufloat(3.42, 0.06)

In [9]: SM/Ta_H_Mars
Out[9]: 6.918128654970761+/-0.7023787314431282

In [10]: 1/(SM/Ta_H_Mars)
Out[10]: 0.1445477599323753+/-0.014675539776395893

In [11]: np.mean([SM/Ta_H_Mars, SU/Ta_H])
---------------------------------------------------------------------------
NameError                                 Traceback (most recent call last)
<ipython-input-11-da764d1ecbdf> in <module>
----> 1 np.mean([SM/Ta_H_Mars, SU/Ta_H])

NameError: name 'np' is not defined

In [12]: import numpy as np

In [13]: np.mean([SM/Ta_H_Mars, SU/Ta_H])
Out[13]: 6.7190643274853805+/-0.5307402477633616

In [14]: 1/np.mean([SM/Ta_H_Mars, SU/Ta_H])
Out[14]: 0.14883024648377663+/-0.011756131217610174

In [15]: 

In [15]: 

In [15]: print("Pol.V")
Pol.V

In [16]: Ta_Uranus_V = ufloat(1.01, 0.06)

In [17]: Ta_Mars_V = ufloat(3.50, 0.07)

In [18]: SU
Out[18]: 6.52+/-0.652

In [19]: SM
Out[19]: 23.66+/-2.366

In [20]: SM/Ta_Mars_V
Out[20]: 6.76+/-0.6893874382377445

In [21]: SU/Ta_Uranus_V
Out[21]: 6.455445544554455+/-0.7508620003682754

In [22]: STa_MV = SM/Ta_Mars_V

In [23]: STa_UV = SU/Ta_Uranus_V

In [24]: STa_MV
Out[24]: 6.76+/-0.6893874382377445

In [25]: STa_UV
Out[25]: 6.455445544554455+/-0.7508620003682754

In [26]: np.mean([STa_MV, STa_UV])
Out[26]: 6.607722772277228+/-0.5096687119092774

In [27]: 1/np.mean([STa_MV, STa_UV])
Out[27]: 0.15133806826695437+/-0.011673050001441681

In [28]: 1/SU/Ta_Uranus_V
Out[28]: 0.15185567636518252+/-0.01766301894050104

In [29]: 1/SM/Ta_Mars_V
Out[29]: 0.012075836251660426+/-0.0012314984938033534

In [30]: 1/(SM/Ta_Mars_V)
Out[30]: 0.14792899408284024+/-0.015085856549091082

In [31]: 1/(SU/Ta_Uranus_V)
Out[31]: 0.15490797546012272+/-0.018018045621205115

In [32]: np.mean([6.92, 6.52])
Out[32]: 6.72

In [33]: np.mean([6.45, 6.76])
Out[33]: 6.605


