Changeset 106 for palm/trunk/SOURCE/production_e.f90

Timestamp:

Aug 16, 2007 2:30:26 PM (17 years ago)

Author:

raasch

Message:

preliminary update of bugfixes and extensions for non-cyclic BCs

File:

• palm/trunk/SOURCE/production_e.f90 (modified) (3 diffs)

palm/trunk/SOURCE/production_e.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
-! 
+! Bugfix: wrong sign removed from the buoyancy production term in the case
+! use_reference = .T.,
+! u_0 and v_0 are calculated for nxr+1, nyn+1 also (otherwise these values are
+! not available in case of non-cyclic boundary conditions)
 !
 ! Former revisions:
@@ -386 +389 @@
                    DO  j = nys, nyn
                       DO  k = nzb_diff_s_inner(j,i), nzt_diff
-                         tend(k,j,i) = tend(k,j,i) +                  &
+                         tend(k,j,i) = tend(k,j,i) -                  &
                                        kh(k,j,i) * g / pt_reference * &
                                        ( pt(k+1,j,i) - pt(k-1,j,i) ) * dd2zu(k)
@@ -976 +979 @@
 !--       (otherwise the timestep may be significantly reduced by large
 !--       surface winds).
+!--       Upper bounds are nxr+1 and nyn+1 because otherwise these values are
+!--       not available in case of non-cyclic boundary conditions.
 !--       WARNING: the exact analytical solution would require the determination
 !--                of the eddy diffusivity by km = u* * kappa * zp / phi_m.
           !$OMP PARALLEL DO PRIVATE( ku, kv )
-          DO  i = nxl, nxr
-             DO  j = nys, nyn
+          DO  i = nxl, nxr+1
+             DO  j = nys, nyn+1
 
                 ku = nzb_u_inner(j,i)+1