Changeset 3176 for palm/trunk/SOURCE/urban_surface_mod.f90
- Timestamp:
- Jul 26, 2018 5:12:48 PM (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/urban_surface_mod.f90
r3152 r3176 28 28 ! ----------------- 29 29 ! $Id$ 30 ! Bugfix, update virtual potential surface temparture, else heat fluxes on 31 ! roofs might become unphysical 32 ! 33 ! 3152 2018-07-19 13:26:52Z suehring 30 34 ! Initialize q_surface, which might be used in surface_layer_fluxes 31 35 ! … … 7438 7442 surf_usm_h%tt_surface_green_m(m) 7439 7443 ! 7440 !-- Store surface temperature 7444 !-- Store surface temperature on pt_surface. Further, in case humidity is used 7445 !-- store also vpt_surface, which is, due to the lack of moisture on roofs simply 7446 !-- assumed to be the surface temperature. 7441 7447 surf_usm_h%pt_surface(m) = ( surf_usm_h%frac(ind_veg_wall,m) * t_surf_h_p(m) & 7442 7448 + surf_usm_h%frac(ind_wat_win,m) * t_surf_window_h_p(m) & 7443 7449 + surf_usm_h%frac(ind_pav_green,m) * t_surf_green_h_p(m) ) & 7444 7450 / exn(k) 7451 7452 IF ( humidity ) surf_usm_h%vpt_surface(m) = & 7453 surf_usm_h%pt_surface(m) 7445 7454 7446 7455 !-- calculate true tendency … … 7637 7646 surf_usm_v(l)%tt_surface_green_m(m) 7638 7647 ! 7639 !-- Store surface temperature 7648 !-- Store surface temperature. Further, in case humidity is used 7649 !-- store also vpt_surface, which is, due to the lack of moisture on roofs simply 7650 !-- assumed to be the surface temperature. 7640 7651 surf_usm_v(l)%pt_surface(m) = ( surf_usm_v(l)%frac(ind_veg_wall,m) * t_surf_v_p(l)%t(m) & 7641 7652 + surf_usm_v(l)%frac(ind_wat_win,m) * t_surf_window_v_p(l)%t(m) & 7642 7653 + surf_usm_v(l)%frac(ind_pav_green,m) * t_surf_green_v_p(l)%t(m) ) & 7643 7654 / exn(k) 7655 7656 IF ( humidity ) surf_usm_v(l)%vpt_surface(m) = & 7657 surf_usm_v(l)%pt_surface(m) 7644 7658 7645 7659 !-- calculate true tendency
Note: See TracChangeset
for help on using the changeset viewer.