PRO dustem_write_isrf_lv,file,st

Ncomments=6
c=strarr(Ncomments)

c(0)='# DUSTEM: exciting radiation field featuring'
c(1)='# Mathis ISRF'
c(2)='# Blackbody with     T=   0.0000E+00'
c(3)='# dilution factor wdil=   1.0000E+00'
c(4)='# Nbr of points'
c(5)='# wave (microns), 4*pi*Inu (erg/cm2/s/Hz)'

openw,unit,file,/get_lun

FOR i=0,Ncomments-1 DO BEGIN
  printf,unit,c(i)
ENDFOR

N_waves=n_elements(st)
printf,unit,n_waves

FOR i=0L,N_waves-1 DO BEGIN
  printf,unit,st(i).lambisrf,st(i).isrf
ENDFOR

close,unit
free_lun,unit

;stop

END