--- trunk/dyn3d/leapfrog.f90	2013/11/15 18:45:49	76
+++ trunk/dyn3d/leapfrog.f	2014/03/11 15:09:02	88
@@ -73,7 +73,7 @@
     REAL massem1(iim + 1, jjm + 1, llm)
 
     ! Tendances dynamiques
-    REAL dv((iim + 1) * jjm, llm), dudyn((iim + 1) * (jjm + 1), llm)
+    REAL dv((iim + 1) * jjm, llm), dudyn(iim + 1, jjm + 1, llm)
     REAL dteta(iim + 1, jjm + 1, llm)
     real dp((iim + 1) * (jjm + 1))
 
@@ -136,7 +136,7 @@
        CALL geopot(teta, pk, pks, phis, phi)
        CALL caldyn(itau, ucov, vcov, teta, ps, masse, pk, pkf, phis, phi, &
             dudyn, dv, dteta, dp, w, pbaru, pbarv, time_0, &
-            conser=MOD(itau, iconser)==0)
+            conser = MOD(itau, iconser) == 0)
 
        CALL caladvtrac(q, pbaru, pbarv, p3d, masse, teta, pk)
 
@@ -158,7 +158,7 @@
           CALL geopot(teta, pk, pks, phis, phi)
           CALL caldyn(itau + 1, ucov, vcov, teta, ps, masse, pk, pkf, phis, &
                phi, dudyn, dv, dteta, dp, w, pbaru, pbarv, time_0, &
-               conser=.false.)
+               conser = .false.)
 
           ! integrations dynamique et traceurs:
           CALL integrd(vcovm1, ucovm1, tetam1, psm1, massem1, dv, dudyn, &
@@ -214,10 +214,6 @@
              teta(:, jjm + 1, l) = SUM(aire_2d(:iim, jjm+1) &
                   * teta(:iim, jjm + 1, l)) / apols
           END forall
-
-          ps(:, 1) = SUM(aire_2d(:iim, 1) * ps(:iim, 1)) / apoln
-          ps(:, jjm + 1) = SUM(aire_2d(:iim, jjm+1) * ps(:iim, jjm + 1)) &
-               / apols
        END IF
 
        IF (MOD(itau + 1, iperiod) == 0) THEN