;
pro prob08d
;
defsysv, '!micrometer_by_meter', 1d6 / 1d
;
sav_filename = 'prob08a_060_.sav'
restore, filename = sav_filename, /verbose
;
ps_filename  = 'prob08d.ps'
;
xsize_inch            =  8.5
ysize_inch            = 11.0
width_inch            =  2.5
width_normal          = width_inch / xsize_inch
height_inch           =  2.5
height_normal         = height_inch / ysize_inch
height_spacer_inch    =  1.0
height_spacer_normal  = height_spacer_inch / ysize_inch
width_spacer_inch     =  1.0
width_spacer_normal   = width_spacer_inch / xsize_inch
yoffset_inch          =  1.0
yoffset_normal        = yoffset_inch / ysize_inch
xoffset_inch          =  1.5
xoffset_normal        = xoffset_inch / xsize_inch
;
set_plot, 'ps'
!p.font    =    0
device, /times, /inches, /color, font_size = 15, filename = ps_filename, $
  yoffset = 0, xoffset = 0., xsize = xsize_inch, ysize = ysize_inch
;
ix                    = 1
iy                    = 1
x_coord_min_normal    = xoffset_normal + (ix-1)*(width_spacer_normal + width_normal)
x_coord_max_normal    = x_coord_min_normal + width_normal
y_coord_min_normal    = 1.0 - yoffset_normal - iy*height_normal - (iy-1)*height_spacer_normal
y_coord_max_normal    = y_coord_min_normal + height_normal
!p.position           = [x_coord_min_normal,y_coord_min_normal,x_coord_max_normal,y_coord_max_normal]
;
;..plot terms in the kohler equation
;
yrange = [0.0001, 10.0]
yticks = 5
yminor = 9
ystyle = 1
ytitle = 'Terms in Kohler Equation'
;
xrange = [0.01, 10.0]
xticks = 3
xminor = 9
xstyle = 1
xtitle = 'Wet Radius, !9m!7m'
;
plot, /xlog, /ylog, /noerase, /nodata, fltarr(1), ystyle = ystyle, yrange = yrange, yticks = yticks, yminor = yminor, $
  xstyle = xstyle, xrange = xrange, xticks = xticks, xminor = xminor, ytitle = ytitle, $
  xtitle = xtitle, xticklen = -0.02, yticklen = -0.02
oplot, !micrometer_by_meter*ar, kohler_b/ar^3
oplot, !micrometer_by_meter*ar, kohler_a/ar
xyouts, x_coord_max_normal - 0.15, y_coord_max_normal - 0.02, 'r!Dd!N = ' + string(!micrometer_by_meter*rd, format='(f5.3)') + ' !9m!7m', /normal
;
device, /close
;
stop
;
end