/[lmdze]/trunk/dyn3d/caldyn0.f
ViewVC logotype

Diff of /trunk/dyn3d/caldyn0.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 23 by guez, Mon Dec 14 15:25:16 2009 UTC revision 46 by guez, Mon May 16 14:52:30 2011 UTC
# Line 16  contains Line 16  contains
16      USE paramet_m, ONLY : iip1, ip1jm, ip1jmp1, jjp1, llmp1      USE paramet_m, ONLY : iip1, ip1jm, ip1jmp1, jjp1, llmp1
17      USE comvert, ONLY : ap, bp      USE comvert, ONLY : ap, bp
18      USE comgeom, ONLY : airesurg      USE comgeom, ONLY : airesurg
     USE pression_m, ONLY : pression  
19    
20      !   Arguments:      !   Arguments:
21      REAL, INTENT (IN) :: vcov(ip1jm, llm), ucov(ip1jmp1, llm)      REAL, INTENT (IN) :: vcov(ip1jm, llm), ucov(ip1jmp1, llm)
# Line 31  contains Line 30  contains
30      !   Local:      !   Local:
31    
32      REAL :: p(ip1jmp1, llmp1)      REAL :: p(ip1jmp1, llmp1)
33      REAL :: massebx(ip1jmp1, llm), masseby(ip1jm, llm), psexbarxy(ip1jm)      REAL :: massebx(ip1jmp1, llm), masseby(ip1jm, llm)
34      REAL :: vorpot(ip1jm, llm)      REAL :: vorpot(ip1jm, llm)
35      REAL :: w(ip1jmp1, llm), ecin(ip1jmp1, llm), convm(ip1jmp1, llm)      REAL :: w(ip1jmp1, llm), ecin(ip1jmp1, llm), convm(ip1jmp1, llm)
36      REAL :: bern(ip1jmp1, llm)      REAL :: bern(ip1jmp1, llm)
37      REAL :: massebxy(ip1jm, llm), dp(ip1jmp1)      REAL :: massebxy(ip1jm, llm), dp(ip1jmp1)
38    
39      INTEGER :: ij      INTEGER :: ij, l
40    
41      !-----------------------------------------------------------------------      !-----------------------------------------------------------------------
42    
# Line 46  contains Line 45  contains
45      !   Calcul des tendances dynamiques:      !   Calcul des tendances dynamiques:
46    
47      CALL covcont(llm, ucov, vcov, ucont, vcont)      CALL covcont(llm, ucov, vcov, ucont, vcont)
48      CALL pression(ip1jmp1, ap, bp, ps, p)      forall (l = 1: llm + 1) p(:, l) = ap(l) + bp(l) * ps
     CALL psextbar(ps, psexbarxy)  
49      CALL massdair(p, masse)      CALL massdair(p, masse)
50      CALL massbar(masse, massebx, masseby)      CALL massbar(masse, massebx, masseby)
51      CALL massbarxy(masse, massebxy)      CALL massbarxy(masse, massebxy)

Legend:
Removed from v.23  
changed lines
  Added in v.46

  ViewVC Help
Powered by ViewVC 1.1.21