!!****f* source/Simulation/Simulation_initSpecies
!!
!! NAME
!!
!!  Simulation_initSpecies
!!
!! SYNOPSIS
!!
!!  Simulation_initSpecies()
!!
!! DESCRIPTION
!!
!!  This routine will initialize the species and species values needed for a 
!!  given setup.   The user should add the 
!!  implementation of this routine to the setups directory of a simulation 
!!  that needs to use the multispecies capabilities of the code.
!!
!!  There two general purpose implementations available in the code, one which sets standard  
!!  isotope properties for the nuclear burning source terms, and another one for the 
!!  Ionization source term.
!!
!!  This routine is called from Multispecies_init, and is called BEFORE
!!  the call to Simulation_init.  
!!
!! SEE ALSO
!!  Multispecies_init
!!  Simulation/SimulationComposition/Simulation_initSpecies
!!
!!***

subroutine Simulation_initSpecies()
  use Multispecies_interface, ONLY : Multispecies_setProperty

  implicit none
#include "Flash.h"
#include "Multispecies.h"
#include "Eos.h"

  call Multispecies_setProperty(CHAM_SPEC, A, 4.)
  call Multispecies_setProperty(CHAM_SPEC, Z, 2.)
  call Multispecies_setProperty(CHAM_SPEC, MS_ZMIN, 0.02)

  call Multispecies_setProperty(PLAS_SPEC, A, 5.7)
  call Multispecies_setProperty(PLAS_SPEC, Z, 3.5)
  call Multispecies_setProperty(PLAS_SPEC, MS_ZMIN, 0.02)

  call Multispecies_setProperty(GOLD_SPEC, A, 26.0)
  call Multispecies_setProperty(GOLD_SPEC, Z, 13.0)
  call Multispecies_setProperty(GOLD_SPEC, MS_ZMIN, 0.02)

  call Multispecies_setProperty(PART_SPEC, A, 26.0)
  call Multispecies_setProperty(PART_SPEC, Z, 13.0)
  call Multispecies_setProperty(PART_SPEC, MS_ZMIN, 0.02)

  call Multispecies_setProperty(TARG_SPEC, A, 5.7)
  call Multispecies_setProperty(TARG_SPEC, Z, 3.5)
  call Multispecies_setProperty(TARG_SPEC, MS_ZMIN, 0.02)

end subroutine Simulation_initSpecies