# Path: DiscEvolution/constants.py
# G = 1.
# AU = 1.496e13


import numpy as np
from scipy import optimize
from .constants import AU, sig_H2, m_H, yr, Msun
    @property
    def T(self):
        """Temperature"""
        return self._eos.T

    @property
    def mu(self):
        return self._eos.mu
    
    @property
    def H(self):
        """Scale-height"""
        return self._eos.H

    @property
    def h(self):
        """Aspect ratio"""
        return self.H/self.R

    @property
    def H_edge(self):
        """Scale-height at cell edge"""
        return self._eos.H_edge

    @property
    def P(self):
        return self.midplane_gas_density * self.cs**2

    @property
    def midplane_gas_density(self):

# Path: DiscEvolution/constants.py
# G = 1.
# AU = 1.496e13


# Path: DiscEvolution/disc.py
import numpy as np
from scipy import optimize
from .constants import AU, sig_H2, m_H, yr, Msun

    @property
    def T(self):
        """Temperature"""
        return self._eos.T

    @property
    def mu(self):
        return self._eos.mu
    
    @property
    def H(self):
        """Scale-height"""
        return self._eos.H

    @property
    def h(self):
        """Aspect ratio"""
        return self.H/self.R

    @property
    def H_edge(self):
        """Scale-height at cell edge"""
        return self._eos.H_edge

    @property
    def P(self):
        return self.midplane_gas_density * self.cs**2

    @property
    def midplane_gas_density(self):

        return self.Sigma_G / (np.sqrt(2*np.pi) * self.H * AU)