######################################################################## # 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