--- trunk/dyn3d/leapfrog.f90 2014/02/05 17:51:07 78 +++ trunk/dyn3d/leapfrog.f 2014/03/11 15:09:02 88 @@ -73,7 +73,7 @@ REAL massem1(iim + 1, jjm + 1, llm) ! Tendances dynamiques - REAL dv((iim + 1) * jjm, llm), dudyn((iim + 1) * (jjm + 1), llm) + REAL dv((iim + 1) * jjm, llm), dudyn(iim + 1, jjm + 1, llm) REAL dteta(iim + 1, jjm + 1, llm) real dp((iim + 1) * (jjm + 1)) @@ -214,10 +214,6 @@ teta(:, jjm + 1, l) = SUM(aire_2d(:iim, jjm+1) & * teta(:iim, jjm + 1, l)) / apols END forall - - ps(:, 1) = SUM(aire_2d(:iim, 1) * ps(:iim, 1)) / apoln - ps(:, jjm + 1) = SUM(aire_2d(:iim, jjm+1) * ps(:iim, jjm + 1)) & - / apols END IF IF (MOD(itau + 1, iperiod) == 0) THEN