--- trunk/libf/phylmd/radiornpb.f90	2008/08/05 13:31:32	17
+++ trunk/phylmd/radiornpb.f	2018/02/05 10:39:38	254
@@ -16,9 +16,9 @@
     ! Le pas de temps "pdtphys" est supposé beaucoup plus petit que la
     ! constante de temps de décroissance.
 
-    use dimens_m, only: llm
-    use dimphy, only: klon, nbtr
-    use numer_rec, only: assert
+    use dimens_m, only: llm, nqmx
+    use dimphy, only: klon
+    use nr_util, only: assert
 
     REAL, intent(in):: tr_seri(:, :, :), pdtphys, tautr(:)
     real radiornpb(klon, llm, 2)
@@ -28,8 +28,8 @@
 
     !-----------------------------------------------
 
-    call assert(shape(tr_seri) == (/klon, llm, nbtr/), "radiornpb tr_seri")
-    call assert(size(tautr) == nbtr, "radiornpb tautr")
+    call assert(shape(tr_seri) == (/klon, llm, nqmx - 2/), "radiornpb tr_seri")
+    call assert(size(tautr) == nqmx - 2, "radiornpb tautr")
 
     DO it = 1, 2
        IF (tautr(it) > 0.) THEN