mwdust 1.6__cp312-cp312-macosx_10_13_x86_64.whl

mwdust/util/extCurves/apj398709t6_ascii.txt

@@ -0,0 +1,51 @@
+Table 6
+F99 Reddening in Different Bandpasses
+
+Bandpass	lambda_eff	R_V	Bandpass	lambda_eff	R_V	
+	2.1	3.1	4.1	5.1		2.1	3.1	4.1	5.1	
+Landolt U	3508.2	5.614	4.334	3.773	3.460	WFPC2 F300W	3087.6	6.777	4.902	4.127	3.710	
+Landolt B	4329.0	4.355	3.626	3.290	3.096	WFPC2 F450W	4587.0	4.014	3.410	3.132	2.971	
+Landolt V	5421.7	2.953	2.742	2.645	2.589	WFPC2 F555W	5439.4	2.976	2.755	2.653	2.594	
+Landolt R	6427.8	2.124	2.169	2.189	2.201	WFPC2 F606W	5984.8	2.469	2.415	2.389	2.375	
+Landolt I	8048.4	1.410	1.505	1.548	1.573	WFPC2 F702W	6887.9	1.850	1.948	1.994	2.020	
+CTIO U	3733.9	5.170	4.107	3.628	3.355	WFPC2 F814W	7940.0	1.452	1.549	1.594	1.620	
+CTIO B	4308.9	4.382	3.641	3.300	3.104	WFC3 F105W	10438.9	0.981	0.969	0.964	0.961	
+CTIO V	5516.6	2.857	2.682	2.600	2.553	WFC3 F110W	11169.7	0.907	0.881	0.870	0.863	
+CTIO R	6520.2	2.055	2.119	2.149	2.166	WFC3 F125W	12335.5	0.778	0.726	0.701	0.687	
+CTIO I	8006.9	1.420	1.516	1.561	1.587	WFC3 F140W	13692.3	0.672	0.613	0.586	0.570	
+UKIRT J	12482.9	0.764	0.709	0.684	0.669	WFC3 F160W	15258.3	0.570	0.512	0.485	0.469	
+UKIRT H	16588.4	0.502	0.449	0.425	0.411	WFC3 F200LP	5515.2	3.457	2.958	2.743	2.625	
+UKIRT K	21897.7	0.331	0.302	0.288	0.280	WFC3 F218W	2248.3	12.405	7.760	5.956	5.027	
+UKIRT L"	37772.5	0.159	0.153	0.150	0.148	WFC3 F225W	2394.0	10.907	6.989	5.458	4.666	
+Gunn g	5200.0	3.225	2.914	2.770	2.687	WFC3 F275W	2742.5	7.986	5.487	4.488	3.963	
+Gunn r	6628.5	1.959	2.055	2.099	2.125	WFC3 F300X	2934.5	7.437	5.228	4.331	3.854	
+Gunn i	7898.6	1.454	1.555	1.601	1.628	WFC3 F336W	3366.4	5.835	4.453	3.853	3.519	
+Gunn z	9050.1	1.188	1.234	1.255	1.267	WFC3 F350LP	5877.1	2.876	2.624	2.509	2.443	
+Spinrad R_S	6927.3	1.810	1.921	1.972	2.002	WFC3 F390W	3994.8	4.803	3.896	3.481	3.244	
+Stromgren u	3510.0	5.539	4.305	3.759	3.452	WFC3 F438W	4335.3	4.347	3.623	3.288	3.095	
+Stromgren b	4670.5	3.916	3.350	3.089	2.938	WFC3 F475W	4785.0	3.755	3.248	3.013	2.878	
+Stromgren v	4119.2	4.619	3.793	3.411	3.191	WFC3 F475X	4969.7	3.548	3.116	2.917	2.803	
+Stromgren beta	4861.3	3.655	3.183	2.966	2.840	WFC3 F555W	5302.8	3.135	2.855	2.726	2.652	
+Stromgren y	5478.9	2.862	2.686	2.605	2.557	WFC3 F600LP	7362.3	1.688	1.781	1.824	1.849	
+SDSS u	3586.8	5.419	4.239	3.715	3.419	WFC3 F606W	5868.5	2.581	2.488	2.445	2.421	
+SDSS g	4716.7	3.843	3.303	3.054	2.910	WFC3 F625W	6225.8	2.230	2.259	2.273	2.281	
+SDSS r	6165.1	2.255	2.285	2.300	2.308	WFC3 F775W	7630.9	1.533	1.643	1.694	1.724	
+SDSS i	7475.9	1.583	1.698	1.751	1.782	WFC3 F814W	7983.1	1.441	1.536	1.580	1.605	
+SDSS z	8922.9	1.211	1.263	1.286	1.300	WFC3 F850LP	9149.7	1.168	1.208	1.226	1.237	
+DSS-II g	4620.6	3.970	3.381	3.110	2.954	ACS clear	6211.1	2.612	2.436	2.356	2.309	
+DSS-II r	6545.5	1.991	2.088	2.133	2.159	ACS F435W	4348.3	4.330	3.610	3.278	3.087	
+DSS-II i	8111.0	1.396	1.487	1.530	1.554	ACS F475W	4760.3	3.787	3.268	3.028	2.890	
+PS1 g	4876.7	3.634	3.172	2.958	2.835	ACS F550M	5581.0	2.754	2.620	2.558	2.522	
+PS1 r	6200.1	2.241	2.271	2.284	2.292	ACS F555W	5361.3	3.031	2.792	2.682	2.618	
+PS1 i	7520.8	1.568	1.682	1.734	1.765	ACS F606W	5901.0	2.555	2.471	2.431	2.409	
+PS1 z	8665.3	1.258	1.322	1.352	1.369	ACS F625W	6298.1	2.171	2.219	2.241	2.254	
+PS1 y	9706.3	1.074	1.087	1.094	1.097	ACS F775W	7673.5	1.520	1.629	1.679	1.708	
+PS1 w	6240.8	2.425	2.341	2.302	2.280	ACS F814W	8012.4	1.432	1.526	1.569	1.594	
+LSST u	3693.2	5.243	4.145	3.652	3.373	ACS F850LP	9007.5	1.196	1.243	1.265	1.277	
+LSST g	4797.3	3.739	3.237	3.006	2.872	DES g	4796.6	3.739	3.237	3.006	2.872	
+LSST r	6195.8	2.245	2.273	2.286	2.294	DES r	6382.6	2.113	2.176	2.205	2.221	
+LSST i	7515.3	1.571	1.684	1.737	1.767	DES i	7769.0	1.490	1.595	1.644	1.672	
+LSST z	8664.4	1.259	1.323	1.353	1.370	DES z	9108.2	1.175	1.217	1.236	1.247	
+LSST y	9710.3	1.075	1.088	1.094	1.098	DES Y	9850.4	1.051	1.058	1.061	1.063	
+
+Notes. Ab/E(B - V)SFD in different bandpasses b, evaluated according to an F99 reddening law with normalization N = 0.78 and RV = 2.1, 3.1, 4.1, and 5.1, using a 7000 K source spectrum. The column lambdaeff gives the throughput-weighted mean wavelength in the bandpass. When used with RV = 3.1, these give the coefficients to use with E(B - V)SFD to get reddenings consistent with the results of this work and Schlafly et al. (2010). The values for other RV are provided only for convenience.

