Changeset 4230 for palm/trunk/SOURCE/urban_surface_mod.f90

Timestamp:

Sep 11, 2019 1:58:14 PM (5 years ago)

Author:

suehring

Message:

Several bugfixes: profile initialization of chemical species in restart runs; Runge-Kutta tendency array not initialized in chemistry model in restart runs; fixed determination of time indices for chemical emissions (introduced with commit -4227); Update chemistry profiles in offline nesting; initialize canopy resistances for greened building walls (even if green fraction is zero)

File:

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

palm/trunk/SOURCE/urban_surface_mod.f90

@@ -28 +28 @@
 ! -----------------
 ! $Id$
+! Bugfix, initialize canopy resistance. Even if no green fraction is set, 
+! r_canopy must be initialized for output purposes. 
+! 
+! 4227 2019-09-10 18:04:34Z gronemeier
 ! implement new palm_date_time_mod
 ! 
@@ -3555 +3559 @@
 !--    Map values onto horizontal elemements
        DO  m = 1, surf_usm_h%ns
-             surf_usm_h%r_canopy_min(m)     = 200.0_wp !< min_canopy_resistance
-             surf_usm_h%g_d(m)              = 0.0_wp   !< canopy_resistance_coefficient
+          surf_usm_h%r_canopy(m)     = 200.0_wp !< canopy_resistance
+          surf_usm_h%r_canopy_min(m) = 200.0_wp !< min_canopy_resistance
+          surf_usm_h%g_d(m)          = 0.0_wp   !< canopy_resistance_coefficient
        ENDDO
 !
@@ -3563 +3568 @@
        DO  l = 0, 3
           DO  m = 1, surf_usm_v(l)%ns
-                surf_usm_v(l)%r_canopy_min(m)     = 200.0_wp !< min_canopy_resistance
-                surf_usm_v(l)%g_d(m)              = 0.0_wp   !< canopy_resistance_coefficient
+             surf_usm_v(l)%r_canopy(m)     = 200.0_wp !< canopy_resistance
+             surf_usm_v(l)%r_canopy_min(m) = 200.0_wp !< min_canopy_resistance
+             surf_usm_v(l)%g_d(m)          = 0.0_wp   !< canopy_resistance_coefficient
           ENDDO
        ENDDO
@@ -7585 +7591 @@
 !--           vegetated surfaces and bare soils.
               m_liq_max = m_max_depth * ( surf_usm_h%lai(m) )
+
               surf_usm_h%c_liq(m) = MIN( 1.0_wp, ( m_liq_usm_h%var_usm_1d(m) / m_liq_max )**0.67 )
 !
@@ -7987 +7994 @@
  
 
-              IF ( surf_usm_v(l)%frac(1,m) > 0.0_wp ) THEN
+              IF ( surf_usm_v(l)%frac(ind_pav_green,m) > 0.0_wp ) THEN
 !
 !--             Adapted from LSM: