--- trunk/dyn3d/laplacien.f90 2013/11/15 18:45:49 76 +++ trunk/Sources/dyn3d/laplacien.f 2016/03/11 18:47:26 178 @@ -10,21 +10,21 @@ ! P. Le Van ! Calcul de div(grad) de teta. - USE dimens_m, ONLY: llm + use dimens_m, only: iim, jjm use divergf_m, only: divergf - use filtreg_m, only: filtreg + use filtreg_scal_m, only: filtreg_scal use grad_m, only: grad - USE paramet_m, ONLY: ip1jm, ip1jmp1, jjp1 + USE paramet_m, ONLY: ip1jm, ip1jmp1 INTEGER, intent(in):: klevel - REAL, intent(inout):: teta(ip1jmp1, klevel) + REAL, intent(inout):: teta(iim + 1, jjm + 1, klevel) ! Variables locales: - REAL ghy(ip1jm, llm), ghx(ip1jmp1, llm) + REAL ghy(ip1jm, klevel), ghx(ip1jmp1, klevel) !----------------------------------------------------------------- - CALL filtreg(teta, jjp1, klevel, 2, 1, .TRUE.) + CALL filtreg_scal(teta, direct = .true., intensive = .true.) CALL grad(klevel, teta, ghx, ghy) CALL divergf(klevel, ghx, ghy, teta)