Using the radius parameter with file

tomasstolker · Feb 8, 2020 · 483be7b · 483be7b
1 parent bce87d3
commit 483be7b
Show file tree

Hide file tree

Showing 3 changed files with 5 additions and 4 deletions.
diff --git a/diskmap/__init__.py b/diskmap/__init__.py
@@ -5,4 +5,4 @@
 __version__ = '0.0.3'
 __maintainer__ = 'Tomas Stolker'
 __email__ = '[email protected]'
-__status__ = 'Development'
+__status__ = 'Development'4
diff --git a/diskmap/diskmap.py b/diskmap/diskmap.py
@@ -28,7 +28,8 @@ def __init__(self,
         pixscale : float
             Pixel scale (arcsec).
         inclination : float
-            Inclination of the disk (deg).
+            Inclination of the disk (deg). Include a minus sign to exchange the near and far side
+            in the mapping of the disk.
         pos_angle : float
             Position angle of the disk (deg). Defined in counterclockwise direction with respect
             to the vertical axis (i.e. east of north).
@@ -124,7 +125,7 @@ def power_law_height(x_power, a_power, b_power, c_power):
             data = np.loadtxt(filename)
 
             # midplane radius [au]
-            disk_radius = np.linspace(radius[0], radius[-1], 100)
+            disk_radius = np.linspace(radius[0], radius[1], radius[2])
 
             # disk height [au]
             height_interp = interp1d(data[:, 0], data[:, 1])

diff --git a/setup.py b/setup.py
@@ -12,7 +12,7 @@
 
 setup(
     name='diskmap',
-    version='0.0.3',
+    version='0.0.4',
     description='Scattered light mapping of protoplanetary disks',
     long_description=open('README.rst').read(),
     long_description_content_type='text/x-rst',