source: LMDZ.3.3/trunk/libf/dyn3d/caldyn.F @ 346

Last change on this file since 346 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: 3.8 KB
Line 
1      SUBROUTINE caldyn
2     $ (itau,ucov,vcov,teta,ps,masse,pk,pkf,phis ,
3     $  phi,conser,du,dv,dteta,dp,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      LOGICAL conser
33
34      INTEGER itau
35      REAL vcov(ip1jm,llm),ucov(ip1jmp1,llm),teta(ip1jmp1,llm)
36      REAL ps(ip1jmp1),phis(ip1jmp1)
37      REAL pk(iip1,jjp1,llm),pkf(ip1jmp1,llm)
38      REAL vcont(ip1jm,llm),ucont(ip1jmp1,llm)
39      REAL phi(ip1jmp1,llm),masse(ip1jmp1,llm)
40      REAL dv(ip1jm,llm),du(ip1jmp1,llm)
41      REAL dteta(ip1jmp1,llm),dp(ip1jmp1)
42      REAL pbaru(ip1jmp1,llm),pbarv(ip1jm,llm)
43      REAL time
44
45c   Local:
46c   ------
47
48      REAL ang(ip1jmp1,llm),p(ip1jmp1,llmp1)
49      REAL massebx(ip1jmp1,llm),masseby(ip1jm,llm),psexbarxy(ip1jm)
50      REAL vorpot(ip1jm,llm)
51      REAL w(ip1jmp1,llm),ecin(ip1jmp1,llm),convm(ip1jmp1,llm)
52      REAL bern(ip1jmp1,llm)
53      REAL massebxy(ip1jm,llm)
54   
55
56      INTEGER   ij,l
57      EXTERNAL  advect,bernoui,convmas,covcont,dteta1,dudv1,dudv2,
58     *          enercin,flumass,pbar,tourpot,vitvert,sortvarc,traceur,
59     *          pression,psextbar,massdair
60
61c-----------------------------------------------------------------------
62c   Calcul des tendances dynamiques:
63c   --------------------------------
64
65      CALL covcont  ( llm    , ucov    , vcov , ucont, vcont        )
66      CALL pression ( ip1jmp1, ap      , bp   ,  ps  , p            )
67      CALL psextbar (   ps   , psexbarxy                            )
68      CALL massdair (    p   , masse                                )
69      CALL massbar  (   masse, massebx , masseby                    )
70      call massbarxy(   masse, massebxy                             )
71      CALL flumass  ( massebx, masseby , vcont, ucont ,pbaru, pbarv )
72      CALL dteta1   (   teta , pbaru   , pbarv, dteta               )
73      CALL convmas  (   pbaru, pbarv   , convm                      )
74
75      DO ij =1, ip1jmp1
76         dp( ij ) = convm( ij,1 ) / airesurg( ij )
77      ENDDO
78
79      CALL vitvert ( convm  , w                                  )
80      CALL tourpot ( vcov   , ucov  , massebxy  , vorpot         )
81      CALL dudv1   ( vorpot , pbaru , pbarv     , du     , dv    )
82      CALL enercin ( vcov   , ucov  , vcont     , ucont  , ecin  )
83      CALL bernoui ( ip1jmp1, llm   , phi       , ecin   , bern  )
84      CALL dudv2   ( teta   , pkf   , bern      , du     , dv    )
85
86
87      DO l=1,llm
88         DO ij=1,ip1jmp1
89            ang(ij,l) = ucov(ij,l) + constang(ij)
90      ENDDO
91      ENDDO
92
93
94      CALL advect( ang, vcov, teta, w, massebx, masseby, du, dv,dteta )
95
96C  WARNING probleme de peridocite de dv sur les PC/linux. Pb d'arrondi
97C          probablement. Observe sur le code compile avec pgf90 3.0-1
98
99      DO l = 1, llm
100         DO ij = 1, ip1jm, iip1
101           IF( dv(ij,l).NE.dv(ij+iim,l) )  THEN
102c         PRINT *,'!!!ATTENTION!!! probleme de periodicite sur vcov', 
103c    ,   ' dans caldyn'
104c         PRINT *,' l,  ij = ', l, ij, ij+iim,dv(ij+iim,l),dv(ij,l)
105          dv(ij+iim,l) = dv(ij,l)
106          endif
107         enddo
108      enddo
109c-----------------------------------------------------------------------
110c   Sorties eventuelles des variables de controle:
111c   ----------------------------------------------
112
113      IF( conser )  THEN
114        CALL sortvarc
115     $ ( itau,ucov,teta,ps,masse,pk,phis,vorpot,phi,bern,dp,time,vcov )
116
117      ENDIF
118
119      RETURN
120      END
Note: See TracBrowser for help on using the repository browser.