Index: /branches/eam_branches/20090715/psModules/src/objects/models/pmModel_PS1_V1.c
===================================================================
--- /branches/eam_branches/20090715/psModules/src/objects/models/pmModel_PS1_V1.c	(revision 25677)
+++ /branches/eam_branches/20090715/psModules/src/objects/models/pmModel_PS1_V1.c	(revision 25678)
@@ -323,8 +323,7 @@
 
     // choose a z value guaranteed to be beyond our limit
-    float z0 = pow((1.0 / limit), (1.0 / ALPHA));
-    float z1 = (1.0 / limit) / PAR[PM_PAR_7];
-    z1 = PS_MAX (z0, z1);
-    z0 = 0.0;
+    float z0 = 0.0;
+    float z1 = pow((1.0 / limit), (1.0 / ALPHA));
+    if (PAR[PM_PAR_7] < 0.0) z1 *= 2.0;
 
     // perform a type of bisection to find the value