pro parameter_space

  ;these ofs's let you shift the image on the screen
  zoom=5
  xofs=-200
  yofs=0

  ;Load the data
  if n_elements(trial) eq 0 then trial=76
  if n_elements(user_name) eq 0 then user_name="jeppesen"
  if n_elements(orbit_num) eq 0 then orbit_num=0
;  restore,get_level_1c_filename(trial_num=trial,orbit_num,user_name=user_name)
;  print,get_level_1c_filename(trial_num=trial,orbit_num,user_name=user_name)
;  restore,get_debug_filename("tile",trial_num=trial,orbit_num,user_name=user_name)
;  print,get_debug_filename("tile",trial_num=trial,orbit_num,user_name=user_name)

  ;Plot the image
  wset,0
  s=size(*(level_1c_retrievals.scale_height_ratio),/dimensions)
  tvscl,/nan,rebin(*(level_1c_retrievals.albedo),s[0]*zoom,s[1]*zoom,/sample),xofs*zoom,yofs*zoom

  cursor,x,y,/device,/up

  x=x/zoom-xofs
  y=y/zoom-yofs
  print,x,y

  ;X marks the spot
  plots,([x-5,x+5]+xofs)*zoom+zoom/2,([y-5,y+5]+yofs)*zoom+zoom/2,/device,color=255
  plots,([x-5,x+5]+xofs)*zoom+zoom/2,([y+5,y-5]+yofs)*zoom+zoom/2,/device,color=255

  w=where(finite(tiles.rad[x,y,*]),nw)
  ssa=tiles.ssa[x,y,w]
  sva=tiles.sva[x,y,w]
  sza=tiles.sza[x,y,w]
  rad=tiles.rad[x,y,w]

  tiles=0
  str_free,level_1c_retrievals

  for i=0,10 do begin
    print,i
    parameter_space_step,a,rad,ssa,sva,sza
    print,n_elements(a)
    print,(size(a,/dimensions))[0]
    save,a,filename="parameter_space.sav"
  end
end