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

Diff of /trunk/dyn3d/leapfrog.f90

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

revision 55 by guez, Mon Dec 12 13:25:01 2011 UTC revision 67 by guez, Tue Oct 2 15:50:56 2012 UTC
# Line 6  contains Line 6  contains
6    
7    SUBROUTINE leapfrog(ucov, vcov, teta, ps, masse, phis, q, time_0)    SUBROUTINE leapfrog(ucov, vcov, teta, ps, masse, phis, q, time_0)
8    
9      ! From dyn3d/leapfrog.F, version 1.6, 2005/04/13 08:58:34      ! From dyn3d/leapfrog.F, version 1.6, 2005/04/13 08:58:34 revision 616
10      ! Authors: P. Le Van, L. Fairhead, F. Hourdin      ! Authors: P. Le Van, L. Fairhead, F. Hourdin
11      ! Matsuno-leapfrog scheme.      ! Matsuno-leapfrog scheme.
12    
# Line 15  contains Line 15  contains
15      use caladvtrac_m, only: caladvtrac      use caladvtrac_m, only: caladvtrac
16      use caldyn_m, only: caldyn      use caldyn_m, only: caldyn
17      USE calfis_m, ONLY: calfis      USE calfis_m, ONLY: calfis
     USE com_io_dyn, ONLY: histaveid  
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 disvert_m, 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           periodav           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 31  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
36        use writedynav_m, only: writedynav
37    
38      ! Variables dynamiques:      ! Variables dynamiques:
39      REAL, intent(inout):: ucov(:, :, :) ! (iim + 1, jjm + 1, llm) vent covariant      REAL, intent(inout):: ucov(:, :, :) ! (iim + 1, jjm + 1, llm) vent covariant
# Line 44  contains Line 43  contains
43      ! potential temperature      ! potential temperature
44    
45      REAL, intent(inout):: ps(:, :) ! (iim + 1, jjm + 1) pression au sol, en Pa      REAL, intent(inout):: ps(:, :) ! (iim + 1, jjm + 1) pression au sol, en Pa
46      REAL masse((iim + 1) * (jjm + 1), llm) ! masse d'air      REAL masse(:, :, :) ! (iim + 1, jjm + 1, llm) masse d'air
47      REAL phis((iim + 1) * (jjm + 1)) ! geopotentiel au sol      REAL phis(:, :) ! (iim + 1, jjm + 1) geopotentiel au sol
48    
49      REAL, intent(inout):: q(:, :, :, :) ! (iim + 1, jjm + 1, llm, nqmx)      REAL, intent(inout):: q(:, :, :, :) ! (iim + 1, jjm + 1, llm, nqmx)
50      ! mass fractions of advected fields      ! mass fractions of advected fields
# Line 58  contains Line 57  contains
57    
58      REAL pks((iim + 1) * (jjm + 1)) ! exner au sol      REAL pks((iim + 1) * (jjm + 1)) ! exner au sol
59      REAL pk(iim + 1, jjm + 1, llm) ! exner au milieu des couches      REAL pk(iim + 1, jjm + 1, llm) ! exner au milieu des couches
60      REAL pkf((iim + 1) * (jjm + 1), llm) ! exner filt.au milieu des couches      REAL pkf(iim + 1, jjm + 1, llm) ! exner filtré au milieu des couches
61      REAL phi(iim + 1, jjm + 1, llm) ! geopotential      REAL phi(iim + 1, jjm + 1, llm) ! geopotential
62      REAL w((iim + 1) * (jjm + 1), llm) ! vitesse verticale      REAL w((iim + 1) * (jjm + 1), llm) ! vitesse verticale
63    
# Line 66  contains Line 65  contains
65      ! Flux de masse :      ! Flux de masse :
66      REAL pbaru((iim + 1) * (jjm + 1), llm), pbarv((iim + 1) * jjm, llm)      REAL pbaru((iim + 1) * (jjm + 1), llm), pbarv((iim + 1) * jjm, llm)
67    
68      ! variables dynamiques au pas - 1      ! Variables dynamiques au pas - 1
69      REAL vcovm1(iim + 1, jjm, llm), ucovm1(iim + 1, jjm + 1, llm)      REAL vcovm1(iim + 1, jjm, llm), ucovm1(iim + 1, jjm + 1, llm)
70      REAL tetam1(iim + 1, jjm + 1, llm), psm1(iim + 1, jjm + 1)      REAL tetam1(iim + 1, jjm + 1, llm), psm1(iim + 1, jjm + 1)
71      REAL massem1((iim + 1) * (jjm + 1), llm)      REAL massem1(iim + 1, jjm + 1, llm)
72    
73      ! tendances dynamiques      ! Tendances dynamiques
74      REAL dv((iim + 1) * jjm, llm), dudyn((iim + 1) * (jjm + 1), llm)      REAL dv((iim + 1) * jjm, llm), dudyn((iim + 1) * (jjm + 1), llm)
75      REAL dteta(iim + 1, jjm + 1, llm), dq((iim + 1) * (jjm + 1), llm, nqmx)      REAL dteta(iim + 1, jjm + 1, llm), dq((iim + 1) * (jjm + 1), llm, nqmx)
76      real dp((iim + 1) * (jjm + 1))      real dp((iim + 1) * (jjm + 1))
77    
78      ! tendances de la dissipation      ! Tendances de la dissipation :
79      REAL dvdis(iim + 1, jjm, llm), dudis(iim + 1, jjm + 1, llm)      REAL dvdis(iim + 1, jjm, llm), dudis(iim + 1, jjm + 1, llm)
80      REAL dtetadis(iim + 1, jjm + 1, llm)      REAL dtetadis(iim + 1, jjm + 1, llm)
81    
82      ! tendances physiques      ! Tendances physiques
83      REAL dvfi((iim + 1) * jjm, llm), dufi((iim + 1) * (jjm + 1), llm)      REAL dvfi((iim + 1) * jjm, llm), dufi((iim + 1) * (jjm + 1), llm)
84      REAL dtetafi(iim + 1, jjm + 1, llm), dqfi((iim + 1) * (jjm + 1), llm, nqmx)      REAL dtetafi(iim + 1, jjm + 1, llm), dqfi((iim + 1) * (jjm + 1), llm, nqmx)
85      real dpfi((iim + 1) * (jjm + 1))      real dpfi((iim + 1) * (jjm + 1))
86    
87      ! variables pour le fichier histoire      ! Variables pour le fichier histoire
88    
89      INTEGER itau ! index of the time step of the dynamics, starts at 0      INTEGER itau ! index of the time step of the dynamics, starts at 0
90      INTEGER itaufin      INTEGER itaufin
91      REAL time ! time of day, as a fraction of day length      REAL time ! time of day, as a fraction of day length
92      real finvmaold((iim + 1) * (jjm + 1), llm)      real finvmaold(iim + 1, jjm + 1, llm)
93      INTEGER l      INTEGER l
94      REAL rdayvrai, rdaym_ini      REAL rdayvrai, rdaym_ini
95    
96      ! Variables test conservation energie      ! Variables test conservation energie
97      REAL ecin(iim + 1, jjm + 1, llm), ecin0(iim + 1, jjm + 1, llm)      REAL ecin(iim + 1, jjm + 1, llm), ecin0(iim + 1, jjm + 1, llm)
98    
     REAL dtetaecdt(iim + 1, jjm + 1, llm)  
     ! tendance de la température potentielle due à la transformation  
     ! d'énergie cinétique en énergie thermique par la dissipation  
   
