--- trunk/libf/dyn3d/leapfrog.f90	2013/07/08 18:12:18	71
+++ trunk/dyn3d/leapfrog.f	2014/03/12 21:16:36	90
@@ -59,7 +59,7 @@
 
     REAL pks(iim + 1, jjm + 1) ! exner au sol
     REAL pk(iim + 1, jjm + 1, llm) ! exner au milieu des couches
-    REAL pkf(iim + 1, jjm + 1, llm) ! exner filtré au milieu des couches
+    REAL pkf(iim + 1, jjm + 1, llm) ! exner filtr\'e au milieu des couches
     REAL phi(iim + 1, jjm + 1, llm) ! geopotential
     REAL w((iim + 1) * (jjm + 1), llm) ! vitesse verticale
 
@@ -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))
 
@@ -95,7 +95,7 @@
     INTEGER l
     REAL rdayvrai, rdaym_ini
 
-    ! Variables test conservation énergie
+    ! Variables test conservation \'energie
     REAL ecin(iim + 1, jjm + 1, llm), ecin0(iim + 1, jjm + 1, llm)
 
     REAL vcont((iim + 1) * jjm, llm), ucont((iim + 1) * (jjm + 1), llm)
@@ -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)
 
@@ -144,7 +144,7 @@
        IF (offline) CALL fluxstokenc(pbaru, pbarv, masse, teta, phi, phis, &
             dtvr, itau)
 
-       ! Intégrations dynamique et traceurs:
+       ! Int\'egrations dynamique et traceurs:
        CALL integrd(vcovm1, ucovm1, tetam1, psm1, massem1, dv, dudyn, dteta, &
             dp, vcov, ucov, teta, q(:, :, :, :2), ps, masse, finvmaold, dt, &
             leapf)
@@ -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, &
@@ -178,8 +178,7 @@
           IF (time > 1.) time = time - 1.
 
           CALL calfis(rdayvrai, time, ucov, vcov, teta, q, ps, pk, phis, phi, &
-               dudyn, dv, w, dufi, dvfi, dtetafi, dqfi, dpfi, &
-               lafin = itau + 1 == itaufin)
+               w, dufi, dvfi, dtetafi, dqfi, dpfi, lafin = itau + 1 == itaufin)
 
           ! Ajout des tendances physiques:
           CALL addfi(ucov, vcov, teta, q, ps, dufi, dvfi, dtetafi, dqfi, dpfi)
@@ -189,9 +188,9 @@
        CALL exner_hyb(ps, p3d, pks, pk, pkf)
 
        IF (MOD(itau + 1, idissip) == 0) THEN
-          ! Dissipation horizontale et verticale des petites échelles
+          ! Dissipation horizontale et verticale des petites \'echelles
 
-          ! calcul de l'énergie cinétique avant dissipation
+          ! calcul de l'\'energie cin\'etique avant dissipation
           call covcont(llm, ucov, vcov, ucont, vcont)
           call enercin(vcov, ucov, vcont, ucont, ecin0)
 
@@ -200,28 +199,24 @@
           ucov = ucov + dudis
           vcov = vcov + dvdis
 
-          ! On ajoute la tendance due à la transformation énergie
-          ! cinétique en énergie thermique par la dissipation
+          ! On ajoute la tendance due \`a la transformation \'energie
+          ! cin\'etique en \'energie thermique par la dissipation
           call covcont(llm, ucov, vcov, ucont, vcont)
           call enercin(vcov, ucov, vcont, ucont, ecin)
           dtetadis = dtetadis + (ecin0 - ecin) / pk
           teta = teta + dtetadis
 
-          ! Calcul de la valeur moyenne aux pôles :
+          ! Calcul de la valeur moyenne aux p\^oles :
           forall (l = 1: llm)
              teta(:, 1, l) = SUM(aire_2d(:iim, 1) * teta(:iim, 1, l)) &
                   / apoln
              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
-          ! Écriture du fichier histoire moyenne:
+          ! \'Ecriture du fichier histoire moyenne:
           CALL writedynav(vcov, ucov, teta, pk, phi, q, masse, ps, phis, &
                time = itau + 1)
           call bilan_dyn(ps, masse, pk, pbaru, pbarv, teta, phi, ucov, vcov, &