#
#
#

include 'almasensitivity.g'

doit := function() {

  mysens := almasensitivity();
  mysens.usejanskies();


  mysens.setinstrument(nants=64, diam=12, baseline=12);

 
  elev := [20,30,40,50,60,70,80]
  sens1 := 0 * elev;
  sens2 := 0 * elev;
  for (ii in ind(elev)) {
     airmass := 1/sin(elev[ii]/ (180/pi) );

     freq := 680; 
     mysens.setobservation(npol=2, bandwidth=4, freq=freq, dt=60);

#     mysens.setatmos(pwv=2.4, airmass=airmass);
#     sens1[ii] := as_float( split( mysens.calcsensitivity1('obs') )[1] );

     mysens.setatmos(pwv=0.55, airmass=airmass);
     sens2[ii] :=  as_float( split( mysens.calcsensitivity1('obs') )[1] );
  }
  for (ii in ind(elev)) {
     print elev[ii],  sens2[ii]
  }
  mysens.done();
  return T;
}