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:   character (len= 5), save :: short_spec_names(nspec)
    character (len= 5), save :: short_aux_names(naux)
  
    character (len=10), save :: reac_names(nrat)
  
    real(kind=dp_t), save :: aion(nspec), zion(nspec), ebin(nspec)
  
    logical, save :: network_initialized = .false.
  
negative.code: !
  !  network_species_index  -- return the index of the species given its name
  !
  
  module network
  
  