HTCM example setup

This page is part of the Human Thermal Comfort Module (HTCM) documentation. It contains a documentation about a test setup for human thermal comfort.
For an overview of all HTCM-related pages, see the HTCM main page.

HTCM can easily be tested by making small modifications to ​the generic "test_urban" setup. The modifications required are indicated below.

HTCM can be enabled by adding the indicated lines to the test_urban_p3d file located in the INPUT folder.

!-------------------------------------------------------------------------------
!-- INITIALIZATION PARAMETER NAMELIST
!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/initialization_parameters
!-------------------------------------------------------------------------------
 &initialization_parameters

!
!-- Grid
!------------------------------------------------------------------------------
    nx                   = 19,
    ny                   = 19,
    nz                   = 60,

    dx                   = 2.0,
    dy                   = 2.0,
    dz                   = 2.0,

!
!-- Numerics
!------------------------------------------------------------------------------
    fft_method           = 'fftw',

!
!-- Mode
!------------------------------------------------------------------------------     
    humidity             = .T.,

!
!-- Initialization
!------------------------------------------------------------------------------
    initializing_actions = 'set_constant_profiles', !'inifor',
    
    pt_surface           = 293.15,
    ug_surface           = 1.0,
    vg_surface           = 0.0,
    
    day_of_year_init     = 202, 
    time_utc_init        = 43200.0,

!
!-- Wall/soil spinup
!------------------------------------------------------------------------------    
    spinup_time         = 86400.0,
    spinup_pt_mean      = 293.15,          
    spinup_pt_amplitude = 15.0,           
    dt_spinup           = 120.0,
    data_output_during_spinup = .T.,

!
!-- Topography
!------------------------------------------------------------------------------    
    topography           = 'read_from_file',

!
!-- Physics
!------------------------------------------------------------------------------
    longitude            = 13.4,
    latitude             = 52.5,
      
 / ! end of initialization_parameters namelist


!-------------------------------------------------------------------------------
!-- RUNTIME PARAMETER NAMELIST
!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/d3par
!-------------------------------------------------------------------------------
 &runtime_parameters

!
!-- Run steering
!------------------------------------------------------------------------------
    end_time              = 86400.0,
    create_disturbances   = .T.,

!
!-- Run-control/timeseries output settings
!------------------------------------------------------------------------------
    dt_run_control        =  60.0,
    dt_dots               =  60.0,

!
!-- Profile output settings
!------------------------------------------------------------------------------
    skip_time_dopr        =   0.0,
    dt_dopr               = 3600.0,
    averaging_interval_pr = 3600.0,
    dt_averaging_input_pr =   0.0,

    data_output_pr        = '#u',  'u*2',  'wu',  'w*u*',  'w"u"',
                            '#v',  'v*2',  'wv',  'w*v*',  'w"v"',
                            'w',   'w*2',
                            '#pt', 'pt*2', 'wpt', 'w*pt*', 'w"pt"',
                            'e', 'e*', '#km', '#l',
                            '#t_soil', '#m_soil',
!
!-- 2D/3D output settings
!------------------------------------------------------------------------------
    do3d_at_begin      = .T., 
    do2d_at_begin      = .T.,

    dt_data_output     = 3600.0, 
    averaging_interval = 3600.0,
    dt_averaging_input = 0.0,

    section_xy         = 0,

    data_output = 'u',  'u_av',
                  'v',  'v_av',
                  'w',  'w_av',
                  'pt', 'pt_av',
                  'e',  'e_av',
                  'p',  'p_av',                  
                  'm_soil', 'm_soil_av', 
                  't_soil', 't_soil_av',
                  'shf*_xy', 'shf*_xy_av',
                  'ghf*_xy', 'ghf*_xy_av',
                  'u*_xy', 'z0*_xy',
                  't*_xy',
                  'm_liq*_xy', 
                  'r_s*_xy', 'r_a*_xy',
                  'rad_net*_xy', 
                  'qsws*_xy', 'qsws*_xy_av',
                  'tsurf*_xy', 'tsurf*_xy_av',
                  'pcm_lad', 'pcm_heatrate',
                  'kc_PM10', 'kc_PM25',
                  'usm_t_surf_roof', 
                  'usm_t_wall_roof', 'usm_t_wall_north', 'usm_t_wall_east', 'usm_t_wall_south', 'usm_t_wall_west',
                  'usm_t_green_roof', 'usm_t_green_north', 'usm_t_green_east', 'usm_t_green_south', 'usm_t_green_west',
                  'usm_t_window_roof', 'usm_t_window_north', 'usm_t_window_east', 'usm_t_window_south', 'usm_t_window_west',
