/[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 95 by guez, Wed Apr 2 12:59:54 2014 UTC revision 96 by guez, Fri Apr 4 11:30:34 2014 UTC
# Line 165  contains Line 165  contains
165                 finvmaold, dtvr, leapf=.false.)                 finvmaold, dtvr, leapf=.false.)
166         end if         end if
167    
168           forall (l = 1: llm + 1) p3d(:, :, l) = ap(l) + bp(l) * ps
169           CALL exner_hyb(ps, p3d, pks, pk, pkf)
170    
171         IF (MOD(itau + 1, iphysiq) == 0 .AND. iflag_phys /= 0) THEN         IF (MOD(itau + 1, iphysiq) == 0 .AND. iflag_phys /= 0) THEN
172            ! Calcul des tendances physiques:            ! Calcul des tendances physiques:
173    
           forall (l = 1: llm + 1) p3d(:, :, l) = ap(l) + bp(l) * ps  
           CALL exner_hyb(ps, p3d, pks, pk)  
   
174            rdaym_ini = itau * dtvr / daysec            rdaym_ini = itau * dtvr / daysec
175            rdayvrai = rdaym_ini + day_ini            rdayvrai = rdaym_ini + day_ini
176            time = REAL(mod(itau, day_step)) / day_step + time_0            time = REAL(mod(itau, day_step)) / day_step + time_0
# Line 183  contains Line 183  contains
183            CALL addfi(ucov, vcov, teta, q, dufi, dvfi, dtetafi, dqfi)            CALL addfi(ucov, vcov, teta, q, dufi, dvfi, dtetafi, dqfi)
184         ENDIF         ENDIF
185    
        forall (l = 1: llm + 1) p3d(:, :, l) = ap(l) + bp(l) * ps  
        CALL exner_hyb(ps, p3d, pks, pk, pkf)  
   
186         IF (MOD(itau + 1, idissip) == 0) THEN         IF (MOD(itau + 1, idissip) == 0) THEN
187            ! Dissipation horizontale et verticale des petites \'echelles            ! Dissipation horizontale et verticale des petites \'echelles
188    
Legend:
Removed from v.95  
changed lines
  Added in v.96
  ViewVC Help
Powered by ViewVC 1.1.21