pro correlate_1c,user_name,trial=trial,orbit_num
  common curvefit_trace,rad,input_c,input_sigma,input_a,input_r,coeff_record
  coeff_record=0
  input_c=39
  input_sigma=0.608d
  input_r=30
  input_a=25

  device,decomposed=0
  loadct,39

  ;Load the data
  if n_elements(trial) eq 0 then trial=84
  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)
;  restore,get_debug_filename("cloud_initial_guess",trial_num=trial,orbit_num,user_name=user_name)
;  print,get_debug_filename("cloud_initial_guess",trial_num=trial,orbit_num,user_name=user_name)

  !p.multi=[0,4,4,0,0]

  range=[[30,60],[0,1],[0,60],[0,100]]
  zmin=min(/nan,(*(level_1c_retrievals.zenith_angle)),max=zmax)
  names=tag_names(level_1c_retrievals)
  k=[6,8,4,10]
  for i=0,3 do begin
    for j=0,3 do begin
      if i eq j then begin
        plot,(*(level_1c_retrievals.zenith_angle)),(*(level_1c_retrievals.(k[i]))),xtitle="Zenith angle",ytitle=names[k[i]],psym=3,charsize=2,yrange=reform(range[*,i]),/nodata
        plots,(*(level_1c_retrievals.zenith_angle)),(*(level_1c_retrievals.(k[i])))>range[0,i]<range[1,i],psym=3,color=linterp(zmin,0,zmax,254,(*(level_1c_retrievals.zenith_angle)))
      end else begin
        plot,(*(level_1c_retrievals.(k[i]))),(*(level_1c_retrievals.(k[j]))),xtitle=names[k[i]],ytitle=names[k[j]],xrange=reform(range[*,i]),yrange=reform(range[*,j]),psym=3,charsize=2,/nodata
        plots,(*(level_1c_retrievals.(k[i])))>range[0,i]<range[1,i],(*(level_1c_retrievals.(k[j])))>range[0,j]<range[1,j],psym=3,color=linterp(zmin,0,zmax,254,(*(level_1c_retrievals.zenith_angle)))
      end
    end
  end
  saveimage,"matrix.png",/png

  !p.multi=[0,1,1]

  for i=0,3 do begin
    for j=0,3 do begin
      erase
      if i eq j then begin
        plot,(*(level_1c_retrievals.zenith_angle)),(*(level_1c_retrievals.(k[i]))),xtitle="Zenith angle",ytitle=names[k[i]],psym=3,charsize=2,yrange=reform(range[*,i]),/nodata
        plots,(*(level_1c_retrievals.zenith_angle)),(*(level_1c_retrievals.(k[i])))>range[0,i]<range[1,i],psym=3,color=linterp(zmin,0,zmax,254,(*(level_1c_retrievals.zenith_angle)))
      end else begin
        plot,(*(level_1c_retrievals.(k[i]))),(*(level_1c_retrievals.(k[j]))),xtitle=names[k[i]],ytitle=names[k[j]],xrange=reform(range[*,i]),yrange=reform(range[*,j]),psym=3,charsize=2,/nodata
        plots,(*(level_1c_retrievals.(k[i])))>range[0,i]<range[1,i],(*(level_1c_retrievals.(k[j])))>range[0,j]<range[1,j],psym=3,color=linterp(zmin,0,zmax,254,(*(level_1c_retrievals.zenith_angle)))
      end
      saveimage,string(i,j,format='(%"matrix%1d_%1d.png")'),/png
    end
  end
end