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

Timestamp:

Mar 7, 2007 8:38:00 AM (17 years ago)

Author:

raasch

Message:

preliminary version, several changes to be explained later

File:

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

palm/trunk/SOURCE/diffusion_v.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
-! 
+! wall functions now include diabatic conditions, call of routine wall_fluxes
 !
 ! Former revisions:
@@ -55 +55 @@
 
        INTEGER ::  i, j, k
-       REAL    ::  kmxm_x, kmxm_y, kmxp_x, kmxp_y, kmzm, kmzp, vsus
+       REAL    ::  kmxm_x, kmxm_y, kmxp_x, kmxp_y, kmzm, kmzp
        REAL    ::  ddzu(1:nzt+1), ddzw(1:nzt+1), km_damp_x(nxl-1:nxr+1)
        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
@@ -100 +101 @@
 !--          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)
-                   vsus   = kappa * v(k,j,i) / LOG( 0.5 * dx / z0(j,i))
-                   vsus   = -vsus * ABS( vsus )
                    kmxp_x = 0.25 * &
                             ( km(k,j,i)+km(k,j,i+1)+km(k,j-1,i)+km(k,j-1,i+1) )
@@ -132 +137 @@
                                 + kmxm_y * ( u(k,j,i) - u(k,j-1,i) ) * ddy     &
                                                   )                            &
-                                     + wall_v(j,i) * vsus                      &
+                                     + wall_v(j,i) * vsus(k)                   &
                                    ) * ddx
                 ENDDO
@@ -204 +209 @@
 
        INTEGER ::  i, j, k
-       REAL    ::  kmxm_x, kmxm_y, kmxp_x, kmxp_y, kmzm, kmzp, vsus
+       REAL    ::  kmxm_x, kmxm_y, kmxp_x, kmxp_y, kmzm, kmzp
        REAL    ::  ddzu(1:nzt+1), ddzw(1:nzt+1), km_damp_x(nxl-1:nxr+1)
        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
@@ -245 +251 @@
 !--    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)
-             vsus   = kappa * v(k,j,i) / LOG( 0.5 * dx / z0(j,i))
-             vsus   = -vsus * ABS( vsus )
              kmxp_x = 0.25 * &
                       ( km(k,j,i)+km(k,j,i+1)+km(k,j-1,i)+km(k,j-1,i+1) )
@@ -277 +287 @@
                                 + kmxm_y * ( u(k,j,i) - u(k,j-1,i) ) * ddy     &
                                                   )                            &
-                                     + wall_v(j,i) * vsus                      &
+                                     + wall_v(j,i) * vsus(k)                   &
                                    ) * ddx
           ENDDO