source: LMDZ.3.3/trunk/libf/dyn3d/angumom.F @ 604

Last change on this file since 604 was 2, checked in by lmdz, 25 years ago

Initial revision
  • Property svn:eol-style set to native
  • Property svn:executable set to *
  • Property svn:keywords set to Author Date Id Revision
File size: 843 bytes
Line 
1      SUBROUTINE angumom(imx,lmx,u,ang)
2      IMPLICIT NONE
3
4#include "dimensions.h"
5#include "paramet.h"
6#include "comconst.h"
7#include "comgeom.h"
8
9      INTEGER imx,lmx
10      REAL u(imx,jjp1,lmx)
11      REAL ang(imx,jjp1,lmx)
12
13      INTEGER i,j,l
14      REAL  radomeg,zcosphi(jjp1),omegcosp(jjp1)
15      REAL zrad2,c1,angsol
16
17c-----------------------------------------------------------------------
18c   moment cinetique:
19c   -----------------
20
21
22      radomeg=omeg*rad
23      DO j=2,jjm
24         zcosphi(j)=COS( rlatu(j ))
25         omegcosp(j)=radomeg*zcosphi(j)
26      ENDDO
27
28      zrad2=rad*rad
29      angsol=2./3.*zrad2*omeg
30      c1=rad/angsol
31
32      DO l = 1  ,lmx
33         DO j =2,jjm
34            DO i=1,imx
35               ang(i,j,l) =c1*(u(i,j,l)+omegcosp(j))*zcosphi(j)
36            ENDDO
37         ENDDO
38      ENDDO
39
40      RETURN
41      END
Note: See TracBrowser for help on using the repository browser.