Changeset 9529 for trunk/psphot/src/models/pmModel_TGAUSS.c

Timestamp:

Oct 12, 2006, 4:36:54 PM (20 years ago)

Author:

rhl

Message:

Renamed PM_PAR_FLUX to PM_PAR_I0

File:

• trunk/psphot/src/models/pmModel_TGAUSS.c (modified) (4 diffs)

trunk/psphot/src/models/pmModel_TGAUSS.c

@@ -35 +35 @@
 
     psF32 r  = 1.0 / (er + PAR[PM_PAR_7]*z + pow(z, TG_S));    // (1/R)
-    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) {
-        psF32 t = PAR[PM_PAR_FLUX]*r*r; // df/dR
+        psF32 t = PAR[PM_PAR_I0]*r*r;   // df/dR
         psF32 q = t*(TRF*er + PAR[PM_PAR_7] + TG_S*pow(z, dTG_S));  // (df/dR)(dR/dz)
 
@@ -127 +127 @@
     norm *= 0.01;
     
-    psF64 Flux = PAR[PM_PAR_FLUX] * Area * norm;
+    psF64 Flux = PAR[PM_PAR_I0] * Area * norm;
 
     return(Flux);
@@ -140 +140 @@
 
     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);
 
     // if Sx == Sy, sigma = Sx == Sy
     psF64 sigma = hypot (1.0 / PAR[PM_PAR_SXX], 1.0 / PAR[PM_PAR_SYY]) / sqrt(2.0);
     psF64 dz = 1.0 / (2.0 * sigma*sigma);
-    psF64 limit = flux / PAR[PM_PAR_FLUX];
+    psF64 limit = flux / PAR[PM_PAR_I0];
 
     // we can do this much better with intelligent choices here
@@ -193 +193 @@
     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);
 
     if (status) return true;