--- trunk/Sources/filtrez/filtreg_hemisph.f	2015/06/05 19:21:08	141
+++ trunk/Sources/filtrez/filtreg_hemisph.f	2015/06/05 19:31:37	142
@@ -17,14 +17,12 @@
 
     !-----------------------------------------------------------------
 
-    DO l = 1, size(champ, 3)
-       DO j = 1, size(champ, 2)
-          champ(:iim, j, l) = champ(:iim, j, l) * sdd
-          champ(:iim, j, l) = (champ(:iim, j, l) &
-               + matmul(matri(:, :, j), champ(:iim, j, l))) / sdd
-          champ(iim + 1, j, l) = champ(1, j, l)
-       END DO
-    END DO
+    forall (j = 1:size(champ, 2), l = 1:size(champ, 3))
+       champ(:iim, j, l) = champ(:iim, j, l) &
+            + matmul(matri(:, :, j), champ(:iim, j, l) * sdd) / sdd
+    END forall
+
+    champ(iim + 1, :, :) = champ(1, :, :)
 
   end subroutine filtreg_hemisph