--- trunk/libf/dyn3d/inter_barxy.f90	2008/02/27 13:16:39	3
+++ trunk/libf/dyn3d/inter_barxy.f90	2010/03/03 13:23:49	24
@@ -16,8 +16,7 @@
 
     ! Author: P. Le Van
 
-    use nrutil, only: assert_eq, assert
-
+    use numer_rec, only: assert_eq, assert
     use dimens_m, only: iim, jjm
     use comgeom, only: aire_2d, apoln, apols
 
@@ -115,7 +114,7 @@
     !      ( L'indice 1 correspond a l'interface mailLE 1 / maille 2)
     !      ( Les abscisses sont exprimées en degres)
 
-    use nrutil, only: assert_eq
+    use numer_rec, only: assert_eq
 
     IMPLICIT NONE
 
@@ -300,7 +299,7 @@
     ! Version unidimensionnelle, en latitude.
     ! L'indice 1 correspond à l'interface maille 1 -- maille 2.
 
-    use nrutil, only: assert
+    use numer_rec, only: assert
 
     IMPLICIT NONE
 
@@ -349,7 +348,7 @@
           dym        = 0.
           jmod       = jmod + 1
        ELSE
-          ! yjmod(jmod) == yjdat(jdat)
+          ! {yjmod(jmod) == yjdat(jdat)}
           dy         = yjmod(jmod) - y0
           dym        = dym + dy
           inter_bary(jmod) = (inter_bary(jmod) + dy * fdat(jdat)) / dym
@@ -377,7 +376,7 @@
     ! reverses their order.
     ! Finally, the procedure adds 90° as the last value of the array.
 
-    use nrutil, only: assert_eq
+    use numer_rec, only: assert_eq
     use comconst, only: pi
 
     IMPLICIT NONE