--- trunk/Sources/dyn3d/fluxstokenc.f 2016/04/14 15:15:56 190 +++ trunk/Sources/dyn3d/fluxstokenc.f 2016/12/07 16:44:53 208 @@ -8,11 +8,12 @@ ! Author: F. Hourdin + USE comgeom, ONLY: aire + USE dimens_m, ONLY: jjm, llm + use groupe_m, only: groupe USE histwrite_m, ONLY: histwrite use initfluxsto_m, only: initfluxsto - USE dimens_m, ONLY: jjm, llm USE paramet_m, ONLY: iip1, ijmllm, ijp1llm, ip1jm, ip1jmp1 - USE comgeom, ONLY: aire USE tracstoke, ONLY: istdyn, istphy REAL, intent(in):: pbaru(ip1jmp1, llm), pbarv(ip1jm, llm) @@ -75,19 +76,18 @@ ! normalisation DO l = 1, llm DO ij = 1, ip1jmp1 - pbaruc(ij, l) = pbaruc(ij, l)/float(istdyn) - tetac(ij, l) = tetac(ij, l)/float(istdyn) - phic(ij, l) = phic(ij, l)/float(istdyn) + pbaruc(ij, l) = pbaruc(ij, l)/real(istdyn) + tetac(ij, l) = tetac(ij, l)/real(istdyn) + phic(ij, l) = phic(ij, l)/real(istdyn) END DO DO ij = 1, ip1jm - pbarvc(ij, l) = pbarvc(ij, l)/float(istdyn) + pbarvc(ij, l) = pbarvc(ij, l)/real(istdyn) END DO END DO ! Traitement des flux de masse avant advection. ! 1. Calcul de w ! 2. Groupement des mailles pr\`es du p\^ole. - CALL groupe(pbaruc, pbarvc, pbarug, pbarvg, wg) CALL histwrite(fluxid, 'masse', itau, massem)