--- trunk/libf/dyn3d/leapfrog.f90	2009/07/31 15:18:47	22
+++ trunk/libf/dyn3d/leapfrog.f90	2010/03/03 13:23:49	24
@@ -1,34 +1,31 @@
 module leapfrog_m
 
-  ! This module is clean: no C preprocessor directive, no include line.
-
   IMPLICIT NONE
 
 contains
 
-  SUBROUTINE leapfrog(ucov, vcov, teta, ps, masse, phis, nq, q, time_0)
+  SUBROUTINE leapfrog(ucov, vcov, teta, ps, masse, phis, q, time_0)
 
     ! From dyn3d/leapfrog.F, version 1.6 2005/04/13 08:58:34
 
-    ! Version du 10/01/98, avec coordonnees verticales hybrides, avec
-    ! nouveaux operat. dissipation * (gradiv2, divgrad2, nxgraro2)
+    ! Version du 10/01/98, avec coordonnées verticales hybrides, avec
+    ! nouveaux opérateurs dissipation "*" (gradiv2, divgrad2, nxgraro2)
 
-    ! Auteur: P. Le Van /L. Fairhead/F.Hourdin
-    ! Objet:
-    ! GCM LMD nouvelle grille
+    ! Auteurs : P. Le Van, L. Fairhead, F. Hourdin
+    ! Objet: nouvelle grille
 
-    ! ... Dans inigeom, nouveaux calculs pour les elongations cu, cv
-    ! et possibilite d'appeler une fonction f(y) a derivee tangente
-    ! hyperbolique a la place de la fonction a derivee sinusoidale.
+    ! Dans "inigeom", nouveaux calculs pour les élongations cu, cv
+    ! et possibilité d'appeler une fonction f(y) à dérivée tangente
+    ! hyperbolique à la place de la fonction à dérivée sinusoïdale.
 
-    ! ... Possibilité de choisir le schéma pour l'advection de
-    ! q, en modifiant iadv dans "traceur.def" (10/02) .
+    ! Possibilité de choisir le schéma pour l'advection de
+    ! q, en modifiant iadv dans "traceur.def".
 
-    ! Pour Van-Leer + Vapeur d'eau saturee, iadv(1)=4. (F.Codron, 10/99)
+    ! Pour Van-Leer + vapeur d'eau saturée, iadv(1)=4.
     ! Pour Van-Leer iadv=10 
 
-    use dimens_m, only: iim, jjm, llm, nqmx
-    use paramet_m, only: ip1jmp1, ip1jm, ijmllm, ijp1llm, jjp1, iip1, iip2
+    use dimens_m, only: iim, llm, nqmx
+    use paramet_m, only: ip1jmp1, ip1jm, ijp1llm, jjp1, iip1
     use comconst, only: dtvr, daysec, dtphys
     use comvert, only: ap, bp
     use conf_gcm_m, only: day_step, iconser, idissip, iphysiq, iperiod, nday, &
@@ -46,8 +43,6 @@
     use pression_m, only: pression
     use pressure_var, only: p3d
 
-    integer, intent(in):: nq
-
     ! Variables dynamiques:
     REAL vcov(ip1jm, llm), ucov(ip1jmp1, llm) ! vents covariants
     REAL teta(ip1jmp1, llm) ! temperature potentielle 
@@ -56,7 +51,7 @@
     REAL masse(ip1jmp1, llm) ! masse d'air
     REAL phis(ip1jmp1) ! geopotentiel au sol
 
-    REAL time_0
+    REAL, intent(in):: time_0
 
     ! Variables local to the procedure:
 
@@ -93,7 +88,6 @@
     REAL tppn(iim), tpps(iim), tpn, tps
 
     INTEGER itau ! index of the time step of the dynamics, starts at 0
-    integer itaufinp1
     INTEGER iday ! jour julien
     REAL time ! time of day, as a fraction of day length
 
@@ -104,9 +98,8 @@
     INTEGER ij, l
 
     REAL rdayvrai, rdaym_ini
-    LOGICAL:: callinigrads = .true.
 
-    !+jld variables test conservation energie
+    ! Variables test conservation energie
     REAL ecin(ip1jmp1, llm), ecin0(ip1jmp1, llm)
     ! Tendance de la temp. potentiel d (theta) / d t due a la 
     ! tansformation d'energie cinetique en energie thermique
@@ -117,7 +110,6 @@
     INTEGER:: ip_ebil_dyn = 0 ! PRINT level for energy conserv. diag.
 
     logical:: dissip_conservative = .true.
-    LOGICAL:: prem = .true.
     logical forward, leapf, apphys, conser, apdiss
 
     !---------------------------------------------------
@@ -125,8 +117,6 @@
     print *, "Call sequence information: leapfrog"
 
     itaufin = nday * day_step
-    itaufinp1 = itaufin + 1
-
     itau = 0
     iday = day_ini
     time = time_0
@@ -135,21 +125,21 @@
        iday = iday + 1
     ENDIF
 
+    dq = 0.
     ! On initialise la pression et la fonction d'Exner :
-    dq=0.
     CALL pression(ip1jmp1, ap, bp, ps, p3d)
     CALL exner_hyb(ps, p3d, pks, pk, pkf)
 
     ! Debut de l'integration temporelle:
     outer_loop:do
-       if (ok_guide.and.(itaufin - itau - 1) * dtvr > 21600) then
+       if (ok_guide.and.(itaufin - itau - 1) * dtvr > 21600.) then
           call guide(itau, ucov, vcov, teta, q, masse, ps)
        else
           IF (prt_level > 9) print *, &
-               'Attention : on ne guide pas les 6 dernieres heures.'
+               'Attention : on ne guide pas les 6 dernières heures.'
        endif
 
-       CALL SCOPY(ijmllm, vcov, 1, vcovm1, 1)
+       CALL SCOPY(ip1jm * llm, vcov, 1, vcovm1, 1)
        CALL SCOPY(ijp1llm, ucov, 1, ucovm1, 1)
        CALL SCOPY(ijp1llm, teta, 1, tetam1, 1)
        CALL SCOPY(ijp1llm, masse, 1, massem1, 1)
@@ -217,7 +207,7 @@
                      teta, q(:, :, 1), q(:, :, 2))
              ENDIF
 
-             CALL calfis(nq, lafin, rdayvrai, time, ucov, vcov, teta, q, &
+             CALL calfis(nqmx, lafin, rdayvrai, time, ucov, vcov, teta, q, &
                   masse, ps, pk, phis, phi, du, dv, dteta, dq, w, &
                   dufi, dvfi, dtetafi, dqfi, dpfi)
 
@@ -304,7 +294,7 @@
              ENDIF
           ENDIF
 
-          IF (itau == itaufinp1) exit outer_loop
+          IF (itau == itaufin + 1) exit outer_loop
 
           ! ecriture du fichier histoire moyenne:
 
@@ -318,7 +308,7 @@
           ENDIF
 
           IF (itau == itaufin) THEN
-             CALL dynredem1("restart.nc", 0., vcov, ucov, teta, q, masse, ps)
+             CALL dynredem1("restart.nc", vcov, ucov, teta, q, masse, ps)
              CLOSE(99)
           ENDIF