Changeset 2895 for palm

Timestamp:

Mar 15, 2018 10:26:12 AM (7 years ago)

Author:

hellstea

Message:

Change in the nest initialization (pmci_interp_tril_all). Bottom wall BC is no longer overwritten

File:

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

palm/trunk/SOURCE/pmc_interface_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Change in the nest initialization (pmci_interp_tril_all). Bottom wall BC is no
+! longer overwritten.
+! 
+! 2868 2018-03-09 13:25:09Z hellstea
 ! Local conditional Neumann conditions for one-way coupling removed.  
 ! 
@@ -3569 +3573 @@
        INTEGER(iwp) ::  k_wall   !<
        INTEGER(iwp) ::  k1       !<
+       INTEGER(iwp) ::  kb       !<
        INTEGER(iwp) ::  kbc      !<
        INTEGER(iwp) ::  l        !<
@@ -3627 +3632 @@
        DO  i = ib, ie
           DO  j = jb, je
-             DO  k = nzb, nzt + 1
+!
+!--          Determine the vertical index of the first node above the
+!--          topography top at grid point (j,i) in order to not overwrite 
+!--          the bottom BC.
+             kb = get_topography_top_index_ji( j, i, TRIM ( var ) ) + 1
+             DO  k = kb, nzt + 1
                 l = ic(i)
                 m = jc(j)