Commit 6cfdcb1c4c5614f3e73e1b3d5a72b703c3f53614

Authored by Annie Hughes
1 parent ba74c712
Exists in master

updated SED_1

Data/EXAMPLE_OBSDATA/example_SED_1.xcat
1   -\Written by write_xcat.pro on Tue Jul 26 14:17:56 2022
2   -\Written by dustem_old2new_sed_format.pro on Tue Jul 26 14:17:56 2022
3   -\Sample SED values for test of Dustem (see file dustem_fit_sed_readme.txt)
4   -\Wavelengths are mixed on purpose, to show that it does not affect processing
5   -\Dustem-Wrap SED
6   -\instru: Instrument name
  1 +\DustemWrap SED
  2 +\Comment lines start with \
  3 +\Undefined/missing value is -32768.0
  4 +\An example SED in .xcat format for illustrating how to use DustEMWrap
  5 +\The wavelengths are disordered on purpose, to show that the ordering of the input data does not affect the fitting
  6 +\=====================================================
  7 +\instru: Instrument name [for user convenience]
7 8 \filter: Instrument filter
8   -\wave: Filter reference wavelength [mic]
  9 +\wave: Filter reference wavelength [micron, for user convenience]
9 10 \StokesI: SED intensity [MJy/sr]
10 11 \StokesQ: Q Stokes parameter intensity [MJy/sr]
11 12 \StokesU: U Stokes parameter intensity [MJy/sr]
12   -\largeP: Polarized intensity P [MJy/sr]
13   -\smallp: polarization fraction []
14   -\psi: polarization angle [deg]
  13 +\largeP: Polarized intensity P [MJy/sr, for user convenience]
  14 +\smallp: polarization fraction [0.01 = 1%, for user convenience]
  15 +\psi: polarization angle [deg, for user convenience]
15 16 \sigmaII: variance on Stokes I [(MJy/sr)^2]
16 17 \sigmaQQ: variance on Stokes Q [(MJy/sr)^2]
17 18 \sigmaUU: variance on Stokes U [(MJy/sr)^2]
18 19 \sigmaIQ: covariance on Stokes I and Q [(MJy/sr)^2]
19 20 \sigmaIU: covariance on Stokes I and U [(MJy/sr)^2]
20 21 \sigmaQU: covariance on Stokes Q and U [(MJy/sr)^2]
21   -\sigma_largep: variance on Polarized Intensity (P) [(MJy/sr)^2]
22   -\sigma_smallp: variance on polarisation fraction (p) []
23   -\sigma_psi: variance on polarisation fraction (psi) [deg^2]
24   -\CAUTION: p, P and psi values are NOT debiased
  22 +\sigma_largep: variance on Polarized Intensity (P) [(MJy/sr)^2, for user convenience]
  23 +\sigma_smallp: variance on polarisation fraction (p) [for user convenience]
  24 +\sigma_psi: variance on polarisation fraction (psi) [deg^2, for user convenience]
  25 +\=====================================================
25 26 |INSTRU |FILTER |WAVE |STOKESI |STOKESQ |STOKESU |LARGEP |SMALLP |PSI |SIGMAII |SIGMAQQ |SIGMAUU |SIGMAIQ |SIGMAIU |SIGMAQU |SIGMA_LARGEP |SIGMA_SMALLP |SIGMA_PSI |
26 27 |char |char |real |real |real |real |real |real |real |real |real |real |real |real |real |real |real |real |
27 28 |NULL |NULL |NULL |NULL |NULL |NULL |NULL |NULL |NULL |NULL |NULL |NULL |NULL |NULL |NULL |NULL |NULL |NULL |
  29 +SPIRE SPIRE1 250.0 1.8 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 2.90004e-02 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0
  30 +SPIRE SPIRE2 350.0 0.9 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 2.90004e-02 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0
  31 +SPIRE SPIRE3 500 0.4 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 5.90004e-03 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0
  32 +NIKA NIKA21 1200 0.15 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 1.90004e-03 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0
  33 +NIKA NIKA22 2000 0.2 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 5.10004e-04 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0
