--- trunk/libf/dyn3d/leapfrog.f90 2012/01/30 12:54:02 57 +++ trunk/libf/dyn3d/leapfrog.f90 2012/09/20 09:57:03 65 @@ -57,7 +57,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 filt.au milieu des couches + REAL pkf(iim + 1, jjm + 1, llm) ! exner filtré au milieu des couches REAL phi(iim + 1, jjm + 1, llm) ! geopotential REAL w((iim + 1) * (jjm + 1), llm) ! vitesse verticale @@ -106,7 +106,7 @@ call assert(shape(ucov) == (/iim + 1, jjm + 1, llm/), "leapfrog") itaufin = nday * day_step - ! "day_step" is a multiple of "iperiod", therefore "itaufin" is one too + ! "day_step" is a multiple of "iperiod", therefore so is "itaufin". dq = 0. @@ -129,7 +129,7 @@ massem1 = masse psm1 = ps finvmaold = masse - CALL filtreg(finvmaold, jjm + 1, llm, - 2, 2, .TRUE., 1) + CALL filtreg(finvmaold, jjm + 1, llm, - 2, 2, .TRUE.) end if ! Calcul des tendances dynamiques: @@ -224,8 +224,8 @@ IF (MOD(itau + 1, iperiod) == 0) THEN ! Écriture du fichier histoire moyenne: - CALL writedynav(nqmx, itau + 1, vcov, ucov, teta, pk, phi, q, & - masse, ps, phis) + 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, & q(:, :, :, 1)) ENDIF