99      REAL vcont((iim + 1) * jjm, llm), ucont((iim + 1) * (jjm + 1), llm)      REAL vcont((iim + 1) * jjm, llm), ucont((iim + 1) * (jjm + 1), llm)
100      logical leapf      logical leapf
101      real dt      real dt
# Line 111  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 134  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 185  contains Line 180  contains
180    
181            CALL calfis(rdayvrai, time, ucov, vcov, teta, q, masse, ps, pk, &            CALL calfis(rdayvrai, time, ucov, vcov, teta, q, masse, ps, pk, &
182                 phis, phi, dudyn, dv, dq, w, dufi, dvfi, dtetafi, dqfi, dpfi, &                 phis, phi, dudyn, dv, dq, w, dufi, dvfi, dtetafi, dqfi, dpfi, &
183                 lafin=itau+1==itaufin)                 lafin = itau + 1 == itaufin)
184    
185            ! ajout des tendances physiques:            ! ajout des tendances physiques:
186            CALL addfi(nqmx, dtphys, ucov, vcov, teta, q, ps, dufi, dvfi, &            CALL addfi(nqmx, dtphys, ucov, vcov, teta, q, ps, dufi, dvfi, &
# Line 211  contains Line 206  contains
206            ! cinétique en énergie thermique par la dissipation            ! cinétique en énergie thermique par la dissipation
207            call covcont(llm, ucov, vcov, ucont, vcont)            call covcont(llm, ucov, vcov, ucont, vcont)
208            call enercin(vcov, ucov, vcont, ucont, ecin)            call enercin(vcov, ucov, vcont, ucont, ecin)
209            dtetaecdt= (ecin0 - ecin) / pk            dtetadis = dtetadis + (ecin0 - ecin) / pk
           dtetadis = dtetadis + dtetaecdt  
210            teta = teta + dtetadis            teta = teta + dtetadis
211    
212            ! Calcul de la valeur moyenne aux pôles :            ! Calcul de la valeur moyenne aux pôles :
# Line 230  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(histaveid, nqmx, itau + 1, vcov, ucov, teta, pk, &            CALL writedynav(vcov, ucov, teta, pk, phi, q, masse, ps, phis, &
228                 phi, q, 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))
                dt_cum = dtvr * day_step * periodav)  
231         ENDIF         ENDIF
232      end do time_integration      end do time_integration
233    
234      CALL dynredem1("restart.nc", vcov, ucov, teta, q, masse, ps, &      CALL dynredem1("restart.nc", vcov, ucov, teta, q, masse, ps, &
235           itau=itau_dyn+itaufin)           itau = itau_dyn + itaufin)
236    
237      ! Calcul des tendances dynamiques:      ! Calcul des tendances dynamiques:
238      CALL geopot((iim + 1) * (jjm + 1), teta, pk, pks, phis, phi)      CALL geopot((iim + 1) * (jjm + 1), teta, pk, pks, phis, phi)
239      CALL caldyn(itaufin, ucov, vcov, teta, ps, masse, pk, pkf, phis, phi, &      CALL caldyn(itaufin, ucov, vcov, teta, ps, masse, pk, pkf, phis, phi, &
240           dudyn, dv, dteta, dp, w, pbaru, pbarv, time_0, &           dudyn, dv, dteta, dp, w, pbaru, pbarv, time_0, &
241           conser=MOD(itaufin, iconser)==0)           conser = MOD(itaufin, iconser) == 0)
242    
243    END SUBROUTINE leapfrog    END SUBROUTINE leapfrog
244    
245  end module leapfrog_m  end module leapfrog_m
Legend:
Removed from v.55  
changed lines
  Added in v.67
  ViewVC Help
Powered by ViewVC 1.1.21