Changeset 19 for palm/trunk/SOURCE/diffusion_e.f90

Timestamp:

Feb 23, 2007 4:53:48 AM (18 years ago)

Author:

raasch

Message:

preliminary version of modified boundary conditions at top

File:

• palm/trunk/SOURCE/diffusion_e.f90 (modified) (8 diffs)

palm/trunk/SOURCE/diffusion_e.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
-! 
+! Calculation extended for gridpoint nzt
 !
 ! Former revisions:
@@ -56 +56 @@
        REAL, DIMENSION(:,:), POINTER   ::  rif
        REAL, DIMENSION(:,:,:), POINTER ::  e, km, theta
-       REAL, DIMENSION(nzb+1:nzt-1,nys:nyn) ::  dissipation, l, ll
+       REAL, DIMENSION(nzb+1:nzt,nys:nyn) ::  dissipation, l, ll
  
 
@@ -72 +72 @@
              ENDIF
 
-             DO  k = nzb_s_inner(j,i)+1, nzt-1
+             DO  k = nzb_s_inner(j,i)+1, nzt
 !
 !--             Calculate the mixing length (for dissipation)
@@ -102 +102 @@
 !--       Calculate the tendency terms
           DO  j = nys, nyn
-             DO  k = nzb_s_inner(j,i)+1, nzt-1
+             DO  k = nzb_s_inner(j,i)+1, nzt
 
                  dissipation(k,j) = ( 0.19 + 0.74 * l(k,j) / ll(k,j) ) * &
@@ -130 +130 @@
           IF ( use_sgs_for_particles )  THEN
              DO  j = nys, nyn
-                DO  k = nzb_s_inner(j,i)+1, nzt-1
+                DO  k = nzb_s_inner(j,i)+1, nzt
                    diss(k,j,i) = dissipation(k,j)
                 ENDDO
@@ -171 +171 @@
        REAL, DIMENSION(:,:), POINTER   ::  rif
        REAL, DIMENSION(:,:,:), POINTER ::  e, km, theta
-       REAL, DIMENSION(nzb+1:nzt-1)    ::  dissipation, l, ll
+       REAL, DIMENSION(nzb+1:nzt)    ::  dissipation, l, ll
 
 
@@ -187 +187 @@
 !
 !--    Calculate the mixing length (for dissipation)
-       DO  k = nzb_s_inner(j,i)+1, nzt-1
+       DO  k = nzb_s_inner(j,i)+1, nzt
           dpt_dz = ( theta(k+1,j,i) - theta(k-1,j,i) ) * dd2zu(k)
           IF ( dpt_dz > 0.0 ) THEN
@@ -234 +234 @@
 !--    Store dissipation if needed for calculating the sgs particle velocities
        IF ( use_sgs_for_particles )  THEN
-          DO  k = nzb_s_inner(j,i)+1, nzt-1
+          DO  k = nzb_s_inner(j,i)+1, nzt
              diss(k,j,i) = dissipation(k)
           ENDDO