--- trunk/dyn3d/laplacien.f 2014/08/29 13:00:05 103 +++ trunk/dyn3d/laplacien.f 2014/09/11 15:09:15 107 @@ -10,20 +10,21 @@ ! P. Le Van ! Calcul de div(grad) de teta. + use dimens_m, only: iim, jjm use divergf_m, only: divergf use filtreg_m, only: filtreg use grad_m, only: grad USE paramet_m, ONLY: ip1jm, ip1jmp1, jjp1 INTEGER, intent(in):: klevel - REAL, intent(inout):: teta(ip1jmp1, klevel) + REAL, intent(inout):: teta(iim + 1, jjm + 1, klevel) ! Variables locales: REAL ghy(ip1jm, klevel), ghx(ip1jmp1, klevel) !----------------------------------------------------------------- - CALL filtreg(teta, jjp1, klevel, 2, 1, .TRUE.) + CALL filtreg(teta, direct = .true., intensive = .true.) CALL grad(klevel, teta, ghx, ghy) CALL divergf(klevel, ghx, ghy, teta)