Changeset 4644 for palm


Ignore:
Timestamp:
Aug 22, 2020 3:05:46 PM (2 months ago)
Author:
pavelkrc
Message:

Bugfix: remove invalid reordering of vertical surfaces in RTM

File:
1 edited

Legend:

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

    r4643 r4644  
    2828! -----------------
    2929! $Id$
     30! Bugfix: remove invalid reordering of vertical surfaces in RTM
     31!
     32! 4643 2020-08-17 09:19:42Z pavelkrc
    3033! - Bugfix: typo in index of surf%* arrays
    3134!
     
    58315834     IMPLICIT NONE
    58325835
    5833      INTEGER(iwp)                      ::  i, j, k, kk, d, refstep, m, mm, l, ll
     5836     INTEGER(iwp)                      ::  i, j, k, kk, d, refstep, m, mm, l
    58345837     INTEGER(iwp)                      ::  isurf, isurfsrc, isvf, icsf, ipcgb
    58355838     INTEGER(iwp)                      ::  imrt, imrtf
    58365839     INTEGER(iwp)                      ::  isd                !< solar direction number
    58375840     INTEGER(iwp)                      ::  pc_box_dimshift    !< transform for best accuracy
    5838      INTEGER(iwp), DIMENSION(0:3)      ::  reorder = (/ 1, 0, 3, 2 /)
    58395841
    58405842     REAL(wp), DIMENSION(3,3)          ::  mrot               !< grid rotation matrix (zyx)
     
    59625964     pabs_surf_lwdifl = 0.0_wp
    59635965     pabs_pc_lwdifl   = 0.0_wp
    5964 
     5966!
    59655967!--  Set up thermal radiation from surfaces
    5966 !--  emiss_surf is defined only for surfaces for which energy balance is calculated
    5967 !--  Workaround: reorder surface data type back on 1D array including all surfaces,
    5968 !--  which implies to reorder horizontal and vertical surfaces
    5969 !
    59705968!--  Horizontal walls
    59715969     mm = 1
     
    60046002     DO  i = nxl, nxr
    60056003        DO  j = nys, nyn
    6006            DO  ll = 0, 3
    6007               l = reorder(ll)
     6004           DO  l = 0, 3
    60086005!
    60096006!--           urban
Note: See TracChangeset for help on using the changeset viewer.