mwdust/util/extCurves/extinction.tbl

@@ -0,0 +1,44 @@
+\Coordinates = 180. 0. gal (   86.404890000    28.936160000 equ J2000)
+\E(B-V)_SFD_1998 =     1.418 (mag)
+\SandF: Schlafly and Finkbeiner 2011 (ApJ 737, 103)
+\SFD: Schlegel et al. 1998 (ApJ 500, 525)
+\ Filter_name
+\ ___ Name of filter
+\ LamEff
+\ ___ Effective wavelength of filter (in microns) as defined in SandF
+\ A_over_E_B_V_SandF
+\ ___ Ratio of SandF band extinction to SandF E(B-V)
+\ A_SandF
+\ ___ Band extinction (mags) for SandF
+\ A_over_E_B_V_SFD
+\ ___ Ratio of SFD band extinction to SandF E(B-V)
+\ A_SFD
+\ ___ Band extinction (mags) for SFD
+|Filter_name         |LamEff      |A_over_E_B_V_SandF      |A_SandF                 |A_over_E_B_V_SFD        |A_SFD                   |
+|char                |float       |float                   |float                   |float                   |float                   |
+|                    |microns     |                        |mags                    |                        |mags                    |
+ CTIO U               0.3734       4.107                    5.824                    4.968                     7.045                    
+ CTIO B               0.4309       3.641                    5.164                    4.325                     6.134                    
+ CTIO V               0.5517       2.682                    3.804                    3.240                     4.595                    
+ CTIO R               0.6520       2.119                    3.005                    2.634                     3.735                    
+ CTIO I               0.8007       1.516                    2.150                    1.962                     2.782                    
+ DSS-II g             0.4621       3.381                    4.795                    3.907                     5.541                    
+ DSS-II r             0.6546       2.088                    2.961                    2.649                     3.757                    
+ DSS-II i             0.8111       1.487                    2.109                    1.893                     2.685                    
+ SDSS u               0.3587       4.239                    6.012                    5.155                     7.311                    
+ SDSS g               0.4717       3.303                    4.684                    3.793                     5.379                    
+ SDSS r               0.6165       2.285                    3.241                    2.751                     3.901                    
+ SDSS i               0.7476       1.698                    2.408                    2.086                     2.958                    
+ SDSS z               0.8923       1.263                    1.791                    1.479                     2.097                    
+ UKIRT J              1.248        0.709                    1.005                    0.902                     1.279                    
+ UKIRT H              1.659        0.449                    0.637                    0.576                     0.817                    
+ UKIRT K              2.190        0.302                    0.428                    0.367                     0.520                    
+ 2MASS J              1.23         0.723                    1.025                    0.937                     1.329                    
+ 2MASS H              1.64         0.460                    0.652                    0.591                     0.838                    
+ 2MASS Ks             2.16         0.310                    0.440                    0.382                     0.542                    
+ IRAC-1               3.52         0.178                    0.252                    0.220                     0.312                    
+ IRAC-2               4.46         0.148                    0.210                    0.183                     0.260                    
+ IRAC-3               5.66         0.130                    0.184                    0.162                     0.230                    
+ IRAC-4               7.68         0.122                    0.173                    0.151                     0.214                    
+ WISE-1               3.32         0.189                    0.268                    0.234                     0.332                    
+ WISE-2               4.57         0.146                    0.207                    0.180                     0.255                    

