fortran careful with implicit...

PIV.F90

  DO IPIV = 1, NPIV
     WHERE (VALID(:,:,IPIV)==0)
        UMOY = UMOY + U(:,:,IPIV)
        VMOY = VMOY + V(:,:,IPIV)
        UU = UU + U(:,:,IPIV)**2
        VV = VV + V(:,:,IPIV)**2
        UV = UV + U(:,:,IPIV)*V(:,:,IPIV)
        NVALID = NVALID + 1
     END WHERE
  END DO
  WHERE (NVALID /= 0)
     UMOY = UMOY/NVALID
     VMOY = VMOY/NVALID
     UU = UU/NVALID
     VV = VV/NVALID
     UV = UV/NVALID
  END WHERE
  UU=UU-UMOY**2
  VV=VV-VMOY**2
  UV=UV-UMOY*VMOY

!!$  DO J = 1, NY
!!$     DO I = 1, NX
!!$        DO IPIV = 1, NPIV
!!$           NVALID(I,J)=0
!!$           IF (VALID(I,J,IPIV) == 0) THEN
!!$              NVALID(I,J) = NVALID(I,J) + 1
!!$              UMOY(I,J) = UMOY(I,J) + U(I,J,IPIV)
!!$              VMOY(I,J) = VMOY(I,J) + V(I,J,IPIV)
!!$              UU(I,J) = UU(I,J) + U(I,J,IPIV)**2
!!$              VV(I,J) = VV(I,J) + V(I,J,IPIV)**2
!!$              UV(I,J) = UV(I,J) + U(I,J,IPIV)*V(I,J,IPIV)
!!$           ENDIF
!!$        END DO
!!$     END DO
!!$  END DO
!!$  DO J = 1, NY
!!$     DO I = 1, NX
!!$        IF (NVALID(I,J) /= 0) THEN
!!$           UMOY(I,J) = UMOY(I,J)/NVALID(I,J)
!!$           VMOY(I,J) = VMOY(I,J)/NVALID(I,J)
!!$           UV(I,J) = UV(I,J) + U(I,J,IPIV)*V(I,J,IPIV)
!!$        ENDIF
!!$     END DO
!!$  END DO
!!$  UU=UU-UMOY**2
!!$  VV=VV-VMOY**2
!!$  UV=UV-UMOY*VMOY
!!$  
!!$  DO I=1,NX
!!$     DO J=1,NY
!!$        NVALID(I,J)=0
!!$        DO IPIV=1,NPIV
!!$           IF (VALID(I,J,IPIV).EQ.0) THEN
!!$              NVALID(I,J)=NVALID(I,J)+1
!!$              UMOY(I,J)=(UMOY(I,J)*REAL(NVALID(I,J)-1)+U(I,J,IPIV))&
!!$                   /REAL(NVALID(I,J))
!!$              VMOY(I,J)=(VMOY(I,J)*REAL(NVALID(I,J)-1)+V(I,J,IPIV))&
!!$                   /REAL(NVALID(I,J))
!!$              UU(I,J)=(UU(I,J)*REAL(NVALID(I,J)-1)+U(I,J,IPIV)**2)&
!!$                   /REAL(NVALID(I,J))
!!$              VV(I,J)=(VV(I,J)*REAL(NVALID(I,J)-1)+V(I,J,IPIV)**2)&
!!$                   /REAL(NVALID(I,J))
!!$              UV(I,J)=(UV(I,J)*REAL(NVALID(I,J)-1)+U(I,J,IPIV)*V(I,J,IPIV))&
!!$                   /REAL(NVALID(I,J))
!!$           ENDIF
!!$        END DO
!!$     END DO
!!$  END DO
!!$  UU=UU-UMOY**2
!!$  VV=VV-VMOY**2
!!$  UV=UV-UMOY*VMOY