/[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 73 by guez, Fri Nov 15 17:48:30 2013 UTC
# Line 12  contains Line 12  contains
12      ! Objet: calcul des tendances dynamiques      ! Objet: calcul des tendances dynamiques
13      ! Modif 04/93 F.Forget      ! Modif 04/93 F.Forget
14    
15        USE comgeom, ONLY : airesurg
16      USE dimens_m, ONLY : llm      USE dimens_m, ONLY : llm
17        USE disvert_m, ONLY : ap, bp
18        use massdair_m, only: massdair
19      USE paramet_m, ONLY : iip1, ip1jm, ip1jmp1, jjp1, llmp1      USE paramet_m, ONLY : iip1, ip1jm, ip1jmp1, jjp1, llmp1
     USE comvert, ONLY : ap, bp  
     USE comgeom, ONLY : airesurg  
     USE pression_m, ONLY : pression  
20    
21      !   Arguments:      !   Arguments:
22      REAL, INTENT (IN) :: vcov(ip1jm, llm), ucov(ip1jmp1, llm)      REAL, INTENT (IN) :: vcov(ip1jm, llm), ucov(ip1jmp1, llm)
23      REAL :: teta(ip1jmp1, llm)      REAL, INTENT(IN):: teta(ip1jmp1, llm)
24      REAL, INTENT (IN) :: ps(ip1jmp1)      REAL, INTENT (IN) :: ps(ip1jmp1)
25      REAL, INTENT (IN) :: phis(ip1jmp1)      REAL, INTENT (IN) :: phis(ip1jmp1)
26      REAL, INTENT (IN) :: pk(iip1, jjp1, llm)      REAL, INTENT (IN) :: pk(iip1, jjp1, llm)
# Line 31  contains Line 31  contains
31      !   Local:      !   Local:
32    
33      REAL :: p(ip1jmp1, llmp1)      REAL :: p(ip1jmp1, llmp1)
34      REAL :: massebx(ip1jmp1, llm), masseby(ip1jm, llm), psexbarxy(ip1jm)      REAL :: massebx(ip1jmp1, llm), masseby(ip1jm, llm)
35      REAL :: vorpot(ip1jm, llm)      REAL :: vorpot(ip1jm, llm)
36      REAL :: w(ip1jmp1, llm), ecin(ip1jmp1, llm), convm(ip1jmp1, llm)      REAL :: w(ip1jmp1, llm), ecin(ip1jmp1, llm), convm(ip1jmp1, llm)
37      REAL :: bern(ip1jmp1, llm)      REAL :: bern(ip1jmp1, llm)
38      REAL :: massebxy(ip1jm, llm), dp(ip1jmp1)      REAL :: massebxy(ip1jm, llm), dp(ip1jmp1)
39    
40      INTEGER :: ij      INTEGER :: ij, l
41    
42      !-----------------------------------------------------------------------      !-----------------------------------------------------------------------
43    
# Line 46  contains Line 46  contains
46      !   Calcul des tendances dynamiques:      !   Calcul des tendances dynamiques:
47    
48      CALL covcont(llm, ucov, vcov, ucont, vcont)      CALL covcont(llm, ucov, vcov, ucont, vcont)
49      CALL pression(ip1jmp1, ap, bp, ps, p)      forall (l = 1: llm + 1) p(:, l) = ap(l) + bp(l) * ps
     CALL psextbar(ps, psexbarxy)  
50      CALL massdair(p, masse)      CALL massdair(p, masse)
51      CALL massbar(masse, massebx, masseby)      CALL massbar(masse, massebx, masseby)
52      CALL massbarxy(masse, massebxy)      CALL massbarxy(masse, massebxy)
Legend:
Removed from v.23  
changed lines
  Added in v.73
  ViewVC Help
Powered by ViewVC 1.1.21