!-- HTCM
                  'bio_perct*_xy',     ! Perceived Temperature based on instantaneous input conditions
                  'bio_pet*_xy',       ! Physiologically Equivalent Temperature, instantaneous
                  'bio_utci*_xy',      ! Universal Thermal Climate Index, instantaneous
                  'bio_perct*_xy_av',  ! Perceived Temperature based on time-averaged input parameters
                  'bio_pet*_xy_av',    ! PET with time-averaged input
                  'bio_utci*_xy_av',   ! UTCI with time-averaged input
                  'bio_mrt',           ! The mean radiant temperature used as input for the thermal comfort calculations (3D)
                  'bio_mrt_av',        ! Mean radiant temperature, time-averaged, 3D
                  'bio_mrt_xy',        ! Mean radiant temperature, instantaneous, 2D
                  'bio_mrt_xy_av',     ! Mean radiant temperature, time-averaged, 2D

 / ! end of d3par namelist


!-------------------------------------------------------------------------------
!-- RADIATION MODEL PARAMETER NAMELIST
!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/radpar
!-------------------------------------------------------------------------------
 &radiation_parameters

    radiation_scheme =  'clear-sky',
    albedo_type = 5,
    constant_albedo = .F.,

    dt_radiation = 60.0,
   
    split_diffusion_radiation = .T.,
    surface_reflections = .T.,
    nrefsteps = 2,
    mrt_factors = .F.,

 /


!-------------------------------------------------------------------------------
!-- LAND SURFACE MODEL PARAMETER NAMELIST
!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/lsmpar
!------------------------------------------------------------------------------- 
 &land_surface_parameters 

    constant_roughness     = .T.,
    aero_resist_kray       = .T., 
    
    vegetation_type        = 2, 
    soil_type              = 3,
    conserve_water_content = .T.,
    
    root_fraction          = 0.49, 0.26, 0.25, 0.0, 0.0, 0.0, 0.0, 0.0,
    soil_temperature       = 293.5, 293.6, 293.1 293.1, 293.1, 293.1 293.1, 293.1, 
    deep_soil_temperature  = 293.1,
    
 / ! end of lsm_par namelist


!-------------------------------------------------------------------------------
!-- URBAN SURFACE MODEL PARAMETER NAMELIST
!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/usmpar
!-------------------------------------------------------------------------------
 &urban_surface_parameters
  
    usm_material_model     = .T.,
    
 / ! end of urban_surface_par namelist
 

!-------------------------------------------------------------------------------
!-- PLANT CANOPY MODEL PARAMETER NAMELIST
!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/canpar
!-------------------------------------------------------------------------------
 &plant_canopy_parameters
 
    canopy_mode       = 'read_from_file_3d',
    canopy_drag_coeff = 0.3,
      
 / ! end of canopy_par namelist
 

!-------------------------------------------------------------------------------
!-- CHEMISTRY MODEL PARAMETER NAMELIST
!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/chempar
!-------------------------------------------------------------------------------
 &chemistry_parameters

    chem_gasphase_on           = .F.,
    
    surface_csflux_name        = 'PM10','PM25',
    surface_csflux             = 0.75, 0.375,
    cs_name                    = 'PM10','PM25',
    cs_surface                 = 2.0, 1.0,
    bc_cs_b                    = 'neumann',
    bc_cs_t                    = 'neumann',
    emiss_factor_main = 1.667,
    emiss_factor_side = 0.334,

    main_street_id = 12,
    side_street_id =  8,
    max_street_id = 19,

 / ! end of chemistry_par namelist

!-------------------------------------------------------------------------------
!-- BIOMETEOROLOGY MODULE PARAMETER NAMELIST
!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/biomet
!-------------------------------------------------------------------------------
 &biometeorology_parameters
    thermal_comfort = .T.,       ! This parameter is optional: It will automatically be set if omitted but thermal comfort output specified in data_output.
 / ! end of biometeorology_parameters namelist
 
!-------------------------------------------------------------------------------

return to HTCM overview