Changeset 1845 for palm/trunk/SOURCE/init_3d_model.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/init_3d_model.f90 (modified) (7 diffs)

palm/trunk/SOURCE/init_3d_model.f90

@@ -19 +19 @@
 ! Current revisions:
 ! ------------------
-! 
+! nzb_2d replaced by nzb_u|v_inner
 ! 
 ! Former revisions:
@@ -1488 +1488 @@
           IF ( nxr == nx )  THEN
              DO  j = nys, nyn
-                DO  k = nzb_2d(j,nx)+1, nzt
+                DO  k = nzb_u_inner(j,nx)+1, nzt
                    volume_flow_initial_l(1) = volume_flow_initial_l(1) +       &
                                               u_init(k) * dzw(k)
@@ -1498 +1498 @@
           IF ( nyn == ny )  THEN
              DO  i = nxl, nxr
-                DO  k = nzb_2d(ny,i)+1, nzt  
+                DO  k = nzb_v_inner(ny,i)+1, nzt
                    volume_flow_initial_l(2) = volume_flow_initial_l(2) + &
                                               v_init(k) * dzw(k)
@@ -1524 +1524 @@
           IF ( nxr == nx )  THEN
              DO  j = nys, nyn
-                DO  k = nzb_2d(j,nx)+1, nzt
+                DO  k = nzb_u_inner(j,nx)+1, nzt
                    volume_flow_initial_l(1) = volume_flow_initial_l(1) +       &
                                               hom_sum(k,1,0) * dzw(k)
@@ -1534 +1534 @@
           IF ( nyn == ny )  THEN
              DO  i = nxl, nxr
-                DO  k = nzb_2d(ny,i)+1, nzt  
+                DO  k = nzb_v_inner(ny,i)+1, nzt
                    volume_flow_initial_l(2) = volume_flow_initial_l(2) +       &
                                               hom_sum(k,2,0) * dzw(k)
@@ -1560 +1560 @@
           IF ( nxr == nx )  THEN
              DO  j = nys, nyn
-                DO  k = nzb_2d(j,nx)+1, nzt
+                DO  k = nzb_u_inner(j,nx)+1, nzt
                    volume_flow_initial_l(1) = volume_flow_initial_l(1) + &
                                               u(k,j,nx) * dzw(k)
@@ -1570 +1570 @@
           IF ( nyn == ny )  THEN
              DO  i = nxl, nxr
-                DO  k = nzb_2d(ny,i)+1, nzt  
+                DO  k = nzb_v_inner(ny,i)+1, nzt
                    volume_flow_initial_l(2) = volume_flow_initial_l(2) +       &
                                               v(k,ny,i) * dzw(k)