Changeset 3051 for palm/trunk/SOURCE/land_surface_model_mod.f90

Timestamp:

May 30, 2018 5:43:55 PM (6 years ago)

Author:

suehring

Message:

Speed-up NetCDF input; Revise NetCDF-input routines and remove input via io-blocks; Temporarily revoke renaming of input variables in dynamic driver; More detailed error messages created; Bugfix in mapping 3D buildings; Bugfix in land-surface model at pavement surfaces; Bugfix in initialization with inifor

File:

• palm/trunk/SOURCE/land_surface_model_mod.f90 (modified) (4 diffs)

palm/trunk/SOURCE/land_surface_model_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Bugfix in surface-element loops for pavement surfaces 
+! 
+! 3045 2018-05-28 07:55:41Z Giersch
 ! Error messages revised
 ! 
@@ -3788 +3791 @@
                    surf_lsm_v(l)%emissivity(ind_pav_green,m)  =                &
                                               pavement_pars(ind_p_emis,st)    
-
-                   DO  k = nzb_soil, surf_lsm_h%nzt_pavement(m)
+                                              
+                   DO  k = nzb_soil, surf_lsm_v(l)%nzt_pavement(m)
                       surf_lsm_v(l)%lambda_h_def(k,m)    =                     &
                                     pavement_subsurface_pars_1(k,pavement_type)                       
@@ -4410 +4413 @@
                 DO  k = nzb_soil, nzt_soil
                    IF ( surf_lsm_v(l)%pavement_surface(m)  .AND.               &
-                        k <= surf_lsm_h%nzt_pavement(m) )  THEN
+                        k <= surf_lsm_v(l)%nzt_pavement(m) )  THEN
                       surf_lsm_v(l)%root_fr(k,m) = 0.0_wp
                    ELSE
@@ -4692 +4695 @@
           ALLOCATE ( surf_lsm_v(l)%water_surface(1:surf_lsm_v(l)%ns)       )
 
-          surf_lsm_v(l)%water_surface     = .FALSE.
-          surf_lsm_v(l)%pavement_surface  = .FALSE.
-          surf_lsm_v(l)%vegetation_surface   = .FALSE. 
+          surf_lsm_v(l)%water_surface       = .FALSE.
+          surf_lsm_v(l)%pavement_surface    = .FALSE.
+          surf_lsm_v(l)%vegetation_surface  = .FALSE.