--- trunk/libf/phylmd/fisrtilp.f90	2009/07/31 15:18:47	22
+++ trunk/libf/phylmd/fisrtilp.f90	2012/04/20 14:58:43	61
@@ -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_95, only: nr_erf
 
   IMPLICIT NONE
 
@@ -24,8 +24,8 @@
   REAL, INTENT (IN) :: & ! intervalle du temps (s)                
        dtime
   REAL, INTENT (IN) :: paprs(klon,klev+1) ! pression a inter-couche   
-  REAL, INTENT (IN) :: pplay(klon,klev) ! pression au milieu de couche
-  REAL t(klon,klev) ! temperature (K)                               
+  REAL, INTENT (IN):: pplay(klon,klev) ! pression au milieu de couche
+  REAL, INTENT (IN):: t(klon,klev) ! temperature (K)
   REAL q(klon,klev) ! humidite specifique (kg/kg)                   
   REAL d_t(klon,klev) ! incrementation de la temperature (K)        
   REAL d_q(klon,klev) ! incrementation de la vapeur d'eau           
@@ -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)