Ignore:
Timestamp:
Mar 4, 2020 7:20:35 PM (19 months ago)
Author:
suehring
Message:

Change order of dimension in surface arrays %frac, %emissivity and %albedo to allow for better vectorization in the radiation interactions; Set back turbulent length scale to 8 x grid spacing in the parametrized mode for the synthetic turbulence generator (was accidentally changed in last commit)

File:
1 edited

Legend:

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

    r4402 r4441  
    2121! Current revisions:
    2222! -----------------
    23 !
     23! Change order of dimension in surface array %frac to allow for better
     24! vectorization.
    2425!
    2526! Former revisions:
     
    12251226                                          ( buildings(nb)%num_facades_per_building_h +              &
    12261227                                            buildings(nb)%num_facades_per_building_v )                !< [m2] area of total facade
    1227              window_area_per_facade       = surf_usm_h%frac(ind_wat_win,m)  * facade_element_area     !< [m2] window area per facade element
     1228             window_area_per_facade       = surf_usm_h%frac(m,ind_wat_win)  * facade_element_area     !< [m2] window area per facade element
    12281229
    12291230             buildings(nb)%net_floor_area = buildings(nb)%vol_tot / ( buildings(nb)%height_storey )
     
    12751276             near_facade_temperature = surf_usm_h%pt_10cm(m)
    12761277             indoor_wall_window_temperature =                                  &
    1277                   surf_usm_h%frac(ind_veg_wall,m) * t_wall_h(nzt_wall,m)       &
    1278                 + surf_usm_h%frac(ind_wat_win,m)  * t_window_h(nzt_wall,m)
     1278                  surf_usm_h%frac(m,ind_veg_wall) * t_wall_h(nzt_wall,m)       &
     1279                + surf_usm_h%frac(m,ind_wat_win)  * t_window_h(nzt_wall,m)
    12791280!
    12801281!--          Solar thermal gains. If net_sw_in larger than sun-protection
     
    14351436                                          ( buildings(nb)%num_facades_per_building_h +              &
    14361437                                            buildings(nb)%num_facades_per_building_v )                !< [m2] area of total facade
    1437              window_area_per_facade       = surf_usm_v(l)%frac(ind_wat_win,m)  * facade_element_area  !< [m2] window area per facade element
     1438             window_area_per_facade       = surf_usm_v(l)%frac(m,ind_wat_win)  * facade_element_area  !< [m2] window area per facade element
    14381439
    14391440             buildings(nb)%net_floor_area = buildings(nb)%vol_tot / ( buildings(nb)%height_storey )
     
    14861487             near_facade_temperature = surf_usm_v(l)%pt_10cm(m)
    14871488             indoor_wall_window_temperature =                                                       &
    1488                   surf_usm_v(l)%frac(ind_veg_wall,m) * t_wall_v(l)%t(nzt_wall,m)                    &
    1489                 + surf_usm_v(l)%frac(ind_wat_win,m)  * t_window_v(l)%t(nzt_wall,m)
     1489                  surf_usm_v(l)%frac(m,ind_veg_wall) * t_wall_v(l)%t(nzt_wall,m)                    &
     1490                + surf_usm_v(l)%frac(m,ind_wat_win)  * t_window_v(l)%t(nzt_wall,m)
    14901491!
    14911492!--          Solar thermal gains. If net_sw_in larger than sun-protection
Note: See TracChangeset for help on using the changeset viewer.