/[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 56 by guez, Tue Jan 10 19:02:02 2012 UTC revision 64 by guez, Wed Aug 29 14:47:17 2012 UTC
# Line 18  contains Line 18  contains
18      USE comconst, ONLY: daysec, dtphys, dtvr      USE comconst, ONLY: daysec, dtphys, dtvr
19      USE comgeom, ONLY: aire_2d, apoln, apols      USE comgeom, ONLY: aire_2d, apoln, apols
20      USE comvert, ONLY: ap, bp      USE comvert, ONLY: ap, bp
21      USE conf_gcm_m, ONLY: day_step, iconser, iperiod, iphysiq, nday, offline      USE conf_gcm_m, ONLY: day_step, iconser, iperiod, iphysiq, nday, offline, &
22             iflag_phys, ok_guide
23      USE dimens_m, ONLY: iim, jjm, llm, nqmx      USE dimens_m, ONLY: iim, jjm, llm, nqmx
24      use dissip_m, only: dissip      use dissip_m, only: dissip
25      USE dynetat0_m, ONLY: day_ini      USE dynetat0_m, ONLY: day_ini
# Line 29  contains Line 30  contains
30      USE guide_m, ONLY: guide      USE guide_m, ONLY: guide
31      use inidissip_m, only: idissip      use inidissip_m, only: idissip
32      use integrd_m, only: integrd      use integrd_m, only: integrd
     USE logic, ONLY: iflag_phys, ok_guide  
33      use nr_util, only: assert      use nr_util, only: assert
34      USE pressure_var, ONLY: p3d      USE pressure_var, ONLY: p3d
35      USE temps, ONLY: itau_dyn      USE temps, ONLY: itau_dyn
# Line 106  contains Line 106  contains
106      call assert(shape(ucov) == (/iim + 1, jjm + 1, llm/), "leapfrog")      call assert(shape(ucov) == (/iim + 1, jjm + 1, llm/), "leapfrog")
107    
108      itaufin = nday * day_step      itaufin = nday * day_step
109      ! "day_step" is a multiple of "iperiod", therefore "itaufin" is one too      ! "day_step" is a multiple of "iperiod", therefore so is "itaufin".
110    
111      dq = 0.      dq = 0.
112    
# Line 129  contains Line 129  contains
129            massem1 = masse            massem1 = masse
130            psm1 = ps            psm1 = ps
131            finvmaold = masse            finvmaold = masse
132            CALL filtreg(finvmaold, jjm + 1, llm, - 2, 2, .TRUE., 1)            CALL filtreg(finvmaold, jjm + 1, llm, - 2, 2, .TRUE.)
133         end if         end if
134    
135         ! Calcul des tendances dynamiques:         ! Calcul des tendances dynamiques:
# Line 224  contains Line 224  contains
224    
225         IF (MOD(itau + 1, iperiod) == 0) THEN         IF (MOD(itau + 1, iperiod) == 0) THEN
226            ! Écriture du fichier histoire moyenne:            ! Écriture du fichier histoire moyenne:
227            CALL writedynav(nqmx, itau + 1, vcov, ucov, teta, pk, phi, q, &            CALL writedynav(vcov, ucov, teta, pk, phi, q, masse, ps, phis, &
228                 masse, ps, phis)                 time = itau + 1)
229            call bilan_dyn(ps, masse, pk, pbaru, pbarv, teta, phi, ucov, vcov, &            call bilan_dyn(ps, masse, pk, pbaru, pbarv, teta, phi, ucov, vcov, &
230                 q(:, :, :, 1), dt_app = dtvr * iperiod)                 q(:, :, :, 1))
231         ENDIF         ENDIF
232      end do time_integration      end do time_integration
233    
Legend:
Removed from v.56  
changed lines
  Added in v.64
  ViewVC Help
Powered by ViewVC 1.1.21