Ignore:
Timestamp:
Jun 15, 2017 4:46:30 PM (4 years ago)
Author:
suehring
Message:

Bugfix in determination topography-top index

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/urban_surface_mod.f90

    r2269 r2287  
    2626! -----------------
    2727! $Id$
     28! Bugfix in determination topography-top index
     29!
     30! 2269 2017-06-09 11:57:32Z suehring
    2831! Enable restart runs with different number of PEs
    2932! Bugfixes nopointer branch
     
    37873790                i = pcbl(ix, ipcgb)
    37883791                k = pcbl(iz, ipcgb)
    3789                 kk = k - MAXLOC(                                               &
     3792!
     3793!--             Following expression equals former kk = k - nzb_s_inner(j,i)
     3794                kk = k - ( MAXLOC(                                             &
    37903795                                MERGE( 1, 0,                                   &
    37913796                                       BTEST( wall_flags_0(:,j,i), 12 )        &
    37923797                                     ), DIM = 1                                &
    3793                                ) - 1 ! kk = k - nzb_s_inner(j,i)  !- lad arrays are defined flat
     3798                               ) - 1                                           &
     3799                         )  !- lad arrays are defined flat
    37943800                pc_heating_rate(kk, j, i) = (pcbinsw(ipcgb) + pcbinlw(ipcgb)) &
    37953801                    * pchf_prep(k) * pt(k, j, i) !-- = dT/dt
Note: See TracChangeset for help on using the changeset viewer.