Const.py
#!/usr/bin/env python
#
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
#
#<LicenseText>
#
# CitcomS.py by Eh Tan, Eun-seo Choi, and Pururav Thoutireddy.
# Copyright (C) 2002-2005, California Institute of Technology.
#
# This program is free software; you can redistribute it and/or modify
# it under the terms of the GNU General Public License as published by
# the Free Software Foundation; either version 2 of the License, or
# (at your option) any later version.
#
# This program is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
# GNU General Public License for more details.
#
# You should have received a copy of the GNU General Public License
# along with this program; if not, write to the Free Software
# Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
#
#</LicenseText>
#
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
#
from CitcomComponent import CitcomComponent
class Const(CitcomComponent):
def __init__(self, name="const", facility="const"):
CitcomComponent.__init__(self, name, facility)
return
def setProperties(self, stream):
from CitcomSLib import Const_set_properties
Const_set_properties(self.all_variables, self.inventory, stream)
return
class Inventory(CitcomComponent.Inventory):
import pyre.inventory
from pyre.units.length import m
from pyre.units.length import km
from pyre.units.mass import kg
from pyre.units.time import s
from pyre.units.temperature import K
from pyre.units.pressure import Pa
from pyre.units.energy import J
#ref_density = pyre.inventory.float("ref_density", default=3500.0*kg/m**3)
#thermdiff = pyre.inventory.float("thermdiff", default=1.0e-6*m**2/s)
#gravacc = pyre.inventory.float("gravacc", default=10.0*m/s**2)
#thermexp = pyre.inventory.float("thermexp", default=3.0e-5/K)
#ref_visc = pyre.inventory.float("ref_visc", default=1.0e21*Pa*s)
#heatcapacity = pyre.inventory.float("heatcapacity", default=1250.0*J/kg/K)
#water_density = pyre.inventory.float("water_density", default=0.0*kg/m**3)
#depth_lith = pyre.inventory.float("depth_lith", default=89e3*m)
#depth_410 = pyre.inventory.float("depth_410", default=410e3*m)
#depth_660 = pyre.inventory.float("depth_660", default=660e3*m)
#depth_d_double_prime = pyre.inventory.float("depth_d_double_prime", default=2691e3*m)
#depth_cmb = pyre.inventory.float("depth_cmb", default=2891e3*m)
# everything in SI units
radius = pyre.inventory.float("radius", default=6371.0e3)
density = pyre.inventory.float("density", default=3340.0)
thermdiff = pyre.inventory.float("thermdiff", default=1.0e-6)
gravacc = pyre.inventory.float("gravacc", default=9.81)
thermexp = pyre.inventory.float("thermexp", default=3.0e-5)
refvisc = pyre.inventory.float("refvisc", default=1.0e21)
cp = pyre.inventory.float("cp", default=1200.0)
density_above = pyre.inventory.float("density_above", default=1030.0)
density_below = pyre.inventory.float("density_below", default=6600.0)
z_lith = pyre.inventory.float("z_lith", default=0.014)
z_410 = pyre.inventory.float("z_410", default=0.06435)
z_lmantle = pyre.inventory.float("z_lmantle", default=0.105)
z_cmb = pyre.inventory.float("z_cmb", default=0.439)
# version
__id__ = "$Id$"
# End of file