diff -u a/Xspec/src/XSFunctions/sumdem.cxx b/Xspec/src/XSFunctions/sumdem.cxx
--- a/Xspec/src/XSFunctions/sumdem.cxx 2019-05-17 14:35:25.000000000 -0500
+++ b/Xspec/src/XSFunctions/sumdem.cxx 2022-02-16 10:42:36.946827375 -0600
@@ -14,7 +14,7 @@
RealArray& fluxArray, RealArray& fluxErrArray);
-void sumdem(int itype, int swtch, float* ear, int ne, float* abun,
+extern "C" void sumdem(int itype, int swtch, float* ear, int ne, float* abun,
float dens, float z, int ninputt, float* inputt, float* dem,
int ifl, bool qtherm, float velocity, float* photar, float* photer,
int* status);
diff -u a/Xspec/src/XSFunctions/xsmaug.cxx b/Xspec/src/XSFunctions/xsmaug.cxx
--- a/Xspec/src/XSFunctions/xsmaug.cxx 2019-05-17 14:35:25.000000000 -0500
+++ b/Xspec/src/XSFunctions/xsmaug.cxx 2022-02-16 10:42:28.949802306 -0600
@@ -49,10 +49,10 @@
// REAL ear(0:*), abun(*), inputt(*), dem(*), photar(*)
// REAL dens, z
-void sumdem_(int& itype, int& flag, float* ear, int& ne,
- float* abun, float& dens, float& z, int& ninputt,
- float* inputt, float* dem, int& ifl, int& qtherm,
- float& velocity, float* photar, int& status);
+void sumdem(int itype, int flag, float* ear, int ne,
+ float* abun, float dens, float z, int ninputt,
+ float* inputt, float* dem, int ifl, bool qtherm,
+ float velocity, float* photar, float* photer, int* status);
}
/*---------------------------------------------------------------------------*/
@@ -325,7 +325,7 @@
const double MIN2RAD = 3437.75;
int i, ie, iel;
int mesh2, meshpt;
- int itype, mode, no_el, status;
+ int itype, mode, no_el, status=0;
double angfac, da, ei, evol, inner, H0, L0, outer, q0, z, zfac;
double a1, a2, Ab_shell, ei_shell, h1, h2, r1, r2, rcutoff, t1, t2, w1, w2, w12;
double elden, norm;
@@ -487,10 +487,10 @@
static float ONE(1.);
static int one(1);
float fz (z);
- static int NOT(0);
+ static bool NOT(false);
static float ZERO(0.0);
- sumdem_(itype, mode, fear, ne, abund, ONE , fz, one, kT_shell, dem, ifl,
- NOT, ZERO, phoshell, status);
+ sumdem(itype, mode, fear, ne, abund, ONE , fz, one, kT_shell, dem, ifl,
+ NOT, ZERO, phoshell, 0, &status);
/* ADD THE SHELL CONTRIBUTION TO THE OUTPUT SPECTRUM,
* CORRECTING IT FOR <angfac>, DISTANCE AND REDSHIFT,
