;
; K-NGAO particulars
;
D_tele = 10.
subapsAcross = 44
pattern = 'K-NGAO 5'; 'K-NGAO 10'
model = 'KPAO model'
seed = 5
n_guidestars = 5; 10
radius = 1.0*arcmin/2.; 2.5*arcmin/2. ; 2.0*arcmin/2. ; 1.5*arcmin/2.
starList = 1
zenithAngle = 45*degrees
r0 = 18*cm
field_angle = radius
n_scienceStars = 1
science_angle = fltarr(n_scienceStars,2)
science_angle[0,*] = [0,0]
;science_angle[1,*] = [.5,0]
;science_angle[2,*] = [1.,0]
;science_angle[3,*] = [.5,.5]
;science_angle[4,*] = [0,.5]
;science_angle[5,*] = [0,1.]
science_angle *= field_angle
science_lambda = 2.2*microns
platescale = 0.02*arcsec
n_iterations = 5
n_realizations = 10
controlGain = 1.0
warm_restart = 1
AO_mode = 'MCAO'
dm_heights = [0, 10]*km
dm_weights = [[1,.8,.6,.34,.09,-.13,-.27], $
              [0,.2,.4,.66,.91,1.13,1.27]]
;
;  case 1: 1.5 arcmin diameter
;  case 2: 2 arcmin diameter
;  case 3: 2.5 arcmin diameter
;
;  - muck with GUI -
TomographySphericalWave_widget_realize, widget_id
forceColdRestart = 1
;
;  - program script -
TomographySphericalWave_Widget_DoSimulation,{top:widget_id}
;
field_angle = radius
n_scienceStars = 49; 25; 6
science_angle = fltarr(n_scienceStars,2)
  n = sqrt(n_scienceStars)
  sax = (ones(n) ## findgen(n))/float(n-1) - 0.5
  say = transpose(sax)
  sax = reform(sax,n_scienceStars)
  say = reform(say,n_scienceStars)
science_angle = [[sax],[say]]
science_angle *= field_angle
wavelengths = [0.9]*microns; [1.2,1.6,2.2]*microns; [0.5,0.7,0.9]*microns; [1.2,1.6,2.2]*microns
;
postProcess,wavelengths,science_angle
;postProcess,/display
end