########################################################################
# Program code
########################################################################
# Import general numpy and plotting stuff
from numpy import *
import matplotlib.pyplot as plt
# Import the line profile function and the emitter data class from
# DESPOTIC
from despotic.lineProfLTE import lineProfLTE
from despotic.emitterData import emitterData
# Read the data files for HCN and N2H+
HCN=emitterData('hcn')
N2Hp=emitterData('n2h+')
# Core radius; set to 0.02 pc
R = 0.02*3.09e18
# Abundances and absolute number densities of HCN and N2H+, in cm^-3;
# abundances are taken from the models of Lee et al. (2004, ApJ, 617,
# 360).
xHCN = 2e-8
xN2Hp = 2e-9
nH = 3e6
nHCN = nH*xHCN
nN2Hp = nH*xN2Hp
# Define a temperature profile so that we can get a nice p Cygni
########################################################################
# Program code
########################################################################
# Import general numpy and plotting stuff
from numpy import *
import matplotlib.pyplot as plt
# Import the line profile function and the emitter data class from
# DESPOTIC
from despotic.lineProfLTE import lineProfLTE
from despotic.emitterData import emitterData
# Read the data files for CS and C^34S
cs = emitterData('cs')
c34s = emitterData('c34s')
# Core radius; set to 0.02 pc
R = 0.02 * 3.09e18
# Number densities of CS and C^34S, in cm^-3; the abundance ratio is
# set to the terrestrial value
ncs = 1e-1
nc34s = ncs / 22.0
# Define a temperature profile so that we can get a nice p Cygni
# profile. The choice here is somewhat arbitrary. For simplicity we'll
# just take 8 K at large radii, added to a Gaussian component rising
# to 20 K at small radii. This function takes as an argument the
########################################################################
# Program code
########################################################################
# Import general numpy and plotting stuff
from numpy import *
import matplotlib.pyplot as plt
# Import the line profile function and the emitter data class from
# DESPOTIC
from despotic.lineProfLTE import lineProfLTE
from despotic.emitterData import emitterData
# Read the data files for CS and C^34S
cs=emitterData('cs')
c34s=emitterData('c34s')
# Core radius; set to 0.02 pc
R = 0.02*3.09e18
# Number densities of CS and C^34S, in cm^-3; the abundance ratio is
# set to the terrestrial value
ncs = 1e-1
nc34s = ncs / 22.0
# Define a temperature profile so that we can get a nice p Cygni
# profile. The choice here is somewhat arbitrary. For simplicity we'll
# just take 8 K at large radii, added to a Gaussian component rising
# to 20 K at small radii. This function takes as an argument the
# radius r normalized to the core radius, and returns the temperature
########################################################################
# Program code
########################################################################
# Import general numpy and plotting stuff
from numpy import *
import matplotlib.pyplot as plt
# Import the line profile function and the emitter data class from
# DESPOTIC
from despotic.lineProfLTE import lineProfLTE
from despotic.emitterData import emitterData
# Read the data files for HCN and N2H+
HCN = emitterData("hcn")
N2Hp = emitterData("n2h+")
# Core radius; set to 0.02 pc
R = 0.02 * 3.09e18
# Abundances and absolute number densities of HCN and N2H+, in cm^-3;
# abundances are taken from the models of Lee et al. (2004, ApJ, 617,
# 360).
xHCN = 2e-8
xN2Hp = 2e-9
nH = 3e6
nHCN = nH * xHCN
nN2Hp = nH * xN2Hp
# Define a temperature profile so that we can get a nice p Cygni
