/*outsum sum outward continuum beams */
void outsum(double *outtot, 
  double *outin, 
  double *outout);