Changeset 1393 for trunk/LMDZ.MARS/libf/phymars

Timestamp:

Mar 9, 2015, 4:02:47 PM (10 years ago)

Author:

aslmd

Message:

bug fix pbl_parameter in cases of specific Ri

File:

• trunk/LMDZ.MARS/libf/phymars/pbl_parameters.F (modified) (7 diffs)

trunk/LMDZ.MARS/libf/phymars/pbl_parameters.F

@@ -218 +218 @@
 !------------------------------------------------------------------------
 
-
-! u* theta* computation:
+! u* theta* computation
 
       DO ig=1,ngrid
@@ -225 +224 @@
            ustar(ig)=0.
            tstar(ig)=0.
-
          ELSE
            ustar(ig)=sqrt(pcdv(ig))
@@ -231 +229 @@
            tstar(ig)=-pcdh(ig)*(pts(ig)-ph(ig,1))
      &        /sqrt(pcdv(ig))
-
-         ENDIF
-      ENDDO
+         ENDIF
+      ENDDO
+
+! Interpolation:
 
       DO ig=1,ngrid
@@ -239 +238 @@
            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 +250 @@
      & (pz0tcomp(ig)))/
      &(karman*fh(ig)))
-
+          ENDIF
         ENDIF
 
@@ -271 +275 @@
       ENDDO   !of n=1,n_out
 
+
 !------------------------------------------------------------------------
 !------------------------------------------------------------------------
@@ -310 +315 @@
          pbl_teta(ig) = ph(ig,ii(1))
       ENDDO
+
 
 !------------------------------------------------------------------------