Changeset 102 for palm/trunk/SOURCE/diffusion_u.f90

Timestamp:

Jul 27, 2007 9:09:17 AM (17 years ago)

Author:

raasch

Message:

preliminary version for coupled runs

File:

• palm/trunk/SOURCE/diffusion_u.f90 (modified) (9 diffs)

palm/trunk/SOURCE/diffusion_u.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
-! 
+! Momentumflux at top (uswst) included as boundary condition
 ! 
 ! Former revisions:
@@ -54 +54 @@
 ! Call for all grid points
 !------------------------------------------------------------------------------!
-    SUBROUTINE diffusion_u( ddzu, ddzw, km, km_damp_y, tend, u, usws, v, w )
+    SUBROUTINE diffusion_u( ddzu, ddzw, km, km_damp_y, tend, u, usws, uswst, &
+                            v, w )
 
        USE control_parameters
@@ -66 +67 @@
        REAL    ::  ddzu(1:nzt+1), ddzw(1:nzt+1), km_damp_y(nys-1:nyn+1)
        REAL    ::  tend(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1)
-       REAL, DIMENSION(:,:),   POINTER ::  usws
+       REAL, DIMENSION(:,:),   POINTER ::  usws, uswst
        REAL, DIMENSION(:,:,:), POINTER ::  km, u, v, w
        REAL, DIMENSION(nzb:nzt+1,nys:nyn,nxl:nxr) ::  usvs
@@ -155 +156 @@
 !--          Compute vertical diffusion. In case of simulating a Prandtl layer,
 !--          index k starts at nzb_u_inner+2.
-             DO  k = nzb_diff_u(j,i), nzt
+             DO  k = nzb_diff_u(j,i), nzt_diff
 !
 !--             Interpolate eddy diffusivities on staggered gridpoints
@@ -195 +196 @@
                       & + ( kmzp * ( w(k,j,i)   - w(k,j,i-1)   ) * ddx       &
                       &   ) * ddzw(k)                                        &
-                      & + ( kmzp * ( u(k+1,j,i) - u(k,j,i)   ) * ddzu(k+1)   &
+                      & + ( kmzp * ( u(k+1,j,i) - u(k,j,i)     ) * ddzu(k+1) &
                       &   + usws(j,i)                                        &
+                      &   ) * ddzw(k)
+             ENDIF
+
+!
+!--          Vertical diffusion at the first gridpoint below the top boundary,
+!--          if the momentum flux at the top is prescribed by the user
+             IF ( use_top_fluxes )  THEN
+                k = nzt
+!
+!--             Interpolate eddy diffusivities on staggered gridpoints
+                kmzp = 0.25 * &
+                       ( km(k,j,i)+km(k+1,j,i)+km(k,j,i-1)+km(k+1,j,i-1) )
+                kmzm = 0.25 * &
+                       ( km(k,j,i)+km(k-1,j,i)+km(k,j,i-1)+km(k-1,j,i-1) )
+
+                tend(k,j,i) = tend(k,j,i)                                    &
+                      & - ( kmzm * ( w(k-1,j,i) - w(k-1,j,i-1) ) * ddx       &
+                      &   ) * ddzw(k)                                        &
+                      & + ( -uswst(j,i)                                      &
+                      &   - kmzm * ( u(k,j,i)   - u(k-1,j,i)   ) * ddzu(k)   &
                       &   ) * ddzw(k)
              ENDIF
@@ -210 +231 @@
 !------------------------------------------------------------------------------!
     SUBROUTINE diffusion_u_ij( i, j, ddzu, ddzw, km, km_damp_y, tend, u, usws, &
-                               v, w )
+                               uswst, v, w )
 
        USE control_parameters
@@ -223 +244 @@
        REAL    ::  tend(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1)
        REAL, DIMENSION(nzb:nzt+1)      ::  usvs
-       REAL, DIMENSION(:,:),   POINTER ::  usws
+       REAL, DIMENSION(:,:),   POINTER ::  usws, uswst
        REAL, DIMENSION(:,:,:), POINTER ::  km, u, v, w
 
@@ -303 +324 @@
 !--    Compute vertical diffusion. In case of simulating a Prandtl layer,
 !--    index k starts at nzb_u_inner+2.
-       DO  k = nzb_diff_u(j,i), nzt
+       DO  k = nzb_diff_u(j,i), nzt_diff
 !
 !--       Interpolate eddy diffusivities on staggered gridpoints
@@ -339 +360 @@
                       & + ( kmzp * ( w(k,j,i)   - w(k,j,i-1)   ) * ddx       &
                       &   ) * ddzw(k)                                        &
-                      & + ( kmzp * ( u(k+1,j,i) - u(k,j,i)   ) * ddzu(k+1)   &
+                      & + ( kmzp * ( u(k+1,j,i) - u(k,j,i)     ) * ddzu(k+1) &
                       &   + usws(j,i)                                        &
                       &   ) * ddzw(k)
        ENDIF
 
+!
+!--    Vertical diffusion at the first gridpoint below the top boundary,
+!--    if the momentum flux at the top is prescribed by the user
+       IF ( use_top_fluxes )  THEN
+          k = nzt
+!
+!--       Interpolate eddy diffusivities on staggered gridpoints
+          kmzp = 0.25 * &
+                 ( km(k,j,i)+km(k+1,j,i)+km(k,j,i-1)+km(k+1,j,i-1) )
+          kmzm = 0.25 * &
+                 ( km(k,j,i)+km(k-1,j,i)+km(k,j,i-1)+km(k-1,j,i-1) )
+
+          tend(k,j,i) = tend(k,j,i)                                          &
+                      & - ( kmzm * ( w(k-1,j,i) - w(k-1,j,i-1) ) * ddx       &
+                      &   ) * ddzw(k)                                        &
+                      & + ( -uswst(j,i)                                      &
+                      &   - kmzm * ( u(k,j,i)   - u(k-1,j,i)   ) * ddzu(k)   &
+                      &   ) * ddzw(k)
+       ENDIF
+
     END SUBROUTINE diffusion_u_ij