Ignore:
Timestamp:
Jul 26, 2018 5:12:48 PM (6 years ago)
Author:
suehring
Message:

Major bugfix in calculation of ol and ts at building roofs; bugfix in restart data for surface elements; in 2D data output, mask latent heat flux at urban-type surfaces

File:
1 edited

Legend:

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

    r3152 r3176  
    2828! -----------------
    2929! $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
    3034! Initialize q_surface, which might be used in surface_layer_fluxes
    3135!
     
    74387442                           surf_usm_h%tt_surface_green_m(m)
    74397443!
    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.
    74417447           surf_usm_h%pt_surface(m) = ( surf_usm_h%frac(ind_veg_wall,m) * t_surf_h_p(m)   &
    74427448                               + surf_usm_h%frac(ind_wat_win,m) * t_surf_window_h_p(m)   &
    74437449                               + surf_usm_h%frac(ind_pav_green,m) * t_surf_green_h_p(m) )  &
    74447450                               / exn(k)
     7451                               
     7452           IF ( humidity )  surf_usm_h%vpt_surface(m) =                        &
     7453                                                   surf_usm_h%pt_surface(m)
    74457454                       
    74467455!--        calculate true tendency
     
    76377646                                surf_usm_v(l)%tt_surface_green_m(m)
    76387647!
    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.     
    76407651              surf_usm_v(l)%pt_surface(m) =  ( surf_usm_v(l)%frac(ind_veg_wall,m) * t_surf_v_p(l)%t(m)  &
    76417652                                      + surf_usm_v(l)%frac(ind_wat_win,m) * t_surf_window_v_p(l)%t(m)  &
    76427653                                      + surf_usm_v(l)%frac(ind_pav_green,m) * t_surf_green_v_p(l)%t(m) ) &
    76437654                                      / exn(k)
     7655                                     
     7656              IF ( humidity )  surf_usm_v(l)%vpt_surface(m) =                  &
     7657                                                    surf_usm_v(l)%pt_surface(m)
    76447658           
    76457659!--           calculate true tendency
Note: See TracChangeset for help on using the changeset viewer.