From cff405d5f14195f219c07e6ff1dd802b29fb2b8e Mon Sep 17 00:00:00 2001
From: Annie Hughes <annie.hughes@irap.omp.eu>
Date: Wed, 29 Mar 2023 10:41:27 +0200
Subject: [PATCH]  improved header and doc help

---
 src/idl/dustem_planck_function.pro | 14 ++++++++++++--
 1 file changed, 12 insertions(+), 2 deletions(-)

diff --git a/src/idl/dustem_planck_function.pro b/src/idl/dustem_planck_function.pro
index 78a111d..70bb778 100755
--- a/src/idl/dustem_planck_function.pro
+++ b/src/idl/dustem_planck_function.pro
@@ -1,4 +1,4 @@
-function dustem_planck_function,temp,wave
+function dustem_planck_function,temp,wave,help=help
 
 ;+
 ; NAME:
@@ -20,6 +20,9 @@ function dustem_planck_function,temp,wave
 ; OUTPUT:
 ;	BBFLUX - Scalar or vector giving the Planck function at the specified
 ;		wavelength points.
+; ACCEPTED KEYWORDS
+;    help      = If set, print this help
+;
 ; COMMON BLOCKS:
 ;    None
 ;
@@ -37,7 +40,13 @@ function dustem_planck_function,temp,wave
 ;    See http://dustemwrap.irap.omp.eu/ for FAQ and help.  
 ;-
 
- On_error,2
+; On_error,2
+
+ IF keyword_set(help) THEN BEGIN
+  doc_library,'dustem_planck_function'
+  bbflux=-1.
+  goto,the_end
+ENDIF
 
  if ( N_elements(wave) LT 1 ) then begin
      print,'Syntax - bbflux = dustem_planck_function(temp,wave)'
@@ -68,6 +77,7 @@ function dustem_planck_function,temp,wave
   fact=1.06177e+06
   bbflux=bbflux*w*w*fact ;MJy/sr
 
+the_end:
 ;  return, bbflux*1.E-8              ; Convert to ergs/cm2/s/A
   RETURN,bbflux
 
--
libgit2 0.21.2