;
;jeff snider, univeristy of wyoming, march 2006
;
pro plot_droplets, title1, title2, nconc, r, specific_volume, icolor
;
common  constants, gascon, rho_h2o, mw_nh42so4, rho_nh42so4, sigma_0,      $
        sigma_t, dsigmadmolality_nh42so4, dsigmadmolality_nacl, mw_h2o,    $
        mw_nacl, rho_nacl, tkmelt, c_nh42so4, a_nh42so4, c_nacl, a_nacl,   $
        cp_air_o, cp_h2o_o, mw_air, gravity, p0, alpha, beta_local,        $
        gascon_h2o, cw_h2o_o, tk_o, lv_o, ew_o, tk_aerosol, salt_type,     $
        epsilon, c_twomey, k_twomey, aerosol_type, right_tail,             $
        gascon_air, tc_base, tk_base, p_base, h_start, hmax,               $
        r0, dt, rmin, rmax, coef, epsilon_aerosol, specific_volume_meas,   $
        nchan, mw_salt, rho_salt, a_salt, c_salt, specific_volume_base,    $
        dsigmadmolality_salt, rho_insoluble, volume_insoluble_to_soluble,  $
        mixrat_tot_1,r_1,geo_sigma_1,mixrat_tot_2,r_2,geo_sigma_2,         $
        tk_start,p_start,sratio_start,mixrat_tot_3,r_3,geo_sigma_3
;
min_y = 1.0
max_y = 100000.
;
min_x = 0.001
max_x = 10.
;
!p.position = [0.15,0.55,0.45,0.85]
;
;..aerosol initialization has largest particles at index = 0, smallest at index = nchan
;
r1 = r(1:nchan)
r2 = r(0:nchan-1)
dmixingratiodlogr = nconc * specific_volume / (alog10 (r2/r1))
;
if icolor eq 1 then begin
 plot, /xlog, /ylog, /noerase, /noclip, /nodata, fltarr(1), psym = 0, xrange = [min_x,max_x], $
  xtitle='Radius, !7l!3m', ytitle='dN/dlog!I10!Nr, cm!E-3!N',yrange = [min_y, max_y], thick = 2, color = 0
 for i = 0, nchan - 2 do begin
  oplot, [1.e6*r2(i), 1.e6*r1(i)], [1.e-6*dmixingratiodlogr(i), 1.e-6*dmixingratiodlogr(i)], color = icolor
  oplot, [1.e6*r1(i), 1.e6*r1(i)], [1.e-6*dmixingratiodlogr(i), 1.e-6*dmixingratiodlogr(i+1)], color = icolor
 endfor
 xyouts, 0.950, 0.05, title2, size=0.4, /normal, ORIENTATION=90.
 xyouts, 0.965, 0.05, title1, size=0.4, /normal, ORIENTATION=90.
endif
;
if icolor eq 2 then begin
 for i = 0, nchan - 2 do begin
  oplot, [1.e6*r2(i), 1.e6*r1(i)], [1.e-6*dmixingratiodlogr(i), 1.e-6*dmixingratiodlogr(i)], color = icolor
  oplot, [1.e6*r1(i), 1.e6*r1(i)], [1.e-6*dmixingratiodlogr(i), 1.e-6*dmixingratiodlogr(i+1)], color = icolor
 endfor
 xyouts, 0.980, 0.05, title1, size=0.4, /normal, ORIENTATION=90.
endif
;
if icolor eq 3 then begin
 plots, [r0*1.e6,r0*1.e6], [min_y,max_y], color = 0
 for i = 0, nchan - 2 do begin
  oplot, [1.e6*r2(i), 1.e6*r1(i)], [1.e-6*dmixingratiodlogr(i), 1.e-6*dmixingratiodlogr(i)], color = icolor
  oplot, [1.e6*r1(i), 1.e6*r1(i)], [1.e-6*dmixingratiodlogr(i), 1.e-6*dmixingratiodlogr(i+1)], color = icolor
 endfor
 xyouts, 0.995, 0.05, title1, size=0.4, /normal, ORIENTATION=90.
endif
;
if icolor eq 4 then begin
 for i = 0, nchan - 2 do begin
  oplot, [1.e6*r2(i), 1.e6*r1(i)], [1.e-6*dmixingratiodlogr(i), 1.e-6*dmixingratiodlogr(i)], color = icolor
  oplot, [1.e6*r1(i), 1.e6*r1(i)], [1.e-6*dmixingratiodlogr(i), 1.e-6*dmixingratiodlogr(i+1)], color = icolor
 endfor
 xyouts, 1.010, 0.05, title1, size=0.4, /normal, ORIENTATION=90.
endif
;
return
;
end