Changeset 2287 for palm/trunk

Timestamp:

Jun 15, 2017 4:46:30 PM (7 years ago)

Author:

suehring

Message:

Bugfix in determination topography-top index

File:

• palm/trunk/SOURCE/urban_surface_mod.f90 (modified) (2 diffs)

palm/trunk/SOURCE/urban_surface_mod.f90

@@ -26 +26 @@
 ! -----------------
 ! $Id$
+! Bugfix in determination topography-top index
+! 
+! 2269 2017-06-09 11:57:32Z suehring
 ! Enable restart runs with different number of PEs
 ! Bugfixes nopointer branch
@@ -3787 +3790 @@
                 i = pcbl(ix, ipcgb)
                 k = pcbl(iz, ipcgb)
-                kk = k - MAXLOC(                                               &
+!
+!--             Following expression equals former kk = k - nzb_s_inner(j,i) 
+                kk = k - ( MAXLOC(                                             &
                                 MERGE( 1, 0,                                   &
                                        BTEST( wall_flags_0(:,j,i), 12 )        &
                                      ), DIM = 1                                &
-                               ) - 1 ! kk = k - nzb_s_inner(j,i)  !- lad arrays are defined flat
+                               ) - 1                                           &
+                         )  !- lad arrays are defined flat
                 pc_heating_rate(kk, j, i) = (pcbinsw(ipcgb) + pcbinlw(ipcgb)) &
                     * pchf_prep(k) * pt(k, j, i) !-- = dT/dt