;
pro prob08a
;
defsysv, '!nanometer_by_meter', 1d9 / 1
;
rd           = 0.060d-6  
sav_filename = 'prob08a_' + string(nearest_integer(!nanometer_by_meter*rd), format='(i03)') + '_.sav'
;
phycon
;
defsysv, '!kohler_tk', 278.15d                     ; K, used in Kelvin term
defsysv, '!kohler_sigma', 0.072d                   ; J/m^2, vapor-liquid surface energy
defsysv, '!kohler_a', 2.d * !kohler_sigma / (!rhol*!Rv*!kohler_tk)
;
defsysv, '!kohler_rhos', 2165.d                    ; kg/m^3, density of sodium chloride
defsysv, '!kohler_ms', 0.05844d                    ; kg/mol, molecular weight of sodium chloride
defsysv, '!kohler_i', 2.d                          ; Van't Hoff factor for sodium chloride
defsysv, '!kohler_mw', 0.018d                      ; kg/mol, molecular weight of water
;
;..create an array of solution droplet radius starting at r = 0.07x10^-6.  Extend to r = 10x10^-6 m
;
x      = 0.07d-6 ; careful here
flag   = 0
while x le 10e-6 do begin
  if flag eq 0 then begin
    ar   = [x]
    flag = 1
  endif else begin
    ar   = [ar, x]
  endelse
  x   = x * 1.05
endwhile
print, 'Here: ', ar[0], ar[-1]
;
;..evaluate Kohler saturation
;      
defsysv, '!kohler_b', !kohler_i * !kohler_mw * !kohler_rhos * rd^3 / (!rhol * !kohler_ms)
as_kohler = 1.d + !kohler_a / ar - !kohler_b / ar^3
kohler_b  = !kohler_b
kohler_a  = !kohler_a
save, filename = sav_filename, ar, rd, as_kohler, kohler_b, kohler_a
;
end