Changeset 2232 for palm/trunk/SOURCE/advec_v_pw.f90

Timestamp:

May 30, 2017 5:47:52 PM (8 years ago)

Author:

suehring

Message:

Adjustments according new topography and surface-modelling concept implemented

File:

• palm/trunk/SOURCE/advec_v_pw.f90 (modified) (7 diffs)

palm/trunk/SOURCE/advec_v_pw.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! topography representation via flags
 ! 
 ! Former revisions:
@@ -96 +96 @@
 
        USE indices,                                                            &
-           ONLY:  nxl, nxr, nyn, nysv, nzb_v_inner, nzt
+           ONLY:  nxl, nxr, nyn, nysv, nzb, nzt, wall_flags_0
 
        USE kinds
@@ -115 +115 @@
        DO  i = nxl, nxr
           DO  j = nysv, nyn
-             DO  k = nzb_v_inner(j,i)+1, nzt
+             DO  k = nzb+1, nzt
                 tend(k,j,i) = tend(k,j,i) - 0.25_wp * (                        &
                          ( v(k,j,i+1) * ( u(k,j-1,i+1) + u(k,j,i+1) - gu )     &
@@ -124 +124 @@
                          - v(k-1,j,i) * ( w(k-1,j-1,i) + w(k-1,j,i) ) )        &
                                                                   * ddzw(k)    &
-                                                      )
+                                                      )                        &
+                                      * MERGE( 1.0_wp, 0.0_wp,                 &
+                                               BTEST( wall_flags_0(k,j,i), 2 ) )
              ENDDO
           ENDDO
@@ -149 +151 @@
 
        USE indices,                                                            &
-           ONLY:  nzb_v_inner, nzt
+           ONLY:  nzb, nzt, wall_flags_0
 
        USE kinds
@@ -166 +168 @@
        gu = 2.0_wp * u_gtrans
        gv = 2.0_wp * v_gtrans
-       DO  k = nzb_v_inner(j,i)+1, nzt
+       DO  k = nzb+1, nzt
           tend(k,j,i) = tend(k,j,i) - 0.25_wp * (                              &
                          ( v(k,j,i+1) * ( u(k,j-1,i+1) + u(k,j,i+1) - gu )     &
@@ -175 +177 @@
                          - v(k-1,j,i) * ( w(k-1,j-1,i) + w(k-1,j,i) ) )        &
                                                                   * ddzw(k)    &
-                                                )
+                                                )                              &
+                                      * MERGE( 1.0_wp, 0.0_wp,                 &
+                                               BTEST( wall_flags_0(k,j,i), 2 ) )
        ENDDO