--- trunk/libf/dyn3d/leapfrog.f90	2008/07/25 19:59:34	13
+++ trunk/libf/dyn3d/leapfrog.f90	2009/12/14 15:25:16	23
@@ -6,7 +6,7 @@
 
 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
 
@@ -27,8 +27,8 @@
     ! Pour Van-Leer + Vapeur d'eau saturee, iadv(1)=4. (F.Codron, 10/99)
     ! 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, ijmllm, 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 +46,6 @@
     use pression_m, only: pression
     use pressure_var, only: p3d
 
-    integer nq
-
     ! Variables dynamiques:
     REAL vcov(ip1jm, llm), ucov(ip1jmp1, llm) ! vents covariants
     REAL teta(ip1jmp1, llm) ! temperature potentielle 
@@ -92,9 +90,10 @@
 
     REAL tppn(iim), tpps(iim), tpn, tps
 
-    INTEGER itau, itaufinp1
+    INTEGER itau ! index of the time step of the dynamics, starts at 0
+    integer itaufinp1
     INTEGER iday ! jour julien
-    REAL time ! Heure de la journee en fraction d'1 jour
+    REAL time ! time of day, as a fraction of day length
 
     REAL SSUM
     real finvmaold(ip1jmp1, llm)
@@ -103,7 +102,6 @@
     INTEGER ij, l
 
     REAL rdayvrai, rdaym_ini
-    LOGICAL:: callinigrads = .true.
 
     !+jld variables test conservation energie
     REAL ecin(ip1jmp1, llm), ecin0(ip1jmp1, llm)
@@ -116,7 +114,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
 
     !---------------------------------------------------
@@ -216,7 +213,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)
 
@@ -317,7 +314,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
 
@@ -335,7 +332,7 @@
                 dt = 2. * dtvr
              END IF
           ELSE
-             ! ...... pas leapfrog .....
+             ! pas leapfrog
              leapf = .TRUE.
              dt = 2. * dtvr
           END IF