Changeset 4678

Timestamp:

Sep 14, 2020 10:23:14 AM (4 years ago)

Author:

pavelkrc

Message:

Fix zdirs indices for full3d in raytrace_2d

File:

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

palm/trunk/SOURCE/radiation_model_mod.f90

@@ -28 +28 @@
 ! -----------------
 ! $Id$
+! Fix zdirs indices for full3d in raytrace_2d
+!
+! 4671 2020-09-09 20:27:58Z pavelkrc
 ! Radiative transfer model RTM version 4.1
 ! - Implementation of downward facing USM and LSM surfaces
@@ -8880 +8883 @@
 !
 !--               Previously mixed rays need to be checked whether they are obstructed
-                  DO  iray = lowest_mixed_ray, lowest_free_ray-1, -1
+                  DO  iray = lowest_mixed_ray, lowest_free_ray+1, -1
                      IF ( zdirs(iray) > horz_mixed_entry )  EXIT
                      IF ( target_surfl(iray) >= 0 )  CYCLE
@@ -8939 +8942 @@
 !
 !--               Previously mixed rays need to be checked whether they are obstructed
-                  DO  iray = lowest_mixed_ray, lowest_free_ray-1, -1
+                  DO  iray = lowest_mixed_ray, lowest_free_ray+1, -1
                      IF ( zdirs(iray) > MAX(horz_mixed_entry, horz_mixed_exit) )  EXIT
                      IF ( target_surfl(iray) >= 0 )  CYCLE