/[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 94 by guez, Wed Mar 26 17:18:58 2014 UTC revision 95 by guez, Wed Apr 2 12:59:54 2014 UTC
# Line 169  contains Line 169  contains
169            ! Calcul des tendances physiques:            ! Calcul des tendances physiques:
170    
171            forall (l = 1: llm + 1) p3d(:, :, l) = ap(l) + bp(l) * ps            forall (l = 1: llm + 1) p3d(:, :, l) = ap(l) + bp(l) * ps
172            CALL exner_hyb(ps, p3d, pks, pk, pkf)            CALL exner_hyb(ps, p3d, pks, pk)
173    
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)
Legend:
Removed from v.94  
changed lines
  Added in v.95
  ViewVC Help
Powered by ViewVC 1.1.21