/[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 47 by guez, Fri Jul 1 15:00:48 2011 UTC revision 56 by guez, Tue Jan 10 19:02:02 2012 UTC
# 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 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
          periodav  
22      USE dimens_m, ONLY: iim, jjm, llm, nqmx      USE dimens_m, ONLY: iim, jjm, llm, nqmx
23      use dissip_m, only: dissip      use dissip_m, only: dissip
24      USE dynetat0_m, ONLY: day_ini      USE dynetat0_m, ONLY: day_ini
# Line 32  contains Line 30  contains
30      use inidissip_m, only: idissip      use inidissip_m, only: idissip
31      use integrd_m, only: integrd      use integrd_m, only: integrd
32      USE logic, ONLY: iflag_phys, ok_guide      USE logic, ONLY: iflag_phys, ok_guide
33      USE paramet_m, ONLY: ip1jmp1      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(ip1jmp1, llm) ! vent covariant      REAL, intent(inout):: ucov(:, :, :) ! (iim + 1, jjm + 1, llm) vent covariant
40      REAL, intent(inout):: vcov((iim + 1) * jjm, llm) ! vent covariant      REAL, intent(inout):: vcov(:, :, :) ! (iim + 1, jjm, llm) ! vent covariant
41    
42      REAL, intent(inout):: teta(:, :, :) ! (iim + 1, jjm + 1, llm)      REAL, intent(inout):: teta(:, :, :) ! (iim + 1, jjm + 1, llm)
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(ip1jmp1, llm) ! masse d'air      REAL masse((iim + 1) * (jjm + 1), llm) ! masse d'air
47      REAL phis(ip1jmp1) ! 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 56  contains Line 55  contains
55    
56      ! Variables dynamiques:      ! Variables dynamiques:
57    
58      REAL pks(ip1jmp1) ! 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(ip1jmp1, llm) ! exner filt.au milieu des couches      REAL pkf((iim + 1) * (jjm + 1), llm) ! exner filt.au milieu des couches
61      REAL phi(iim + 1, jjm + 1, llm) ! geopotential      REAL phi(iim + 1, jjm + 1, llm) ! geopotential
62      REAL w(ip1jmp1, llm) ! vitesse verticale      REAL w((iim + 1) * (jjm + 1), llm) ! vitesse verticale
63    
64      ! variables dynamiques intermediaire pour le transport      ! Variables dynamiques intermediaire pour le transport
65      REAL pbaru(ip1jmp1, llm), pbarv((iim + 1) * jjm, llm) !flux de masse      ! Flux de masse :
66        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(ip1jmp1, 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(ip1jmp1, llm)      REAL massem1((iim + 1) * (jjm + 1), llm)
72    
73      ! tendances dynamiques      ! Tendances dynamiques
74      REAL dv((iim + 1) * jjm, llm), dudyn(ip1jmp1, llm)      REAL dv((iim + 1) * jjm, llm), dudyn((iim + 1) * (jjm + 1), llm)
75      REAL dteta(iim + 1, jjm + 1, llm), dq(ip1jmp1, llm, nqmx), dp(ip1jmp1)      REAL dteta(iim + 1, jjm + 1, llm), dq((iim + 1) * (jjm + 1), llm, nqmx)
76        real dp((iim + 1) * (jjm + 1))
77    
78      ! tendances de la dissipation      ! Tendances de la dissipation :
79      REAL dvdis((iim + 1) * jjm, llm), dudis(ip1jmp1, 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(ip1jmp1, llm)      REAL dvfi((iim + 1) * jjm, llm), dufi((iim + 1) * (jjm + 1), llm)
84      REAL dtetafi(iim + 1, jjm + 1, llm), dqfi(ip1jmp1, llm, nqmx), dpfi(ip1jmp1)      REAL dtetafi(iim + 1, jjm + 1, llm), dqfi((iim + 1) * (jjm + 1), llm, nqmx)
85        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(ip1jmp1, 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    
99      REAL dtetaecdt(iim + 1, jjm + 1, llm)      REAL vcont((iim + 1) * jjm, llm), ucont((iim + 1) * (jjm + 1), llm)
     ! tendance de la température potentielle due à la tansformation  
     ! d'énergie cinétique en énergie thermique créée par la dissipation  
   
     REAL vcont((iim + 1) * jjm, llm), ucont(ip1jmp1, llm)  
100      logical leapf      logical leapf
101      real dt      real dt
102    
103      !---------------------------------------------------      !---------------------------------------------------
104    
105      print *, "Call sequence information: leapfrog"      print *, "Call sequence information: leapfrog"
106        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 "itaufin" is one too
# Line 134  contains Line 133  contains
133         end if         end if
134    
135         ! Calcul des tendances dynamiques:         ! Calcul des tendances dynamiques:
136         CALL geopot(ip1jmp1, teta, pk, pks, phis, phi)         CALL geopot((iim + 1) * (jjm + 1), teta, pk, pks, phis, phi)
137         CALL caldyn(itau, ucov, vcov, teta, ps, masse, pk, pkf, phis, phi, &         CALL caldyn(itau, ucov, vcov, teta, ps, masse, pk, pkf, phis, phi, &
138              dudyn, dv, dteta, dp, w, pbaru, pbarv, time_0, &              dudyn, dv, dteta, dp, w, pbaru, pbarv, time_0, &
139              conser=MOD(itau, iconser)==0)              conser=MOD(itau, iconser)==0)
# Line 157  contains Line 156  contains
156            CALL exner_hyb(ps, p3d, pks, pk, pkf)            CALL exner_hyb(ps, p3d, pks, pk, pkf)
157    
158            ! Calcul des tendances dynamiques:            ! Calcul des tendances dynamiques:
159            CALL geopot(ip1jmp1, teta, pk, pks, phis, phi)            CALL geopot((iim + 1) * (jjm + 1), teta, pk, pks, phis, phi)
160            CALL caldyn(itau + 1, ucov, vcov, teta, ps, masse, pk, pkf, phis, &            CALL caldyn(itau + 1, ucov, vcov, teta, ps, masse, pk, pkf, phis, &
161                 phi, dudyn, dv, dteta, dp, w, pbaru, pbarv, time_0, &                 phi, dudyn, dv, dteta, dp, w, pbaru, pbarv, time_0, &
162                 conser=.false.)                 conser=.false.)
# Line 192  contains Line 191  contains
191         CALL exner_hyb(ps, p3d, pks, pk, pkf)         CALL exner_hyb(ps, p3d, pks, pk, pkf)
192    
193         IF (MOD(itau + 1, idissip) == 0) THEN         IF (MOD(itau + 1, idissip) == 0) THEN
194            ! dissipation horizontale et verticale des petites echelles:            ! Dissipation horizontale et verticale des petites échelles
195    
196            ! calcul de l'energie cinetique avant dissipation            ! calcul de l'énergie cinétique avant dissipation
197            call covcont(llm, ucov, vcov, ucont, vcont)            call covcont(llm, ucov, vcov, ucont, vcont)
198            call enercin(vcov, ucov, vcont, ucont, ecin0)            call enercin(vcov, ucov, vcont, ucont, ecin0)
199    
200            ! dissipation            ! dissipation
201            CALL dissip(vcov, ucov, teta, p3d, dvdis, dudis, dtetadis)            CALL dissip(vcov, ucov, teta, p3d, dvdis, dudis, dtetadis)
202            ucov=ucov + dudis            ucov = ucov + dudis
203            vcov=vcov + dvdis            vcov = vcov + dvdis
204    
205            ! On rajoute la tendance due à la transformation Ec -> E            ! On ajoute la tendance due à la transformation énergie
206            ! thermique créée lors de 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
210            dtetadis=dtetadis + dtetaecdt            teta = teta + dtetadis
           teta=teta + dtetadis  
211    
212            ! Calcul de la valeur moyenne aux pôles :            ! Calcul de la valeur moyenne aux pôles :
213            forall (l = 1: llm)            forall (l = 1: llm)
# Line 226  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(nqmx, itau + 1, vcov, ucov, teta, pk, phi, q, &
228                 phi, q, masse, ps, phis)                 masse, ps, phis)
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_app = dtvr * iperiod)
                dt_cum = dtvr * day_step * periodav)  
231         ENDIF         ENDIF
232      end do time_integration      end do time_integration
233    
# Line 238  contains Line 235  contains
235           itau=itau_dyn+itaufin)           itau=itau_dyn+itaufin)
236    
237      ! Calcul des tendances dynamiques:      ! Calcul des tendances dynamiques:
238      CALL geopot(ip1jmp1, 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.47  
changed lines
  Added in v.56
  ViewVC Help
Powered by ViewVC 1.1.21