#
#
iband1 := 1
iband2 := 10
#
#  note:  300 GHz: 0.00477 degrees beam
#             so... 0.05 degrees source size?  
#             10 x 10, or 20 x 20 = 400  source size?
#
sourcesize := 0.05
el := 60
#
#
#
include 'sdatmosim.g'
include 'almatau.g'
#
anoise1 := [1:10];  anoiseatm1 := [1:10];
bnoise1 := [1:10];  bnoiseatm1 := [1:10];

#
#  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 (iband in [iband1:iband2]) {
  pwv := mytautool.getpwv(tau225good[iband]);

  anoise1[iband]    := results.bestsigmanoatm;
  anoiseatm1[iband] := results.bestsigma;

  pwv := mytautool.getpwv(tau225bad[iband]);

  bnoise1[iband]    := results.bestsigmanoatm;
  bnoiseatm1[iband] := results.bestsigma;

}

for (iband in [iband1:iband2]) {
  print 'A:', iband, anoise1[iband],anoiseatm1[iband]
  print 'B:', iband, bnoise1[iband],bnoiseatm1[iband]
}