replicated ext format procedures

Ilyes Choubani
1 parent 09e61d64
Showing 2 changed files with 119 additions and 0 deletions Show diff stats
src/idl/dustem_initialize_ext.pro
src/idl/dustem_old2new_ext_format.pro
@@ -0,0 +1,90 @@
+FUNCTION dustem_initialize_ext,Next,comments=comments,help=help
+
+;+
+; NAME:
+;    dustem_initialize_ext
+; PURPOSE:
+;    Initializes an empty Dustem EXT structure
+; CATEGORY:
+;    DUSTEM Wrapper
+; CALLING SEQUENCE:
+;    st=dustem_initialize_ext(Next[,/help])
+; INPUTS:
+;    Next: Number of elements/filters in the EXT structure
+; OPTIONAL INPUT PARAMETERS:
+;    None
+; OUTPUTS:
+;    st = EXT structure
+;    The structure contains the following tags
+;      instru  : name of the instrument. Actually not used by DustemWrap.
+;      filter  : name of the filter. Actually defines the wavelength, the flux convention, etc ...
+;      wave    : wavelength in microns. Actually not used by DustemWrap. [microns]
+;      EXT_I   : Total Extinction Cross-section (a.k.a. Ext_I) [cm^2/H]
+;      EXT_Q   : Q Extinction Cross-section [cm^2/H]
+;      EXT_U   : U Extintion Cross-section  [cm^2/H]
+;      EXT_P   : P Extinction Cross-section [cm^2/H] (computed from EXT_Q,EXT_U or directly from VG's models)
+;      EXT_smallp  : Smallp Extinction Cross-section (Extinction polarization fraction) []
+;      psi     : polarization angle [deg] (computed from EXT_Q,EXT_U)
+;      sigextII : = error^2 (for completness on sigmas) [cm^2/H^2]
+;      sigextQQ : variance on Q [cm^2/H^2]
+;      sigextUU : variance on U [cm^2/H^2]
+;      sigextIQ : IQ covariance [cm^2/H^2]
+;      sigextIU : IU covariance [cm^2/H^2]
+;      sigextQU : QU covariance [cm^2/H^2]
+;      sigextP : uncertainty on largep [cm^2/H]
+;      sigextsmallp = uncertainty on smallp [] 
+;      sigextpsi  : uncertainty on psi [deg]
+; OPTIONAL OUTPUT PARAMETERS:
+;    comments : string comments about the EXT SED
+; ACCEPTED KEY-WORDS:
+;    help                  = if set, print this help
+; COMMON BLOCKS:
+;    None
+; SIDE EFFECTS:
+;    None
+; RESTRICTIONS:
+;    None
+; PROCEDURE:
+;    Straightforward
+; EXAMPLES
+;    st=dustem_initialize_ext(10)
+;    help,st
+; MODIFICATION HISTORY:
+;    Written by JPB 
+;-
+
+IF keyword_set(help) THEN BEGIN
+  doc_library,'dustem_initialize_ext'
+  sed=0.
+  goto,the_end
+ENDIF
+
+iv=la_undef(4)
+one_st={instru:'',filter:'',wave:iv,EXT_I:iv,EXT_Q:iv,EXT_U:iv,EXT_P:iv,EXT_smallp:iv,psi:iv,sigextII:iv,sigextQQ:iv,sigextUU:iv,sigextIQ:iv,sigextIU:iv,sigextQU:iv,sigextP:iv,sigextsmallp:iv,sigext_psi:iv}
+
+comments=['Dustem-Wrap EXT']
+comments=[comments,'instru: Instrument name']
+comments=[comments,'filter: Instrument filter']
+comments=[comments,'wave: Filter reference wavelength [mic]']
+comments=[comments,'EXT_I: EXT intensity [cm^2/H]']
+comments=[comments,'EXT_Q: Q Stokes parameter extinction cross-section [cm^2/H]']
+comments=[comments,'EXT_U: U Stokes parameter extinction cross-section [cm^2/H]']
+comments=[comments,'EXT_P: Extinction Polarized Intensity [cm^2/H]']
+comments=[comments,'EXT_smallp: Extinction Polarized Fraction []']
+comments=[comments,'psi: polarization angle [deg]']
+comments=[comments,'sigextII: variance on EXT_I [(cm^2/H)^2]']
+comments=[comments,'sigextQQ: variance on EXT_Q [(cm^2/H)^2]']
+comments=[comments,'sigextUU: variance on EXT_U [(cm^2/H)^2]']
+comments=[comments,'sigextIQ: covariance on EXT_I and EXT_Q [(cm^2/H)^2]']
+comments=[comments,'sigextIU: covariance on EXT_I and EXT_U [(cm^2/H)^2]']
+comments=[comments,'sigextQU: covariance on EXT_Q and EXT_U [(cm^2/Hs)^2]']
+comments=[comments,'sigextP: variance on Extinction Polarized Intensity (P) [(cm^2/H)^2]']
+comments=[comments,'sigextsmallp: variance on Extinction Polarized fraction (EXT_P) []']
+comments=[comments,'sigextpsi: variance on polarisation fraction (psi) [deg^2]']
+
+sed=replicate(one_st,Next)
+
+the_end:
+RETURN,sed
+
+END
@@ -0,0 +1,29 @@
+PRO dustem_old2new_ext_format
+
+
+model='THEMIS'
+dustem_init,mode=model
+
+;dir=!dustem_wrap_soft_dir+'/Data/SEDs/'
+dir='/Users/ilyeschoubani/Downloads/test_wrap_NY/'
+
+file=dir+'Serkowski_HD21.xcat'
+
+str='cp '+file+' '+file+'.old'
+spawn,str
+
+old_ext=read_xcat(file,/silent,comments=old_comments)
+Next=n_elements(old_ext)
+
+new_ext=dustem_initialize_ext(Next,comments=new_comments,help=help)
+comments=['Written by dustem_old2new_ext_format.pro on '+systime(0),old_comments,new_comments]
+
+new_ext.instru=old_ext.instru
+new_ext.filter=old_ext.filter
+new_ext.wave=old_ext.wave
+new_ext.EXT_I=old_ext.spec
+new_ext.sigextII=(old_ext.error)^2
+
+write_xcat,new_ext,file,comments=comments
+
+END