Changeset 4692


Ignore:
Timestamp:
Sep 22, 2020 5:17:52 PM (4 years ago)
Author:
maronga
Message:

minor bugfix in urban surface model, further update of output for urban_environment_salsa test case

Location:
palm/trunk
Files:
2 edited

Legend:

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

    r4687 r4692  
    2727! -----------------
    2828! $Id$
     29! Bugfix for previous revision
     30!
     31! 4687 2020-09-21 19:40:16Z maronga
    2932! Optimized code structure for treatment of inner wall and window heat flux
    3033!
     
    52005203!--    heat diffusion (innermost wall temperature as boundary conditions)
    52015204       IF ( .NOT. indoor_model ) THEN
    5202           surf%iwghf_eb(m) = - surf%lambda_h(nzt_wall-1,m) * ( t_wall%val(nzt_wall,m) &
    5203                              - t_wall%val(nzt_wall-1,m) ) * surf%ddz_wall(nzt_wall,m)
    5204           ENDIF
    5205 
    5206           DO  kw = nzb_wall+1, nzt_wall-1
    5207              wtend(kw) = ( 1.0_wp / surf%rho_c_wall(kw,m) )                             &
    5208                          * ( surf%lambda_h(kw,m) * wall_mod(kw)                         &
    5209                              * ( t_wall%val(kw+1,m) - t_wall%val(kw,m) )                &
    5210                              * surf%ddz_wall(kw+1,m)                                    &
    5211                              - surf%lambda_h(kw-1,m) * wall_mod(kw-1)                   &
    5212                              * ( t_wall%val(kw,m) - t_wall%val(kw-1,m) )                &
    5213                              * surf%ddz_wall(kw,m)                                      &
    5214                            ) * surf%ddz_wall_stag(kw,m)
    5215           ENDDO
    5216           wtend(nzt_wall) = ( 1.0_wp / surf%rho_c_wall(nzt_wall,m) )                    &
    5217                             * ( -surf%lambda_h(nzt_wall-1,m) * wall_mod(nzt_wall-1)     &
    5218                                 * ( t_wall%val(nzt_wall,m) - t_wall%val(nzt_wall-1,m) ) &
    5219                                 * surf%ddz_wall(nzt_wall,m)                            &
    5220                                 + surf%iwghf_eb(m)                                      &
    5221                               ) * surf%ddz_wall_stag(nzt_wall,m)
     5205          surf%iwghf_eb(m) = - surf%lambda_h(nzt_wall-1,m) * wall_mod(kw-1)          &
     5206                             * ( t_wall%val(nzt_wall,m) - t_wall%val(nzt_wall-1,m) )  &
     5207                             * surf%ddz_wall(nzt_wall,m)
     5208       ENDIF
     5209
     5210       DO  kw = nzb_wall+1, nzt_wall-1
     5211          wtend(kw) = ( 1.0_wp / surf%rho_c_wall(kw,m) )                         &
     5212                      * ( surf%lambda_h(kw,m) * wall_mod(kw)                     &
     5213                      * ( t_wall%val(kw+1,m) - t_wall%val(kw,m) )                &
     5214                      * surf%ddz_wall(kw+1,m)                                    &
     5215                      - surf%lambda_h(kw-1,m) * wall_mod(kw-1)                   &
     5216                      * ( t_wall%val(kw,m) - t_wall%val(kw-1,m) )                &
     5217                      * surf%ddz_wall(kw,m)                                      &
     5218                      ) * surf%ddz_wall_stag(kw,m)
     5219       ENDDO
     5220       wtend(nzt_wall) = ( 1.0_wp / surf%rho_c_wall(nzt_wall,m) )                    &
     5221                         * ( -surf%lambda_h(nzt_wall-1,m) * wall_mod(nzt_wall-1)    &
     5222                             * ( t_wall%val(nzt_wall,m) - t_wall%val(nzt_wall-1,m) ) &
     5223                              * surf%ddz_wall(nzt_wall,m) + surf%iwghf_eb(m)         &
     5224                           ) * surf%ddz_wall_stag(nzt_wall,m)
    52225225
    52235226       t_wall_p%val(nzb_wall:nzt_wall,m) = t_wall%val(nzb_wall:nzt_wall,m) + dt_3d      &
Note: See TracChangeset for help on using the changeset viewer.