mwdust/util/extCurves.py

@@ -0,0 +1,65 @@
+import re
+import os, os.path
+import numpy
+
+#A / E(B-v)
+avebv=  {}
+avebvsf= {}
+def _read_extCurves():
+    """Read the extinction curves files in extCurves"""
+    extFile= open(os.path.join(os.path.dirname(os.path.realpath(__file__)),
+                               'extCurves','extinction.tbl'),'r')
+    global avebv
+    global avebvsf
+    for line in extFile:
+        if line[0] == '\\': continue
+        if line[0] == '|': continue
+        vals= re.split('\s\s+',line)
+        avebv[vals[0].strip()]= float(vals[4])
+        avebvsf[vals[0].strip()]= float(vals[2])    
+    # Add filters from Schlafly & Finkbeiner
+    extFile= open(os.path.join(os.path.dirname(os.path.realpath(__file__)),
+                               'extCurves','apj398709t6_ascii.txt'),'r')
+    cnt= 0
+    for line in extFile:
+        cnt+= 1
+        if cnt < 6: continue
+        if cnt > 49: break
+        vals= line.split()
+        # Each line has 2 filters
+        filter1= '%s %s' % (vals[0],vals[1])
+        if not filter1 in avebv:
+            avebv[filter1]= numpy.nan
+            avebvsf[filter1]= float(vals[4])
+        filter2= '%s %s' % (vals[7],vals[8])
+        if not filter2 in avebv:
+            avebv[filter2]= numpy.nan
+            avebvsf[filter2]= float(vals[11])
+
+_read_extCurves()
+
+def aebv(filter,sf10=True):
+    """
+    NAME:
+       aebv
+    PURPOSE:
+       return A_filter / E(B-V), necessary to turn SFD E(B-V) into total extinctions
+    INPUT:
+       filter - filter to use (e.g., '2MASS Ks')
+       sf10= (True) if True, use the values from Schlafly & Finkbeiner 2010, which use an updated extinction law, source spectrum, and recalibrated SFD map
+    OUTPUT:
+       A_filter / E(B-V)
+    HISTORY:
+       2013-11-24 - Written - Bovy (IAS)
+    """
+    if filter.lower() == 'e(b-v)':
+        return 0.884
+    if sf10:
+        if not filter in avebvsf:
+            raise ValueError("Requested filter is not supported")
+        else:
+            return avebvsf[filter]
+    else:
+        if not filter in avebv:
+            raise ValueError("Requested filter is not supported")
+        return avebv[filter]

