This is software written in IDL for calculating the zodiacal light predicted
from the DIRBE team interplanetary dust model (Kelsall et al. 1998, ApJ, 508, 
44).  This model was obtained by fitting the time variation of DIRBE all-sky 
observations over 10 months with a parameterized model of the dust cloud.  
Thus it is constrained over solar elongation angles ranging from 64 to 124 
degrees. 

The Kelsall_zodi_model.tar file contains

  Readme   - this file
  colcorr.tab - ascii table of color correction factors
  confint.pro
  get_zmod.pro 
  mk_zdata.pro
  str_len.pro
  zkernel6p8.pro
  zl_coeffs_str.pro
  zmodel_init.pro
  zpars.xdr - IDL saveset of model parameters
  zutils.pro


The top level routine is the function get_zmod.pro.


; Function GET_ZMOD
;
; Function to compute the Kelsall et al. 1998 (ApJ,508,44) ZODI model 
; intensity for a given set of LOS, Wavelength, and Time.
;
; Written By:   BA Franz
; Date:         23 January 1995
;
; Usage:
;     zodi = GET_ZMOD(lambda,day,lon,lat)
;
; Inputs:
;     lambda - nominal DIRBE wavelength(s) in microns (any combination
;               of 1.25, 2.2, 3.5, 4.9, 12, 25, 60, 100, 140, 240) 
;     day    - 1990 day number(s), where 1.0 = 1 Jan 1990
;     lon    - ecliptic longitude(s)
;     lat    - ecliptic latitude(s)
;
;     note: lambda, day, lon, lat can be any mixture of scalars or arrays,
;           but all arrays must be of same length.
;
; Outputs:
;
;     zodi - scalar or array of zodi model intensity in MJy sr-1
;            For comparison with DIRBE data, this is 'quoted' intensity
;            at the nominal wavelength(s) for an assumed source
;            spectrum nu*F_nu = constant
;
; Optional Keyword Inputs:
;
;     zpar=zpar - array of zodi model params.  If not specified, the
;                 default will be restored from zpars.xdr
;     no_colcorr=no_colcorr - if set, the result will be actual
;                  intensity at the specified wavelength(s) instead of
;                  quoted intensity for a spectrum nu*F_nu = constant