/[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 91 by guez, Wed Mar 26 17:18:58 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, pkf)  
   
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
177            IF (time > 1.) time = time - 1.            IF (time > 1.) time = time - 1.
178    
179            CALL calfis(rdayvrai, time, ucov, vcov, teta, q, ps, pk, phis, phi, &            CALL calfis(rdayvrai, time, ucov, vcov, teta, q, pk, phis, phi, w, &
180                 w, dufi, dvfi, dtetafi, dqfi, lafin = itau + 1 == itaufin)                 dufi, dvfi, dtetafi, dqfi, lafin = itau + 1 == itaufin)
181    
182            ! Ajout des tendances physiques:            ! Ajout des tendances physiques:
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.91  
changed lines
  Added in v.96
  ViewVC Help
Powered by ViewVC 1.1.21