Changeset 1845 for palm/trunk/SOURCE/pres.f90

Timestamp:

Apr 8, 2016 8:29:13 AM (8 years ago)

Author:

raasch

Message:

nzb_2d replaced by nzb_..._inner, Kessler precipitation stored on surface grid point

File:

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

palm/trunk/SOURCE/pres.f90

@@ -19 +19 @@
 ! Current revisions:
 ! ------------------
-! 
+! nzb_2d replace by nzb_u|v_inner
 ! 
 ! Former revisions:
@@ -119 +119 @@
         ONLY:  nbgp, ngp_2dh_outer, nx, nxl, nxlg, nxl_mg, nxr, nxrg, nxr_mg,  &
                ny, nys, nysg, nys_mg, nyn, nyng, nyn_mg, nzb, nzb_s_inner,     &
-               nzb_u_inner, nzb_v_inner, nzb_w_inner, nzb_2d, nzt, nzt_mg,     &
+               nzb_u_inner, nzb_v_inner, nzb_w_inner, nzt, nzt_mg,             &
                rflags_s_inner
 
@@ -212 +212 @@
 !
 !--       Sum up the volume flow through the south/north boundary
-          DO  k = nzb_2d(j,i)+1, nzt
+          DO  k = nzb_u_inner(j,i)+1, nzt
              volume_flow_l(1) = volume_flow_l(1) + u(k,j,i) * dzw(k)
           ENDDO
@@ -228 +228 @@
 
        DO  j = nysg, nyng
-          DO  k = nzb_2d(j,i)+1, nzt
+          DO  k = nzb_u_inner(j,i)+1, nzt
              u(k,j,i) = u(k,j,i) + volume_flow_offset(1)
           ENDDO
@@ -251 +251 @@
 !
 !--       Sum up the volume flow through the south/north boundary
-          DO  k = nzb_2d(j,i)+1, nzt
+          DO  k = nzb_v_inner(j,i)+1, nzt
              volume_flow_l(2) = volume_flow_l(2) + v(k,j,i) * dzw(k)
           ENDDO
@@ -648 +648 @@
        DO  j = nys, nyn
           !$OMP CRITICAL
-          DO  k = nzb_2d(j,nx) + 1, nzt
+          DO  k = nzb_u_inner(j,nx) + 1, nzt
              volume_flow_l(1) = volume_flow_l(1) + u(k,j,nx) * dzw(k)
           ENDDO
@@ -664 +664 @@
        DO  i = nxl, nxr
           !$OMP CRITICAL
-          DO  k = nzb_2d(ny,i) + 1, nzt
+          DO  k = nzb_v_inner(ny,i) + 1, nzt
              volume_flow_l(2) = volume_flow_l(2) + v(k,ny,i) * dzw(k)
            ENDDO
@@ -695 +695 @@
                 u(k,j,i) = u(k,j,i) + volume_flow_offset(1)
              ENDDO
-             DO k = nzb_v_inner(j,i) + 1, nzt
+             DO  k = nzb_v_inner(j,i) + 1, nzt
                 v(k,j,i) = v(k,j,i) + volume_flow_offset(2)
              ENDDO