FUNCTION BB_FLUX,TEMP,WAVE
;******************************************************************************
;+
; NAME:
;	BB_FLUX
; PURPOSE:
;       gives the blackbody flux (NOT the "astrophysical flux") in 
;       ergs/cm2/s/sr from allen astrophysical quantities section 44
; CALLING SEQUENCE:
;       RESULT = BB_FLUX(TEMPERATURE,WAVELENGTH)
; INPUTS:
;	TEMP = tempreature
;	WAVE = wavelength  - must be in microns
; OPTIONAL INPUT PARAMETERS:
;	NONE
; OUTPUTS:
;       blackbody flux
; COMMON BLOCKS:
;	NONE.
; SIDE EFFECTS:
;	NONE.
; PROCEDURE:
;	Straightforward.
; RESTRICTIONS
;       Works on V 1 and V 2 of IDL
; MODIFICATION HISTORY:
;	Written, Stephen Merkowitz 1989
;-
;******************************************************************************
  TEMP = DOUBLE(TEMP)
   W  = (WAVE > 0.001)/1D4      ; microns to cm
   C1 = 1.19106D-5
   C2 = 1.438786D0
RETURN, C1/((W^5)*(EXP(C2/(W*TEMP))-1.0D0)) > 0.0D0
END