Index: trunk/psphot/src/models/pmModel_SGAUSS.c
===================================================================
--- trunk/psphot/src/models/pmModel_SGAUSS.c	(revision 8882)
+++ trunk/psphot/src/models/pmModel_SGAUSS.c	(revision 9529)
@@ -37,9 +37,9 @@
     psF32 p  = pow(z, PAR[PM_PAR_7] - 1.0);
     psF32 r  = 1.0 / (1 + z*p + pr*pr3);
-    psF32 f  = PAR[PM_PAR_FLUX]*r + PAR[PM_PAR_SKY];
+    psF32 f  = PAR[PM_PAR_I0]*r + PAR[PM_PAR_SKY];
 
     if (deriv != NULL) {
         // note difference from a pure gaussian: q = params->data.F32[1]*r
-        psF32 t = PAR[PM_PAR_FLUX]*r*r;
+        psF32 t = PAR[PM_PAR_I0]*r*r;
 	psF32 q = t*(PAR[PM_PAR_7]*p + 4*PAR[PM_PAR_8]*pr3);
 
@@ -245,5 +245,5 @@
     norm *= 0.01;
     
-    psF64 Flux = PAR[PM_PAR_FLUX] * Area * norm;
+    psF64 Flux = PAR[PM_PAR_I0] * Area * norm;
 
     return(Flux);
@@ -262,6 +262,6 @@
 
     if (flux <= 0) return (1.0);
-    if (PAR[PM_PAR_FLUX] <= 0) return (1.0);
-    if (flux >= PAR[PM_PAR_FLUX]) return (1.0);
+    if (PAR[PM_PAR_I0] <= 0) return (1.0);
+    if (flux >= PAR[PM_PAR_I0]) return (1.0);
 
     // convert Sx,Sy,Sxy to major/minor axes
@@ -272,5 +272,5 @@
     axes = EllipseShapeToAxes (shape);
     psF64 dr = 1.0 / axes.major;
-    psF64 limit = flux / PAR[PM_PAR_FLUX];
+    psF64 limit = flux / PAR[PM_PAR_I0];
 
     // XXX : we can do this faster with an intelligent starting choice
@@ -329,6 +329,6 @@
     status = true;
     status &= (dP < 0.5);
-    status &= (PAR[PM_PAR_FLUX] > 0);
-    status &= ((dPAR[PM_PAR_FLUX]/PAR[PM_PAR_FLUX]) < 0.5);
+    status &= (PAR[PM_PAR_I0] > 0);
+    status &= ((dPAR[PM_PAR_I0]/PAR[PM_PAR_I0]) < 0.5);
     status &= (fabs(PAR[PM_PAR_8]) < 0.5);
     status &= (dPAR[PM_PAR_8] < 0.1);