Re: [FLASH-USERS] About abundances . . .

From: Salvatore Orlando <orlando@oapa23.astropa.unipa.it>
Date: Thu May 22 2003 - 07:59:46 CDT

Claudio Dalla Vecchia wrote:
> All this, supposing H is neutral hydrogen and the fraction of electrons
> takes account of HII. There is no HII in the routines - at least, I have
> not seen it.

Sorry!
I forgot to answer the above question: the ioniz module assumes that the
hydrogen is *fully* ionized.
 
Finally, in the routine init_mat.F90 it is assumed that N_e = 1.213 N_H
following Cox (2000 ``Allen's Astrophysical Quantities'', 4th ed.~
Publisher: New York: AIP Press; Springer, 2000. Edited by Arthur N.
Cox). If you wish a different ratio N_e/N_H, you need to modify
init_mat.F90 (line 45).

By the way, I have noticed a bug on line 45 of init_mat.F90:

Line
      integer, parameter :: ep_ratio = 1.213

should be

      real, parameter :: ep_ratio = 1.213

Best wishes
Salvatore 

-- 
  I N A F   Istituto Nazionale di Astrofisica
  Osservatorio Astronomico di Palermo "G.S. Vaiana"
  Piazza del Parlamento 1, I-90134 Palermo, Italy
  Tel. +39-91-233 423, Fax. +39-91-233 444
  Web http://www.astropa.unipa.it/~orlando
Received on Thu May 22 08:00:04 2003

This archive was generated by hypermail 2.1.8 : Thu Aug 31 2006 - 21:20:48 CDT