28 34 IRAC IRAC1 3.60000 0.00130775 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 7.90004e-08 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0
29 35 IRAC IRAC2 4.50000 0.00123874 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 4.18501e-08 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0
30 36 IRAC IRAC4 8.00000 0.0504011 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 7.57350e-06 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0
... ... @@ -36,3 +42,4 @@ MIPS MIPS2 70.0000 0.861524 -32768.0 -32768.0 -32768.0 -32768.0 -32
36 42 IRAS IRAS4 100.000 1.46069 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 0.00677211 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0
37 43 MIPS MIPS3 160.000 2.12215 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 0.0104897 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0
38 44 IRAC IRAC3 5.80000 0.0154945 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 7.54617e-07 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0
  45 +\=====================================================
... ...
Data/EXAMPLE_OBSDATA/example_SED_1b.xcat deleted
... ... @@ -1,45 +0,0 @@
1   -\DustemWrap SED
2   -\Comment lines start with \
3   -\Undefined/missing value is -32768.0
4   -\An example SED in .xcat format for illustrating how to use DustEMWrap
5   -\The wavelengths are disordered on purpose, to show that the ordering of the input data does not affect the fitting
6   -\=====================================================
7   -\instru: Instrument name [for user convenience]
8   -\filter: Instrument filter
9   -\wave: Filter reference wavelength [micron, for user convenience]
10   -\StokesI: SED intensity [MJy/sr]
11   -\StokesQ: Q Stokes parameter intensity [MJy/sr]
12   -\StokesU: U Stokes parameter intensity [MJy/sr]
13   -\largeP: Polarized intensity P [MJy/sr, for user convenience]
14   -\smallp: polarization fraction [0.01 = 1%, for user convenience]
15   -\psi: polarization angle [deg, for user convenience]
16   -\sigmaII: variance on Stokes I [(MJy/sr)^2]
17   -\sigmaQQ: variance on Stokes Q [(MJy/sr)^2]
18   -\sigmaUU: variance on Stokes U [(MJy/sr)^2]
19   -\sigmaIQ: covariance on Stokes I and Q [(MJy/sr)^2]
20   -\sigmaIU: covariance on Stokes I and U [(MJy/sr)^2]
21   -\sigmaQU: covariance on Stokes Q and U [(MJy/sr)^2]
22   -\sigma_largep: variance on Polarized Intensity (P) [(MJy/sr)^2, for user convenience]
23   -\sigma_smallp: variance on polarisation fraction (p) [for user convenience]
24   -\sigma_psi: variance on polarisation fraction (psi) [deg^2, for user convenience]
25   -\=====================================================
26   -|INSTRU |FILTER |WAVE |STOKESI |STOKESQ |STOKESU |LARGEP |SMALLP |PSI |SIGMAII |SIGMAQQ |SIGMAUU |SIGMAIQ |SIGMAIU |SIGMAQU |SIGMA_LARGEP |SIGMA_SMALLP |SIGMA_PSI |
27   -|char |char |real |real |real |real |real |real |real |real |real |real |real |real |real |real |real |real |
28   -|NULL |NULL |NULL |NULL |NULL |NULL |NULL |NULL |NULL |NULL |NULL |NULL |NULL |NULL |NULL |NULL |NULL |NULL |
29   -SPIRE SPIRE1 250.0 1.8 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 2.90004e-02 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0
30   -SPIRE SPIRE2 350.0 0.9 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 2.90004e-02 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0
31   -SPIRE SPIRE3 500 0.4 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 5.90004e-03 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0
32   -NIKA NIKA21 1200 0.15 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 1.90004e-03 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0
33   -NIKA NIKA22 2000 0.2 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 5.10004e-04 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0
34   -IRAC IRAC1 3.60000 0.00130775 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 7.90004e-08 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0
35   -IRAC IRAC2 4.50000 0.00123874 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 4.18501e-08 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0
36   -IRAC IRAC4 8.00000 0.0504011 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 7.57350e-06 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0
37   -IRAS IRAS1 12.0000 0.0344639 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 3.87873e-06 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0
38   -MIPS MIPS1 24.0000 0.0373517 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 8.64360e-06 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0
39   -IRAS IRAS3 60.0000 0.732637 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 0.00243502 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0
40   -IRAS IRAS2 25.0000 0.0557090 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 3.15081e-05 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0
41   -MIPS MIPS2 70.0000 0.861524 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 0.00355327 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0
42   -IRAS IRAS4 100.000 1.46069 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 0.00677211 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0
43   -MIPS MIPS3 160.000 2.12215 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 0.0104897 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0
44   -IRAC IRAC3 5.80000 0.0154945 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 7.54617e-07 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0 -32768.0
45   -\=====================================================