Index: /trunk/psastro/src/psastro.h
===================================================================
--- /trunk/psastro/src/psastro.h	(revision 17555)
+++ /trunk/psastro/src/psastro.h	(revision 17556)
@@ -23,4 +23,5 @@
     double mPeak;			// mag of peak bin 
     int nPeak;				// # of stars in peak bin
+    int sPeak;				// sum of stars to peak bin
 } pmLumFunc;
 
Index: /trunk/psastro/src/psastroLuminosityFunction.c
===================================================================
--- /trunk/psastro/src/psastroLuminosityFunction.c	(revision 17555)
+++ /trunk/psastro/src/psastroLuminosityFunction.c	(revision 17556)
@@ -22,4 +22,5 @@
   func->mPeak = 0.0;
   func->nPeak = 0;
+  func->sPeak = 0;
 
   return func;
@@ -62,8 +63,10 @@
   int iPeak = 0;
   int nPeak = 0;
+  int sPeak = 0;
   for (int i = 0; i < nMags->n; i++) {
     if (nMags->data.F32[i] < nPeak) continue;
     iPeak = i;
     nPeak = nMags->data.F32[i];
+    sPeak += nPeak;
   }
 
@@ -122,4 +125,5 @@
   lumFunc->mPeak = mMin + (iPeak + 0.5)*dMag;
   lumFunc->nPeak = nPeak;
+  lumFunc->sPeak = sPeak;
 
   psFree (lnMag);