Index: trunk/LMDZ.COMMON/libf/dyn3d/addfi.F
===================================================================
--- trunk/LMDZ.COMMON/libf/dyn3d/addfi.F	(revision 1209)
+++ trunk/LMDZ.COMMON/libf/dyn3d/addfi.F	(revision 1238)
@@ -103,5 +103,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
+!***********************
 
       DO k = 1,llm
Index: trunk/LMDZ.COMMON/libf/dyn3d/inidissip.F90
===================================================================
--- trunk/LMDZ.COMMON/libf/dyn3d/inidissip.F90	(revision 1209)
+++ trunk/LMDZ.COMMON/libf/dyn3d/inidissip.F90	(revision 1238)
@@ -268,5 +268,6 @@
     DO l=1,llm
      zz      = 1. - preff/presnivs(l)
-     zvert(l)= dissip_fac_mid -( dissip_fac_mid-1.)/( 1.+zz*zz )
+!     zvert(l)= dissip_fac_mid -( dissip_fac_mid-1.)/( 1.+zz*zz )
+     zvert(l)= fac_mid -( fac_mid-1.)/( 1.+zz*zz )
      
      zvert(l)= zvert(l)*(1.0+((fac_up/fac_mid-1)*    &
@@ -275,5 +276,7 @@
                 ))
     ENDDO
-
+    write(lunout,*) "inidissip: vert_prof_disip=1, scaleheight=",scaleheight
+    write(lunout,*) "           fac_mid=",fac_mid,", fac_up=",fac_up
+    
    else
      write(lunout,*) 'wrong value for vert_prof_dissip:',vert_prof_dissip