Ignore:
Timestamp:
Sep 11, 2019 1:58:14 PM (22 months 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:
1 edited

Legend:

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

    r4227 r4230  
    2828! -----------------
    2929! $Id$
     30! Bugfix, initialize canopy resistance. Even if no green fraction is set,
     31! r_canopy must be initialized for output purposes.
     32!
     33! 4227 2019-09-10 18:04:34Z gronemeier
    3034! implement new palm_date_time_mod
    3135!
     
    35553559!--    Map values onto horizontal elemements
    35563560       DO  m = 1, surf_usm_h%ns
    3557              surf_usm_h%r_canopy_min(m)     = 200.0_wp !< min_canopy_resistance
    3558              surf_usm_h%g_d(m)              = 0.0_wp   !< canopy_resistance_coefficient
     3561          surf_usm_h%r_canopy(m)     = 200.0_wp !< canopy_resistance
     3562          surf_usm_h%r_canopy_min(m) = 200.0_wp !< min_canopy_resistance
     3563          surf_usm_h%g_d(m)          = 0.0_wp   !< canopy_resistance_coefficient
    35593564       ENDDO
    35603565!
     
    35633568       DO  l = 0, 3
    35643569          DO  m = 1, surf_usm_v(l)%ns
    3565                 surf_usm_v(l)%r_canopy_min(m)     = 200.0_wp !< min_canopy_resistance
    3566                 surf_usm_v(l)%g_d(m)              = 0.0_wp   !< canopy_resistance_coefficient
     3570             surf_usm_v(l)%r_canopy(m)     = 200.0_wp !< canopy_resistance
     3571             surf_usm_v(l)%r_canopy_min(m) = 200.0_wp !< min_canopy_resistance
     3572             surf_usm_v(l)%g_d(m)          = 0.0_wp   !< canopy_resistance_coefficient
    35673573          ENDDO
    35683574       ENDDO
     
    75857591!--           vegetated surfaces and bare soils.
    75867592              m_liq_max = m_max_depth * ( surf_usm_h%lai(m) )
     7593
    75877594              surf_usm_h%c_liq(m) = MIN( 1.0_wp, ( m_liq_usm_h%var_usm_1d(m) / m_liq_max )**0.67 )
    75887595!
     
    79877994 
    79887995
    7989               IF ( surf_usm_v(l)%frac(1,m) > 0.0_wp ) THEN
     7996              IF ( surf_usm_v(l)%frac(ind_pav_green,m) > 0.0_wp ) THEN
    79907997!
    79917998!--             Adapted from LSM:
Note: See TracChangeset for help on using the changeset viewer.