dustem_define_isrf_stellar_modifyer_variable.pro 702 Bytes
FUNCTION dustem_define_isrf_stellar_modifyer_variable,stellar_values=stellar_values

one_st={radius:0.,temperature:0.,distance:0.,amplitude:0.}

IF not keyword_set(stellar_values) THEN BEGIN
  ;The defaults is one sun at 1au, but with zero amplitude
  Nstar=1
  st=replicate(one_st,Nstar)
  one_pc=3.08e16
  one_au=1.49e11
  st[0].radius=1.          ;Rsun
  st[0].temperature=5770.    ;K
  st[0].distance=one_au/one_pc         ;pc
  st[0].amplitude=0.         ;--
  status=1
ENDIF ELSE BEGIN
  Nstars=n_elements(stellar_values)
  ;st=replicate(one_st,Nstar)
  st=stellar_values
  status=2
ENDELSE

defsysv,'!dustem_isrf_star_add',st
message,'defined !dustem_isrf_star_add',/continue

RETURN,status

END