Extended IDL Help

This page was created by the IDL library routine mk_html_help. For more information on this routine, refer to the IDL Online Help Navigator or type: 

     ? mk_html_help

at the IDL command line prompt.

Last modified: Sat Sep 23 18:13:26 2006.

List of Routines

Routine Descriptions

OCLC_EY92_AAREADME

 [Next Routine] [List of Routines] 
 NAME:
   oclc_ey92_AAREADME
 PURPOSE
   Print the file 'oclc_ey92_AAREADME.pro'
 CATEGORY:
  Occultation lightcurve (oclc)
 CALLING SEQUENCE:
  oclc_ey92_AAREADME

 PROCEDURE

 README_oclc_ey92
 
 These routines implement the model of Elliot, J. L. and L. A. Young 1992.
 Analysis of stellar occultation data for planetary atmospheres. I - Model
 fitting, with application to Pluto. Astronomical Journal 103, 991-1015.
 
 The routines are called oclc_ey92_*.  oclc stands for OCcultation
 LightCurve.  EY92 refers to the Elliot and Young 1992 paper cited above.
 Other occultation lightcurves can be implemented, and will be prefixed
 oclc__.pro.
 
 ---------------------------------------------------------------------------
 VARIABLES
 ---------------------------------------------------------------------------
 
 a
 exponent for non-constant mean moleculr weight mu(r) = mu0 * (r/r0)^-a
 Unitless
 
 b
 Exponent for non-constant temperature, T(r) = T * (r/r0)^b
 Unitless.
 
 aser       A(a, b, delta) in EY92
 Series in delta for use in calculating theta
 Unitless.
 
 bser       B(a,b,delta) in EY92
 Series in delta for use in calculating dtheta
 Unitless.
 
 cser       C(delta_tau) in EY92
 Series in delta_tau for use in calculating optical depth
 Unitless.
 
 distobs       D in EY92
 Distance between target and observer
 cm
 
 dtheta     d theta/d r in EY92
 Derivative of the bending angle with respect to planet radius,
 positive for an isothermal atmosphere.
 radian/cm
 
 htau1
 scale height of absorbtion coefficient
 cm
 
 kappa
 linear absorption coefficient
 kappa = kappa1 exp( - (r-r1)/(htau1 * (r/r1) ) )     EY92 3.23
 1/cm (or cm^2/cm^3, cross section per particle * particle/volume)
 
 kappa1
 linear absorption coefficient at r1 (top of haze)
 1/cm (or cm^2/cm^3, cross section per particle * particle/volume)
 
 lam        lambda in EY92
 Ratio of radius (r) to scale height (H).  Since H = kT r^2/(mu m_adu M G),
 this is the same as the ratio of kinetic energy (kT) to potential
 energy (mu m_adu M G/r).
 unitless
 
 lam0        lambda0 in EY92
 Reference energy ratio. Energy ratio, lambda, evaluated at the
 reference radius, r0.
 unitless
 
 nu
 Refractivity. 
 Unitless
 
 nu0 
 Reference refractivity. Refractivity at r0
 Unitless
 
 phi
 normalized stellar flux.
 unitless.
 
 phiref
 normalized stellar flux for just refraction.
 unitless.
 
 r
 planet radius (distance from planet center).
 cm
 
 r0
 reference planet radius (distance from planet center).
 cm
 
 rh
 half-light planet radius
 cm
 
 r1
 planet radius of top of haze layer
 cm
 
 r2
 planet radius where line-of-sight optical depth = 1
 cm
 
 rtau_ref
 planet radius for reference tauobs, tauobs_ref
 (normally tauobs_ref = 1, whence rtau_ref = r2)
 
 rsurf
 planet radius of surface.
 flux is cut off (no diffraction) at r < rsurf
 cm
 
 rho
 shadow radius (distance from shadow center), negative for
 far-limb contribution.
 cm.
 
 rhoh
 half-light shadow radius
 cm
 
 rho_min
 minimum observer shadow radius
 cm
 
 t
 temperature
 K
 
 t0
 temperature at reference radius, r0
 K
 
 theta
 bending angle, negative for an isothermal atmosphere.
 radians.
 
 tauobs
 line-of-sight optical depth
 unitless.
 
 tauobs_ref
 value of tauobs at planet radius rtau_ref
 
 v
 shadow velocity
 cm/s
 
 vperph
 perpendicular shadow velocity
 cm/s
 
 
 ---------------------------------------------------------------------------
 FUNCTIONS
 ---------------------------------------------------------------------------
 
 aser   = oclc_ey92_aseries(a, b, delta)
 bser   = oclc_ey92_bseries(a, b, delta)
 cser   = oclc_ey92_cseries(delta_tau)
 dtheta = oclc_ey92_dtheta(r0, theta0,lam0,a,b, r)
 phi    = oclc_ey92_phi_of_r(r0,nu0,lam0,a,b,r1,kappa1,htau1,distobs,rsurf, r)
 lam    = oclc_ey92_lam(r0, lam0, a, b, r)
 nu     = oclc_ey92_nu(r0, nu0,lam0,a,b, r)
 rho    = oclc_ey92_r_of_rho(r0, nu0,lam0,a,b, distobs, r)
 rho    = oclc_ey92_rho_of_r(r0, nu0,lam0,a,b, distobs, r)
 tauobs = oclc_ey92_tauobs(r1, r2, htau1, r)
 theta  = oclc_ey92_theta(r0,nu0,lam0,a,b, r)

(See oclc_ey92_AAREADME.pro)

OCLC_EY92_ASERIES

 [Previous Routine] [Next Routine] [List of Routines] 
 NAME:
  oclc_ey92_aseries
 PURPOSE: (one line)
  return power series in 1/lam used by EY92
 DESCRIPTION:
  return power series in 1/lam used by EY92
 CATEGORY:
  Occultation lightcurve (oclc)
 CALLING SEQUENCE:
  aser = oclc_ey92_aseries(a,b,order,delta)
 INPUTS:
  a - exponent for non-constant mean moleculr weight
      mu(r) = mu0 * (r/r0)^-a
  b - exponent for non-constant temperature
      T(r) = T * (r/r0)^b
  order - order of delta in the series (default 4)
  delta - 1/lambda (may be an array)
 OPTIONAL INPUT PARAMETERS:
 KEYWORD INPUT PARAMETERS:
  none
 KEYWORD OUTPUT PARAMETERS:
  none
 OUTPUTS:
  A(delta,a,b) 
  Series in delta for use in calculating theta
  Unitless.
 COMMON BLOCKS:
  None
 SIDE EFFECTS:
 RESTRICTIONS:
  None
 PROCEDURE:
  Elliot and Young 1992, AJ 103, pp 991; eg A2
 MODIFICATION HISTORY:
  Written 2006 Jun 28, Leslie Young

(See oclc_ey92_aseries.pro)

