/[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 212 by guez, Thu Jan 12 12:31:31 2017 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 219  contains Line 219  contains
219         DO j = 2, jjm         DO j = 2, jjm
220            ig0 = 1 + (j - 2) * iim            ig0 = 1 + (j - 2) * iim
221            DO i = 1, iim - 1            DO i = 1, iim - 1
222               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)) &
223                    * cu_2d(i, j)                    * cu_2d(i, j)
224            ENDDO            ENDDO
225            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 234  contains
234         DO j = 2, jjm - 1         DO j = 2, jjm - 1
235            ig0 = 1 + (j - 2) * iim            ig0 = 1 + (j - 2) * iim
236            DO i = 1, iim            DO i = 1, iim
237               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)) &
238                    * cv_2d(i, j)                    * cv_2d(i, j)
239            ENDDO            ENDDO
240            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.212
  ViewVC Help
Powered by ViewVC 1.1.21