Changeset 3261


Ignore:
Timestamp:
Sep 19, 2018 12:32:39 PM (3 years ago)
Author:
suehring
Message:

Bugfix in raytrace_2d calls

File:
1 edited

Legend:

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

    r3248 r3261  
    2828! -----------------
    2929! $Id$
     30! Bugfix in raytrace_2d calls
     31!
     32! 3248 2018-09-14 09:42:06Z sward
    3033! Minor formating changes
    3134!
     
    58565859               ENDIF
    58575860               yxdir = (/ COS(azmid), SIN(azmid) /)
    5858                CALL raytrace_2d(ta, yxdir, zdirs,     &
    5859                                    -999, -999._wp, vffrac, .FALSE., .FALSE., &
    5860                                    win_lad, horizon, ztransp)
     5861               CALL raytrace_2d(ta, yxdir, zdirs,                              &
     5862                               surfstart(myid) + isurflt, facearea(td),        &
     5863                               vffrac, .TRUE., .FALSE., win_lad, horizon,       &
     5864                               ztransp)
     5865               
    58615866
    58625867               azen = pi/2 - ATAN(horizon)
     
    59845989              azmid = az0 + (REAL(iaz, wp) - .5_wp) * azs
    59855990              yxdir = (/ COS(azmid), SIN(azmid) /)
    5986               CALL raytrace_2d(ta, yxdir, zdirs,                               &
    5987                                surfstart(myid) + isurflt, facearea(td),        &
    5988                                vffrac, .TRUE., .TRUE., win_lad, horizon,       &
    5989                                ztransp)
     5991              CALL raytrace_2d(ta, yxdir, zdirs,     &
     5992                                   -999, -999._wp, vffrac, .FALSE., .TRUE., &
     5993                                   win_lad, horizon, ztransp)
    59905994
    59915995              !--Save direct solar transparency
     
    66916695         !--Assert that we have space allocated for CSFs
    66926696         !--
    6693          maxboxes = (ntrack + MAX(origin(1) - nzub, nzpt - origin(1))) * SIZE(zdirs, 1)
     6697         maxboxes = (ntrack + MAX(origin(1) - nzub, nzut - origin(1))) * SIZE(zdirs, 1)
    66946698!          IF ( ncsfl + maxboxes > ncsfla )  THEN
    66956699! !--         use this code for growing by fixed exponential increments (equivalent to case where ncsfl always increases by 1)
Note: See TracChangeset for help on using the changeset viewer.