--- trunk/libf/dyn3d/leapfrog.f90 2012/04/20 14:58:43 61 +++ trunk/libf/dyn3d/leapfrog.f90 2012/09/20 13:00:41 66 @@ -17,7 +17,7 @@ USE calfis_m, ONLY: calfis USE comconst, ONLY: daysec, dtphys, dtvr USE comgeom, ONLY: aire_2d, apoln, apols - USE comvert, ONLY: ap, bp + USE disvert_m, ONLY: ap, bp USE conf_gcm_m, ONLY: day_step, iconser, iperiod, iphysiq, nday, offline, & iflag_phys, ok_guide USE dimens_m, ONLY: iim, jjm, llm, nqmx @@ -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: