SUBROUTINE flumass3d(ucont,vcont,pext,pbaru,pbarv,w,dp)
c%W%	%G%
c=======================================================================
c
c   Author:    Frederic Hourdin      original: 17/02/93
c   -------
c
c   Subject:
c   ------
c   calcul du flux de masse dans les trois directions
c
c   Method:
c   --------
c
c   Interface:
c   ----------
c
c      Input:
c      ------
c
c      Output:
c      -------
c
c=======================================================================
      IMPLICIT NONE
c-----------------------------------------------------------------------
c   Declararations:
c   ---------------

#include "dimensions.h"
#include "paramet.h"
#include "comconst.h"
#include "comvert.h"
#include "comgeom.h"

c   Arguments:
c   ----------


      REAL pext(ip1jmp1)
      REAL vcont(ip1jm,llm),ucont(ip1jmp1,llm)
      REAL pbaru(ip1jmp1,llm),pbarv(ip1jm,llm)
      REAL dp(ip1jmp1)
      REAL w(ip1jmp1,llm)

c   Local:
c   ------

      REAL pbarx(ip1jmp1),pbary(ip1jm),pbarxy(ip1jm)
      REAL convm(ip1jmp1,llm)

c-----------------------------------------------------------------------

      CALL pbar(pext,pbarx,pbary,pbarxy)
      CALL flumass(pbarx,pbary,vcont,ucont,pbaru,pbarv)
      CALL convmas(pbaru,pbarv,convm)
      CALL SCOPY(ip1jmp1,convm,1,dp,1)
      CALL vitvert(convm,w)

      RETURN
      END