--- trunk/libf/phylmd/fisrtilp.f90	2010/06/08 15:37:21	35
+++ trunk/libf/phylmd/fisrtilp.f90	2010/12/02 17:11:04	36
@@ -15,7 +15,7 @@
   USE yoethf
   USE fcttre
   USE comfisrtilp
-  use numer_rec, only: erf
+  use numer_rec, only: nr_erf
 
   IMPLICIT NONE
 
@@ -338,10 +338,10 @@
               zpdf_b(i) = zpdf_k(i)/(2.*sqrt(2.))
               zpdf_e1(i) = zpdf_a(i) - zpdf_b(i)
               zpdf_e1(i) = sign(min(abs(zpdf_e1(i)),5.),zpdf_e1(i))
-              zpdf_e1(i) = 1. - erf(zpdf_e1(i))
+              zpdf_e1(i) = 1. - nr_erf(zpdf_e1(i))
               zpdf_e2(i) = zpdf_a(i) + zpdf_b(i)
               zpdf_e2(i) = sign(min(abs(zpdf_e2(i)),5.),zpdf_e2(i))
-              zpdf_e2(i) = 1. - erf(zpdf_e2(i))
+              zpdf_e2(i) = 1. - nr_erf(zpdf_e2(i))
               IF (zpdf_e1(i)<1.E-10) THEN
                  rneb(i,k) = 0.
                  zqn(i) = zqs(i)