Changeset 3597 for palm/trunk/SOURCE/indoor_model_mod.f90

Timestamp:

Dec 4, 2018 8:40:18 AM (5 years ago)

Author:

maronga

Message:

revised calculation of near surface air potential temperature

File:

• palm/trunk/SOURCE/indoor_model_mod.f90 (modified) (5 diffs)

palm/trunk/SOURCE/indoor_model_mod.f90

@@ -26 +26 @@
 ! -----------------
 ! $Id$
+! Renamed t_surf_10cm to pt_10cm
+! 
+! 3593 2018-12-03 13:51:13Z kanani
 ! Replace degree symbol by degree_C
 !
@@ -318 +321 @@
 !     i_global             -->  net_sw_in                         !global radiation [W/m2]
 !     theta_e              -->  pt(k,j,i)                         !undisturbed outside temperature, 1. PALM volume, for windows
-!     theta_sup = theta_f  -->  surf_usm_h%t_surf_10cm(m)
-!                               surf_usm_v(l)%t_surf_10cm(m)   !Air temperature, facade near (10cm) air temperature from 1. Palm volume
+!     theta_sup = theta_f  -->  surf_usm_h%pt_10cm(m)
+!                               surf_usm_v(l)%pt_10cm(m)   !Air temperature, facade near (10cm) air temperature from 1. Palm volume
 !     theta_node           -->  t_wall_h(nzt_wall,m)
 !                               t_wall_v(l)%t(nzt_wall,m)         !Temperature of innermost wall layer, for opaque wall
@@ -844 +847 @@
 
     USE urban_surface_mod,                                                     &
-        ONLY:  nzt_wall, t_surf_10cm_h, t_surf_10cm_v, t_wall_h, t_wall_v,     &
-               t_window_h, t_window_v, building_type
+        ONLY:  nzt_wall, t_wall_h, t_wall_v, t_window_h, t_window_v,           &
+               building_type
 
 
@@ -1000 +1003 @@
              j = surf_usm_h%j(m)
              k = surf_usm_h%k(m)
-             near_facade_temperature = t_surf_10cm_h(m)
+             near_facade_temperature = surf_usm_h%pt_10cm(m)
              indoor_wall_window_temperature =                                  &
                   surf_usm_h%frac(ind_veg_wall,m) * t_wall_h(nzt_wall,m)       &
@@ -1192 +1195 @@
              j = surf_usm_v(l)%j(m)
              k = surf_usm_v(l)%k(m)
-             near_facade_temperature = t_surf_10cm_v(l)%t(m)
+             near_facade_temperature = surf_usm_v(l)%pt_10cm(m)
              indoor_wall_window_temperature =                                    &
                   surf_usm_v(l)%frac(ind_veg_wall,m) * t_wall_v(l)%t(nzt_wall,m) &