#
#
iband := 7;
#
el := 60
#
#
#
include 'sdatmosim.g'
include 'almatau.g'
#
#  345 GHz:  lambd/2D = 0.002076 deg
#
sizes := [0.005, 0.01, 0.02, 0.04, 0.06, 0.08, 0.10];

anoise1 := 0*sizes;  anoiseatm1 := 0*sizes; aslew1 := 0*sizes; adt1 := 0*sizes;
bnoise1 := 0*sizes;  bnoiseatm1 := 0*sizes; bslew1 := 0*sizes; bdt1 := 0*sizes;

#
#  integrate through the various observing conditions
#


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();
#
#
for (isize in ind(sizes)) {
  pwv := mytautool.getpwv(tau225good[iband]);
  print 'HHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHH', isize, pwv, tau225good[iband]
  results := optimalvslew(freq=freq[iband], sourcesize=sizes[isize], 
                       pwv=pwv, stiphase=phagood[iband], el=el );

  print results;

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

#  pwv := mytautool.getpwv(tau225bad[iband]);
#  print 'IIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIII', isize, pwv, tau225good[iband]
#  results := optimalvslew(freq=freq[iband], sourcesize=sizes[isize], 
#                       pwv=pwv, stiphase=phabad[iband], el=el );
#
#  bnoise1[isize]    := results.bestsigmanoatm;
#  bnoiseatm1[isize] := results.bestsigma;
#  bslew1[isize]     := results.bestvslew;
#  bdt1[isize]       := results.bestdt;

}

print 'We did a bunch of different source sizes'
for (isize in ind(sizes)) {
  print 'A:', isize, sizes[isize], anoise1[isize],anoiseatm1[isize] ,aslew1[isize], adt1[isize]       
#  print 'B:', isize, sizes[isize], bnoise1[isize],bnoiseatm1[isize], bslew1[isize], bdt1[isize]       
}