anchor.diff:      return
     end function network_species_index
   
  -
  -  function network_reaction_index(name)
  -    
  -    character(len=*) :: name
  -    integer :: network_reaction_index, n
  -
  -    network_reaction_index = -1
  -
  -    do n = 1, nrat
  -       if (name == reac_names(n)) then
  -          network_reaction_index = n
  -          exit
  -       endif
  -    enddo
  -
  -    return
  -  end function network_reaction_index
  -
  -
     subroutine network_finalize()
   
     end subroutine network_finalize
  
positive.code:     ! done initializing
      network_initialized = .true.
  
      ! rpar is VODE's way of passing information into the RHS and                
      ! jacobian routines.  Here we initialize some indices to make               
      ! sense of what is stored in the rpar() array.                              
      call init_rpar_indices(nrat, nspec)
  
    end subroutine network_init
  
    
  
negative.code:     short_spec_names(ihe4) = "He4"
      short_spec_names(ih1) = "H1"    
  
      ! set the species properties
      aion(ic12) = TWELVE
      aion(io14) = 14.0_dp_t
  