mwdust/util/fortranfile.py

@@ -0,0 +1,260 @@
+# Copyright 2008-2010 Neil Martinsen-Burrell
+#
+# Permission is hereby granted, free of charge, to any person obtaining a copy
+# of this software and associated documentation files (the "Software"), to deal
+# in the Software without restriction, including without limitation the rights
+# to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+# copies of the Software, and to permit persons to whom the Software is
+# furnished to do so, subject to the following conditions:
+
+# The above copyright notice and this permission notice shall be included in
+# all copies or substantial portions of the Software.
+
+# THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+# IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+# FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+# AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+# LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+# OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
+# THE SOFTWARE.
+
+# Copyright 2018 Jo Bovy
+#
+# Made small changes to allow this module to be used in Python 3, but only
+# to be able to read the files used in the mwdust package. Keep the same 
+# license as above
+"""Defines a file-derived class to read/write Fortran unformatted files.
+
+The assumption is that a Fortran unformatted file is being written by
+the Fortran runtime as a sequence of records.  Each record consists of
+an integer (of the default size [usually 32 or 64 bits]) giving the
+length of the following data in bytes, then the data itself, then the
+same integer as before.
+
+Examples
+--------
+
+To use the default endian and precision settings, one can just do::
+    >>> f = FortranFile('filename')
+    >>> x = f.readReals()
+
+One can read arrays with varying precisions::
+    >>> f = FortranFile('filename')
+    >>> x = f.readInts('h')
+    >>> y = f.readInts('q')
+    >>> z = f.readReals('f')
+Where the format codes are those used by Python's struct module.
+
+One can change the default endian-ness and header precision::
+    >>> f = FortranFile('filename', endian='>', header_prec='l')
+for a file with little-endian data whose record headers are long
+integers.
+"""
+
+__docformat__ = "restructuredtext en"
+import sys
+_PY3= sys.version > '3'
+if _PY3:
+    from io import FileIO
+    file= FileIO
+import numpy
+
+class FortranFile(file):
+
+    """File with methods for dealing with fortran unformatted data files"""
+
+    def _get_header_length(self):
+        return numpy.dtype(self._header_prec).itemsize
+    _header_length = property(fget=_get_header_length)
+
+    def _set_endian(self,c):
+        """Set endian to big (c='>') or little (c='<') or native (c='=')
+
+        :Parameters:
+          `c` : string
+            The endian-ness to use when reading from this file.
+        """
+        if c in '<>@=':
+            if c == '@':
+                c = '='
+            self._endian = c
+        else:
+            raise ValueError('Cannot set endian-ness')
+    def _get_endian(self):
+        return self._endian
+    ENDIAN = property(fset=_set_endian,
+                      fget=_get_endian,
+                      doc="Possible endian values are '<', '>', '@', '='"
+                     )
+
+    def _set_header_prec(self, prec):
+        if prec in 'hilq':
+            self._header_prec = prec
+        else:
+            raise ValueError('Cannot set header precision')
+    def _get_header_prec(self):
+        return self._header_prec
+    HEADER_PREC = property(fset=_set_header_prec,
+                           fget=_get_header_prec,
+                           doc="Possible header precisions are 'h', 'i', 'l', 'q'"
+                          )
+
+    def __init__(self, fname, endian='@', header_prec='i', *args, **kwargs):
+        """Open a Fortran unformatted file for writing.
+        
+        Parameters
+        ----------
+        endian : character, optional
+            Specify the endian-ness of the file.  Possible values are
+            '>', '<', '@' and '='.  See the documentation of Python's
+            struct module for their meanings.  The deafult is '>' (native
+            byte order)
+        header_prec : character, optional
+            Specify the precision used for the record headers.  Possible
+            values are 'h', 'i', 'l' and 'q' with their meanings from
+            Python's struct module.  The default is 'i' (the system's
+            default integer).
+
+        """
+        file.__init__(self, fname, *args, **kwargs)
+        self.ENDIAN = endian
+        self.HEADER_PREC = header_prec
+
+    def _read_exactly(self, num_bytes):
+        """Read in exactly num_bytes, raising an error if it can't be done."""
+        if _PY3:
+            data = b''
+        else:
+            data = ''
+        while True:
+            l = len(data)
+            if l == num_bytes:
+                return data
+            else:
+                read_data = self.read(num_bytes - l)
+            if read_data == '':
+                raise IOError('Could not read enough data.'
+                              '  Wanted %d bytes, got %d.' % (num_bytes, l))
+            data += read_data
+
+    def _read_check(self):
+        return numpy.frombuffer(
+            self._read_exactly(self._header_length),
+            dtype=self.ENDIAN + self.HEADER_PREC,
+        )[0]
+
+    def _write_check(self, number_of_bytes):
+        """Write the header for the given number of bytes"""
+        self.write(numpy.array(number_of_bytes, 
+                               dtype=self.ENDIAN+self.HEADER_PREC,).tostring()
+                  )
+
+    def readRecord(self):
+        """Read a single fortran record"""
+        l = self._read_check()
+        data_str = self._read_exactly(l)
+        check_size = self._read_check()
+        if check_size != l:
+            raise IOError('Error reading record from data file')
+        return data_str
+
+    def writeRecord(self,s):
+        """Write a record with the given bytes.
+
+        Parameters
+        ----------
+        s : the string to write
+
+        """
+        length_bytes = len(s)
+        self._write_check(length_bytes)
+        self.write(s)
+        self._write_check(length_bytes)
+
+    def readString(self):
+        """Read a string."""
+        return self.readRecord()
+
+    def writeString(self,s):
+        """Write a string
+
+        Parameters
+        ----------
+        s : the string to write
+        
+        """
+        self.writeRecord(s)
+
+    _real_precisions = 'df'
+
+    def readReals(self, prec='f'):
+        """Read in an array of real numbers.
+        
+        Parameters
+        ----------
+        prec : character, optional
+            Specify the precision of the array using character codes from
+            Python's struct module.  Possible values are 'd' and 'f'.
+            
+        """
+        
+        _numpy_precisions = {'d': numpy.float64,
+                             'f': numpy.float32
+                            }
+
+        if prec not in self._real_precisions:
+            raise ValueError('Not an appropriate precision')
+            
+        data_str = self.readRecord()
+        return numpy.frombuffer(data_str, dtype=self.ENDIAN + prec)
+
+    def writeReals(self, reals, prec='f'):
+        """Write an array of floats in given precision
+
+        Parameters
+        ----------
+        reals : array
+            Data to write
+        prec` : string
+            Character code for the precision to use in writing
+        """
+        if prec not in self._real_precisions:
+            raise ValueError('Not an appropriate precision')
+        
+        nums = numpy.array(reals, dtype=self.ENDIAN+prec)
+        self.writeRecord(nums.tostring())
+    
+    _int_precisions = 'hilq'
+
+    def readInts(self, prec='i'):
+        """Read an array of integers.
+        
+        Parameters
+        ----------
+        prec : character, optional
+            Specify the precision of the data to be read using 
+            character codes from Python's struct module.  Possible
+            values are 'h', 'i', 'l' and 'q'
+            
+        """
+        if prec not in self._int_precisions:
+            raise ValueError('Not an appropriate precision')
+            
+        data_str = self.readRecord()
+        return numpy.frombuffer(data_str, dtype=self.ENDIAN + prec)
+
+    def writeInts(self, ints, prec='i'):
+        """Write an array of integers in given precision
+
+        Parameters
+        ----------
+        reals : array
+            Data to write
+        prec : string
+            Character code for the precision to use in writing
+        """
+        if prec not in self._int_precisions:
+            raise ValueError('Not an appropriate precision')
+        
+        nums = numpy.array(ints, dtype=self.ENDIAN+prec)
+        self.writeRecord(nums.tostring())