Index: trunk/LMDZ.MARS/libf/dyn3d/gcm.F
===================================================================
--- trunk/LMDZ.MARS/libf/dyn3d/gcm.F	(revision 795)
+++ trunk/LMDZ.MARS/libf/dyn3d/gcm.F	(revision 798)
@@ -189,8 +189,6 @@
 c
 
-      CALL pression ( ip1jmp1, ap, bp, ps, p       )
 
       call dump2d(iip1,jjp1,ps,'PRESSION SURFACE')
-      CALL exner_hyb( ip1jmp1, ps, p,beta, pks, pk, pkf )
 
 c
@@ -249,4 +247,9 @@
    1  CONTINUE
 c
+c TN 09/2012. To ensure "1+1=2" in dynamical core :
+c update atmospheric pressure IN the main loop
+      CALL pression ( ip1jmp1, ap, bp, ps, p       )
+      CALL exner_hyb( ip1jmp1, ps, p,beta, pks, pk, pkf )
+
       IF( MOD( itau, 10* day_step ).EQ.0 )  THEN
         CALL test_period ( ucov,vcov,teta,q,p,phis )
