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

Diff of /trunk/dyn3d/leapfrog.f90

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

trunk/Sources/dyn3d/leapfrog.f revision 224 by guez, Fri Apr 28 13:40:59 2017 UTC trunk/dyn3d/leapfrog.f revision 261 by guez, Wed Mar 7 13:33:15 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    
80      ! Tendances de la dissipation :      ! Tendances de la dissipation :
81      REAL dvdis(iim + 1, jjm, llm), dudis(iim + 1, jjm + 1, llm)      REAL dvdis(iim + 1, jjm, llm), dudis(iim + 1, jjm + 1, llm)
# 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 204  contains Line 204  contains
204    
205         IF (MOD(itau + 1, iperiod) == 0) THEN         IF (MOD(itau + 1, iperiod) == 0) THEN
206            ! \'Ecriture du fichier histoire moyenne:            ! \'Ecriture du fichier histoire moyenne:
207            CALL writedynav(vcov, ucov, teta, pk, phi, q, masse, ps, phis, &            CALL writedynav(vcov, ucov, teta, pk, phi, q, masse, ps, &
208                 time = itau + 1)                 time = itau + 1)
209            call bilan_dyn(ps, masse, pk, pbaru, pbarv, teta, phi, ucov, vcov, &            call bilan_dyn(ps, masse, pk, pbaru, pbarv, teta, phi, ucov, vcov, &
210                 q(:, :, :, 1))                 q(:, :, :, 1))
# Line 212  contains Line 212  contains
212    
213         IF (MOD(itau + 1, iecri * day_step) == 0) THEN         IF (MOD(itau + 1, iecri * day_step) == 0) THEN
214            CALL geopot(teta, pk, pks, phis, phi)            CALL geopot(teta, pk, pks, phis, phi)
215            CALL writehist(itau, vcov, ucov, teta, phi, masse, ps)            CALL writehist(vcov, ucov, teta, pk, phi, q, masse, ps, itau)
216         END IF         END IF
217      end do time_integration      end do time_integration
218    
# 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.224  
changed lines
  Added in v.261
  ViewVC Help
Powered by ViewVC 1.1.21