#
#
sourcesize := 0.05
el := 60
#
#
#
include 'sdatmosim.g'
include 'almatau.g'
#

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



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];
#
mytautool := almatau();
#
#
iband := 9;  #  ie, 680 GHZ

gerr := [0.00001, 0.0000215, 0.0000464, 0.0001, 0.000215, 0.000464, 0.001];
anoise1 := 0*gerr;  anoiseatm1 := 0*gerr; aslew1 := 0*gerr; adt1 := 0*gerr;


for (ii in ind(gerr) ) {

  pwv := mytautool.getpwv(tau225good[iband]);
  print '***************************'
  print '***  gainerr[', ii,']=', gerr[ii], ' ****'
  print '***************************'
  results := optimalvslew(freq=freq[iband], sourcesize=sourcesize, 
                       pwv=pwv, stiphase=phagood[iband], el=el, gainerr=gerr[ii]);

  print results;

  anoise1[ii]    := results.bestsigmanoatm;
  anoiseatm1[ii] := results.bestsigma;
  aslew1[ii]     := results.bestvslew;
  adt1[ii]       := results.bestdt;

}


print 'ii  gerr   noise  noise+atm  slew  dt';
for (ii in ind(gerr) ) {
   print ii, gerr[ii], huh(anoise1[ii], 5), huh(anoiseatm1[ii], 5),
	huh(aslew1[ii]), huh(adt1[ii], 4);

}