source: LMDZ.3.3/branches/LF/libf/dyn3d/caldyn0.F

Last change on this file 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: 2.7 KB
Line 
1      SUBROUTINE caldyn0
2     $ (itau,ucov,vcov,teta,ps,masse,pk,phis ,
3     $  phi,w,pbaru,pbarv,time )
4
5      IMPLICIT NONE
6
7c=======================================================================
8c
9c  Auteur :  P. Le Van
10c
11c   Objet:
12c   ------
13c
14c   Calcul des tendances dynamiques.
15c
16c Modif 04/93 F.Forget
17c=======================================================================
18
19c-----------------------------------------------------------------------
20c   0. Declarations:
21c   ----------------
22
23#include "dimensions.h"
24#include "paramet.h"
25#include "comconst.h"
26#include "comvert.h"
27#include "comgeom.h"
28
29c   Arguments:
30c   ----------
31
32      INTEGER itau
33      REAL vcov(ip1jm,llm),ucov(ip1jmp1,llm),teta(ip1jmp1,llm)
34      REAL ps(ip1jmp1),phis(ip1jmp1)
35      REAL pk(iip1,jjp1,llm)
36      REAL vcont(ip1jm,llm),ucont(ip1jmp1,llm)
37      REAL phi(ip1jmp1,llm),masse(ip1jmp1,llm)
38      REAL pbaru(ip1jmp1,llm),pbarv(ip1jm,llm)
39      REAL time
40
41c   Local:
42c   ------
43
44      REAL ang(ip1jmp1,llm),p(ip1jmp1,llmp1)
45      REAL massebx(ip1jmp1,llm),masseby(ip1jm,llm),psexbarxy(ip1jm)
46      REAL vorpot(ip1jm,llm)
47      REAL w(ip1jmp1,llm),ecin(ip1jmp1,llm),convm(ip1jmp1,llm)
48      REAL bern(ip1jmp1,llm)
49      REAL massebxy(ip1jm,llm), dp(ip1jmp1)
50   
51
52      INTEGER   ij,l
53      EXTERNAL  advect,bernoui,convmas,covcont ,
54     *          enercin,flumass,pbar,tourpot,vitvert,sortvarc,traceur,
55     *          pression,psextbar,massdair
56
57c-----------------------------------------------------------------------
58c   Calcul des tendances dynamiques:
59c   --------------------------------
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
Note: See TracBrowser for help on using the repository browser.