#!/usr/bin/python
# Support for python2
from __future__ import print_function
# Importing Numpy (math, arrays, etc...)
import numpy as np
# Importing Matplotlib (plotting)
import matplotlib.pyplot as plt

# Avogadros number
NA = 6.02214076e23
KB = 1.380650524e-23

def calc_reduced_T(Ta, eps):
    """ Calculate reduced temperature
    """
    Tstar = np.zeros_like(Ta)
    Tstar = Ta/eps
    return Tstar

def calc_reduced_rho(rhoa, sigma):
    """ Calculate reduced density
    """
    rhoStar = np.zeros_like(rhoa)
    rhoStar = sigma**3*NA*rhoa
    return rhoStar

def calc_real_T(Tstar, eps):
    """ Calculate temperature from reduced temperature
    """
    Ta = np.zeros_like(Tstar)
    Ta = Tstar*eps
    return Ta

def calc_real_rho(rhoStar, sigma):
    """ Calculate density from reduced density
    """
    rhoa = np.zeros_like(rhoStar)
    rhoa = rhoStar/(sigma**3*NA)
    return rhoa

def calc_reduced_entropy(s):
    """ Calclate reduced entropy
    """
    sStar = np.zeros_like(s)
    sStar = s/(NA*KB)
    return sStar

def calc_reduced_P(Pa, eps, sigma):
    """ Calculate reduced pressure
    """
    Pstar = np.zeros_like(Pa)
    Pstar = Pa*sigma**3/eps/KB
    return Pstar

def calc_real_P(Pstar, eps, sigma):
    """ Calculate reduced pressure
    """
    Pa = np.zeros_like(Pstar)
    Pa = Pstar*eps*KB/sigma**3
    return Pa

def calc_reduced_heat_capacity(C):
    """ Calclate reduced heat capcity
    """
    CStar = np.zeros_like(C)
    CStar = C/(NA*KB)
    return CStar

#!/usr/bin/python
# Support for python2
from __future__ import print_function
# Importing Numpy (math, arrays, etc...)
import numpy as np
# Importing Matplotlib (plotting)
import matplotlib.pyplot as plt

# Avogadros number
NA = 6.02214076e23
KB = 1.380650524e-23

def calc_reduced_T(Ta, eps):
    """ Calculate reduced temperature
    """
    Tstar = np.zeros_like(Ta)
    Tstar = Ta/eps
    return Tstar

def calc_reduced_rho(rhoa, sigma):
    """ Calculate reduced density
    """
    rhoStar = np.zeros_like(rhoa)
    rhoStar = sigma**3*NA*rhoa
    return rhoStar

def calc_real_T(Tstar, eps):
    """ Calculate temperature from reduced temperature
    """
    Ta = np.zeros_like(Tstar)
    Ta = Tstar*eps
    return Ta

def calc_real_rho(rhoStar, sigma):
    """ Calculate density from reduced density
    """
    rhoa = np.zeros_like(rhoStar)
    rhoa = rhoStar/(sigma**3*NA)
    return rhoa

def calc_reduced_entropy(s):
    """ Calclate reduced entropy
    """
    sStar = np.zeros_like(s)
    sStar = s/(NA*KB)
    return sStar

def calc_reduced_P(Pa, eps, sigma):
    """ Calculate reduced pressure
    """
    Pstar = np.zeros_like(Pa)
    Pstar = Pa*sigma**3/eps/KB
    return Pstar

def calc_real_P(Pstar, eps, sigma):
    """ Calculate reduced pressure
    """
    Pa = np.zeros_like(Pstar)
    Pa = Pstar*eps*KB/sigma**3
    return Pa

def calc_reduced_heat_capacity(C):
    """ Calclate reduced heat capcity
    """
    CStar = np.zeros_like(C)
    CStar = C/(NA*KB)
    return CStar

#!/usr/bin/python # Support for python2 # Importing Numpy (math, arrays, etc...) # Importing Matplotlib (plotting) # Avogadros number Calculate reduced temperature Calculate reduced density Calculate temperature from reduced temperature Calculate density from reduced density Calclate reduced entropy Calculate reduced pressure Calculate reduced pressure Calclate reduced heat capcity