/[lmdze]/trunk/dyn3d/calfis.f
ViewVC logotype

Diff of /trunk/dyn3d/calfis.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 162 by guez, Fri Jul 24 16:54:30 2015 UTC revision 252 by guez, Mon Jan 22 15:02:56 2018 UTC
# Line 52  contains Line 52  contains
52      REAL, intent(in):: q(:, :, :, :) ! (iim + 1, jjm + 1, llm, nqmx)      REAL, intent(in):: q(:, :, :, :) ! (iim + 1, jjm + 1, llm, nqmx)
53      ! mass fractions of advected fields      ! mass fractions of advected fields
54    
55      REAL, intent(in):: p3d(:, :, :) ! (iim + 1, jjm + 1, llm+1)      REAL, intent(in):: p3d(:, :, :) ! (iim + 1, jjm + 1, llm + 1)
56      ! pressure at layer interfaces, in Pa      ! pressure at layer interfaces, in Pa
57      ! ("p3d(i, j, l)" is at longitude "rlonv(i)", latitude "rlatu(j)",      ! ("p3d(i, j, l)" is at longitude "rlonv(i)", latitude "rlatu(j)",
58      ! for interface "l")      ! for interface "l")
# Line 181  contains Line 181  contains
181    
182      forall(l = 1: llm) v(:, l) = pack(zvfi(:, :, l), dyn_phy)      forall(l = 1: llm) v(:, l) = pack(zvfi(:, :, l), dyn_phy)
183    
     ! Appel de la physique :  
184      CALL physiq(lafin, dayvrai, time, paprs, play, pphi, pphis, u, v, t, qx, &      CALL physiq(lafin, dayvrai, time, paprs, play, pphi, pphis, u, v, t, qx, &
185           omega, d_u, d_v, d_t, d_qx)           omega, d_u, d_v, d_t, d_qx)
186    
# Line 219  contains Line 218  contains
218         DO j = 2, jjm         DO j = 2, jjm
219            ig0 = 1 + (j - 2) * iim            ig0 = 1 + (j - 2) * iim
220            DO i = 1, iim - 1            DO i = 1, iim - 1
221               dufi(i, j, l) = 0.5 * (d_u(ig0 + i, l) + d_u(ig0 + i+1, l)) &               dufi(i, j, l) = 0.5 * (d_u(ig0 + i, l) + d_u(ig0 + i + 1, l)) &
222                    * cu_2d(i, j)                    * cu_2d(i, j)
223            ENDDO            ENDDO
224            dufi(iim, j, l) = 0.5 * (d_u(ig0 + 1, l) + d_u(ig0 + iim, l)) &            dufi(iim, j, l) = 0.5 * (d_u(ig0 + 1, l) + d_u(ig0 + iim, l)) &
# Line 234  contains Line 233  contains
233         DO j = 2, jjm - 1         DO j = 2, jjm - 1
234            ig0 = 1 + (j - 2) * iim            ig0 = 1 + (j - 2) * iim
235            DO i = 1, iim            DO i = 1, iim
236               dvfi(i, j, l) = 0.5 * (d_v(ig0 + i, l) + d_v(ig0 + i+iim, l)) &               dvfi(i, j, l) = 0.5 * (d_v(ig0 + i, l) + d_v(ig0 + i + iim, l)) &
237                    * cv_2d(i, j)                    * cv_2d(i, j)
238            ENDDO            ENDDO
239            dvfi(iim + 1, j, l) = dvfi(1, j, l)            dvfi(iim + 1, j, l) = dvfi(1, j, l)
Legend:
Removed from v.162  
changed lines
  Added in v.252
  ViewVC Help
Powered by ViewVC 1.1.21