Index: trunk/LMDZ.MARS/libf/dyn3d/addfi.F
===================================================================
--- trunk/LMDZ.MARS/libf/dyn3d/addfi.F	(revision 1232)
+++ trunk/LMDZ.MARS/libf/dyn3d/addfi.F	(revision 1238)
@@ -104,5 +104,12 @@
        ENDDO
       ENDDO
-c
+!***********************
+! Correction on teta due to surface pressure changes
+      DO k = 1,llm
+        DO j = 1,ip1jmp1
+           pteta(j,k)= pteta(j,k)*(1+pdpfi(j)*pdt/pps(j))**kappa
+        ENDDO
+      ENDDO
+!***********************
 
 ! increment covariant zonal wind