Index: trunk/LMDZ.MARS/libf/phymars/physiq.F
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/physiq.F	(revision 268)
+++ trunk/LMDZ.MARS/libf/phymars/physiq.F	(revision 269)
@@ -1160,4 +1160,13 @@
       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,6 +1582,5 @@
          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,8 +1691,6 @@
          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))