Index: trunk/libf/dyn3d/disvert_noterre.F
===================================================================
--- trunk/libf/dyn3d/disvert_noterre.F	(revision 109)
+++ trunk/libf/dyn3d/disvert_noterre.F	(revision 110)
@@ -199,16 +199,16 @@
 c     Specifique.  Ce choix est spécifié ici ET dans exner_hyb.F
 
-      DO l = 1, llm
+      DO l = 1, llm-1
        aps(l) =  0.5 *( ap(l) +ap(l+1)) 
        bps(l) =  0.5 *( bp(l) +bp(l+1)) 
       ENDDO
      
-      if (hybrid) then
+c     if (hybrid) then
          aps(llm) = aps(llm-1)**2 / aps(llm-2) 
          bps(llm) = 0.5*(bp(llm) + bp(llm+1))
-      else
-         bps(llm) = bps(llm-1)**2 / bps(llm-2) 
-         aps(llm) = 0.
-      end if
+c     else
+c        bps(llm) = bps(llm-1)**2 / bps(llm-2) 
+c        aps(llm) = 0.
+c     end if
 
       PRINT *,' BPs '
Index: trunk/libf/dyn3d/disvert_terre.F90
===================================================================
--- trunk/libf/dyn3d/disvert_terre.F90	(revision 109)
+++ trunk/libf/dyn3d/disvert_terre.F90	(revision 110)
@@ -140,3 +140,3 @@
   write(lunout, *) presnivs
 
-END SUBROUTINE disvert
+END SUBROUTINE disvert_terre
Index: trunk/libf/dyn3d/inidissip.F
===================================================================
--- trunk/libf/dyn3d/inidissip.F	(revision 109)
+++ trunk/libf/dyn3d/inidissip.F	(revision 110)
@@ -185,5 +185,5 @@
 
          write(*,*) 'Dissipation : '
-         write(*,*) 'Multiplication de la dissipation en altitude :',
+         write(*,*) 'Multiplication de la dissipation en altitude :'
          write(*,*) '  dissip_fac_mid =', dissip_fac_mid
 
Index: trunk/libf/dyn3d/leapfrog.F
===================================================================
--- trunk/libf/dyn3d/leapfrog.F	(revision 109)
+++ trunk/libf/dyn3d/leapfrog.F	(revision 110)
@@ -470,5 +470,6 @@
 c      -------------------
          IF (ok_strato) THEN
-           CALL top_bound( vcov,ucov,teta,masse,dutop,dvtop,dtetatop)
+           CALL top_bound( vcov,ucov,teta,phi,masse,
+     $                     dutop,dvtop,dtetatop)
 c dqtop=0, dptop=0
            CALL addfi( dtphys, leapf, forward   ,
@@ -516,5 +517,6 @@
         ! Sponge layer (if any)
         IF (ok_strato) THEN
-          CALL top_bound(vcov,ucov,teta,masse,dutop,dvtop,dtetatop)
+          CALL top_bound(vcov,ucov,teta,phi,
+     $                   masse,dutop,dvtop,dtetatop)
 c dqtop=0, dptop=0
           CALL addfi( dtvr, leapf, forward   ,
Index: trunk/libf/dyn3d/top_bound.F
===================================================================
--- trunk/libf/dyn3d/top_bound.F	(revision 109)
+++ trunk/libf/dyn3d/top_bound.F	(revision 110)
@@ -1,3 +1,3 @@
-      SUBROUTINE top_bound( vcov,ucov,teta,masse, du,dv,dh )
+      SUBROUTINE top_bound( vcov,ucov,teta,phi,masse, du,dv,dh )
       IMPLICIT NONE
 c
@@ -35,4 +35,5 @@
 
       REAL ucov(iip1,jjp1,llm),vcov(iip1,jjm,llm),teta(iip1,jjp1,llm)
+      REAL phi(iip1,jjp1,llm)                  ! geopotentiel
       REAL masse(iip1,jjp1,llm)
       REAL dv(iip1,jjm,llm),du(iip1,jjp1,llm),dh(iip1,jjp1,llm)
@@ -77,11 +78,9 @@
          print*,'TOP_BOUND mode',mode_top_bound
          print*,'Coeffs pour la couche eponge a l equateur'
-         print*,'p (Pa)  z(km)  tau (s)   dt*rdamp'
+         print*,'p (Pa)  z(km)  tau (s)'
          do l=1,llm
            if (rdamp(l).ne.0.) then
             zkm        = phi(iip1/2,jjp1/2,l)/(1000*g)
-          print*,presnivs(l),zkm,
-     .          1./rdamp(l),
-     .          dt*rdamp(l)
+          print*,presnivs(l),zkm,1./rdamp(l)
            endif
          enddo