Changeset 3261 for palm/trunk/SOURCE

Timestamp:

Sep 19, 2018 12:32:39 PM (6 years ago)

Author:

suehring

Message:

Bugfix in raytrace_2d calls

File:

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

palm/trunk/SOURCE/radiation_model_mod.f90

@@ -28 +28 @@
 ! -----------------
 ! $Id$
+! Bugfix in raytrace_2d calls
+! 
+! 3248 2018-09-14 09:42:06Z sward
 ! Minor formating changes
 ! 
@@ -5856 +5859 @@
                ENDIF
                yxdir = (/ COS(azmid), SIN(azmid) /)
-               CALL raytrace_2d(ta, yxdir, zdirs,     &
-                                   -999, -999._wp, vffrac, .FALSE., .FALSE., &
-                                   win_lad, horizon, ztransp)
+               CALL raytrace_2d(ta, yxdir, zdirs,                              &
+                               surfstart(myid) + isurflt, facearea(td),        &
+                               vffrac, .TRUE., .FALSE., win_lad, horizon,       &
+                               ztransp)
+               
 
                azen = pi/2 - ATAN(horizon)
@@ -5984 +5989 @@
               azmid = az0 + (REAL(iaz, wp) - .5_wp) * azs
               yxdir = (/ COS(azmid), SIN(azmid) /)
-              CALL raytrace_2d(ta, yxdir, zdirs,                               &
-                               surfstart(myid) + isurflt, facearea(td),        &
-                               vffrac, .TRUE., .TRUE., win_lad, horizon,       &
-                               ztransp)
+              CALL raytrace_2d(ta, yxdir, zdirs,     &
+                                   -999, -999._wp, vffrac, .FALSE., .TRUE., &
+                                   win_lad, horizon, ztransp)
 
               !--Save direct solar transparency
@@ -6691 +6695 @@
          !--Assert that we have space allocated for CSFs
          !--
-         maxboxes = (ntrack + MAX(origin(1) - nzub, nzpt - origin(1))) * SIZE(zdirs, 1)
+         maxboxes = (ntrack + MAX(origin(1) - nzub, nzut - origin(1))) * SIZE(zdirs, 1)
 !          IF ( ncsfl + maxboxes > ncsfla )  THEN
 ! !--         use this code for growing by fixed exponential increments (equivalent to case where ncsfl always increases by 1)