/[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 27 by guez, Thu Mar 25 14:29:07 2010 UTC revision 28 by guez, Fri Mar 26 18:33:04 2010 UTC
# Line 27  contains Line 27  contains
27      USE paramet_m, ONLY: iip1, ip1jm, ip1jmp1, jjp1      USE paramet_m, ONLY: iip1, ip1jm, ip1jmp1, jjp1
28      USE pression_m, ONLY: pression      USE pression_m, ONLY: pression
29      USE pressure_var, ONLY: p3d      USE pressure_var, ONLY: p3d
30      USE temps, ONLY: dt, itau_dyn      USE temps, ONLY: itau_dyn
31    
32      ! Variables dynamiques:      ! Variables dynamiques:
33      REAL vcov(ip1jm, llm), ucov(ip1jmp1, llm) ! vents covariants      REAL vcov(ip1jm, llm), ucov(ip1jmp1, llm) ! vents covariants
# Line 91  contains Line 91  contains
91      REAL dtetaecdt(ip1jmp1, llm)      REAL dtetaecdt(ip1jmp1, llm)
92      REAL vcont(ip1jm, llm), ucont(ip1jmp1, llm)      REAL vcont(ip1jm, llm), ucont(ip1jmp1, llm)
93      logical forward, leapf      logical forward, leapf
94        REAL dt
95    
96      !---------------------------------------------------      !---------------------------------------------------
97    
# Line 140  contains Line 141  contains
141            ! integrations dynamique et traceurs:            ! integrations dynamique et traceurs:
142            CALL integrd(2, vcovm1, ucovm1, tetam1, psm1, massem1, dv, du, &            CALL integrd(2, vcovm1, ucovm1, tetam1, psm1, massem1, dv, du, &
143                 dteta, dq, dp, vcov, ucov, teta, q, ps, masse, phis, &                 dteta, dq, dp, vcov, ucov, teta, q, ps, masse, phis, &
144                 finvmaold, leapf)                 finvmaold, leapf, dt)
145    
146            IF (MOD(itau + 1, iphysiq) == 0 .AND. iflag_phys /= 0) THEN            IF (MOD(itau + 1, iphysiq) == 0 .AND. iflag_phys /= 0) THEN
147               ! calcul des tendances physiques:               ! calcul des tendances physiques:
# Line 177  contains Line 178  contains
178               ucov=ucov + dudis               ucov=ucov + dudis
179               vcov=vcov + dvdis               vcov=vcov + dvdis
180    
181               ! On rajoute la tendance due a la transform. Ec -> E               ! On rajoute la tendance due à la transformation Ec -> E
182               ! therm. cree lors de la dissipation               ! thermique créée lors de la dissipation
183               call covcont(llm, ucov, vcov, ucont, vcont)               call covcont(llm, ucov, vcov, ucont, vcont)
184               call enercin(vcov, ucov, vcont, ucont, ecin)               call enercin(vcov, ucov, vcont, ucont, ecin)
185               dtetaecdt= (ecin0 - ecin) / pk               dtetaecdt= (ecin0 - ecin) / pk
186               dtetadis=dtetadis + dtetaecdt               dtetadis=dtetadis + dtetaecdt
187               teta=teta + dtetadis               teta=teta + dtetadis
188    
189               ! Calcul de la valeur moyenne, unique de h aux poles .....               ! Calcul de la valeur moyenne unique de h aux pôles
190               DO l = 1, llm               DO l = 1, llm
191                  DO ij = 1, iim                  DO ij = 1, iim
192                     tppn(ij) = aire(ij) * teta(ij, l)                     tppn(ij) = aire(ij) * teta(ij, l)
Legend:
Removed from v.27  
changed lines
  Added in v.28
  ViewVC Help
Powered by ViewVC 1.1.21