OCLC_EY92_ATMPARAM

 [Previous Routine] [Next Routine] [List of Routines] 
 NAME:
  oclc_ey92_atmparam
 PURPOSE: (one line)
  convert "data parameters" to "atmospheric parameters"
 DESCRIPTION:
  convert "data parameters" to "atmospheric parameters"
 CATEGORY:
  Occultation lightcurve (oclc)
 CALLING SEQUENCE:
  oclc_ey92_atmparam, v, distobs, f, r0, order,
  timeim, timeem, timehiso, a, b,timeh1, timeh2, timehtau2, timemin, 
  timemid, rhomin, nu0, lam0, r1, kappa1, htau1, rhof = rhof, r2 = r2
 INPUTS:
  v - sky-plane velocity, cm/s
  distobs - Observer-target distance, cm
  f - flux level for timeim, timeem (unitless)
  r0 - reference radius (cm)
  order - expansion order for power series
  timeim - time of immersion (s)
  timeem - time of emersion (s)
  timehiso - equivalent isothermal scale height in units of s.  
   timehiso ~= 1/(8 dphi/dt) at halflight)
  a - exponent for molecular weight
  b - exponent for temperature
  timeh1 - time between tim and haze onset (s)
  timeh2 - time between tim and optical-depth = 1
  timehtau2 - scale height of haze at r2 in units of s
  timemin - minimum shadow radius in s (rhomin/v)
 OPTIONAL INPUT PARAMETERS:
  none
 KEYWORD INPUT PARAMETERS:
  none
 KEYWORD OUTPUT PARAMETERS:
  rh - radius at half-light (phi = f)
  r2 - radius of tau = 1
 OUTPUTS:
  timemid - midtiime (s)
  rhomin - minimum rho (cm)
  nu0 - reference refractivity (unitless)
  lam0 - reference energy ratio (unitless)
  r1 - radius of haze onset (cm)
  kappa1 - linear absorption at r1
  Htau1 - scale height of absorption at r1
 COMMON BLOCKS:
  None
 SIDE EFFECTS:
 RESTRICTIONS:
  None
 PROCEDURE:
  Elliot and Young 1992, AJ 103, pp 991; eg 3.18
 MODIFICATION HISTORY:
  Written 2006 Jun 28, Leslie Young

(See oclc_ey92_atmparam.pro)

OCLC_EY92_BSERIES

 [Previous Routine] [Next Routine] [List of Routines] 
 NAME:
  oclc_ey92_bseries
 PURPOSE: (one line)
  return power series in 1/lam used by EY92
 DESCRIPTION:
  return power series in 1/lam used by EY92
 CATEGORY:
  Occultation lightcurve (oclc)
 CALLING SEQUENCE:
  bser = oclc_ey92_bseries(a,b,order,delta)
 INPUTS:
  a - exponent for non-constant mean moleculr weight
      mu(r) = mu0 * (r/r0)^-a
  b - exponent for non-constant temperature
      T(r) = T * (r/r0)^b
  order - order of delta in the series (default 4)
  delta - 1/lambda (may be an array)
 OPTIONAL INPUT PARAMETERS:
 KEYWORD INPUT PARAMETERS:
  none
 KEYWORD OUTPUT PARAMETERS:
  none
 OUTPUTS:
  B(a,b,delta)
  Series in delta for use in calculating dtheta
  Unitless.
 COMMON BLOCKS:
  None
 SIDE EFFECTS:
 RESTRICTIONS:
  None
 PROCEDURE:
  Elliot and Young 1992, AJ 103, pp 991; eg A2
 MODIFICATION HISTORY:
  Written 2006 Jun 28, Leslie Young

(See oclc_ey92_bseries.pro)

OCLC_EY92_CSERIES

 [Previous Routine] [Next Routine] [List of Routines] 
 NAME:
  oclc_ey92_cseries
 PURPOSE: (one line)
  return power series in 1/lam used by EY92
 DESCRIPTION:
  return power series in 1/lam used by EY92
 CATEGORY:
  Occultation lightcurve (oclc)
 CALLING SEQUENCE:
  cser = oclc_ey92_cseries(order,delta_tau)
 INPUTS:
  order - order for series
  delta_tau - H_tau1 r/r_1^2
 OPTIONAL INPUT PARAMETERS:
  none
 KEYWORD INPUT PARAMETERS:
  none
 KEYWORD OUTPUT PARAMETERS:
  none
 OUTPUTS:
  C(delta_tau)
  Series in delta_tau for use in calculating optical depth
  Unitless.
 COMMON BLOCKS:
  None
 SIDE EFFECTS:
 RESTRICTIONS:
  None
 PROCEDURE:
  Elliot and Young 1992, AJ 103, pp 991; eg A7
 MODIFICATION HISTORY:
  Written 2006 Jun 28, Leslie Young

(See oclc_ey92_cseries.pro)

OCLC_EY92_D2T

 [Previous Routine] [Next Routine] [List of Routines] 
 NAME:
  oclc_ey92_d2t
 PURPOSE: (one line)
  return temperature second derivative as function of radius
 DESCRIPTION:
  return temperature second derivative (d^2t/dr^2) as function of radius
 CATEGORY:
  Occultation lightcurve (oclc)
 CALLING SEQUENCE:
  d2t = oclc_ey92_d2t(r0,lam0, b, Mp, mu0, r)
 INPUTS:
  r0 - reference radius (scalar, cm)
  lam0 - reference energy ratio (radius/scale height) (scalar, unitless)
  b - exponent for non-constant temperature  T(r) = T * (r/r0)^b
  Mp - planet mass (g). 1.30497e+25 for Pluto (Buie et al 2006)
  mu0 - reference molecular weight (unitless). 28.01 for Pluto.
 OPTIONAL INPUT PARAMETERS:
  none
 KEYWORD INPUT PARAMETERS:
  none
 KEYWORD OUTPUT PARAMETERS:
  none
 OUTPUTS:
  d2t - second derivative of temperature (K/cm^2) at r
 COMMON BLOCKS:
  None
 SIDE EFFECTS:
 RESTRICTIONS:
  None
 PROCEDURE:
  Elliot and Young 1992, AJ 103, pp 991; eg 3.3
 MODIFICATION HISTORY:
  Written 2006 Jun 28, Leslie Young

(See oclc_ey92_d2t.pro)

OCLC_EY92_DNU

 [Previous Routine] [Next Routine] [List of Routines] 
 NAME:
  oclc_ey92_dnu
 PURPOSE: (one line)
  return dnu (derivative of refractivity) as a function of radius
 DESCRIPTION:
  return dnu as a function of radius
 CATEGORY:
  Occultation lightcurve (oclc)
 CALLING SEQUENCE:
  dnu = oclc_ey92_nu(r0,nu0,lam0,a,b, r)
 INPUTS:
  r0 - reference radius (scalar, cm)
  nu0 - refractivity at r0 (scalar, unitless)
  lam0 - reference energy ratio (radius/scale height) (scalar, unitless)
  a - exponent for non-constant mean moleculr weight
      mu(r) = mu0 * (r/r0)^-a
  b - exponent for non-constant temperature
      T(r) = T * (r/r0)^b
  r - planet radius at which to evaluate nu (scalar or array, cm)
 OPTIONAL INPUT PARAMETERS:
  none
 KEYWORD INPUT PARAMETERS:
  none
 KEYWORD OUTPUT PARAMETERS:
  none
 OUTPUTS:
  nu (refractivity, unitless) at r (planet radius, cm)
 COMMON BLOCKS:
  None
 SIDE EFFECTS:
 RESTRICTIONS:
  None
 PROCEDURE:
  Elliot and Young 1992, AJ 103, pp 991; eg 3.22
 MODIFICATION HISTORY:
  Written 2006 Jun 28, Leslie Young

(See oclc_ey92_dnu.pro)

OCLC_EY92_DT

 [Previous Routine] [Next Routine] [List of Routines] 
 NAME:
  oclc_ey92_dt
 PURPOSE: (one line)
  return temperature gradient as function of radius
 DESCRIPTION:
  return temperature gradient as function of radius
 CATEGORY:
  Occultation lightcurve (oclc)
 CALLING SEQUENCE:
  dt = oclc_ey92_dt(r0,lam0, b, Mp, mu0, r)
 INPUTS:
  r0 - reference radius (scalar, cm)
  lam0 - reference energy ratio (radius/scale height) (scalar, unitless)
  b - exponent for non-constant temperature  T(r) = T * (r/r0)^b
  Mp - planet mass (g). 1.30497e+25 for Pluto (Buie et al 2006)
  mu0 - reference molecular weight (unitless). 28.01 for Pluto.
 OPTIONAL INPUT PARAMETERS:
  none
 KEYWORD INPUT PARAMETERS:
  none
 KEYWORD OUTPUT PARAMETERS:
  none
 OUTPUTS:
  dt - temperature gradient (K/cm) at r
 COMMON BLOCKS:
  None
 SIDE EFFECTS:
 RESTRICTIONS:
  None
 PROCEDURE:
  Elliot and Young 1992, AJ 103, pp 991; eg 3.3
 MODIFICATION HISTORY:
  Written 2006 Jun 28, Leslie Young

