#
#   HEY!!! PHASE IS MULTIPLIED BY THIS NUMBER!
phmult := 3.0
#
#   Band 5's conditions are non-optimal, and we are limited by the
#   atmosphere, how much can FITORDER > 2 help?
#  
#
sourcesize := 0.05
el := 60
#
include 'sdatmosim.g'
include 'almatau.g'



tau225good :=
[0.328999996,0.175999999,0.116999999,0.0839999989,0.0390000008,0.057,0.0450000018,0.0359999985,0.0309999995,0.023]  
phagood := 
[5.19000006,4.11999989,3.26999998,2.11999989,5.09000015,1.97000003,1.71500003,1.44000006,1.00999999,0.689999998]
tau225bad  :=
[0.138999999,0.129999995,0.111000001,0.0890000015,0.0549999997,0.057,0.0480000004,0.0399999991,0.0340000018,0.0270000007]  
phabad  :=
[8.29500008,5.73000002,4.2249999,3.1500001,4.31500006,2.6099999,2.05500007,1.58000004,1.23500001,0.879999995]
freq := [43, 80, 90, 145, 190, 230, 345, 500, 680, 880];
iband := 5
#
phagood := phagood * phmult;
phabad  := phabad  * phmult;
#
mytautool := almatau();
pwv := mytautool.getpwv(tau225good[iband]);
mytautool.done();
#
rec1 := optimalvslew(freq=freq[iband], sourcesize=sourcesize, pwv=pwv, stiphase=phagood[iband], el=el, fitorder=1 );
rec2 := optimalvslew(freq=freq[iband], sourcesize=sourcesize, pwv=pwv, stiphase=phagood[iband], el=el, fitorder=2 );
rec3 := optimalvslew(freq=freq[iband], sourcesize=sourcesize, pwv=pwv, stiphase=phagood[iband], el=el, fitorder=3 );
rec4 := optimalvslew(freq=freq[iband], sourcesize=sourcesize, pwv=pwv, stiphase=phagood[iband], el=el, fitorder=4 );
rec5 := optimalvslew(freq=freq[iband], sourcesize=sourcesize, pwv=pwv, stiphase=phagood[iband], el=el, fitorder=5 );
rec6 := optimalvslew(freq=freq[iband], sourcesize=sourcesize, pwv=pwv, stiphase=phagood[iband], el=el, fitorder=6 );
rec7 := optimalvslew(freq=freq[iband], sourcesize=sourcesize, pwv=pwv, stiphase=phagood[iband], el=el, fitorder=7 );
rec8 := optimalvslew(freq=freq[iband], sourcesize=sourcesize, pwv=pwv, stiphase=phagood[iband], el=el, fitorder=8 );

print 'Look at rec1, rec2, rec3, rec4, rec5, rec6, rec7, rec8'

huh := function(anum, nn=3) {
  aa := 10^nn
  return( as_integer( aa * anum ) / aa );
}


 print huh(rec1.bestvslew), huh(rec1.bestdt), huh(rec1.bestsigma,6), huh(rec1.bestsigmanoatm,6)
 print huh(rec2.bestvslew), huh(rec2.bestdt), huh(rec2.bestsigma,6), huh(rec2.bestsigmanoatm,6)
 print huh(rec3.bestvslew), huh(rec3.bestdt), huh(rec3.bestsigma,6), huh(rec3.bestsigmanoatm,6)
 print huh(rec4.bestvslew), huh(rec4.bestdt), huh(rec4.bestsigma,6), huh(rec4.bestsigmanoatm,6)
 print huh(rec5.bestvslew), huh(rec5.bestdt), huh(rec5.bestsigma,6), huh(rec5.bestsigmanoatm,6)
 print huh(rec6.bestvslew), huh(rec6.bestdt), huh(rec6.bestsigma,6), huh(rec6.bestsigmanoatm,6)
 print huh(rec7.bestvslew), huh(rec7.bestdt), huh(rec7.bestsigma,6), huh(rec7.bestsigmanoatm,6)
 print huh(rec8.bestvslew), huh(rec8.bestdt), huh(rec8.bestsigma,6), huh(rec8.bestsigmanoatm,6)