--- trunk/dyn3d/fluxstokenc.f90 2013/11/15 18:45:49 76 +++ trunk/Sources/dyn3d/fluxstokenc.f 2015/06/10 16:46:46 144 @@ -15,7 +15,7 @@ USE tracstoke, ONLY: istdyn, istphy REAL pbaru(ip1jmp1, llm), pbarv(ip1jm, llm) - REAL masse(ip1jmp1, llm) + REAL, intent(in):: masse(ip1jmp1, llm) real, intent(in):: phi(ip1jmp1, llm) real, intent(in):: teta(ip1jmp1, llm) REAL, intent(in):: phis(ip1jmp1) @@ -37,7 +37,7 @@ !------------------------------------------------------------- IF (itau == 0) THEN - CALL initfluxsto(time_step, istdyn*time_step, istdyn*time_step, nqmx, & + CALL initfluxsto(time_step, istdyn*time_step, istdyn*time_step, & fluxid, fluxvid, fluxdid) CALL histwrite(fluxid, 'phis', 1, phis) CALL histwrite(fluxid, 'aire', 1, aire) @@ -67,9 +67,7 @@ END DO ! selection de la masse instantannee des mailles avant le transport. - IF (itau == 0) THEN - CALL scopy(ip1jmp1*llm, masse, 1, massem, 1) - END IF + IF (itau == 0) massem = masse IF (mod(itau + 1, istdyn) == 0) THEN ! on advecte a ce pas de temps