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

Diff of /trunk/dyn3d/leapfrog.f

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

revision 254 by guez, Mon Feb 5 10:39:38 2018 UTC revision 260 by guez, Tue Mar 6 17:18:33 2018 UTC
# Line 63  contains Line 63  contains
63      REAL phi(iim + 1, jjm + 1, llm) ! geopotential      REAL phi(iim + 1, jjm + 1, llm) ! geopotential
64      REAL w(iim + 1, jjm + 1, llm) ! vitesse verticale      REAL w(iim + 1, jjm + 1, llm) ! vitesse verticale
65    
66      ! Variables dynamiques intermediaire pour le transport      ! Variables dynamiques interm\'ediaires pour le transport
67      ! Flux de masse :      ! Flux de masse :
68      REAL pbaru(iim + 1, jjm + 1, llm), pbarv(iim + 1, jjm, llm)      REAL pbaru(iim + 1, jjm + 1, llm), pbarv(iim + 1, jjm, llm)
69    
# Line 73  contains Line 73  contains
73      REAL massem1(iim + 1, jjm + 1, llm)      REAL massem1(iim + 1, jjm + 1, llm)
74    
75      ! Tendances dynamiques      ! Tendances dynamiques
76      REAL dv((iim + 1) * jjm, llm), dudyn(iim + 1, jjm + 1, llm)      REAL dv((iim + 1) * jjm, llm), du(iim + 1, jjm + 1, llm)
77      REAL dteta(iim + 1, jjm + 1, llm)      REAL dteta(iim + 1, jjm + 1, llm)
78      real dp(iim + 1, jjm + 1)      real dp(iim + 1, jjm + 1)
79    
# Line 133  contains Line 133  contains
133         ! Calcul des tendances dynamiques:         ! Calcul des tendances dynamiques:
134         CALL geopot(teta, pk, pks, phis, phi)         CALL geopot(teta, pk, pks, phis, phi)
135         CALL caldyn(itau, ucov, vcov, teta, ps, masse, pk, pkf, phis, phi, &         CALL caldyn(itau, ucov, vcov, teta, ps, masse, pk, pkf, phis, phi, &
136              dudyn, dv, dteta, dp, w, pbaru, pbarv, &              du, dv, dteta, dp, w, pbaru, pbarv, &
137              conser = MOD(itau, iconser) == 0)              conser = MOD(itau, iconser) == 0)
138    
139         CALL caladvtrac(q, pbaru, pbarv, p3d, masse, teta, pk)         CALL caladvtrac(q, pbaru, pbarv, p3d, masse, teta, pk)
140    
141         ! Int\'egrations dynamique et traceurs:         ! Int\'egrations dynamique et traceurs:
142         CALL integrd(vcovm1, ucovm1, tetam1, psm1, massem1, dv, dudyn, dteta, &         CALL integrd(vcovm1, ucovm1, tetam1, psm1, massem1, dv, du, dteta, &
143              dp, vcov, ucov, teta, q(:, :, :, :2), ps, masse, dt, leapf)              dp, vcov, ucov, teta, q(:, :, :, :2), ps, masse, dt, leapf)
144    
145         forall (l = 1: llm + 1) p3d(:, :, l) = ap(l) + bp(l) * ps         forall (l = 1: llm + 1) p3d(:, :, l) = ap(l) + bp(l) * ps
# Line 152  contains Line 152  contains
152            ! Calcul des tendances dynamiques:            ! Calcul des tendances dynamiques:
153            CALL geopot(teta, pk, pks, phis, phi)            CALL geopot(teta, pk, pks, phis, phi)
154            CALL caldyn(itau + 1, ucov, vcov, teta, ps, masse, pk, pkf, phis, &            CALL caldyn(itau + 1, ucov, vcov, teta, ps, masse, pk, pkf, phis, &
155                 phi, dudyn, dv, dteta, dp, w, pbaru, pbarv, conser = .false.)                 phi, du, dv, dteta, dp, w, pbaru, pbarv, conser = .false.)
156    
157            ! integrations dynamique et traceurs:            ! integrations dynamique et traceurs:
158            CALL integrd(vcovm1, ucovm1, tetam1, psm1, massem1, dv, dudyn, &            CALL integrd(vcovm1, ucovm1, tetam1, psm1, massem1, dv, du, &
159                 dteta, dp, vcov, ucov, teta, q(:, :, :, :2), ps, masse, dtvr, &                 dteta, dp, vcov, ucov, teta, q(:, :, :, :2), ps, masse, dtvr, &
160                 leapf=.false.)                 leapf=.false.)
161    
# Line 221  contains Line 221  contains
221      ! Calcul des tendances dynamiques:      ! Calcul des tendances dynamiques:
222      CALL geopot(teta, pk, pks, phis, phi)      CALL geopot(teta, pk, pks, phis, phi)
223      CALL caldyn(itaufin, ucov, vcov, teta, ps, masse, pk, pkf, phis, phi, &      CALL caldyn(itaufin, ucov, vcov, teta, ps, masse, pk, pkf, phis, phi, &
224           dudyn, dv, dteta, dp, w, pbaru, pbarv, &           du, dv, dteta, dp, w, pbaru, pbarv, &
225           conser = MOD(itaufin, iconser) == 0)           conser = MOD(itaufin, iconser) == 0)
226    
227    END SUBROUTINE leapfrog    END SUBROUTINE leapfrog
Legend:
Removed from v.254  
changed lines
  Added in v.260
  ViewVC Help
Powered by ViewVC 1.1.21