Changeset 4644 for palm/trunk/SOURCE

Timestamp:

Aug 22, 2020 3:05:46 PM (4 years ago)

Author:

pavelkrc

Message:

Bugfix: remove invalid reordering of vertical surfaces in RTM

File:

• palm/trunk/SOURCE/radiation_model_mod.f90 (modified) (4 diffs)

palm/trunk/SOURCE/radiation_model_mod.f90

@@ -28 +28 @@
 ! -----------------
 ! $Id$
+! Bugfix: remove invalid reordering of vertical surfaces in RTM
+!
+! 4643 2020-08-17 09:19:42Z pavelkrc
 ! - Bugfix: typo in index of surf%* arrays
 !
@@ -5831 +5834 @@
      IMPLICIT NONE
 
-     INTEGER(iwp)                      ::  i, j, k, kk, d, refstep, m, mm, l, ll
+     INTEGER(iwp)                      ::  i, j, k, kk, d, refstep, m, mm, l
      INTEGER(iwp)                      ::  isurf, isurfsrc, isvf, icsf, ipcgb
      INTEGER(iwp)                      ::  imrt, imrtf
      INTEGER(iwp)                      ::  isd                !< solar direction number
      INTEGER(iwp)                      ::  pc_box_dimshift    !< transform for best accuracy
-     INTEGER(iwp), DIMENSION(0:3)      ::  reorder = (/ 1, 0, 3, 2 /)
 
      REAL(wp), DIMENSION(3,3)          ::  mrot               !< grid rotation matrix (zyx)
@@ -5962 +5964 @@
      pabs_surf_lwdifl = 0.0_wp
      pabs_pc_lwdifl   = 0.0_wp
-
+!
 !--  Set up thermal radiation from surfaces
-!--  emiss_surf is defined only for surfaces for which energy balance is calculated
-!--  Workaround: reorder surface data type back on 1D array including all surfaces,
-!--  which implies to reorder horizontal and vertical surfaces
-!
 !--  Horizontal walls
      mm = 1
@@ -6004 +6002 @@
      DO  i = nxl, nxr
         DO  j = nys, nyn
-           DO  ll = 0, 3
-              l = reorder(ll)
+           DO  l = 0, 3
 !
 !--           urban