source: LMDZ5/trunk/libf/dyn3d_common/caldyn0.F @ 2225

Last change on this file since 2225 was 2083, checked in by Ehouarn Millour, 10 years ago
  • Minor fix in dyn3dpar/leapfrog_p.F , should call geopot_p and not geopot
  • Added a sanity check in iniacademic
  • Added flag "resetvarc" to trigger a reset of initial values in sortvarc
  • Removed "sortvarc0" since the job can now be done with "resetvarc" and having set flag resertvarc to true.

EM
  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 2.6 KB
Line 
1!
2! $Header$
3!
4      SUBROUTINE caldyn0
5     $ (itau,ucov,vcov,teta,ps,masse,pk,phis ,
6     $  phi,w,pbaru,pbarv,time )
7
8      USE control_mod, ONLY: resetvarc
9      IMPLICIT NONE
10
11c=======================================================================
12c
13c  Auteur :  P. Le Van
14c
15c   Objet:
16c   ------
17c
18c   Calcul des tendances dynamiques.
19c
20c Modif 04/93 F.Forget
21c=======================================================================
22
23c-----------------------------------------------------------------------
24c   0. Declarations:
25c   ----------------
26
27#include "dimensions.h"
28#include "paramet.h"
29#include "comconst.h"
30#include "comvert.h"
31#include "comgeom.h"
32
33c   Arguments:
34c   ----------
35
36      INTEGER itau
37      REAL vcov(ip1jm,llm),ucov(ip1jmp1,llm),teta(ip1jmp1,llm)
38      REAL ps(ip1jmp1),phis(ip1jmp1)
39      REAL pk(iip1,jjp1,llm)
40      REAL vcont(ip1jm,llm),ucont(ip1jmp1,llm)
41      REAL phi(ip1jmp1,llm),masse(ip1jmp1,llm)
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), dp(ip1jmp1)
54   
55
56      INTEGER   ij,l
57
58c-----------------------------------------------------------------------
59c   Calcul des tendances dynamiques:
60c   --------------------------------
61
62      CALL covcont  ( llm    , ucov    , vcov , ucont, vcont        )
63      CALL pression ( ip1jmp1, ap      , bp   ,  ps  , p            )
64      CALL psextbar (   ps   , psexbarxy                            )
65      CALL massdair (    p   , masse                                )
66      CALL massbar  (   masse, massebx , masseby                    )
67      CALL massbarxy(   masse, massebxy                             )
68      CALL flumass  ( massebx, masseby , vcont, ucont ,pbaru, pbarv )
69      CALL convmas  (   pbaru, pbarv   , convm                      )
70
71      DO ij =1, ip1jmp1
72         dp( ij ) = convm( ij,1 ) / airesurg( ij )
73      ENDDO
74
75      CALL vitvert ( convm  , w                                  )
76      CALL tourpot ( vcov   , ucov  , massebxy  , vorpot         )
77      CALL enercin ( vcov   , ucov  , vcont     , ucont  , ecin  )
78      CALL bernoui ( ip1jmp1, llm   , phi       , ecin   , bern  )
79
80      DO l=1,llm
81         DO ij=1,ip1jmp1
82            ang(ij,l) = ucov(ij,l) + constang(ij)
83         ENDDO
84      ENDDO
85
86      resetvarc=.true. ! force a recomputation of initial values in sortvarc
87      CALL sortvarc
88     $ ( itau,ucov,teta,ps,masse,pk,phis,vorpot,phi,bern,dp,time,vcov )
89
90      RETURN
91      END
Note: See TracBrowser for help on using the repository browser.