FUNCTION dustem_create_continuum,key=key,val=val
  
  temp=1000.
  ampl=1.d-2
  
  IF keyword_set(key) THEN BEGIN 
     a=where(key EQ 1,count1)
     b=where(key EQ 2,count2)
     IF count1 NE 0 then temp=(val(a))(0)
     IF count2 NE 0 then ampl=(val(b))(0)
;     IF count1 NE 0 then temp=1d0*10.^((val(a))(0))
;     IF count2 NE 0 then ampl=1d0*10.^((val(b))(0))
;     print, "key=",key
;     print, "val=",val
  ENDIF
  
;  stop
;  IF getenv('DUSTEM_WHICH') EQ 'DESERT' THEN BEGIN
  IF !dustem_which EQ 'DESERT' THEN BEGIN
    lambir=((*!dustem_params).gemissiv.lambir)
  ENDIF ELSE BEGIN
    lambir=((*!dustem_params).lambda.lambda)
  ENDELSE

  norm = max(dustem_planck_function(temp,lambir))
  output = ampl*dustem_planck_function(temp,lambir)/norm
 
  return,output
  
END