(See oclc_ey92_dt.pro)

OCLC_EY92_DTHETA

 [Previous Routine] [Next Routine] [List of Routines] 
 NAME:
  oclc_ey92_dtheta
 PURPOSE: (one line)
  return dtheta (bending angle) as a function of radius
 DESCRIPTION:
  return dtheta as a function of radius
 CATEGORY:
  Occultation lightcurve (oclc)
 CALLING SEQUENCE:
  dtheta = oclc_ey92_dtheta(r0,nu0,lam0,a,b,order, r)
 INPUTS:
  r0 - reference radius (scalar, cm)
  nu0 - refractivity at r0 (scalar, unitless)
  lam0 - reference energy ratio (radius/scale height) (scalar, unitless)
  a - exponent for non-constant mean moleculr weight
      mu(r) = mu0 * (r/r0)^-a
  b - exponent for non-constant temperature
      T(r) = T * (r/r0)^b
  order - order for series expantion
  r - radius at which to evaluat dtheta (scalar or array, cm)
 OPTIONAL INPUT PARAMETERS:
  none
 KEYWORD INPUT PARAMETERS:
  none
 KEYWORD OUTPUT PARAMETERS:
  none
 OUTPUTS:
  dtheta at r
  Derivative of the bending angle with respect to planet radius,
  positive for an isothermal atmosphere.
  radian/cm
 COMMON BLOCKS:
  None
 SIDE EFFECTS:
 RESTRICTIONS:
  None
 PROCEDURE:
  Elliot and Young 1992, AJ 103, pp 991; eg 4.11
 MODIFICATION HISTORY:
  Written 2006 Jun 28, Leslie Young

(See oclc_ey92_dtheta.pro)

OCLC_EY92_FUNC_DATAPARAM

 [Previous Routine] [Next Routine] [List of Routines] 
 NAME:
  oclc_ey92_func_dataparam
 PURPOSE: (one line)
  lightcurve fit of EY92 to flux(time)
 DESCRIPTION:
  lightcurve fit of EY92 to flux(time)
 CATEGORY:
  Occultation lightcurve (oclc)
 CALLING SEQUENCE:
  oclc_ey92_funcdataparam, y, param, z, pder
 INPUTS:
  y = shadow radius (km)
  param = [sb, dsb, sf, $
     timeim, timeem, timehiso, a, b,timeh1, timeh2, timehtau2, timemin]
 OPTIONAL INPUT PARAMETERS:
  none
 KEYWORD INPUT PARAMETERS:
  none
 KEYWORD OUTPUT PARAMETERS:
  none
 OUTPUTS:
  z = lightcurve
 COMMON BLOCKS:
  None
 SIDE EFFECTS:
 RESTRICTIONS:
  None
 PROCEDURE:
  Elliot and Young 1992, AJ 103, pp 991;
 MODIFICATION HISTORY:
  Written 2006 Jun 28, Leslie Young

(See oclc_ey92_funcdataparam.pro)

OCLC_EY92_H

 [Previous Routine] [Next Routine] [List of Routines] 
 NAME:
  oclc_ey92_h
 PURPOSE: (one line)
  return scale height (H) as a function of radius
 DESCRIPTION:
  return scale height (H) as a function of radius
 CATEGORY:
  Occultation lightcurve (oclc)
 CALLING SEQUENCE:
  h = oclc_ey92_h(r0,lam0, a, b, r)
 INPUTS:
  r0 - reference radius (scalar, cm)
  lam0 - reference energy ratio (radius/scale height) (scalar, unitless)
  a - exponent for non-constant mean moleculr weight
      mu(r) = mu0 * (r/r0)^-a
  b - exponent for non-constant temperature
      T(r) = T * (r/r0)^b
  r - radius at which to evaluat H (scalar or array, cm)
 OPTIONAL INPUT PARAMETERS:
  none
 KEYWORD INPUT PARAMETERS:
  none
 KEYWORD OUTPUT PARAMETERS:
  none
 OUTPUTS:
  H = scale height (cm), defined by
  H = k T / (mu m_amu g) or
  - (1/p) dp/dr = 1/H
 COMMON BLOCKS:
  None
 SIDE EFFECTS:
 RESTRICTIONS:
  None
 PROCEDURE:
  Elliot and Young 1992, AJ 103, pp 991; eg 3.21
 MODIFICATION HISTORY:
  Written 2006 Jun 28, Leslie Young

(See oclc_ey92_h.pro)

OCLC_EY92_HN

 [Previous Routine] [Next Routine] [List of Routines] 
 NAME:
  oclc_ey92_hn
 PURPOSE: (one line)
  return scale height for density (Hn) as a function of radius
 DESCRIPTION:
  return scale height for density (Hn) as a function of radius
 CATEGORY:
  Occultation lightcurve (oclc)
 CALLING SEQUENCE:
  hn = oclc_ey92_hn(r0,lam0, a, b, r)
 INPUTS:
  r0 - reference radius (scalar, cm)
  lam0 - reference energy ratio (radius/scale height) (scalar, unitless)
  a - exponent for non-constant mean moleculr weight
      mu(r) = mu0 * (r/r0)^-a
  b - exponent for non-constant temperature
      T(r) = T * (r/r0)^b
  r - radius at which to evaluat H (scalar or array, cm)
 OPTIONAL INPUT PARAMETERS:
  none
 KEYWORD INPUT PARAMETERS:
  none
 KEYWORD OUTPUT PARAMETERS:
  none
 OUTPUTS:
  Hn = scale height (cm), defined by
  H = - (1/n(r) d n(r) / dr)
 COMMON BLOCKS:
  None
 SIDE EFFECTS:
 RESTRICTIONS:
  None
 PROCEDURE:
  Elliot and Young 1992, AJ 103, pp 991; eg 3.22
 MODIFICATION HISTORY:
  Written 2006 Jun 28, Leslie Young

(See oclc_ey92_hn.pro)

OCLC_EY92_KAPPA

 [Previous Routine] [Next Routine] [List of Routines] 
 NAME:
  oclc_ey92_kappa
 PURPOSE: (one line)
  return kappa as a function of radius
 DESCRIPTION:
  return kappa (absorption coeffient, cm^-1) as a function of radius
 CATEGORY:
  Occultation lightcurve (oclc)
 CALLING SEQUENCE:
  kappa = oclc_ey92_kappa(r1, kappa1, htau1, r)
 INPUTS:
  r1 - top of haze layer
  kappa1 - linear absorption coefficient at r1
  htau1 - scale height of absorbtion coefficient
  r - radius at which to evaluat kappa (scalar or array, cm)
 OPTIONAL INPUT PARAMETERS:
  none
 KEYWORD INPUT PARAMETERS:
  none
 KEYWORD OUTPUT PARAMETERS:
  none
 OUTPUTS:
  kappa (absorption coeffient, cm^-1) at r (planet radius, cm)
 COMMON BLOCKS:
  None
 SIDE EFFECTS:
 RESTRICTIONS:
  None
 PROCEDURE:
  Elliot and Young 1992, AJ 103, pp 991; eg 3.23
 MODIFICATION HISTORY:
  Written 2006 Aug 17, Leslie Young

