--- trunk/libf/phylmd/fisrtilp.f90	2009/07/31 15:18:47	22
+++ trunk/libf/phylmd/fisrtilp.f90	2011/01/06 17:52:19	38
@@ -11,11 +11,11 @@
   USE dimens_m
   USE dimphy
   USE tracstoke
-  USE yomcst
-  USE yoethf
+  USE suphec_m
+  USE yoethf_m
   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)