/[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 161 by guez, Fri Jul 24 14:27:59 2015 UTC revision 162 by guez, Fri Jul 24 16:54:30 2015 UTC
# Line 35  contains Line 35  contains
35      use inidissip_m, only: idissip      use inidissip_m, only: idissip
36      use integrd_m, only: integrd      use integrd_m, only: integrd
37      use nr_util, only: assert      use nr_util, only: assert
     USE pressure_var, ONLY: p3d  
38      USE temps, ONLY: itau_dyn      USE temps, ONLY: itau_dyn
39      use writedynav_m, only: writedynav      use writedynav_m, only: writedynav
40      use writehist_m, only: writehist      use writehist_m, only: writehist
# Line 98  contains Line 97  contains
97      logical leapf      logical leapf
98      real dt ! time step, in s      real dt ! time step, in s
99    
100        REAL p3d(iim + 1, jjm + 1, llm+1) ! pressure at layer interfaces, in Pa
101        ! ("p3d(i, j, l)" is at longitude "rlonv(i)", latitude "rlatu(j)",
102        ! for interface "l")
103    
104      !---------------------------------------------------      !---------------------------------------------------
105    
106      print *, "Call sequence information: leapfrog"      print *, "Call sequence information: leapfrog"
# Line 167  contains Line 170  contains
170         end if         end if
171    
172         IF (MOD(itau + 1, iphysiq) == 0 .AND. iflag_phys /= 0) THEN         IF (MOD(itau + 1, iphysiq) == 0 .AND. iflag_phys /= 0) THEN
173            CALL calfis(ucov, vcov, teta, q, pk, phis, phi, w, dufi, dvfi, &            CALL calfis(ucov, vcov, teta, q, p3d, pk, phis, phi, w, dufi, dvfi, &
174                 dtetafi, dqfi, dayvrai = itau / day_step + day_ini, &                 dtetafi, dqfi, dayvrai = itau / day_step + day_ini, &
175                 time = REAL(mod(itau, day_step)) / day_step, &                 time = REAL(mod(itau, day_step)) / day_step, &
176                 lafin = itau + 1 == itaufin)                 lafin = itau + 1 == itaufin)
Legend:
Removed from v.161  
changed lines
  Added in v.162
  ViewVC Help
Powered by ViewVC 1.1.21