(See oclc_ey92_kappa.pro)

OCLC_EY92_KAPPA1_OF_R2

 [Previous Routine] [Next Routine] [List of Routines] 
 NAME:
  oclc_ey92_kappa1_of_r2
 PURPOSE: (one line)
  return kappa1 (linear absorption coefficient at r1) s.t. tauobs(r2) = 1
 DESCRIPTION:
  return kappa1 (linear absorption coefficient at r1) s.t. tauobs(r2) = 1
 CATEGORY:
  Occultation lightcurve (oclc)
 CALLING SEQUENCE:
  kappa1 = oclc_ey92_kappa1_of_r2(r1, htau1, r2)
 INPUTS:
  r1 - top of haze layer
  htau1 - scale height of absorbtion coefficient
  r2 - planet radius of unit optical depth (scalar or array, cm)
 OPTIONAL INPUT PARAMETERS:
  none
 KEYWORD INPUT PARAMETERS:
  none
 KEYWORD OUTPUT PARAMETERS:
  none
 OUTPUTS:
  kappa1 - linear absorption coefficient at r1
 COMMON BLOCKS:
  None
 SIDE EFFECTS:
 RESTRICTIONS:
  None
 PROCEDURE:
  Elliot and Young 1992, AJ 103, pp 991; eg 4.20, 4.21
 MODIFICATION HISTORY:
  Written 2006 Jun 30, Leslie Young
  2006 Sep 14 LAY fixed documentation

(See oclc_ey92_kappa1_of_r2.pro)

OCLC_EY92_KAPPA1_OF_TAUOBS

 [Previous Routine] [Next Routine] [List of Routines] 
 NAME:
  oclc_ey92_kappa1_of_tauobs
 PURPOSE: (one line)
  return kappa1 (linear absorption coefficient at r1) s.t. tauobs(rtau_ref) = tauobs
 DESCRIPTION:
  return kappa1 (linear absorption coefficient at r1) s.t. tauobs(rtau_ref) = tauobs
 CATEGORY:
  Occultation lightcurve (oclc)
 CALLING SEQUENCE:
  kappa1 = oclc_ey92_kappa1_of_tauobs(r1, htau1, rtau_ref, tauobs_ref)
 INPUTS:
  r1 - top of haze layer
  htau1 - scale height of absorbtion coefficient
  rtau_ref - radius at which tauobs = tauobs_ref (scalar or array, cm)
  tauobs_ref (line-of-sight optical depth, unitless) at rtau_ref
 OPTIONAL INPUT PARAMETERS:
  none
 KEYWORD INPUT PARAMETERS:
  none
 KEYWORD OUTPUT PARAMETERS:
  none
 OUTPUTS:
  kappa1 - linear absorption coefficient at r1
 COMMON BLOCKS:
  None
 SIDE EFFECTS:
 RESTRICTIONS:
  None
 PROCEDURE:
  Elliot and Young 1992, AJ 103, pp 991; eg 4.20, 4.21
 MODIFICATION HISTORY:
  Written 2006 Jun 30, Leslie Young

(See oclc_ey92_kappa1_of_tauobs.pro)

OCLC_EY92_LAM

 [Previous Routine] [Next Routine] [List of Routines] 
 NAME:
  oclc_ey92_lam
 PURPOSE: (one line)
  return lambda (r/H) as a function of radius
 DESCRIPTION:
  return lambda (r/H) as a function of radius
 CATEGORY:
  Occultation lightcurve (oclc)
 CALLING SEQUENCE:
  lam = oclc_ey92_lam(r0,lam0, a, b, r)
 INPUTS:
  r0 - reference radius (scalar, cm)
  lam0 - reference energy ratio (radius/scale height) (scalar, unitless)
  a - exponent for non-constant mean moleculr weight
      mu(r) = mu0 * (r/r0)^-a
  b - exponent for non-constant temperature
      T(r) = T * (r/r0)^b
  r - radius at which to evaluat lam (scalar or array, cm)
 OPTIONAL INPUT PARAMETERS:
  none
 KEYWORD INPUT PARAMETERS:
  none
 KEYWORD OUTPUT PARAMETERS:
  none
 OUTPUTS:
  lambda = r/H = (gravitational potential energy/kinetic energy)
  Ratio of radius (r) to scale height (H).  Since H = kT r^2/(mu m_adu M G),
  this is the same as the ratio of kinetic energy (kT) to potential
  energy (mu m_adu M G/r).
  unitless
 COMMON BLOCKS:
  None
 SIDE EFFECTS:
 RESTRICTIONS:
  None
 PROCEDURE:
  Elliot and Young 1992, AJ 103, pp 991; eg 3.9
 MODIFICATION HISTORY:
  Written 2006 Jun 28, Leslie Young

(See oclc_ey92_lam.pro)

OCLC_EY92_N

 [Previous Routine] [Next Routine] [List of Routines] 
 NAME:
  oclc_ey92_n
 PURPOSE: (one line)
  return n (number density) as a function of radius
 DESCRIPTION:
  return n as a function of radius
 CATEGORY:
  Occultation lightcurve (oclc)
 CALLING SEQUENCE:
  n = oclc_ey92_n(r0,nu0,lam0,a,b,nustp, r)
 INPUTS:
  r0 - reference radius (scalar, cm)
  nu0 - refractivity at r0 (scalar, unitless)
  lam0 - reference energy ratio (radius/scale height) (scalar, unitless)
  a - exponent for non-constant mean moleculr weight
      mu(r) = mu0 * (r/r0)^-a
  b - exponent for non-constant temperature
      T(r) = T * (r/r0)^b
  r - planet radius at which to evaluate nu (scalar or array, cm)
 OPTIONAL INPUT PARAMETERS:
  none
 KEYWORD INPUT PARAMETERS:
  none
 KEYWORD OUTPUT PARAMETERS:
  none
 OUTPUTS:
  nu (refractivity, unitless) at r (planet radius, cm)
 COMMON BLOCKS:
  None
 SIDE EFFECTS:
 RESTRICTIONS:
  None
 PROCEDURE:
  Elliot and Young 1992, AJ 103, pp 991; eg 3.20, 3.18
 MODIFICATION HISTORY:
  Written 2006 Jun 28, Leslie Young

(See oclc_ey92_n.pro)

OCLC_EY92_N0

 [Previous Routine] [Next Routine] [List of Routines] 
 NAME:
  oclc_ey92_n0
 PURPOSE: (one line)
  return reference number density
 DESCRIPTION:
  return reference number density
 CATEGORY:
  Occultation lightcurve (oclc)
 CALLING SEQUENCE:
  n0 = oclc_ey92_n0(nu0, nustp)
 INPUTS:
  nu0 - reference refractivity (scalar, cm)
  nustp - refractivity at STP (273.15 K, 1.013250e6 microsbar). 2.98e-4 for Pluto.
 OPTIONAL INPUT PARAMETERS:
  none
 KEYWORD INPUT PARAMETERS:
  none
 KEYWORD OUTPUT PARAMETERS:
  none
 OUTPUTS:
  n0 - number density (cm^-3) at r0
 COMMON BLOCKS:
  None
 SIDE EFFECTS:
 RESTRICTIONS:
  None
 PROCEDURE:
  Elliot and Young 1992, AJ 103, pp 991; eg 3.18
 MODIFICATION HISTORY:
  Written 2006 Jun 28, Leslie Young

(See oclc_ey92_n0.pro)

