model_utilities.py | model_utilities.py | |||
---|---|---|---|---|
# -*- coding: utf-8 -*- | # -*- coding: utf-8 -*- | |||
import numpy as np | import numpy as np | |||
from matplotlib import pyplot as plt | from matplotlib import pyplot as plt | |||
from astropy.table import Table | ||||
# LOGGING CONFIGURATION | # LOGGING CONFIGURATION | |||
from .galpak3d import * | from .galpak3d import * | |||
from .galaxy_parameters import GalaxyParameters | from .galaxy_parameters import GalaxyParameters | |||
import logging | import logging | |||
logging.basicConfig(level=logging.INFO) | logging.basicConfig(level=logging.INFO) | |||
logger = logging.getLogger('GalPaK: DiskUtilities:') | logger = logging.getLogger('GalPaK: DiskUtilities:') | |||
################################################################# | ################################################################# | |||
# | # | |||
class ModelExt: | class ModelExt: | |||
##Methods to update galaxy parameters | ##Methods to update galaxy parameters | |||
def set_v22(self, galaxy): | def set_v22(self, galaxy): | |||
""" | """ | |||
Computes velocity at 2.2 Re | Computes velocity at 2.2 Rd | |||
:param galaxy: GalaxyParameters | Table | :param galaxy: GalaxyParameters | Table | |||
:return: return V(2.2 Re) | :return: return V(2.2 Rd) | |||
This Model method will add/update the galaxy property | This Model method will add/update the galaxy property | |||
""" | """ | |||
#V at 2.2 x Rd | #V at 2.2 x Rd | |||
x = 2.2 * galaxy['radius'] / 1.68 | x = 2.2 * galaxy['radius'] / 1.68 | |||
_v22 = self.set_velocity_profile(galaxy, x) | _v22 = self.set_velocity_profile(galaxy, x) | |||
#set v22 for galaxy | #set v22 for galaxy | |||
if isinstance(galaxy, GalaxyParameters): | if isinstance(galaxy, GalaxyParameters): | |||
logger.info("Adding v22 property to GalaxyParameters class") | logger.info("Adding v22 property to GalaxyParameters class") | |||
galaxy.set_v22(_v22) | galaxy.set_v22(_v22) | |||
elif isinstance(galaxy, Table): | elif isinstance(galaxy, Table): | |||
logger.info("Adding v22 property to Table class") | logger.info("Adding v22 property to Table class") | |||
if 'v22' not in galaxy.colnames: | if 'v22' not in galaxy.colnames: | |||
galaxy.add_column(_v22,name='v22') | galaxy.add_column(_v22, name='v22') | |||
else: | else: | |||
galaxy['v22']=_v22 | galaxy['v22']=_v22 | |||
else: | else: | |||
raise Exception("set_v22 must use a GalaxyParameters or a Table ") | raise Exception("set_v22 must use a GalaxyParameters or a Table ") | |||
def set_dvdx(self, galaxy): | def set_dvdx(self, galaxy): | |||
""" | """ | |||
Computes Velocity inner slope at R=0 | Computes Velocity inner slope at R=0 | |||
:param galaxy: GalaxyParameters | Table | :param galaxy: GalaxyParameters | Table | |||
:return: dv/dx(r=0) using the units of GalaxyParameters (km/s / pix ) | :return: dv/dx(r=0) using the units of GalaxyParameters (km/s / pix ) | |||
End of changes. 4 change blocks. | ||||
3 lines changed or deleted | 4 lines changed or added | |||
This html diff was produced by rfcdiff 1.41. The latest version is available from http://tools.ietf.org/tools/rfcdiff/ |