--- trunk/dyn3d/tourpot.f 2014/03/05 14:57:53 82 +++ trunk/dyn3d/tourpot.f 2014/09/11 15:09:15 107 @@ -11,8 +11,8 @@ ! Author: P. Le Van ! Objet : calcul du tourbillon potentiel - USE dimens_m, ONLY: iim, jjm, llm USE comgeom, ONLY: fext_2d + USE dimens_m, ONLY: iim, jjm, llm use filtreg_m, only: filtreg REAL, intent(in):: vcov(:, :, :) ! (iim + 1, jjm, llm) @@ -23,6 +23,7 @@ ! = (Filtre(d(vcov)/dx - d(ucov)/dy) + fext) / massebxy ! Local: + REAL rot(iim + 1, jjm, llm) ! relative vorticity multiplied by cell area, in m2 s-1 @@ -36,7 +37,7 @@ = vcov(i + 1, j, :) - vcov(i, j, :) + ucov(i, j + 1, :) - ucov(i, j, :) rot(iim + 1, :, :) = rot(1, :, :) - CALL filtreg(rot, jjm, llm, 2, 1, .FALSE.) + CALL filtreg(rot, direct = .true., intensive = .true.) forall (l = 1: llm) vorpot(:iim, :, l) & = (rot(:iim, :, l) + fext_2d(:iim, :)) / massebxy(:iim, :, l)