OCLC_EY92_NU

 [Previous Routine] [Next Routine] [List of Routines] 
 NAME:
  oclc_ey92_nu
 PURPOSE: (one line)
  return nu (refractivity) as a function of radius
 DESCRIPTION:
  return nu as a function of radius
 CATEGORY:
  Occultation lightcurve (oclc)
 CALLING SEQUENCE:
  nu = oclc_ey92_nu(r0,nu0,lam0,a,b, r)
 INPUTS:
  r0 - reference radius (scalar, cm)
  nu0 - refractivity at r0 (scalar, unitless)
  lam0 - reference energy ratio (radius/scale height) (scalar, unitless)
  a - exponent for non-constant mean moleculr weight
      mu(r) = mu0 * (r/r0)^-a
  b - exponent for non-constant temperature
      T(r) = T * (r/r0)^b
  r - planet radius at which to evaluate nu (scalar or array, cm)
 OPTIONAL INPUT PARAMETERS:
  none
 KEYWORD INPUT PARAMETERS:
  none
 KEYWORD OUTPUT PARAMETERS:
  none
 OUTPUTS:
  nu (refractivity, unitless) at r (planet radius, cm)
 COMMON BLOCKS:
  None
 SIDE EFFECTS:
 RESTRICTIONS:
  None
 PROCEDURE:
  Elliot and Young 1992, AJ 103, pp 991; eg 3.20
 MODIFICATION HISTORY:
  Written 2006 Jun 28, Leslie Young

(See oclc_ey92_nu.pro)

OCLC_EY92_P

 [Previous Routine] [Next Routine] [List of Routines] 
 NAME:
  oclc_ey92_p
 PURPOSE: (one line)
  return pressure as a function of radius
 DESCRIPTION:
  return pressure as a function of radius
 CATEGORY:
  Occultation lightcurve (oclc)
 CALLING SEQUENCE:
  p = oclc_ey92_p(r0, nu0, lam0, a, b, mp, mu0, nustp, r)
 INPUTS:
  r0 - reference radius (scalar, cm)
  nu0 - refractivity at r0 (scalar, unitless)
  lam0 - reference energy ratio (radius/scale height) (scalar, unitless)
  a - exponent for non-constant mean moleculr weight: mu(r) = mu0 * (r/r0)^-a
  b - exponent for non-constant temperature; T(r) = T * (r/r0)^b
  Mp - planet mass (g). 1.30497e+25 for Pluto (Buie et al 2006)
  mu0 - molecular weight (unitless). 28.01 for Pluto.
  nustp - refractivity at STP (273.15 K, 1.013250e6 microsbar). 2.98e-4 for Pluto.
 OPTIONAL INPUT PARAMETERS:
  none
 KEYWORD INPUT PARAMETERS:
  none
 KEYWORD OUTPUT PARAMETERS:
  none
 OUTPUTS:
  p0 - pressure (microbar) at r0
 COMMON BLOCKS:
  None
 SIDE EFFECTS:
 RESTRICTIONS:
  None
 PROCEDURE:
  Elliot and Young 1992, AJ 103, pp 991; eg 3.6
 MODIFICATION HISTORY:
  Written 2006 Jun 28, Leslie Young

(See oclc_ey92_p.pro)

OCLC_EY92_P0

 [Previous Routine] [Next Routine] [List of Routines] 
 NAME:
  oclc_ey92_p0
 PURPOSE: (one line)
  return reference pressure
 DESCRIPTION:
  return reference pressure
 CATEGORY:
  Occultation lightcurve (oclc)
 CALLING SEQUENCE:
  p0 = oclc_ey92_p0(r0, nu0, lam0, mp, mu0, nustp)
 INPUTS:
  r0 - reference radius (scalar, cm)
  nu0 - refractivity at r0 (scalar, unitless)
  lam0 - reference energy ratio (radius/scale height) (scalar,
         unitless)
  Mp - planet mass (g). 1.30497e+25 for Pluto (Buie et al 2006)
  mu0 - molecular weight (unitless). 28.01 for Pluto.
  nustp - refractivity at STP (273.15 K, 1.013250e6 microsbar). 2.98e-4 for Pluto.
 OPTIONAL INPUT PARAMETERS:
  none
 KEYWORD INPUT PARAMETERS:
  none
 KEYWORD OUTPUT PARAMETERS:
  none
 OUTPUTS:
  p0 - pressure (microbar) at r0
 COMMON BLOCKS:
  None
 SIDE EFFECTS:
 RESTRICTIONS:
  None
 PROCEDURE:
  Elliot and Young 1992, AJ 103, pp 991; eg 3.18
 MODIFICATION HISTORY:
  Written 2006 Jun 28, Leslie Young

(See oclc_ey92_p0.pro)

OCLC_EY92_PHIREF_OF_R

 [Previous Routine] [Next Routine] [List of Routines] 
 NAME:
  oclc_ey92_phiref_of_r
 PURPOSE: (one line)
  return phi from refaction as a function of planet radius
 DESCRIPTION:
  return phi from refaction as a function of planet radius
 CATEGORY:
  Occultation lightcurve (oclc)
 CALLING SEQUENCE:
  phi = oclc_ey92_phiref_of_r(r0, nu0, lam0, a, b, distobs, rsurf, order, r)
 INPUTS:
  r0 - reference radius (scalar, cm)
  nu0 - refractivity at r0 (scalar, unitless)
  lam0 - reference energy ratio (radius/scale height) (scalar,
         unitless)
  a - exponent for non-constant mean moleculr weight
      mu(r) = mu0 * (r/r0)^-a
  b - exponent for non-constant temperature
      T(r) = T * (r/r0)^b
  rsurf - surface radius
  order - order for series expansion
  r - radius at which to evaluat phi (scalar or array, cm)
 OPTIONAL INPUT PARAMETERS:
  none
 KEYWORD INPUT PARAMETERS:
  none
 KEYWORD OUTPUT PARAMETERS:
  none
 OUTPUTS:
  stellar flux, phi, at planet radius, r
 COMMON BLOCKS:
  None
 SIDE EFFECTS:
 RESTRICTIONS:
  None
 PROCEDURE:
  Elliot and Young 1992, AJ 103, pp 991; eg 2.1, 2.2, 2.3
 MODIFICATION HISTORY:
  Written 2006 Jun 28, Leslie Young

(See oclc_ey92_phiref_of_r.pro)

OCLC_EY92_PHI_OF_R

 [Previous Routine] [Next Routine] [List of Routines] 
 NAME:
  oclc_ey92_phi_of_r
 PURPOSE: (one line)
  return phi (normalized stellar flux) as a function of planet radius, r
 DESCRIPTION:
  return phi as a function of planet radius
 CATEGORY:
  Occultation lightcurve (oclc)
 CALLING SEQUENCE:
  phi = oclc_ey92_phi_of_r(r0, nu0, lam0, a, b, r1, kappa1, htau1,
  distobs, rsurf, order, r)
 INPUTS:
  r0 - reference radius (scalar, cm)
  nu0 - refractivity at r0 (scalar, unitless)
  lam0 - reference energy ratio (radius/scale height) (scalar,
         unitless)
  a - exponent for non-constant mean moleculr weight
      mu(r) = mu0 * (r/r0)^-a
  b - exponent for non-constant temperature
      T(r) = T * (r/r0)^b
  r1 - top of haze layer
  kappa1 - linear absorption coefficient at r1
  htau1 - scale height of absorbtion coefficient
  rsurf - surface radius
  order - order for series expansion
  r - radius at which to evaluat phi (scalar or array, cm)
 OPTIONAL INPUT PARAMETERS:
  none
 KEYWORD INPUT PARAMETERS:
  none
 KEYWORD OUTPUT PARAMETERS:
  none
 OUTPUTS:
  stellar flux, phi, at planet radius, r
 COMMON BLOCKS:
  None
 SIDE EFFECTS:
 RESTRICTIONS:
  None
 PROCEDURE:
  Elliot and Young 1992, AJ 103, pp 991; eg 2.1, 2.2, 2.3
 MODIFICATION HISTORY:
  Written 2006 Jun 28, Leslie Young

