SUBROUTINE angumom(imx,lmx,u,ang)
      IMPLICIT NONE

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

      INTEGER imx,lmx
      REAL u(imx,jjp1,lmx)
      REAL ang(imx,jjp1,lmx)

      INTEGER i,j,l
      REAL  radomeg,zcosphi(jjp1),omegcosp(jjp1)
      REAL zrad2,c1,angsol

c-----------------------------------------------------------------------
c   moment cinetique:
c   -----------------


      radomeg=omeg*rad
      DO j=2,jjm
         zcosphi(j)=COS( rlatu(j ))
         omegcosp(j)=radomeg*zcosphi(j)
      ENDDO

      zrad2=rad*rad
      angsol=2./3.*zrad2*omeg
      c1=rad/angsol

      DO l = 1  ,lmx
         DO j =2,jjm
            DO i=1,imx
               ang(i,j,l) =c1*(u(i,j,l)+omegcosp(j))*zcosphi(j)
            ENDDO
         ENDDO
      ENDDO

      RETURN
      END