Changeset 56 for palm/trunk/SOURCE/diffusion_v.f90

Timestamp:

Mar 8, 2007 1:57:07 PM (17 years ago)

Author:

raasch

Message:

further checkin of preliminary changes

File:

• palm/trunk/SOURCE/diffusion_v.f90 (modified) (4 diffs)

palm/trunk/SOURCE/diffusion_v.f90

@@ -32 +32 @@
 !------------------------------------------------------------------------------!
 
+    USE wall_fluxes_mod
+
     PRIVATE
     PUBLIC diffusion_v
@@ -59 +61 @@
        REAL    ::  z0(nys-1:nyn+1,nxl-1:nxr+1)
        REAL    ::  tend(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1)
-       REAL, DIMENSION(nzb:nzt+1)      ::  vsus
        REAL, DIMENSION(:,:),   POINTER ::  vsws
        REAL, DIMENSION(:,:,:), POINTER ::  km, u, v, w
+       REAL, DIMENSION(nzb:nzt+1,nys:nyn+vynp,nxl:nxr) ::  vsus
+
+!
+!--    First calculate horizontal momentum flux v'u' at vertical walls,
+!--    if neccessary
+       IF ( topography /= 'flat' )  THEN
+          CALL wall_fluxes( vsus, 0.0, 1.0, 0.0, 0.0, 0, vynp, nzb_v_inner, &
+                            nzb_v_outer, wall_v )
+       ENDIF
 
        DO  i = nxl, nxr
@@ -101 +111 @@
 !--          Wall functions at the left and right walls, respectively
              IF ( wall_v(j,i) /= 0.0 )  THEN
-
-!
-!--             Calculate the horizontal momentum flux v'u'
-                CALL wall_fluxes( i, j, nzb_v_inner(j,i)+1, nzb_v_outer(j,i), &
-                                  vsus, 0.0, 1.0, 0.0, 0.0 )
 
                 DO  k = nzb_v_inner(j,i)+1, nzb_v_outer(j,i)
@@ -137 +142 @@
                                 + kmxm_y * ( u(k,j,i) - u(k,j-1,i) ) * ddy     &
                                                   )                            &
-                                     + wall_v(j,i) * vsus(k)                   &
+                                     + wall_v(j,i) * vsus(k,j,i)               &
                                    ) * ddx
                 ENDDO