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

Diff of /trunk/libf/dyn3d/leapfrog.f90

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

revision 46 by guez, Mon May 16 14:52:30 2011 UTC revision 47 by guez, Fri Jul 1 15:00:48 2011 UTC
# Line 22  contains Line 22  contains
22      USE conf_gcm_m, ONLY: day_step, iconser, iperiod, iphysiq, nday, offline, &      USE conf_gcm_m, ONLY: day_step, iconser, iperiod, iphysiq, nday, offline, &
23           periodav           periodav
24      USE dimens_m, ONLY: iim, jjm, llm, nqmx      USE dimens_m, ONLY: iim, jjm, llm, nqmx
25        use dissip_m, only: dissip
26      USE dynetat0_m, ONLY: day_ini      USE dynetat0_m, ONLY: day_ini
27      use dynredem1_m, only: dynredem1      use dynredem1_m, only: dynredem1
28      USE exner_hyb_m, ONLY: exner_hyb      USE exner_hyb_m, ONLY: exner_hyb
# Line 58  contains Line 59  contains
59      REAL pks(ip1jmp1) ! exner au sol      REAL pks(ip1jmp1) ! exner au sol
60      REAL pk(iim + 1, jjm + 1, llm) ! exner au milieu des couches      REAL pk(iim + 1, jjm + 1, llm) ! exner au milieu des couches
61      REAL pkf(ip1jmp1, llm) ! exner filt.au milieu des couches      REAL pkf(ip1jmp1, llm) ! exner filt.au milieu des couches
62      REAL phi(ip1jmp1, llm) ! geopotential      REAL phi(iim + 1, jjm + 1, llm) ! geopotential
63      REAL w(ip1jmp1, llm) ! vitesse verticale      REAL w(ip1jmp1, llm) ! vitesse verticale
64    
65      ! variables dynamiques intermediaire pour le transport      ! variables dynamiques intermediaire pour le transport
# Line 70  contains Line 71  contains
71      REAL massem1(ip1jmp1, llm)      REAL massem1(ip1jmp1, llm)
72    
73      ! tendances dynamiques      ! tendances dynamiques
74      REAL dv((iim + 1) * jjm, llm), du(ip1jmp1, llm)      REAL dv((iim + 1) * jjm, llm), dudyn(ip1jmp1, llm)
75      REAL dteta(iim + 1, jjm + 1, llm), dq(ip1jmp1, llm, nqmx), dp(ip1jmp1)      REAL dteta(iim + 1, jjm + 1, llm), dq(ip1jmp1, llm, nqmx), dp(ip1jmp1)
76    
77      ! tendances de la dissipation      ! tendances de la dissipation
# Line 135  contains Line 136  contains
136         ! Calcul des tendances dynamiques:         ! Calcul des tendances dynamiques:
137         CALL geopot(ip1jmp1, teta, pk, pks, phis, phi)         CALL geopot(ip1jmp1, teta, pk, pks, phis, phi)
138         CALL caldyn(itau, ucov, vcov, teta, ps, masse, pk, pkf, phis, phi, &         CALL caldyn(itau, ucov, vcov, teta, ps, masse, pk, pkf, phis, phi, &
139              MOD(itau, iconser) == 0, du, dv, dteta, dp, w, pbaru, pbarv, &              dudyn, dv, dteta, dp, w, pbaru, pbarv, time_0, &
140              time_0)              conser=MOD(itau, iconser)==0)
141    
142         ! Calcul des tendances advection des traceurs (dont l'humidité)         ! Calcul des tendances advection des traceurs (dont l'humidité)
143         CALL caladvtrac(q, pbaru, pbarv, p3d, masse, dq, teta, pk)         CALL caladvtrac(q, pbaru, pbarv, p3d, masse, dq, teta, pk)
# Line 146  contains Line 147  contains
147              dtvr, itau)              dtvr, itau)
148    
149         ! Integrations dynamique et traceurs:         ! Integrations dynamique et traceurs:
150         CALL integrd(vcovm1, ucovm1, tetam1, psm1, massem1, dv, du, dteta, dp, &         CALL integrd(vcovm1, ucovm1, tetam1, psm1, massem1, dv, dudyn, dteta, &
151              vcov, ucov, teta, q(:, :, :, :2), ps, masse, finvmaold, dt, leapf)              dp, vcov, ucov, teta, q(:, :, :, :2), ps, masse, finvmaold, dt, &
152                leapf)
153    
154         if (.not. leapf) then         if (.not. leapf) then
155            ! Matsuno backward            ! Matsuno backward
# Line 157  contains Line 159  contains
159            ! Calcul des tendances dynamiques:            ! Calcul des tendances dynamiques:
160            CALL geopot(ip1jmp1, teta, pk, pks, phis, phi)            CALL geopot(ip1jmp1, teta, pk, pks, phis, phi)
161            CALL caldyn(itau + 1, ucov, vcov, teta, ps, masse, pk, pkf, phis, &            CALL caldyn(itau + 1, ucov, vcov, teta, ps, masse, pk, pkf, phis, &
162                 phi, .false., du, dv, dteta, dp, w, pbaru, pbarv, time_0)                 phi, dudyn, dv, dteta, dp, w, pbaru, pbarv, time_0, &
163                   conser=.false.)
164    
165            ! integrations dynamique et traceurs:            ! integrations dynamique et traceurs:
166            CALL integrd(vcovm1, ucovm1, tetam1, psm1, massem1, dv, du, dteta, &            CALL integrd(vcovm1, ucovm1, tetam1, psm1, massem1, dv, dudyn, &
167                 dp, vcov, ucov, teta, q(:, :, :, :2), ps, masse, finvmaold, &                 dteta, dp, vcov, ucov, teta, q(:, :, :, :2), ps, masse, &
168                 dtvr, leapf=.false.)                 finvmaold, dtvr, leapf=.false.)
169         end if         end if
170    
171         IF (MOD(itau + 1, iphysiq) == 0 .AND. iflag_phys /= 0) THEN         IF (MOD(itau + 1, iphysiq) == 0 .AND. iflag_phys /= 0) THEN
# Line 177  contains Line 180  contains
180            IF (time > 1.) time = time - 1.            IF (time > 1.) time = time - 1.
181    
182            CALL calfis(rdayvrai, time, ucov, vcov, teta, q, masse, ps, pk, &            CALL calfis(rdayvrai, time, ucov, vcov, teta, q, masse, ps, pk, &
183                 phis, phi, du, dv, dq, w, dufi, dvfi, dtetafi, dqfi, dpfi, &                 phis, phi, dudyn, dv, dq, w, dufi, dvfi, dtetafi, dqfi, dpfi, &
184                 lafin=itau+1==itaufin)                 lafin=itau+1==itaufin)
185    
186            ! ajout des tendances physiques:            ! ajout des tendances physiques:
# Line 237  contains Line 240  contains
240      ! Calcul des tendances dynamiques:      ! Calcul des tendances dynamiques:
241      CALL geopot(ip1jmp1, teta, pk, pks, phis, phi)      CALL geopot(ip1jmp1, teta, pk, pks, phis, phi)
242      CALL caldyn(itaufin, ucov, vcov, teta, ps, masse, pk, pkf, phis, phi, &      CALL caldyn(itaufin, ucov, vcov, teta, ps, masse, pk, pkf, phis, phi, &
243           MOD(itaufin, iconser) == 0, du, dv, dteta, dp, w, pbaru, pbarv, &           dudyn, dv, dteta, dp, w, pbaru, pbarv, time_0, &
244           time_0)           conser=MOD(itaufin, iconser)==0)
   
245    END SUBROUTINE leapfrog    END SUBROUTINE leapfrog
246    
247  end module leapfrog_m  end module leapfrog_m
Legend:
Removed from v.46  
changed lines
  Added in v.47
  ViewVC Help
Powered by ViewVC 1.1.21