(See oclc_ey92_phi_of_r.pro)

OCLC_EY92_PHI_OF_RHO

 [Previous Routine] [Next Routine] [List of Routines] 
 NAME:
  oclc_ey92_phi_of_rho
 PURPOSE: (one line)
  return phi (normalized stellar flux) as a function of shadow radius
 DESCRIPTION:
  return phi as a function of shadow radius
 CATEGORY:
  Occultation lightcurve (oclc)
 CALLING SEQUENCE:
  phi = oclc_ey92_phi_of_rho(r0,nu0,lam0,a,b,r1,kappa1,htau1,distobs,rsurf,farside, haze, order, rho)
 INPUTS:
  r0 - reference radius (scalar, cm)
  nu0 - refractivity at r0 (scalar, unitless)
  lam0 - reference energy ratio (radius/scale height) (scalar,
         unitless)
  a - exponent for non-constant mean moleculr weight
      mu(r) = mu0 * (r/r0)^-a
  b - exponent for non-constant temperature
      T(r) = T * (r/r0)^b
  r1 - top of haze layer
  kappa1 - linear absorption coefficient at r1
  htau1 - scale height of absorbtion coefficient
  rsurf - surface radius
  farside - if 1, add farside contribution
  haze - if 1, include haze extinction
  order - order for series expansion
  rho - shadow radius at which to evaluat phi (scalar or array, cm)
  
 OPTIONAL INPUT PARAMETERS:
  none
 KEYWORD INPUT PARAMETERS:
  none
 KEYWORD OUTPUT PARAMETERS:
  r_near - the near-side planet radius (cm) 
 OUTPUTS:
  phi (normalized stellar flux) at shadow radius, rho.
 COMMON BLOCKS:
  None
 SIDE EFFECTS:
 RESTRICTIONS:
  None
 PROCEDURE:
  Elliot and Young 1992, AJ 103, pp 991; eg 2.1, 2.2, 2.3
 MODIFICATION HISTORY:
  Written 2006 Jun 28, Leslie Young
  Sep 9 2007 LAY.  Don't calculate tau if haze = 0

(See oclc_ey92_phi_of_rho.pro)

OCLC_EY92_RHOMIN_OF_TIMEMIN

 [Previous Routine] [Next Routine] [List of Routines] 
 NAME:
  oclc_ey92_rhomin_of_timemin
 PURPOSE: (one line)
  return rho_min = minimim shadow radius
 DESCRIPTION:
  return rho_min = minimim shadow radius
 CATEGORY:
  Occultation lightcurve (oclc)
 CALLING SEQUENCE:
  rho_min = oclc_ey92_rhomin_of_timemin(v, timemin)
 INPUTS:
  v = sky-plane velocity (assuming linear motion), cm/s
  timemin = minimum observer distance in units of s.
 OPTIONAL INPUT PARAMETERS:
  none
 KEYWORD INPUT PARAMETERS:
  none
 KEYWORD OUTPUT PARAMETERS:
 OUTPUTS:
  rhomin = minimum observer distance in units of cm
 COMMON BLOCKS:
  None
 SIDE EFFECTS:
 RESTRICTIONS:
  Admitedly silly to have this as a seperate function
 PROCEDURE:
  Elliot and Young 1992, AJ 103, pp 991; eq 5.8
 MODIFICATION HISTORY:
  Written 2006 Aug 18, Leslie Young

(See oclc_ey92_rhomin_of_timemin.pro)

OCLC_EY92_RHO_OF_R

 [Previous Routine] [Next Routine] [List of Routines] 
 NAME:
  oclc_ey92_rho_of_r
 PURPOSE: (one line)
  return r (planet radius) as a function of rho (shadow radius)
 DESCRIPTION:
  return r as a function of rho
 CATEGORY:
  Occultation lightcurve (oclc)
 CALLING SEQUENCE:
  rho = oclc_ey92_rho_of_r(r0,nu0,lam0,a,b, distobs, order, r)
 INPUTS:
  r0 - reference radius (scalar, cm)
  nu0 - refractivity at r0 (scalar, unitless)
  lam0 - reference energy ratio (radius/scale height) (scalar, unitless)
  a - exponent for non-constant mean moleculr weight
      mu(r) = mu0 * (r/r0)^-a
  b - exponent for non-constant temperature
      T(r) = T * (r/r0)^b
  distobs - planet-observer distance (scalar, cm)
  order - order of expansion for series
  r - planet radius (cm)
 OPTIONAL INPUT PARAMETERS:
  none
 KEYWORD INPUT PARAMETERS:
  none
 KEYWORD OUTPUT PARAMETERS:
  none
 OUTPUTS:
  rho at r
  rho - shadow radius at  r (scalar or array, cm),
        negative for far-limb contribution; COMMON BLOCKS:
 SIDE EFFECTS:
 RESTRICTIONS:
  None
 PROCEDURE:
  Elliot and Young 1992, AJ 103, pp 991; eg 2.2
 MODIFICATION HISTORY:
  Written 2006 Jun 28, Leslie Young

(See oclc_ey92_rho_of_r.pro)

OCLC_EY92_RHO_OF_TIME

 [Previous Routine] [Next Routine] [List of Routines] 
 NAME:
  oclc_ey92_rho_of_time
 PURPOSE: (one line)
  return rho (shadow radius) as a function of time assuming linear motion
 DESCRIPTION:
  return rho as a function of time
 CATEGORY:
  Occultation lightcurve (oclc)
 CALLING SEQUENCE:
  rho = oclc_ey92_rho_of_time(rhomin, timemid, v, time)
 INPUTS:
  rhomin - shadow radius at midtime (cm)
  timemid - mid time (s)
  v - sky-plane velocity, cm/s
  time - time (s) at which to evaluate rho
 OPTIONAL INPUT PARAMETERS:
  none
 KEYWORD INPUT PARAMETERS:
  none
 KEYWORD OUTPUT PARAMETERS:
  none
 OUTPUTS:
  rho - shadow radius 
 COMMON BLOCKS:
  None
 SIDE EFFECTS:
 RESTRICTIONS:
  None
 PROCEDURE:
  Elliot and Young 1992, AJ 103, pp 991; eg 5.1
 MODIFICATION HISTORY:
  Written 2006 Jun 28, Leslie Young

(See oclc_ey92_rho_of_time.pro)

