Ignore:
Timestamp:
Dec 4, 2018 8:40:18 AM (3 years ago)
Author:
maronga
Message:

revised calculation of near surface air potential temperature

File:
1 edited

Legend:

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

    r3593 r3597  
    2626! -----------------
    2727! $Id$
     28! Renamed t_surf_10cm to pt_10cm
     29!
     30! 3593 2018-12-03 13:51:13Z kanani
    2831! Replace degree symbol by degree_C
    2932!
     
    318321!     i_global             -->  net_sw_in                         !global radiation [W/m2]
    319322!     theta_e              -->  pt(k,j,i)                         !undisturbed outside temperature, 1. PALM volume, for windows
    320 !     theta_sup = theta_f  -->  surf_usm_h%t_surf_10cm(m)
    321 !                               surf_usm_v(l)%t_surf_10cm(m)   !Air temperature, facade near (10cm) air temperature from 1. Palm volume
     323!     theta_sup = theta_f  -->  surf_usm_h%pt_10cm(m)
     324!                               surf_usm_v(l)%pt_10cm(m)   !Air temperature, facade near (10cm) air temperature from 1. Palm volume
    322325!     theta_node           -->  t_wall_h(nzt_wall,m)
    323326!                               t_wall_v(l)%t(nzt_wall,m)         !Temperature of innermost wall layer, for opaque wall
     
    844847
    845848    USE urban_surface_mod,                                                     &
    846         ONLY:  nzt_wall, t_surf_10cm_h, t_surf_10cm_v, t_wall_h, t_wall_v,     &
    847                t_window_h, t_window_v, building_type
     849        ONLY:  nzt_wall, t_wall_h, t_wall_v, t_window_h, t_window_v,           &
     850               building_type
    848851
    849852
     
    10001003             j = surf_usm_h%j(m)
    10011004             k = surf_usm_h%k(m)
    1002              near_facade_temperature = t_surf_10cm_h(m)
     1005             near_facade_temperature = surf_usm_h%pt_10cm(m)
    10031006             indoor_wall_window_temperature =                                  &
    10041007                  surf_usm_h%frac(ind_veg_wall,m) * t_wall_h(nzt_wall,m)       &
     
    11921195             j = surf_usm_v(l)%j(m)
    11931196             k = surf_usm_v(l)%k(m)
    1194              near_facade_temperature = t_surf_10cm_v(l)%t(m)
     1197             near_facade_temperature = surf_usm_v(l)%pt_10cm(m)
    11951198             indoor_wall_window_temperature =                                    &
    11961199                  surf_usm_v(l)%frac(ind_veg_wall,m) * t_wall_v(l)%t(nzt_wall,m) &
Note: See TracChangeset for help on using the changeset viewer.