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

Timestamp:

Aug 17, 2011, 11:38:15 AM (13 years ago)

Author:

acolaitis

Message:

Bug correction in call to surflayer, added recomputation of potential temperature zh after incrementation of physical tendancies.

File:

• trunk/LMDZ.MARS/libf/phymars/physiq.F (modified) (3 diffs)

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

@@ -1160 +1160 @@
       ENDDO
 
+      ! Potential Temperature
+
+       DO ig=1,ngridmx
+          DO l=1,nlayermx
+              zh(ig,l) = zt(ig,l)*(zplay(ig,l)/zplev(ig,1))**rcp
+          ENDDO
+       ENDDO
+
+
 c    Compute surface stress : (NB: z0 is a common in surfdat.h)
 c     DO ig=1,ngrid
@@ -1573 +1582 @@
          if (calltherm .and. (z_out .gt. 0.)) then
          call surflayer_interpol(ngrid,nlayer,z0,g,zzlay,zu,zv,wmax_th
-     &              ,tsurf,zt(:,:)*(zplay(:,:)/zplev(:,:))**rcp
-     &              ,z_out,Teta_out,u_out,ustar,tstar)
+     &              ,tsurf,zh,z_out,Teta_out,u_out,ustar,tstar)
 
          zu2(:)=sqrt(zu(:,1)*zu(:,1)+zv(:,1)*zv(:,1))
@@ -1683 +1691 @@
          z_out=0.
          if (calltherm .and. (z_out .gt. 0.)) then
-
          call surflayer_interpol(ngrid,nlayer,z0,g,zzlay,zu,zv,wmax_th
-     &              ,tsurf,zt(:,:)*(zplay(:,:)/zplev(:,:))**rcp
-     &              ,z_out,Teta_out,u_out,ustar,tstar)
+     &              ,tsurf,zh,z_out,Teta_out,u_out,ustar,tstar)
 
          zu2(:)=sqrt(zu(:,1)*zu(:,1)+zv(:,1)*zv(:,1))