OCLC_EY92_RH_OF_RHOH

 [Previous Routine] [Next Routine] [List of Routines] 
 NAME:
  oclc_ey92_rh_of_rhoh
 PURPOSE: (one line)
  return rh (half-light planet radius) as a function of rhoh (half-light shadow radius)
 DESCRIPTION:
  return rh (half-light planet radius) and lam0 as a 
  function of rho (half-light shadow radius) and hiso (
 CATEGORY:
  Occultation lightcurve (oclc)
 CALLING SEQUENCE:
  rh = oclc_ey92_rh_of_rhoh(r0,a,b,distobs,order,rhoh,hiso,nu0,lam0)
 INPUTS:
  r0 - reference radius (scalar, cm)
  a - exponent for non-constant mean moleculr weight
      mu(r) = mu0 * (r/r0)^-a
  b - exponent for non-constant temperature
      T(r) = T * (r/r0)^b
  distobs - planet-observer distance (scalar, cm)
  order - order for expansion of series
  rhoh - shadow radius at which flux = 1/2.
  hiso - equivalent isothermal scale height at rhoh, defined such
         that d phi/d rho = (1 + Hiso/2 rh + ...)/8
 OPTIONAL INPUT PARAMETERS:
  none
 KEYWORD INPUT PARAMETERS:
  none
 KEYWORD OUTPUT PARAMETERS:
  none
 OUTPUTS:
  rh such that oclc_ey92_rho_of_r(r0,nu0,lam0,a,b,distobs,r) = rhoh
          and oclc_ey92_phiref_of_r(r0,nu0,lam0,a,b,distobs,rsurf,order,r) = 0.5
  rh is radius (distance from planet center) in cm.
  nu0 - refractivity at r0 (scalar, unitless)
  lam0 - reference energy ratio (radius/scale height) (scalar, unitless)
 COMMON BLOCKS:
  None
 SIDE EFFECTS:
 RESTRICTIONS:
  None
 PROCEDURE:
  Elliot and Young 1992, AJ 103, pp 991; eq 6.3-6.7
 MODIFICATION HISTORY:
  Written 2006 Jun 28, Leslie Young

(See oclc_ey92_rh_of_rhoh.pro)

OCLC_EY92_R_OF_RHO

 [Previous Routine] [Next Routine] [List of Routines] 
 NAME:
  oclc_ey92_r_of_rho
 PURPOSE: (one line)
  return r (planet radius) as a function of rho (shadow radius)
 DESCRIPTION:
  return r as a function of rho
 CATEGORY:
  Occultation lightcurve (oclc)
 CALLING SEQUENCE:
  r = oclc_ey92_r_of_rho(r0,nu0,lam0,a,b, distobs, order, rho)
 INPUTS:
  r0 - reference radius (scalar, cm)
  nu0 - refractivity at r0 (scalar, unitless)
  lam0 - reference energy ratio (radius/scale height) (scalar, unitless)
  a - exponent for non-constant mean moleculr weight
      mu(r) = mu0 * (r/r0)^-a
  b - exponent for non-constant temperature
      T(r) = T * (r/r0)^b
  distobs - planet-observer distance (scalar, cm)
  order - order for expansion of series
  rho - shadow radius at which to evaluate r (scalar or array, cm)
 OPTIONAL INPUT PARAMETERS:
  none
 KEYWORD INPUT PARAMETERS:
  none
 KEYWORD OUTPUT PARAMETERS:
  none
 OUTPUTS:
  r such that oclc_ey92_rho_of_r(r0,nu0,lam0,a,b,distobs,r) = rho
  r is radius (distance from planet center) in cm.
 COMMON BLOCKS:
  None
 SIDE EFFECTS:
 RESTRICTIONS:
  None
 PROCEDURE:
  Elliot and Young 1992, AJ 103, pp 991; eg 2.2
 MODIFICATION HISTORY:
  Written 2006 Jun 28, Leslie Young

(See oclc_ey92_r_of_rho.pro)

OCLC_EY92_T

 [Previous Routine] [Next Routine] [List of Routines] 
 NAME:
  oclc_ey92_t
 PURPOSE: (one line)
  return temperature as function of radius
 DESCRIPTION:
  return temperature as function of radius
 CATEGORY:
  Occultation lightcurve (oclc)
 CALLING SEQUENCE:
  t = oclc_ey92_t(r0,lam0, b, Mp, mu0, r)
 INPUTS:
  r0 - reference radius (scalar, cm)
  lam0 - reference energy ratio (radius/scale height) (scalar, unitless)
  b - exponent for non-constant temperature  T(r) = T * (r/r0)^b
  Mp - planet mass (g). 1.30497e+25 for Pluto (Buie et al 2006)
  mu0 - reference molecular weight (unitless). 28.01 for Pluto.
 OPTIONAL INPUT PARAMETERS:
  none
 KEYWORD INPUT PARAMETERS:
  none
 KEYWORD OUTPUT PARAMETERS:
  none
 OUTPUTS:
  t - temperature (K) at r
 COMMON BLOCKS:
  None
 SIDE EFFECTS:
 RESTRICTIONS:
  None
 PROCEDURE:
  Elliot and Young 1992, AJ 103, pp 991; eg 3.3
 MODIFICATION HISTORY:
  Written 2006 Jun 28, Leslie Young

(See oclc_ey92_t.pro)

OCLC_EY92_T0

 [Previous Routine] [Next Routine] [List of Routines] 
 NAME:
  oclc_ey92_t0
 PURPOSE: (one line)
  return reference temperature
 DESCRIPTION:
  return reference temperature
 CATEGORY:
  Occultation lightcurve (oclc)
 CALLING SEQUENCE:
  t0 = oclc_ey92_t0(r0,lam0, Mp, mu0)
 INPUTS:
  r0 - reference radius (scalar, cm)
  lam0 - reference energy ratio (radius/scale height) (scalar, unitless)
  Mp - planet mass (g). 1.30497e+25 for Pluto (Buie et al 2006)
  mu0 - molecular weight (unitless). 28.01 for Pluto.
 OPTIONAL INPUT PARAMETERS:
  none
 KEYWORD INPUT PARAMETERS:
  none
 KEYWORD OUTPUT PARAMETERS:
  none
 OUTPUTS:
  t0 - temperature (K) at r0
 COMMON BLOCKS:
  None
 SIDE EFFECTS:
 RESTRICTIONS:
  None
 PROCEDURE:
  Elliot and Young 1992, AJ 103, pp 991; eg 3.10
 MODIFICATION HISTORY:
  Written 2006 Jun 28, Leslie Young

(See oclc_ey92_t0.pro)

OCLC_EY92_TABLE2

 [Previous Routine] [Next Routine] [List of Routines] 
 NAME:
  oclc_ey92_table2
 PURPOSE: (one line)
  reproduce table 2
 DESCRIPTION:
  reproduce table 2
 CATEGORY:
  Occultation lightcurve (oclc)
 CALLING SEQUENCE:
  oclc_ey92_table2, time, rho, r, lam, nu, kappa, tauobs, phi_cyl,
  phi_ref, phi
 INPUTS:
  none
 OPTIONAL INPUT PARAMETERS:
  none
 KEYWORD INPUT PARAMETERS:
  none
 KEYWORD OUTPUT PARAMETERS:
  none
 OUTPUTS:
  time - time
 COMMON BLOCKS:
  None
 SIDE EFFECTS:
 RESTRICTIONS:
  None
 PROCEDURE:
  Elliot and Young 1992, AJ 103, pp 991
 MODIFICATION HISTORY:
  Written 2006 Sep 15, Leslie Young

(See oclc_ey92_table2.pro)

OCLC_EY92_TABLE3A

 [Previous Routine] [Next Routine] [List of Routines] 
 NAME:
  oclc_ey92_table3a
 PURPOSE: (one line)
  reproduce table 3a
 DESCRIPTION:
  reproduce table 3a
 CATEGORY:
  Occultation lightcurve (oclc)
 CALLING SEQUENCE:
  oclc_ey92_table3a, time, rho, r, lam, nu, kappa, tauobs, phi_cyl,
  phi_ref, phi
 INPUTS:
  none
 OPTIONAL INPUT PARAMETERS:
  none
 KEYWORD INPUT PARAMETERS:
  none
 KEYWORD OUTPUT PARAMETERS:
  none
 OUTPUTS:
  time - time
 COMMON BLOCKS:
  None
 SIDE EFFECTS:
 RESTRICTIONS:
  None
 PROCEDURE:
  Elliot and Young 1992, AJ 103, pp 991; eg 3.18
 MODIFICATION HISTORY:
  Written 2006 Jun 28, Leslie Young

(See oclc_ey92_table3a.pro)

OCLC_EY92_TABLE7

 [Previous Routine] [Next Routine] [List of Routines] 
 NAME:
  oclc_ey92_table7
 PURPOSE: (one line)
  reproduce table 7
 DESCRIPTION:
  reproduce table 7 from Elliot and Young, 1992
 CATEGORY:
  Occultation lightcurve (oclc)
 CALLING SEQUENCE:
  oclc_ey92_table7, kao=kao
 INPUTS:
  none
 OPTIONAL INPUT PARAMETERS:
  none
 KEYWORD INPUT PARAMETERS:
  kao : optional filename for kao.utts
 KEYWORD OUTPUT PARAMETERS:
  none
 OUTPUTS:
  time - time
 COMMON BLOCKS:
  None
 SIDE EFFECTS:
 RESTRICTIONS:
  None
 PROCEDURE:
  Elliot and Young 1992, AJ 103, pp 991
 MODIFICATION HISTORY:
  Written 2006 Sep 15, Leslie Young

(See oclc_ey92_table7.pro)

OCLC_EY92_TAUOBS

 [Previous Routine] [Next Routine] [List of Routines] 
 NAME:
  oclc_ey92_tauobs
 PURPOSE: (one line)
  return tauobs (line-of-sight optical depth) as a function of radius
 DESCRIPTION:
  return tauobs as a function of radius
 CATEGORY:
  Occultation lightcurve (oclc)
 CALLING SEQUENCE:
  tauobs = oclc_ey92_tauobs(r1, kappa1, htau1, order, r)
 INPUTS:
  r1 - top of haze layer
  kappa1 - linear absorption coefficient at r1
  htau1 - scale height of absorbtion coefficient
  order - order for series
  r - radius at which to evaluat tauobs (scalar or array, cm)
 OPTIONAL INPUT PARAMETERS:
  none
 KEYWORD INPUT PARAMETERS:
  none
 KEYWORD OUTPUT PARAMETERS:
  none
 OUTPUTS:
  tauobs (line-of-sight optical depth, unitless) at r (planet radius, cm)
 COMMON BLOCKS:
  None
 SIDE EFFECTS:
 RESTRICTIONS:
  None
 PROCEDURE:
  Elliot and Young 1992, AJ 103, pp 991; eg 4.20
 MODIFICATION HISTORY:
  Written 2006 Jun 28, Leslie Young

(See oclc_ey92_tauobs.pro)

OCLC_EY92_TAUVERT

 [Previous Routine] [Next Routine] [List of Routines] 
 NAME:
  oclc_ey92_tauvert
 PURPOSE: (one line)
  return tauvert (vertical optical depth) as a function of radius
 DESCRIPTION:
  return tauvert as a function of radius
 CATEGORY:
  Occultation lightcurve (oclc)
 CALLING SEQUENCE:
  tauvert = oclc_ey92_tauvert(r1, kappa1, htau1, order, r)
 INPUTS:
  r1 - top of haze layer
  kappa1 - linear absorption coefficient at r1
  htau1 - scale height of absorbtion coefficient
  order - order for series
  r - radius at which to evaluate tauvert (scalar or array, cm)
 OPTIONAL INPUT PARAMETERS:
  none
 KEYWORD INPUT PARAMETERS:
  none
 KEYWORD OUTPUT PARAMETERS:
  none
 OUTPUTS:
  tauvert (optical depth from top to r, unitless) at r (planet radius, cm)
 COMMON BLOCKS:
  None
 SIDE EFFECTS:
 RESTRICTIONS:
  None
 PROCEDURE:

           /r1
  tau(r) = |   kappa(r) dr
           /r

  kappa(r) from Elliot and Young 1992, eq 3.23


 MODIFICATION HISTORY:
  Written 2006 Jun 28, Leslie Young

(See oclc_ey92_tauvert.pro)

OCLC_EY92_THETA

 [Previous Routine] [Next Routine] [List of Routines] 
 NAME:
  oclc_ey92_theta
 PURPOSE: (one line)
  return theta (bending angle) as a function of radius
 DESCRIPTION:
  return theta as a function of radius
 CATEGORY:
  Occultation lightcurve (oclc)
 CALLING SEQUENCE:
  theta = oclc_ey92_theta(r0,nu0,lam0,a,b,order,r)
 INPUTS:
  r0 - reference radius (scalar, cm)
  nu0 - refractivity at r0 (scalar, unitless)
  lam0 - reference energy ratio (radius/scale height) (scalar, unitless)
  a - exponent for non-constant mean moleculr weight
      mu(r) = mu0 * (r/r0)^-a
  b - exponent for non-constant temperature
      T(r) = T * (r/r0)^b
  order - order for series expansion
  r - radius at which to evaluat theta (scalar or array, cm)
 OPTIONAL INPUT PARAMETERS:
  none
 KEYWORD INPUT PARAMETERS:
  none
 KEYWORD OUTPUT PARAMETERS:
  none
 OUTPUTS:
  theta at r
  bending angle, negative for an isothermal atmosphere.
  radians.
 COMMON BLOCKS:
  None
 SIDE EFFECTS:
 RESTRICTIONS:
  None
 PROCEDURE:
  Elliot and Young 1992, AJ 103, pp 991; eg 3.20
 MODIFICATION HISTORY:
  Written 2006 Jun 28, Leslie Young

(See oclc_ey92_theta.pro)

OCLC_EY92_TIMEMID

 [Previous Routine] [Next Routine] [List of Routines] 
 NAME:
  oclc_ey92_timemid
 PURPOSE: (one line)
  return mid time of event
 DESCRIPTION:
  return mid time of event
 CATEGORY:
  Occultation lightcurve (oclc)
 CALLING SEQUENCE:
  timemid = oclc_ey92_timemid(timeim, timeem)
 INPUTS:
  v = sky-plane velocity (assuming linear motion), cm/s
  rhomin = minimum observer distance (km)
 OPTIONAL INPUT PARAMETERS:
  none
 KEYWORD INPUT PARAMETERS:
  none
 KEYWORD OUTPUT PARAMETERS:
 OUTPUTS:
  timemin = minimum observer distance in units of s.
 COMMON BLOCKS:
  None
 SIDE EFFECTS:
 RESTRICTIONS:
  Admitedly silly to have this as a seperate function
 PROCEDURE:
  Elliot and Young 1992, AJ 103, pp 991; eq 5.8
 MODIFICATION HISTORY:
  Written 2006 Aug 18, Leslie Young

(See oclc_ey92_timemid.pro)

OCLC_EY92_TIMEMIN_OF_RHOMIN

 [Previous Routine] [List of Routines] 
 NAME:
  oclc_ey92_timemin_of_rhomin
 PURPOSE: (one line)
  return minimum observer distance in units of s.
 DESCRIPTION:
  return minimum observer distance in units of s.
 CATEGORY:
  Occultation lightcurve (oclc)
 CALLING SEQUENCE:
  time_min = oclc_ey92_timemin_of_rhomin(v, rhomin)
 INPUTS:
  v = sky-plane velocity (assuming linear motion), cm/s
  rhomin = minimum observer distance (km)
 OPTIONAL INPUT PARAMETERS:
  none
 KEYWORD INPUT PARAMETERS:
  none
 KEYWORD OUTPUT PARAMETERS:
 OUTPUTS:
  timemin = minimum observer distance in units of s.
 COMMON BLOCKS:
  None
 SIDE EFFECTS:
 RESTRICTIONS:
  Admitedly silly to have this as a seperate function
 PROCEDURE:
  Elliot and Young 1992, AJ 103, pp 991; eq 5.8
 MODIFICATION HISTORY:
  Written 2006 Aug 18, Leslie Young

(See oclc_ey92_timemin_of_rhomin.pro)