/* include for rxn.c */
#include <stdio.h>
#define SIZE 100  /* max line length for a reaction record */
#define SPESIZE 4000 /* max length of species list */
#define MSPR 10  /* max species per reaction */
#define MNR  100 /* max no. reactions */
#define MNS  100 /* max no. species */
typedef struct {
  float k;  /* rate constant */
  int s[MSPR], n; /* list of participating species and count */
  float rate;  /* computed reaction rate */
} RXN; 
RXN R[MNR];  /* reactions */
double *Y, *DY;  /* concentrations */
char *SPECIES;  /* list of species */
int RP_FLAG;  /* reactant vs product flag */
int NR;  /* no. reactions read in */
int NS;  /* no. species used */
extern void bail();
FILE *params, *mech, *Ywrite, *Swrite;
char *parameters, *mechanism, *results, *select, *Species;
char *altlist;