My Project
 All Classes Files Functions Variables Macros
caldyn0.F
Go to the documentation of this file.
1 !
2 ! $Header$
3 !
4  SUBROUTINE caldyn0
5  $ (itau,ucov,vcov,teta,ps,masse,pk,phis ,
6  $ phi,w,pbaru,pbarv,time )
7 
8  IMPLICIT NONE
9 
10 c=======================================================================
11 c
12 c Auteur : P. Le Van
13 c
14 c Objet:
15 c ------
16 c
17 c Calcul des tendances dynamiques.
18 c
19 c Modif 04/93 F.Forget
20 c=======================================================================
21 
22 c-----------------------------------------------------------------------
23 c 0. Declarations:
24 c ----------------
25 
26 #include "dimensions.h"
27 #include "paramet.h"
28 #include "comconst.h"
29 #include "comvert.h"
30 #include "comgeom.h"
31 
32 c Arguments:
33 c ----------
34 
35  INTEGER itau
36  REAL vcov(ip1jm,llm),ucov(ip1jmp1,llm),teta(ip1jmp1,llm)
37  REAL ps(ip1jmp1),phis(ip1jmp1)
38  REAL pk(iip1,jjp1,llm)
39  REAL vcont(ip1jm,llm),ucont(ip1jmp1,llm)
40  REAL phi(ip1jmp1,llm),masse(ip1jmp1,llm)
41  REAL pbaru(ip1jmp1,llm),pbarv(ip1jm,llm)
42  REAL time
43 
44 c Local:
45 c ------
46 
47  REAL ang(ip1jmp1,llm),p(ip1jmp1,llmp1)
48  REAL massebx(ip1jmp1,llm),masseby(ip1jm,llm),psexbarxy(ip1jm)
49  REAL vorpot(ip1jm,llm)
50  REAL w(ip1jmp1,llm),ecin(ip1jmp1,llm),convm(ip1jmp1,llm)
51  REAL bern(ip1jmp1,llm)
52  REAL massebxy(ip1jm,llm), dp(ip1jmp1)
53 
54 
55  INTEGER ij,l
56 
57 c-----------------------------------------------------------------------
58 c Calcul des tendances dynamiques:
59 c --------------------------------
60 
61  CALL covcont( llm , ucov , vcov , ucont, vcont )
62  CALL pression( ip1jmp1, ap , bp , ps , p )
63  CALL psextbar( ps , psexbarxy )
64  CALL massdair( p , masse )
65  CALL massbar( masse, massebx , masseby )
66  CALL massbarxy( masse, massebxy )
67  CALL flumass( massebx, masseby , vcont, ucont ,pbaru, pbarv )
68  CALL convmas( pbaru, pbarv , convm )
69 
70  DO ij =1, ip1jmp1
71  dp( ij ) = convm( ij,1 ) / airesurg( ij )
72  ENDDO
73 
74  CALL vitvert( convm , w )
75  CALL tourpot( vcov , ucov , massebxy , vorpot )
76  CALL enercin( vcov , ucov , vcont , ucont , ecin )
77  CALL bernoui( ip1jmp1, llm , phi , ecin , bern )
78 
79  DO l=1,llm
80  DO ij=1,ip1jmp1
81  ang(ij,l) = ucov(ij,l) + constang(ij)
82  ENDDO
83  ENDDO
84 
85  CALL sortvarc0
86  $ ( itau,ucov,teta,ps,masse,pk,phis,vorpot,phi,bern,dp,time,vcov )
87 
88  RETURN
89  END