Setting Program Defaults

There are some parameters used in the program whose default values need to be set. I will probably build their values into the routines at some unspecified later date. Until then you must include the following lines of code in your application:

double
   cf, 		/* convergence factor */
   accuracy;

double 
   e_surface,  	/* value of energy density at star's surface */
   p_surface,   /* value of pressure at star's surface */
   enthalpy_min; /* minimum value of enthalpy */

  /* set program defaults */
  cf=1.0;
  accuracy=1e-5;    
 
  if(strcmp(eos_type,"tab")==0) {
    e_surface=7.8*C*C*KSCALE;
    p_surface=1.01e8*KSCALE;
    enthalpy_min=1.0/(C*C);
  }
  else{
    e_surface=0.0;
    p_surface=0.0;
    enthalpy_min=0.0;
  }