/[lmdze]/trunk/Sources/dyn3d/leapfrog.f
ViewVC logotype

Diff of /trunk/Sources/dyn3d/leapfrog.f

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

revision 55 by guez, Mon Dec 12 13:25:01 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 35  contains Line 33  contains
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 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
# Line 97  contains Line 96  contains
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 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(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 246  contains Line 239  contains
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.56
  ViewVC Help
Powered by ViewVC 1.1.21