Index: trunk/LMDZ.MARS/libf/phymars/pbl_parameters.F
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/pbl_parameters.F	(revision 1390)
+++ trunk/LMDZ.MARS/libf/phymars/pbl_parameters.F	(revision 1393)
@@ -218,6 +218,5 @@
 !------------------------------------------------------------------------
 
-
-! u* theta* computation:
+! u* theta* computation
 
       DO ig=1,ngrid
@@ -225,5 +224,4 @@
            ustar(ig)=0.
            tstar(ig)=0.
-
          ELSE
            ustar(ig)=sqrt(pcdv(ig))
@@ -231,7 +229,8 @@
            tstar(ig)=-pcdh(ig)*(pts(ig)-ph(ig,1))
      &        /sqrt(pcdv(ig))
-
-         ENDIF
-      ENDDO
+         ENDIF
+      ENDDO
+
+! Interpolation:
 
       DO ig=1,ngrid
@@ -239,5 +238,10 @@
            u_out(ig,n)=0.
            Teta_out(ig,n)=pts(ig)
-         ELSE
+
+        ELSE
+          IF (rib(ig) .ge. ric) THEN ! ustar=tstar=0  (and fm=fh=0)
+           u_out(ig,n)=0
+           Teta_out(ig,n)=pts(ig)
+          ELSE
            u_out(ig,n)= ustar(ig)*log(zout/pz0(ig))/
      &(karman*sqrt(fm(ig)))
@@ -246,5 +250,5 @@
      & (pz0tcomp(ig)))/
      &(karman*fh(ig)))
-
+          ENDIF
         ENDIF
 
@@ -271,4 +275,5 @@
       ENDDO   !of n=1,n_out
 
+
 !------------------------------------------------------------------------
 !------------------------------------------------------------------------
@@ -310,4 +315,5 @@
          pbl_teta(ig) = ph(ig,ii(1))
       ENDDO
+
 
 !------------------------------------------------------------------------