my ($initial_enthalpy, $initial_entropy) = get_initials (\@terminal, \@initial_enths, \@initial_entros);


sub get_initials {

  my ($terminal, $initial_enths, $initial_entros) = @_;
  my ($initial_entropy, $initial_enthalpy);

  if (($terminal[0] =~ /^G/)  || ($terminal[0] =~ /^C/))    {
    $initial_enthalpy += initial_enths[1];
    $initial_entropy += $initial_entros[1];
  }

  return ($initial_enthalpy, $initial_entropy);

}