Index: palm/trunk/SOURCE/land_surface_model_mod.f90
===================================================================
--- palm/trunk/SOURCE/land_surface_model_mod.f90	(revision 3880)
+++ palm/trunk/SOURCE/land_surface_model_mod.f90	(revision 3881)
@@ -25,4 +25,8 @@
 ! -----------------
 ! $Id$
+! Bugfix in level 3 initialization of pavement albedo type and pavement 
+! emissivity
+! 
+! 3868 2019-04-08 11:52:36Z suehring
 ! More strict limitation of roughness length when it is in the order of the 
 ! vertical grid spacing 
@@ -4186,9 +4190,9 @@
                      pavement_pars_f%fill )                                    &
                    surf_lsm_h%albedo_type(ind_pav_green,m) =                   &
-                                              INT( pavement_pars(ind_p_at,st) )
+                                   INT( pavement_pars_f%pars_xy(ind_p_at,j,i) )
                 IF ( pavement_pars_f%pars_xy(ind_p_emis,j,i) /=                &
                      pavement_pars_f%fill )                                    &
                    surf_lsm_h%emissivity(ind_pav_green,m)  =                   &
-                                              pavement_pars(ind_p_emis,st) 
+                                   pavement_pars_f%pars_xy(ind_p_emis,j,i)
              ENDIF 
 
@@ -4246,10 +4250,10 @@
                         pavement_pars_f%fill )                                 &
                       surf_lsm_v(l)%albedo_type(ind_pav_green,m) =             &
-                                            INT( pavement_pars(ind_p_at,st) )
+                                   INT( pavement_pars_f%pars_xy(ind_p_at,j,i) )
 
                    IF ( pavement_pars_f%pars_xy(ind_p_emis,j,i) /=             &
                         pavement_pars_f%fill )                                 &
                       surf_lsm_v(l)%emissivity(ind_pav_green,m)  =             &
-                                            pavement_pars(ind_p_emis,st)  
+                                   pavement_pars_f%pars_xy(ind_p_emis,j,i)  
                 ENDIF 
              ENDDO
Index: palm/trunk/SOURCE/prognostic_equations.f90
===================================================================
--- palm/trunk/SOURCE/prognostic_equations.f90	(revision 3880)
+++ palm/trunk/SOURCE/prognostic_equations.f90	(revision 3881)
@@ -25,4 +25,7 @@
 ! -----------------
 ! $Id$
+! Bugfix in OpenMP directive
+! 
+! 3880 2019-04-08 21:43:02Z knoop
 ! Moved wtm_tendencies to module_interface_actions
 ! 
@@ -589,6 +592,5 @@
 !
 !-- Loop over all prognostic equations
-!-- b, c ang g added for SALSA
-    !$OMP PARALLEL PRIVATE (i,i_omp_start,j,k,loop_start,tn,b,c,g)
+    !$OMP PARALLEL PRIVATE (i,i_omp_start,j,k,loop_start,tn)
 
     !$  tn = omp_get_thread_num()
Index: palm/trunk/SOURCE/radiation_model_mod.f90
===================================================================
--- palm/trunk/SOURCE/radiation_model_mod.f90	(revision 3880)
+++ palm/trunk/SOURCE/radiation_model_mod.f90	(revision 3881)
@@ -28,4 +28,8 @@
 ! -----------------
 ! $Id$
+! Output of albedo and emissivity moved from USM, bugfixes in initialization
+! of albedo
+! 
+! 3861 2019-04-04 06:27:41Z maronga
 ! Bugfix: factor of 4.0 required instead of 3.0 in calculation of rad_lw_out_change_0
 ! 
@@ -1397,5 +1401,6 @@
           unit = 'W/m2'
        ELSE IF ( var(1:7) == 'rtm_svf'  .OR.  var(1:7) == 'rtm_dif'  .OR.                &
-                 var(1:9) == 'rtm_skyvf' .OR. var(1:9) == 'rtm_skyvft' )  THEN
+                 var(1:9) == 'rtm_skyvf' .OR. var(1:9) == 'rtm_skyvft'  .OR.             &
+                 var(1:12) == 'rtm_surfalb_'  .OR.  var(1:13) == 'rtm_surfemis_'  ) THEN
           IF ( .NOT.  radiation ) THEN
                 message_string = 'output of "' // TRIM( var ) // '" require'&
@@ -2443,19 +2448,19 @@
                 DO  ind_type = 0, 2
                    IF ( surf_lsm_h%albedo_type(ind_type,m) == 0 )  THEN
-                      IF ( albedo_pars_f%pars_xy(1,j,i) /= albedo_pars_f%fill )&
+                      IF ( albedo_pars_f%pars_xy(0,j,i) /= albedo_pars_f%fill )&
                          surf_lsm_h%albedo(ind_type,m) =                       &
-                                                albedo_pars_f%pars_xy(1,j,i)
+                                                albedo_pars_f%pars_xy(0,j,i)
                       IF ( albedo_pars_f%pars_xy(1,j,i) /= albedo_pars_f%fill )&
                          surf_lsm_h%aldir(ind_type,m) =                        &
                                                 albedo_pars_f%pars_xy(1,j,i)
+                      IF ( albedo_pars_f%pars_xy(1,j,i) /= albedo_pars_f%fill )&
+                         surf_lsm_h%aldif(ind_type,m) =                        &
+                                                albedo_pars_f%pars_xy(1,j,i)
                       IF ( albedo_pars_f%pars_xy(2,j,i) /= albedo_pars_f%fill )&
-                         surf_lsm_h%aldif(ind_type,m) =                        &
+                         surf_lsm_h%asdir(ind_type,m) =                        &
                                                 albedo_pars_f%pars_xy(2,j,i)
-                      IF ( albedo_pars_f%pars_xy(3,j,i) /= albedo_pars_f%fill )&
-                         surf_lsm_h%asdir(ind_type,m) =                        &
-                                                albedo_pars_f%pars_xy(3,j,i)
-                      IF ( albedo_pars_f%pars_xy(4,j,i) /= albedo_pars_f%fill )&
+                      IF ( albedo_pars_f%pars_xy(2,j,i) /= albedo_pars_f%fill )&
                          surf_lsm_h%asdif(ind_type,m) =                        &
-                                                albedo_pars_f%pars_xy(4,j,i)
+                                                albedo_pars_f%pars_xy(2,j,i)
                    ENDIF
                 ENDDO
@@ -2469,24 +2474,54 @@
                    j = surf_usm_h%j(m)
 !
-!--                Spectral albedos for wall/green/window surfaces
+!--                Broadband albedos for wall/green/window surfaces
                    DO  ind_type = 0, 2
                       IF ( surf_usm_h%albedo_type(ind_type,m) == 0 )  THEN
-                         IF ( albedo_pars_f%pars_xy(1,j,i) /= albedo_pars_f%fill )&
+                         IF ( albedo_pars_f%pars_xy(0,j,i) /= albedo_pars_f%fill )&
                             surf_usm_h%albedo(ind_type,m) =                       &
-                                                albedo_pars_f%pars_xy(1,j,i)
-                         IF ( albedo_pars_f%pars_xy(1,j,i) /= albedo_pars_f%fill )&
-                            surf_usm_h%aldir(ind_type,m) =                        &
-                                                albedo_pars_f%pars_xy(1,j,i)
-                         IF ( albedo_pars_f%pars_xy(2,j,i) /= albedo_pars_f%fill )&
-                            surf_usm_h%aldif(ind_type,m) =                        &
-                                                albedo_pars_f%pars_xy(2,j,i)
-                         IF ( albedo_pars_f%pars_xy(3,j,i) /= albedo_pars_f%fill )&
-                            surf_usm_h%asdir(ind_type,m) =                        &
-                                                albedo_pars_f%pars_xy(3,j,i)
-                         IF ( albedo_pars_f%pars_xy(4,j,i) /= albedo_pars_f%fill )&
-                            surf_usm_h%asdif(ind_type,m) =                        &
-                                                albedo_pars_f%pars_xy(4,j,i)
+                                                albedo_pars_f%pars_xy(0,j,i)
                       ENDIF
                    ENDDO
+!
+!--                Spectral albedos especially for building wall surfaces
+                   IF ( albedo_pars_f%pars_xy(1,j,i) /= albedo_pars_f%fill )  THEN
+                      surf_usm_h%aldir(ind_veg_wall,m) =                       &
+                                                albedo_pars_f%pars_xy(1,j,i)
+                      surf_usm_h%aldif(ind_veg_wall,m) =                       &
+                                                albedo_pars_f%pars_xy(1,j,i)
+                   ENDIF
+                   IF ( albedo_pars_f%pars_xy(2,j,i) /= albedo_pars_f%fill )  THEN
+                      surf_usm_h%asdir(ind_veg_wall,m) =                       &
+                                                albedo_pars_f%pars_xy(2,j,i)
+                      surf_usm_h%asdif(ind_veg_wall,m) =                       &
+                                                albedo_pars_f%pars_xy(2,j,i)
+                   ENDIF
+!
+!--                Spectral albedos especially for building green surfaces
+                   IF ( albedo_pars_f%pars_xy(3,j,i) /= albedo_pars_f%fill )  THEN
+                      surf_usm_h%aldir(ind_pav_green,m) =                      &
+                                                albedo_pars_f%pars_xy(3,j,i)
+                      surf_usm_h%aldif(ind_pav_green,m) =                      &
+                                                albedo_pars_f%pars_xy(3,j,i)
+                   ENDIF
+                   IF ( albedo_pars_f%pars_xy(4,j,i) /= albedo_pars_f%fill )  THEN
+                      surf_usm_h%asdir(ind_pav_green,m) =                      &
+                                                albedo_pars_f%pars_xy(4,j,i)
+                      surf_usm_h%asdif(ind_pav_green,m) =                      &
+                                                albedo_pars_f%pars_xy(4,j,i)
+                   ENDIF
+!
+!--                Spectral albedos especially for building window surfaces
+                   IF ( albedo_pars_f%pars_xy(5,j,i) /= albedo_pars_f%fill )  THEN
+                      surf_usm_h%aldir(ind_wat_win,m) =                        &
+                                                albedo_pars_f%pars_xy(5,j,i)
+                      surf_usm_h%aldif(ind_wat_win,m) =                        &
+                                                albedo_pars_f%pars_xy(5,j,i)
+                   ENDIF
+                   IF ( albedo_pars_f%pars_xy(6,j,i) /= albedo_pars_f%fill )  THEN
+                      surf_usm_h%asdir(ind_wat_win,m) =                        &
+                                                albedo_pars_f%pars_xy(6,j,i)
+                      surf_usm_h%asdif(ind_wat_win,m) =                        &
+                                                albedo_pars_f%pars_xy(6,j,i)
+                   ENDIF
 
                 ENDDO
@@ -2505,24 +2540,24 @@
                    DO  ind_type = 0, 2
                       IF ( surf_lsm_v(l)%albedo_type(ind_type,m) == 0 )  THEN
-                         IF ( albedo_pars_f%pars_xy(1,j+joff,i+ioff) /=        &
+                         IF ( albedo_pars_f%pars_xy(0,j+joff,i+ioff) /=        &
                               albedo_pars_f%fill )                             &
                             surf_lsm_v(l)%albedo(ind_type,m) =                 &
-                                          albedo_pars_f%pars_xy(1,j+joff,i+ioff)
+                                          albedo_pars_f%pars_xy(0,j+joff,i+ioff)
                          IF ( albedo_pars_f%pars_xy(1,j+joff,i+ioff) /=        &
                               albedo_pars_f%fill )                             &
                             surf_lsm_v(l)%aldir(ind_type,m) =                  &
                                           albedo_pars_f%pars_xy(1,j+joff,i+ioff)
+                         IF ( albedo_pars_f%pars_xy(1,j+joff,i+ioff) /=        &
+                              albedo_pars_f%fill )                             &
+                            surf_lsm_v(l)%aldif(ind_type,m) =                  &
+                                          albedo_pars_f%pars_xy(1,j+joff,i+ioff)
                          IF ( albedo_pars_f%pars_xy(2,j+joff,i+ioff) /=        &
                               albedo_pars_f%fill )                             &
-                            surf_lsm_v(l)%aldif(ind_type,m) =                  &
+                            surf_lsm_v(l)%asdir(ind_type,m) =                  &
                                           albedo_pars_f%pars_xy(2,j+joff,i+ioff)
-                         IF ( albedo_pars_f%pars_xy(3,j+joff,i+ioff) /=        &
-                              albedo_pars_f%fill )                             &
-                            surf_lsm_v(l)%asdir(ind_type,m) =                  &
-                                          albedo_pars_f%pars_xy(3,j+joff,i+ioff)
-                         IF ( albedo_pars_f%pars_xy(4,j+joff,i+ioff) /=        &
+                         IF ( albedo_pars_f%pars_xy(2,j+joff,i+ioff) /=        &
                               albedo_pars_f%fill )                             &
                             surf_lsm_v(l)%asdif(ind_type,m) =                  &
-                                          albedo_pars_f%pars_xy(4,j+joff,i+ioff)
+                                          albedo_pars_f%pars_xy(2,j+joff,i+ioff)
                       ENDIF
                    ENDDO
@@ -2539,30 +2574,61 @@
                       j = surf_usm_v(l)%j(m)
 !
-!--                   Spectral albedos for wall/green/window surfaces
+!--                   Broadband albedos for wall/green/window surfaces
                       DO  ind_type = 0, 2
                          IF ( surf_usm_v(l)%albedo_type(ind_type,m) == 0 )  THEN
-                            IF ( albedo_pars_f%pars_xy(1,j+joff,i+ioff) /=        &
-                                 albedo_pars_f%fill )                             &
-                               surf_usm_v(l)%albedo(ind_type,m) =                 &
-                                             albedo_pars_f%pars_xy(1,j+joff,i+ioff)
-                            IF ( albedo_pars_f%pars_xy(1,j+joff,i+ioff) /=        &
-                                 albedo_pars_f%fill )                             &
-                               surf_usm_v(l)%aldir(ind_type,m) =                  &
-                                             albedo_pars_f%pars_xy(1,j+joff,i+ioff)
-                            IF ( albedo_pars_f%pars_xy(2,j+joff,i+ioff) /=        &
-                                 albedo_pars_f%fill )                             &
-                               surf_usm_v(l)%aldif(ind_type,m) =                  &
-                                             albedo_pars_f%pars_xy(2,j+joff,i+ioff)
-                            IF ( albedo_pars_f%pars_xy(3,j+joff,i+ioff) /=        &
-                                 albedo_pars_f%fill )                             &
-                               surf_usm_v(l)%asdir(ind_type,m) =                  &
-                                             albedo_pars_f%pars_xy(3,j+joff,i+ioff)
-                            IF ( albedo_pars_f%pars_xy(4,j+joff,i+ioff) /=        &
-                                 albedo_pars_f%fill )                             &
-                               surf_usm_v(l)%asdif(ind_type,m) =                  &
-                                             albedo_pars_f%pars_xy(4,j+joff,i+ioff)
+                            IF ( albedo_pars_f%pars_xy(0,j+joff,i+ioff) /=     &
+                                 albedo_pars_f%fill )                          &
+                               surf_usm_v(l)%albedo(ind_type,m) =              &
+                                             albedo_pars_f%pars_xy(0,j+joff,i+ioff)
                          ENDIF
                       ENDDO
-
+!
+!--                   Spectral albedos especially for building wall surfaces
+                      IF ( albedo_pars_f%pars_xy(1,j+joff,i+ioff) /=           &
+                           albedo_pars_f%fill )  THEN
+                         surf_usm_v(l)%aldir(ind_veg_wall,m) =                 &
+                                         albedo_pars_f%pars_xy(1,j+joff,i+ioff)
+                         surf_usm_v(l)%aldif(ind_veg_wall,m) =                 &
+                                         albedo_pars_f%pars_xy(1,j+joff,i+ioff)
+                      ENDIF
+                      IF ( albedo_pars_f%pars_xy(2,j+joff,i+ioff) /=           &
+                           albedo_pars_f%fill )  THEN
+                         surf_usm_v(l)%asdir(ind_veg_wall,m) =                 &
+                                         albedo_pars_f%pars_xy(2,j+joff,i+ioff)
+                         surf_usm_v(l)%asdif(ind_veg_wall,m) =                 &
+                                         albedo_pars_f%pars_xy(2,j+joff,i+ioff)
+                      ENDIF
+!                     
+!--                   Spectral albedos especially for building green surfaces
+                      IF ( albedo_pars_f%pars_xy(3,j+joff,i+ioff) /=           &
+                           albedo_pars_f%fill )  THEN
+                         surf_usm_v(l)%aldir(ind_pav_green,m) =                &
+                                         albedo_pars_f%pars_xy(3,j+joff,i+ioff)
+                         surf_usm_v(l)%aldif(ind_pav_green,m) =                &
+                                         albedo_pars_f%pars_xy(3,j+joff,i+ioff)
+                      ENDIF
+                      IF ( albedo_pars_f%pars_xy(4,j+joff,i+ioff) /=           &
+                           albedo_pars_f%fill )  THEN
+                         surf_usm_v(l)%asdir(ind_pav_green,m) =                &
+                                         albedo_pars_f%pars_xy(4,j+joff,i+ioff)
+                         surf_usm_v(l)%asdif(ind_pav_green,m) =                &
+                                         albedo_pars_f%pars_xy(4,j+joff,i+ioff)
+                      ENDIF
+!                     
+!--                   Spectral albedos especially for building window surfaces
+                      IF ( albedo_pars_f%pars_xy(5,j+joff,i+ioff) /=           &
+                           albedo_pars_f%fill )  THEN
+                         surf_usm_v(l)%aldir(ind_wat_win,m) =                  &
+                                         albedo_pars_f%pars_xy(5,j+joff,i+ioff)
+                         surf_usm_v(l)%aldif(ind_wat_win,m) =                  &
+                                         albedo_pars_f%pars_xy(5,j+joff,i+ioff)
+                      ENDIF
+                      IF ( albedo_pars_f%pars_xy(6,j+joff,i+ioff) /=           &
+                           albedo_pars_f%fill )  THEN
+                         surf_usm_v(l)%asdir(ind_wat_win,m) =                  &
+                                         albedo_pars_f%pars_xy(6,j+joff,i+ioff)
+                         surf_usm_v(l)%asdif(ind_wat_win,m) =                  &
+                                         albedo_pars_f%pars_xy(6,j+joff,i+ioff)
+                      ENDIF
                    ENDDO
                 ENDIF
@@ -9809,4 +9875,5 @@
          var(1:7) == 'rtm_svf'  .OR.  var(1:7) == 'rtm_dif'  .OR.                       &
          var(1:9) == 'rtm_skyvf' .OR. var(1:10) == 'rtm_skyvft'  .OR.                   &
+         var(1:12) == 'rtm_surfalb_'  .OR.  var(1:13) == 'rtm_surfemis_'  .OR.          &
          var == 'rtm_mrt'  .OR.  var ==  'rtm_mrt_sw'  .OR.  var == 'rtm_mrt_lw' )  THEN
 
@@ -10511,40 +10578,4 @@
          ENDIF
 
-!--   block of RTM output variables
-!--   variables are intended mainly for debugging and detailed analyse purposes
-      CASE ( 'rtm_skyvf' )
-!--        sky view factor
-         DO isurf = dirstart(ids), dirend(ids)
-            IF ( surfl(id,isurf) == idsint_u  .OR.  surfl(id,isurf) == idsint_l )  THEN
-               local_pf(surfl(ix,isurf),surfl(iy,isurf),surfl(iz,isurf)) = skyvf(isurf)
-            ENDIF
-         ENDDO
-
-      CASE ( 'rtm_skyvft' )
-!--      sky view factor
-         DO isurf = dirstart(ids), dirend(ids)
-            IF ( surfl(id,isurf) == ids )  THEN
-               local_pf(surfl(ix,isurf),surfl(iy,isurf),surfl(iz,isurf)) = skyvft(isurf)
-            ENDIF
-         ENDDO
-
-      CASE ( 'rtm_svf', 'rtm_dif' )
-!--      shape view factors or iradiance factors to selected surface
-         IF ( TRIM(var)=='rtm_svf' )  THEN
-             k = 1
-         ELSE
-             k = 2
-         ENDIF
-         DO isvf = 1, nsvfl
-            isurflt = svfsurf(1, isvf)
-            isurfs = svfsurf(2, isvf)
-
-            IF ( surf(ix,isurfs) == is  .AND.  surf(iy,isurfs) == js  .AND. surf(iz,isurfs) == ks  .AND. &
-                 (surf(id,isurfs) == idsint_u .OR. surfl(id,isurfs) == idsint_l ) ) THEN
-!--            correct source surface
-               local_pf(surfl(ix,isurflt),surfl(iy,isurflt),surfl(iz,isurflt)) = svf(k,isvf)
-            ENDIF
-         ENDDO
-
       CASE ( 'rtm_rad_net' )
 !--     array of complete radiation balance
@@ -10785,7 +10816,65 @@
             ENDIF
          ENDIF
-
-       CASE DEFAULT
-          found = .FALSE.
+!         
+!--   block of RTM output variables
+!--   variables are intended mainly for debugging and detailed analyse purposes
+      CASE ( 'rtm_skyvf' )
+!      
+!--      sky view factor
+         DO isurf = dirstart(ids), dirend(ids)
+            IF ( surfl(id,isurf) == idsint_u  .OR.  surfl(id,isurf) == idsint_l )  THEN
+               local_pf(surfl(ix,isurf),surfl(iy,isurf),surfl(iz,isurf)) = skyvf(isurf)
+            ENDIF
+         ENDDO
+
+      CASE ( 'rtm_skyvft' )
+!
+!--      sky view factor
+         DO isurf = dirstart(ids), dirend(ids)
+            IF ( surfl(id,isurf) == ids )  THEN
+               local_pf(surfl(ix,isurf),surfl(iy,isurf),surfl(iz,isurf)) = skyvft(isurf)
+            ENDIF
+         ENDDO
+
+      CASE ( 'rtm_svf', 'rtm_dif' )
+!
+!--      shape view factors or iradiance factors to selected surface
+         IF ( TRIM(var)=='rtm_svf' )  THEN
+             k = 1
+         ELSE
+             k = 2
+         ENDIF
+         DO isvf = 1, nsvfl
+            isurflt = svfsurf(1, isvf)
+            isurfs = svfsurf(2, isvf)
+
+            IF ( surf(ix,isurfs) == is  .AND.  surf(iy,isurfs) == js  .AND. surf(iz,isurfs) == ks  .AND. &
+                 (surf(id,isurfs) == idsint_u .OR. surfl(id,isurfs) == idsint_l ) ) THEN
+!
+!--            correct source surface
+               local_pf(surfl(ix,isurflt),surfl(iy,isurflt),surfl(iz,isurflt)) = svf(k,isvf)
+            ENDIF
+         ENDDO
+
+      CASE ( 'rtm_surfalb' )
+!
+!--      surface albedo
+         DO isurf = dirstart(ids), dirend(ids)
+            IF ( surfl(id,isurf) == idsint_u  .OR.  surfl(id,isurf) == idsint_l )  THEN
+               local_pf(surfl(ix,isurf),surfl(iy,isurf),surfl(iz,isurf)) = albedo_surf(isurf)
+            ENDIF
+         ENDDO
+
+      CASE ( 'rtm_surfemis' )
+!
+!--      surface emissivity, weighted average
+         DO isurf = dirstart(ids), dirend(ids)
+            IF ( surfl(id,isurf) == idsint_u  .OR.  surfl(id,isurf) == idsint_l )  THEN
+               local_pf(surfl(ix,isurf),surfl(iy,isurf),surfl(iz,isurf)) = emiss_surf(isurf)
+            ENDIF
+         ENDDO
+
+      CASE DEFAULT
+         found = .FALSE.
 
     END SELECT
Index: palm/trunk/SOURCE/surface_data_output_mod.f90
===================================================================
--- palm/trunk/SOURCE/surface_data_output_mod.f90	(revision 3880)
+++ palm/trunk/SOURCE/surface_data_output_mod.f90	(revision 3881)
@@ -25,4 +25,7 @@
 ! -----------------
 ! $Id$
+! Check for zero output timestep (not allowed in parallel NetCDF output mode) 
+! 
+! 3817 2019-03-26 13:53:57Z suehring
 ! Correct output coordinates of vertical surface elements
 ! 
@@ -125,5 +128,6 @@
    USE netcdf_interface,                                                       &
        ONLY:  netcdf_create_att, netcdf_create_dim, netcdf_create_file,        &
-              netcdf_create_global_atts, netcdf_create_var, netcdf_handle_error
+              netcdf_create_global_atts, netcdf_create_var, netcdf_data_format,&
+              netcdf_handle_error
 
    USE pegrid
@@ -4316,4 +4320,26 @@
 #endif
 !
+!--    In case of parallel NetCDF output the output timestep must not be zero.
+!--    This is because the number of requiered output timesteps is 
+!--    pre-calculated, which is not possible with zero output timestep. 
+       IF ( netcdf_data_format > 4 )  THEN
+          IF ( dt_dosurf == 0.0_wp )  THEN
+             message_string = 'dt_dosurf = 0.0 while using a ' //             &
+                              'variable timestep and parallel netCDF4 ' //    &
+                              'is not allowed.'
+             CALL message( 'surface_data_output_check_parameters', 'PA0081',  &
+                           1, 2, 0, 6, 0 )
+          ENDIF
+          
+          IF ( dt_dosurf_av == 0.0_wp )  THEN
+             message_string = 'dt_dosurf_av = 0.0 while using a ' //          &
+                              'variable timestep and parallel netCDF4 ' //    &
+                              'is not allowed.'
+             CALL message( 'surface_data_output_check_parameters', 'PA0081',  &
+                           1, 2, 0, 6, 0 )
+          ENDIF
+       ENDIF
+
+!
 !--   In case of restart runs, check it the number of cores has been changed.
 !--   With surface output this is not allowed.
Index: palm/trunk/SOURCE/surface_layer_fluxes_mod.f90
===================================================================
--- palm/trunk/SOURCE/surface_layer_fluxes_mod.f90	(revision 3880)
+++ palm/trunk/SOURCE/surface_layer_fluxes_mod.f90	(revision 3881)
@@ -26,4 +26,7 @@
 ! -----------------
 ! $Id$
+! Assure that Obukhov length does not become zero
+! 
+! 3834 2019-03-28 15:40:15Z forkel
 ! added USE chem_gasphase_mod 
 ! 
@@ -1188,4 +1191,8 @@
                 CYCLE
              ENDIF
+!
+!--          Assure that Obukhov length does not become zero
+             IF ( ABS( surf%ol(m) ) < 1E-6_wp )                                &
+                surf%ol(m) = SIGN( 10E-6_wp, surf%ol(m) )              
 
           ENDDO
Index: palm/trunk/SOURCE/urban_surface_mod.f90
===================================================================
--- palm/trunk/SOURCE/urban_surface_mod.f90	(revision 3880)
+++ palm/trunk/SOURCE/urban_surface_mod.f90	(revision 3881)
@@ -28,4 +28,9 @@
 ! -----------------
 ! $Id$
+! Revise determination of local ground-floor level height. 
+! Make level 3 initalization conform with Palm-input-data standard
+! Move output of albedo and emissivity to radiation module
+! 
+! 3832 2019-03-28 13:16:58Z raasch
 ! instrumented with openmp directives
 ! 
@@ -564,129 +569,136 @@
     REAL(wp)     ::  roughness_concrete = 0.001_wp   !< roughness length of average concrete surface
 !
-!-- Indices of input attributes for (above) ground floor level 
-    INTEGER(iwp) ::  ind_alb_wall_agfl     = 65   !< index in input list for albedo_type of wall above ground floor level
-    INTEGER(iwp) ::  ind_alb_wall_gfl      = 32   !< index in input list for albedo_type of wall ground floor level
-    INTEGER(iwp) ::  ind_alb_wall_r        = 96   !< index in input list for albedo_type of wall roof
-    INTEGER(iwp) ::  ind_alb_green_agfl    = 83   !< index in input list for albedo_type of green above ground floor level
-    INTEGER(iwp) ::  ind_alb_green_gfl     = 50   !< index in input list for albedo_type of green ground floor level
-    INTEGER(iwp) ::  ind_alb_green_r       = 115  !< index in input list for albedo_type of green roof
-    INTEGER(iwp) ::  ind_alb_win_agfl      = 79   !< index in input list for albedo_type of window fraction
+!-- Indices of input attributes in building_pars for (above) ground floor level 
+    INTEGER(iwp) ::  ind_alb_wall_agfl     = 38   !< index in input list for albedo_type of wall above ground floor level
+    INTEGER(iwp) ::  ind_alb_wall_gfl      = 66   !< index in input list for albedo_type of wall ground floor level
+    INTEGER(iwp) ::  ind_alb_wall_r        = 101  !< index in input list for albedo_type of wall roof
+    INTEGER(iwp) ::  ind_alb_green_agfl    = 39   !< index in input list for albedo_type of green above ground floor level
+    INTEGER(iwp) ::  ind_alb_green_gfl     = 78   !< index in input list for albedo_type of green ground floor level
+    INTEGER(iwp) ::  ind_alb_green_r       = 117  !< index in input list for albedo_type of green roof
+    INTEGER(iwp) ::  ind_alb_win_agfl      = 40   !< index in input list for albedo_type of window fraction above ground floor level
+    INTEGER(iwp) ::  ind_alb_win_gfl       = 77   !< index in input list for albedo_type of window fraction ground floor level
+    INTEGER(iwp) ::  ind_alb_win_r         = 115  !< index in input list for albedo_type of window fraction roof
+    INTEGER(iwp) ::  ind_c_surface         = 45   !< index in input list for heat capacity wall surface
+    INTEGER(iwp) ::  ind_c_surface_green   = 48   !< index in input list for heat capacity green surface
+    INTEGER(iwp) ::  ind_c_surface_win     = 47   !< index in input list for heat capacity window surface
+    INTEGER(iwp) ::  ind_emis_wall_agfl    = 14   !< index in input list for wall emissivity, above ground floor level
+    INTEGER(iwp) ::  ind_emis_wall_gfl     = 32   !< index in input list for wall emissivity, ground floor level
+    INTEGER(iwp) ::  ind_emis_wall_r       = 100  !< index in input list for wall emissivity, roof
+    INTEGER(iwp) ::  ind_emis_green_agfl   = 15   !< index in input list for green emissivity, above ground floor level
+    INTEGER(iwp) ::  ind_emis_green_gfl    = 34   !< index in input list for green emissivity, ground floor level
+    INTEGER(iwp) ::  ind_emis_green_r      = 116  !< index in input list for green emissivity, roof
+    INTEGER(iwp) ::  ind_emis_win_agfl     = 16   !< index in input list for window emissivity, above ground floor level
+    INTEGER(iwp) ::  ind_emis_win_gfl      = 33   !< index in input list for window emissivity, ground floor level
+    INTEGER(iwp) ::  ind_emis_win_r        = 113  !< index in input list for window emissivity, roof
+    INTEGER(iwp) ::  ind_gflh              = 20   !< index in input list for ground floor level height
+    INTEGER(iwp) ::  ind_green_frac_w_agfl = 2    !< index in input list for green fraction on wall, above ground floor level
+    INTEGER(iwp) ::  ind_green_frac_w_gfl  = 23   !< index in input list for green fraction on wall, ground floor level
+    INTEGER(iwp) ::  ind_green_frac_r_agfl = 3    !< index in input list for green fraction on roof, above ground floor level
+    INTEGER(iwp) ::  ind_green_frac_r_gfl  = 24   !< index in input list for green fraction on roof, ground floor level
+    INTEGER(iwp) ::  ind_hc1_agfl          = 6    !< index in input list for heat capacity at first wall layer,
                                                   !< above ground floor level
-    INTEGER(iwp) ::  ind_alb_win_gfl       = 46   !< index in input list for albedo_type of window fraction ground floor level
-    INTEGER(iwp) ::  ind_alb_win_r         = 110  !< index in input list for albedo_type of window fraction roof
-    INTEGER(iwp) ::  ind_emis_wall_agfl    = 64   !< index in input list for wall emissivity, above ground floor level
-    INTEGER(iwp) ::  ind_emis_wall_gfl     = 31   !< index in input list for wall emissivity, ground floor level
-    INTEGER(iwp) ::  ind_emis_wall_r       = 95   !< index in input list for wall emissivity, roof
-    INTEGER(iwp) ::  ind_emis_green_agfl   = 82   !< index in input list for green emissivity, above ground floor level
-    INTEGER(iwp) ::  ind_emis_green_gfl    = 49   !< index in input list for green emissivity, ground floor level
-    INTEGER(iwp) ::  ind_emis_green_r      = 114  !< index in input list for green emissivity, roof
-    INTEGER(iwp) ::  ind_emis_win_agfl     = 77   !< index in input list for window emissivity, above ground floor level
-    INTEGER(iwp) ::  ind_emis_win_gfl      = 44   !< index in input list for window emissivity, ground floor level
-    INTEGER(iwp) ::  ind_emis_win_r        = 108  !< index in input list for window emissivity, roof
-    INTEGER(iwp) ::  ind_green_frac_w_agfl = 80   !< index in input list for green fraction on wall, above ground floor level
-    INTEGER(iwp) ::  ind_green_frac_w_gfl  = 47   !< index in input list for green fraction on wall, ground floor level
-    INTEGER(iwp) ::  ind_green_frac_r_agfl = 112  !< index in input list for green fraction on roof, above ground floor level
-    INTEGER(iwp) ::  ind_green_frac_r_gfl  = 111  !< index in input list for green fraction on roof, ground floor level
-    INTEGER(iwp) ::  ind_hc1_agfl          = 58   !< index in input list for heat capacity at first wall layer,
+    INTEGER(iwp) ::  ind_hc1_gfl           = 26   !< index in input list for heat capacity at first wall layer, ground floor level
+    INTEGER(iwp) ::  ind_hc1_wall_r        = 94   !< index in input list for heat capacity at first wall layer, roof
+    INTEGER(iwp) ::  ind_hc1_win_agfl      = 83   !< index in input list for heat capacity at first window layer,
                                                   !< above ground floor level
-    INTEGER(iwp) ::  ind_hc1_gfl           = 25   !< index in input list for heat capacity at first wall layer, ground floor level
-    INTEGER(iwp) ::  ind_hc1_wall_r        = 89   !< index in input list for heat capacity at first wall layer, roof
-    INTEGER(iwp) ::  ind_hc1_win_agfl      = 71   !< index in input list for heat capacity at first window layer,
+    INTEGER(iwp) ::  ind_hc1_win_gfl       = 71   !< index in input list for heat capacity at first window layer,
+                                                  !< ground floor level
+    INTEGER(iwp) ::  ind_hc1_win_r         = 107  !< index in input list for heat capacity at first window layer, roof
+    INTEGER(iwp) ::  ind_hc2_agfl          = 7    !< index in input list for heat capacity at second wall layer,
                                                   !< above ground floor level
-    INTEGER(iwp) ::  ind_hc1_win_gfl       = 38   !< index in input list for heat capacity at first window layer,
+    INTEGER(iwp) ::  ind_hc2_gfl           = 27   !< index in input list for heat capacity at second wall layer, ground floor level
+    INTEGER(iwp) ::  ind_hc2_wall_r        = 95   !< index in input list for heat capacity at second wall layer, roof
+    INTEGER(iwp) ::  ind_hc2_win_agfl      = 84   !< index in input list for heat capacity at second window layer,
+                                                  !< above ground floor level
+    INTEGER(iwp) ::  ind_hc2_win_gfl       = 72   !< index in input list for heat capacity at second window layer,
                                                   !< ground floor level
-    INTEGER(iwp) ::  ind_hc1_win_r         = 102  !< index in input list for heat capacity at first window layer, roof
-    INTEGER(iwp) ::  ind_hc2_agfl          = 59   !< index in input list for heat capacity at second wall layer,
+    INTEGER(iwp) ::  ind_hc2_win_r         = 108  !< index in input list for heat capacity at second window layer, roof
+    INTEGER(iwp) ::  ind_hc3_agfl          = 8    !< index in input list for heat capacity at third wall layer,
                                                   !< above ground floor level
-    INTEGER(iwp) ::  ind_hc2_gfl           = 26   !< index in input list for heat capacity at second wall layer, ground floor level
-    INTEGER(iwp) ::  ind_hc2_wall_r        = 90   !< index in input list for heat capacity at second wall layer, roof
-    INTEGER(iwp) ::  ind_hc2_win_agfl      = 72   !< index in input list for heat capacity at second window layer,
+    INTEGER(iwp) ::  ind_hc3_gfl           = 28   !< index in input list for heat capacity at third wall layer, ground floor level
+    INTEGER(iwp) ::  ind_hc3_wall_r        = 96   !< index in input list for heat capacity at third wall layer, roof
+    INTEGER(iwp) ::  ind_hc3_win_agfl      = 85   !< index in input list for heat capacity at third window layer,
                                                   !< above ground floor level
-    INTEGER(iwp) ::  ind_hc2_win_gfl       = 39   !< index in input list for heat capacity at second window layer,
+    INTEGER(iwp) ::  ind_hc3_win_gfl       = 73   !< index in input list for heat capacity at third window layer,
                                                   !< ground floor level
-    INTEGER(iwp) ::  ind_hc2_win_r         = 103  !< index in input list for heat capacity at second window layer, roof
-    INTEGER(iwp) ::  ind_hc3_agfl          = 60   !< index in input list for heat capacity at third wall layer,
+    INTEGER(iwp) ::  ind_hc3_win_r         = 109  !< index in input list for heat capacity at third window layer, roof
+    INTEGER(iwp) ::  ind_indoor_target_temp_summer = 12
+    INTEGER(iwp) ::  ind_indoor_target_temp_winter = 13
+    INTEGER(iwp) ::  ind_lai_r_agfl        = 4    !< index in input list for LAI on roof, above ground floor level
+    INTEGER(iwp) ::  ind_lai_r_gfl         = 4  !< index in input list for LAI on roof, ground floor level
+    INTEGER(iwp) ::  ind_lai_w_agfl        = 5    !< index in input list for LAI on wall, above ground floor level
+    INTEGER(iwp) ::  ind_lai_w_gfl         = 25   !< index in input list for LAI on wall, ground floor level
+    INTEGER(iwp) ::  ind_lambda_surf       = 46   !< index in input list for thermal conductivity of wall surface
+    INTEGER(iwp) ::  ind_lambda_surf_green = 50   !< index in input list for thermal conductivity of green surface
+    INTEGER(iwp) ::  ind_lambda_surf_win   = 49   !< index in input list for thermal conductivity of window surface
+    INTEGER(iwp) ::  ind_tc1_agfl          = 9    !< index in input list for thermal conductivity at first wall layer,
                                                   !< above ground floor level
-    INTEGER(iwp) ::  ind_hc3_gfl           = 27   !< index in input list for heat capacity at third wall layer, ground floor level
-    INTEGER(iwp) ::  ind_hc3_wall_r        = 91   !< index in input list for heat capacity at third wall layer, roof
-    INTEGER(iwp) ::  ind_hc3_win_agfl      = 73   !< index in input list for heat capacity at third window layer,
+    INTEGER(iwp) ::  ind_tc1_gfl           = 29   !< index in input list for thermal conductivity at first wall layer,
+                                                  !< ground floor level
+    INTEGER(iwp) ::  ind_tc1_wall_r        = 97   !< index in input list for thermal conductivity at first wall layer, roof
+    INTEGER(iwp) ::  ind_tc1_win_agfl      = 86   !< index in input list for thermal conductivity at first window layer,
                                                   !< above ground floor level
-    INTEGER(iwp) ::  ind_hc3_win_gfl       = 40   !< index in input list for heat capacity at third window layer,
+    INTEGER(iwp) ::  ind_tc1_win_gfl       = 74   !< index in input list for thermal conductivity at first window layer,
                                                   !< ground floor level
-    INTEGER(iwp) ::  ind_hc3_win_r         = 104  !< index in input list for heat capacity at third window layer, roof
-    INTEGER(iwp) ::  ind_gflh              = 17   !< index in input list for ground floor level height
-    INTEGER(iwp) ::  ind_lai_r_agfl        = 113  !< index in input list for LAI on roof, above ground floor level
-    INTEGER(iwp) ::  ind_lai_r_gfl         = 113  !< index in input list for LAI on roof, ground floor level
-    INTEGER(iwp) ::  ind_lai_w_agfl        = 81   !< index in input list for LAI on wall, above ground floor level
-    INTEGER(iwp) ::  ind_lai_w_gfl         = 48   !< index in input list for LAI on wall, ground floor level
-    INTEGER(iwp) ::  ind_tc1_agfl          = 61   !< index in input list for thermal conductivity at first wall layer,
+    INTEGER(iwp) ::  ind_tc1_win_r         = 110  !< index in input list for thermal conductivity at first window layer, roof
+    INTEGER(iwp) ::  ind_tc2_agfl          = 10   !< index in input list for thermal conductivity at second wall layer,
                                                   !< above ground floor level
-    INTEGER(iwp) ::  ind_tc1_gfl           = 28   !< index in input list for thermal conductivity at first wall layer,
+    INTEGER(iwp) ::  ind_tc2_gfl           = 30   !< index in input list for thermal conductivity at second wall layer,
                                                   !< ground floor level
-    INTEGER(iwp) ::  ind_tc1_wall_r        = 92   !< index in input list for thermal conductivity at first wall layer, roof
-    INTEGER(iwp) ::  ind_tc1_win_agfl      = 74   !< index in input list for thermal conductivity at first window layer,
+    INTEGER(iwp) ::  ind_tc2_wall_r        = 98   !< index in input list for thermal conductivity at second wall layer, roof
+    INTEGER(iwp) ::  ind_tc2_win_agfl      = 87   !< index in input list for thermal conductivity at second window layer,
                                                   !< above ground floor level
-    INTEGER(iwp) ::  ind_tc1_win_gfl       = 41   !< index in input list for thermal conductivity at first window layer,
+    INTEGER(iwp) ::  ind_tc2_win_gfl       = 75   !< index in input list for thermal conductivity at second window layer,
                                                   !< ground floor level
-    INTEGER(iwp) ::  ind_tc1_win_r         = 105  !< index in input list for thermal conductivity at first window layer, roof
-    INTEGER(iwp) ::  ind_tc2_agfl          = 62   !< index in input list for thermal conductivity at second wall layer,
+    INTEGER(iwp) ::  ind_tc2_win_r         = 111  !< index in input list for thermal conductivity at second window layer,
+                                                  !< ground floor level
+    INTEGER(iwp) ::  ind_tc3_agfl          = 11   !< index in input list for thermal conductivity at third wall layer,
                                                   !< above ground floor level
-    INTEGER(iwp) ::  ind_tc2_gfl           = 29   !< index in input list for thermal conductivity at second wall layer,
+    INTEGER(iwp) ::  ind_tc3_gfl           = 31   !< index in input list for thermal conductivity at third wall layer,
                                                   !< ground floor level
-    INTEGER(iwp) ::  ind_tc2_wall_r        = 93   !< index in input list for thermal conductivity at second wall layer, roof
-    INTEGER(iwp) ::  ind_tc2_win_agfl      = 75   !< index in input list for thermal conductivity at second window layer,
+    INTEGER(iwp) ::  ind_tc3_wall_r        = 99   !< index in input list for thermal conductivity at third wall layer, roof
+    INTEGER(iwp) ::  ind_tc3_win_agfl      = 88   !< index in input list for thermal conductivity at third window layer,
                                                   !< above ground floor level
-    INTEGER(iwp) ::  ind_tc2_win_gfl       = 42   !< index in input list for thermal conductivity at second window layer,
+    INTEGER(iwp) ::  ind_tc3_win_gfl       = 76   !< index in input list for thermal conductivity at third window layer,
                                                   !< ground floor level
-    INTEGER(iwp) ::  ind_tc2_win_r         = 106  !< index in input list for thermal conductivity at second window layer,
-                                                  !< ground floor level
-    INTEGER(iwp) ::  ind_tc3_agfl          = 63   !< index in input list for thermal conductivity at third wall layer,
-                                                  !< above ground floor level
-    INTEGER(iwp) ::  ind_tc3_gfl           = 30   !< index in input list for thermal conductivity at third wall layer,
-                                                  !< ground floor level
-    INTEGER(iwp) ::  ind_tc3_wall_r        = 94   !< index in input list for thermal conductivity at third wall layer, roof
-    INTEGER(iwp) ::  ind_tc3_win_agfl      = 76   !< index in input list for thermal conductivity at third window layer,
-                                                  !< above ground floor level
-    INTEGER(iwp) ::  ind_tc3_win_gfl       = 43   !< index in input list for thermal conductivity at third window layer,
-                                                  !< ground floor level
-    INTEGER(iwp) ::  ind_tc3_win_r         = 107  !< index in input list for thermal conductivity at third window layer, roof
-    INTEGER(iwp) ::  ind_thick_1_agfl      = 54   !< index for wall layer thickness - 1st layer above ground floor level
-    INTEGER(iwp) ::  ind_thick_1_gfl       = 21   !< index for wall layer thickness - 1st layer ground floor level
-    INTEGER(iwp) ::  ind_thick_1_wall_r    = 85   !< index for wall layer thickness - 1st layer roof
-    INTEGER(iwp) ::  ind_thick_1_win_agfl  = 67   !< index for window layer thickness - 1st layer above ground floor level
-    INTEGER(iwp) ::  ind_thick_1_win_gfl   = 34   !< index for window layer thickness - 1st layer ground floor level
-    INTEGER(iwp) ::  ind_thick_1_win_r     = 98   !< index for window layer thickness - 1st layer roof
-    INTEGER(iwp) ::  ind_thick_2_agfl      = 55   !< index for wall layer thickness - 2nd layer above ground floor level
-    INTEGER(iwp) ::  ind_thick_2_gfl       = 22   !< index for wall layer thickness - 2nd layer ground floor level
-    INTEGER(iwp) ::  ind_thick_2_wall_r    = 86   !< index for wall layer thickness - 2nd layer roof
-    INTEGER(iwp) ::  ind_thick_2_win_agfl  = 68   !< index for window layer thickness - 2nd layer above ground floor level
-    INTEGER(iwp) ::  ind_thick_2_win_gfl   = 35   !< index for window layer thickness - 2nd layer ground floor level
-    INTEGER(iwp) ::  ind_thick_2_win_r     = 99   !< index for window layer thickness - 2nd layer roof
-    INTEGER(iwp) ::  ind_thick_3_agfl      = 56   !< index for wall layer thickness - 3rd layer above ground floor level
-    INTEGER(iwp) ::  ind_thick_3_gfl       = 23   !< index for wall layer thickness - 3rd layer ground floor level
-    INTEGER(iwp) ::  ind_thick_3_wall_r    = 87   !< index for wall layer thickness - 3rd layer roof
-    INTEGER(iwp) ::  ind_thick_3_win_agfl  = 69   !< index for window layer thickness - 3rd layer above ground floor level
-    INTEGER(iwp) ::  ind_thick_3_win_gfl   = 36   !< index for window layer thickness - 3rd layer ground floor level  
-    INTEGER(iwp) ::  ind_thick_3_win_r     = 100  !< index for window layer thickness - 3rd layer roof
-    INTEGER(iwp) ::  ind_thick_4_agfl      = 57   !< index for wall layer thickness - 4th layer above ground floor level
-    INTEGER(iwp) ::  ind_thick_4_gfl       = 24   !< index for wall layer thickness - 4th layer ground floor level
-    INTEGER(iwp) ::  ind_thick_4_wall_r    = 88   !< index for wall layer thickness - 4st layer roof
-    INTEGER(iwp) ::  ind_thick_4_win_agfl  = 70   !< index for window layer thickness - 4th layer above ground floor level
-    INTEGER(iwp) ::  ind_thick_4_win_gfl   = 37   !< index for window layer thickness - 4th layer ground floor level
-    INTEGER(iwp) ::  ind_thick_4_win_r     = 101  !< index for window layer thickness - 4th layer roof
-    INTEGER(iwp) ::  ind_trans_agfl        = 78   !< index in input list for window transmissivity, above ground floor level
-    INTEGER(iwp) ::  ind_trans_gfl         = 45   !< index in input list for window transmissivity, ground floor level
-    INTEGER(iwp) ::  ind_trans_r           = 109  !< index in input list for window transmissivity, roof
-    INTEGER(iwp) ::  ind_wall_frac_agfl    = 53   !< index in input list for wall fraction, above ground floor level
-    INTEGER(iwp) ::  ind_wall_frac_gfl     = 20   !< index in input list for wall fraction, ground floor level
-    INTEGER(iwp) ::  ind_wall_frac_r       = 84   !< index in input list for wall fraction, roof
-    INTEGER(iwp) ::  ind_win_frac_agfl     = 66   !< index in input list for window fraction, above ground floor level
-    INTEGER(iwp) ::  ind_win_frac_gfl      = 33   !< index in input list for window fraction, ground floor level
-    INTEGER(iwp) ::  ind_win_frac_r        = 97   !< index in input list for window fraction, roof
-    INTEGER(iwp) ::  ind_z0_agfl           = 51   !< index in input list for z0, above ground floor level
-    INTEGER(iwp) ::  ind_z0_gfl            = 18   !< index in input list for z0, ground floor level
-    INTEGER(iwp) ::  ind_z0qh_agfl         = 52   !< index in input list for z0h / z0q, above ground floor level
-    INTEGER(iwp) ::  ind_z0qh_gfl          = 19   !< index in input list for z0h / z0q, ground floor level
-    INTEGER(iwp) ::  ind_green_type_roof   = 116  !< index in input list for type of green roof
+    INTEGER(iwp) ::  ind_tc3_win_r         = 112  !< index in input list for thermal conductivity at third window layer, roof
+    INTEGER(iwp) ::  ind_thick_1_agfl      = 41   !< index for wall layer thickness - 1st layer above ground floor level
+    INTEGER(iwp) ::  ind_thick_1_gfl       = 62   !< index for wall layer thickness - 1st layer ground floor level
+    INTEGER(iwp) ::  ind_thick_1_wall_r    = 90   !< index for wall layer thickness - 1st layer roof
+    INTEGER(iwp) ::  ind_thick_1_win_agfl  = 79   !< index for window layer thickness - 1st layer above ground floor level
+    INTEGER(iwp) ::  ind_thick_1_win_gfl   = 67   !< index for window layer thickness - 1st layer ground floor level
+    INTEGER(iwp) ::  ind_thick_1_win_r     = 103  !< index for window layer thickness - 1st layer roof
+    INTEGER(iwp) ::  ind_thick_2_agfl      = 42   !< index for wall layer thickness - 2nd layer above ground floor level
+    INTEGER(iwp) ::  ind_thick_2_gfl       = 63   !< index for wall layer thickness - 2nd layer ground floor level
+    INTEGER(iwp) ::  ind_thick_2_wall_r    = 91   !< index for wall layer thickness - 2nd layer roof
+    INTEGER(iwp) ::  ind_thick_2_win_agfl  = 80   !< index for window layer thickness - 2nd layer above ground floor level
+    INTEGER(iwp) ::  ind_thick_2_win_gfl   = 68   !< index for window layer thickness - 2nd layer ground floor level
+    INTEGER(iwp) ::  ind_thick_2_win_r     = 104  !< index for window layer thickness - 2nd layer roof
+    INTEGER(iwp) ::  ind_thick_3_agfl      = 43   !< index for wall layer thickness - 3rd layer above ground floor level
+    INTEGER(iwp) ::  ind_thick_3_gfl       = 64   !< index for wall layer thickness - 3rd layer ground floor level
+    INTEGER(iwp) ::  ind_thick_3_wall_r    = 92   !< index for wall layer thickness - 3rd layer roof
+    INTEGER(iwp) ::  ind_thick_3_win_agfl  = 81   !< index for window layer thickness - 3rd layer above ground floor level
+    INTEGER(iwp) ::  ind_thick_3_win_gfl   = 69   !< index for window layer thickness - 3rd layer ground floor level  
+    INTEGER(iwp) ::  ind_thick_3_win_r     = 105  !< index for window layer thickness - 3rd layer roof
+    INTEGER(iwp) ::  ind_thick_4_agfl      = 44   !< index for wall layer thickness - 4th layer above ground floor level
+    INTEGER(iwp) ::  ind_thick_4_gfl       = 65   !< index for wall layer thickness - 4th layer ground floor level
+    INTEGER(iwp) ::  ind_thick_4_wall_r    = 93   !< index for wall layer thickness - 4st layer roof
+    INTEGER(iwp) ::  ind_thick_4_win_agfl  = 82   !< index for window layer thickness - 4th layer above ground floor level
+    INTEGER(iwp) ::  ind_thick_4_win_gfl   = 70   !< index for window layer thickness - 4th layer ground floor level
+    INTEGER(iwp) ::  ind_thick_4_win_r     = 106  !< index for window layer thickness - 4th layer roof
+    INTEGER(iwp) ::  ind_trans_agfl        = 17   !< index in input list for window transmissivity, above ground floor level
+    INTEGER(iwp) ::  ind_trans_gfl         = 35   !< index in input list for window transmissivity, ground floor level
+    INTEGER(iwp) ::  ind_trans_r           = 114  !< index in input list for window transmissivity, roof
+    INTEGER(iwp) ::  ind_wall_frac_agfl    = 0    !< index in input list for wall fraction, above ground floor level
+    INTEGER(iwp) ::  ind_wall_frac_gfl     = 21   !< index in input list for wall fraction, ground floor level
+    INTEGER(iwp) ::  ind_wall_frac_r       = 89   !< index in input list for wall fraction, roof
+    INTEGER(iwp) ::  ind_win_frac_agfl     = 1    !< index in input list for window fraction, above ground floor level
+    INTEGER(iwp) ::  ind_win_frac_gfl      = 22   !< index in input list for window fraction, ground floor level
+    INTEGER(iwp) ::  ind_win_frac_r        = 102  !< index in input list for window fraction, roof
+    INTEGER(iwp) ::  ind_z0_agfl           = 18   !< index in input list for z0, above ground floor level
+    INTEGER(iwp) ::  ind_z0_gfl            = 36   !< index in input list for z0, ground floor level
+    INTEGER(iwp) ::  ind_z0qh_agfl         = 19   !< index in input list for z0h / z0q, above ground floor level
+    INTEGER(iwp) ::  ind_z0qh_gfl          = 37   !< index in input list for z0h / z0q, ground floor level
+    INTEGER(iwp) ::  ind_green_type_roof   = 118  !< index in input list for type of green roof
 
 
@@ -705,278 +717,949 @@
                                    'bridges                              '  &  !< type  7
                                                                      /)
-!
-!-- building parameters, 6 different types
-!-- Parameter for urban surface model
-!-- 0 - heat capacity wall surface, 1 - heat capacity of window surface, 2 - heat capacity of green surface
-!-- 3 - thermal conductivity of wall surface, 4 - thermal conductivity of window surface, 
-!-- 5 - thermal conductivty of green surface, 6 - wall fraction ground plate, 
-!-- 7 - 1st wall layer thickness ground plate, 8 - 2nd wall layer thickness ground plate
-!-- 9 - 3rd wall layer thickness ground plate, 10 - 4th wall layer thickness ground plate,
-!-- 11 - heat capacity 1st/2nd wall layer ground plate, 12 - heat capacity 3rd wall layer ground plate
-!-- 13 - heat capacity 4th wall layer ground plate, 14 - thermal conductivity 1st/2nd wall layer ground plate,
-!-- 15 - thermal conductivity 3rd wall layer ground plate, 16 - thermal conductivity 4th wall layer ground plate
-!-- 17 - ground floor level height, 18 - z0 roughness ground floor level, 19 - z0h/z0g roughness heaat/humidity,
-!-- 20 - wall fraction ground floor level, 21 - 1st wall layer thickness ground floor level, 
-!-- 22 - 2nd wall layer thickness ground floor level, 23 - 3rd wall layer thickness ground floor level,
-!-- 24 - 4th wall layer thickness ground floor level, 25 - heat capacity 1st/2nd wall layer ground floor level,
-!-- 26 - heat capacity 3rd wall layer ground floor level, 27 - heat capacity 4th wall layer ground floor level,
-!-- 28 - thermal conductivity 1st/2nd wall layer ground floor level, 
-!-- 29 - thermal conductivity 3rd wall layer ground floor level, 30 - thermal conductivity 4th wall layer ground floor level
-!-- 31 - wall emissivity ground floor level, 32 - wall albedo ground floor level, 33 - window fraction ground floor level,
-!-- 34 - 1st window layer thickness ground floor level, 35 - 2nd window layer thickness ground floor level,
-!-- 36 - 3rd window layer thickness ground floor level, 37 - 4th window layer thickness ground floor level,
-!-- 38 - heat capacity 1st/2nd window layer ground floor level, 39 - heat capacity 3rd window layer ground floor level,
-!-- 40 - heat capacity 4th window layer ground floor level, 
-!-- 41 - thermal conductivity 1st/2nd window layer ground floor level, 
-!-- 42 - thermal conductivity 3rd window layer ground floor level, 
-!-- 43 - thermal conductivity 4th window layer ground floor level, 44 - window emissivity ground floor level,
-!-- 45 - window transmissivity ground floor level, 46 - window albedo ground floor level,
-!-- 47 - green fraction ground floor level, 48 - LAI on wall ground floor level, 49 - green emissivity ground floor level,
-!-- 50 - green albedo ground floor level, 51 - z0 roughness above ground floor level, 
-!-- 52 - z0h/z0g roughness heat/humidity above ground floor level, 53 - wall fraction above ground floor level
-!-- 54 - 1st wall layer thickness above ground floor level, 55 - 2nd wall layer thickness above ground floor level
-!-- 56 - 3rd wall layer thickness above ground floor level, 57 - 4th wall layer thickness above ground floor level 
-!-- 58 - heat capacity 1st/2nd wall layer above ground floor level,
-!-- 59 - heat capacity 3rd wall layer above ground floor level, 
-!-- 60 - heat capacity 4th wall layer above ground floor level,
-!-- 61 - thermal conductivity 1st/2nd wall layer above ground floor level,
-!-- 62 - thermal conductivity 3rd wall layer above ground floor level,
-!-- 63 - thermal conductivity 4th wall layer above ground floor level,
-!-- 64 - wall emissivity above ground floor level, 65 - wall albedo above ground floor level,
-!-- 66 - window fraction above ground floor level, 67 - 1st window layer thickness above ground floor level,
-!-- 68 - 2nd thickness window layer above ground floor level, 69 - 3rd window layer thickness above ground floor level,
-!-- 70 - 4th window layer thickness above ground floor level, 
-!-- 71 - heat capacity 1st/2nd window layer above ground floor level,
-!-- 72 - heat capacity 3rd window layer above ground floor level,
-!-- 73 - heat capacity 4th window layer above ground floor level, 
-!-- 74 - conductivity 1st/2nd window layer above ground floor level,
-!-- 75 - thermal conductivity 3rd window layer above ground floor level, 
-!-- 76 - thermal conductivity 4th window layer above ground floor level, 77 - window emissivity above ground floor level,
-!-- 78 - window transmissivity above ground floor level, 79 - window albedo above ground floor level,
-!-- 80 - green fraction above ground floor level, 81 - LAI on wall above ground floor level,
-!-- 82 - green emissivity above ground floor level, 83 - green albedo above ground floor level,
-!-- 84 - wall fraction roof, 85 - 1st wall layer thickness roof, 86 - 2nd wall layer thickness roof,
-!-- 87 - 3rd wall layer thickness roof, 88 - 4th wall layer thickness roof,
-!-- 89 - heat capacity 1st/2nd wall layer roof, 90 - heat capacity 3rd wall layer roof,
-!-- 91 - heat capacity 4th wall layer roof, 92 - thermal conductivity 1st/2nd wall layer roof,
-!-- 93 - thermal conductivity 3rd wall layer roof, 94 - thermal conductivity 4th wall layer roof,
-!-- 95 - wall emissivity roof, 96 - wall albedo roof, 97 - window fraction roof,
-!-- 98 - window 1st layer thickness roof, 99 - window 2nd layer thickness roof, 100 - window 3rd layer thickness roof,
-!-- 101 - window 4th layer thickness, 102 - heat capacity 1st/2nd window layer roof, 
-!-- 103 - heat capacity 3rd window layer roof, 104 - heat capacity 4th window layer roof,
-!-- 105 - thermal conductivity 1st/2nd window layer roof, 106 - thermal conductivity 3rd window layer roof,
-!-- 107 - thermal conductivity 4th window layer roof, 108 - window emissivity roof, 109 - window transmissivity roof,
-!-- 110 - window albedo roof, 111 - green fraction roof ground floor level, 
-!-- 112 - green fraction roof above ground floor level, 113 - LAI roof, 114 - green emissivity roof,
-!-- 115 - green albedo roof, 116 - green type roof,
-!-- Parameter for indoor model
-!-- 117 - indoor target summer temperature, 118 - indoor target winter temperature, 
-!-- 119 - shading factor, 120 - g-value windows, 121 - u-value windows, 122 - basical airflow without occupancy of the room,
-!-- 123 - additional airflow depend of occupancy of the room, 124 - heat recovery efficiency,
-!-- 125 - dynamic parameter specific effective surface, 126 - dynamic parameter innner heatstorage,
-!-- 127 - ratio internal surface/floor area, 128 - maximal heating capacity, 129 - maximal cooling capacity,
-!-- 130 - additional internal heat gains dependent on occupancy of the room, 
-!-- 131 - basic internal heat gains without occupancy of the room, 132 - storey height, 133 - ceiling construction height
-
-
+
+
+!
+!-- Building facade/wall/green/window properties (partly according to PIDS)
     REAL(wp), DIMENSION(0:133,1:7), PARAMETER :: building_pars = RESHAPE( (/   &
-        20000.0_wp, 20000.0_wp, 20000.0_wp, 23.0_wp, 23.0_wp, 10.0_wp,         & !parameter 0-5
-        1.0_wp, 0.005_wp, 0.01_wp, 0.39_wp, 0.63_wp, 2200000.0_wp,             & !parameter 6-11
-        1400000.0_wp, 1300000.0_wp, 0.35_wp, 0.8_wp, 2.1_wp, 4.0_wp,           & !parameter 12-17
-        0.01_wp, 0.001_wp, 0.75_wp,                                            & !parameter 18-20
-        0.005_wp, 0.01_wp, 0.39_wp, 0.63_wp, 2200000.0_wp,                     & !parameter 21-25
-        1400000.0_wp, 1300000.0_wp, 0.35_wp,                                   & !parameter 26-28                     
-        0.8_wp, 2.1_wp, 0.93_wp,                                               & !parameter 29-31       
-        27.0_wp, 0.25_wp, 0.003_wp, 0.006_wp, 0.012_wp, 0.018_wp,              & !parameter 32-37
-        1736000.0_wp, 1736000.0_wp, 1736000.0_wp,                              & !parameter 38-40
-        0.57_wp, 0.57_wp, 0.57_wp, 0.91_wp,                                    & !parameter 41-44
-        0.75_wp, 27.0_wp, 0.0_wp, 1.5_wp, 0.86_wp,                             & !parameter 45-49
-        5.0_wp, 0.001_wp, 0.0001_wp, 0.7_wp, 0.005_wp,                         & !parameter 50-54
-        0.01_wp, 0.39_wp, 0.63_wp, 2200000.0_wp,                               & !parameter 55-58
-        1400000.0_wp, 1300000.0_wp, 0.35_wp, 0.8_wp,                           & !parameter 59-62
-        2.1_wp, 0.93_wp, 27.0_wp, 0.3_wp,                                      & !parameter 63-66
-        0.003_wp, 0.006_wp, 0.012_wp, 0.018_wp,                                & !parameter 67-70
-        1736000.0_wp, 1736000.0_wp, 1736000.0_wp,                              & !parameter 71-73
-        0.57_wp, 0.57_wp, 0.57_wp, 0.91_wp, 0.75_wp,                           & !parameter 74-78
-        27.0_wp, 0.0_wp, 1.5_wp, 0.86_wp, 5.0_wp, 1.0_wp,                      & !parameter 79-84
-        0.005_wp, 0.01_wp, 0.31_wp, 0.63_wp, 2200000.0_wp, 1400000.0_wp,       & !parameter 85-90
-        1300000.0_wp, 0.35_wp, 0.8_wp, 2.1_wp, 0.93_wp, 27.0_wp, 0.0_wp,       & !parameter 91-97
-        0.003_wp, 0.006_wp, 0.012_wp, 0.018_wp, 1736000.0_wp,                  & !parameter 98-102 
-        1736000.0_wp, 1736000.0_wp, 0.57_wp, 0.57_wp, 0.57_wp,                 & !parameter 103-107
-        0.91_wp, 0.75_wp, 27.0_wp, 0.0_wp, 0.0_wp, 1.5_wp,                     & !parameter 108-113
-        0.86_wp, 5.0_wp, 0.0_wp,                                               & !parameter 114-116
-        299.15_wp, 293.15_wp, 0.8_wp, 0.76_wp, 5.0_wp,                         & !parameter 117-121
-        0.1_wp, 0.5_wp, 0.0_wp, 3.5_wp, 370000.0_wp, 4.5_wp,                   & !parameter 122-127
-        100000.0_wp, 0.0_wp, 3.0_wp, 10.0_wp, 3.0_wp, 0.2_wp,                  & !parameter 128-133- end of type 1
-        20000.0_wp, 20000.0_wp, 20000.0_wp, 23.0_wp, 23.0_wp, 10.0_wp,         & !parameter 0-5
-        1.0_wp, 0.005_wp, 0.01_wp, 0.31_wp, 0.42_wp, 2000000.0_wp,             & !parameter 6-11
-        103000.0_wp, 900000.0_wp, 0.35_wp, 0.38_wp, 0.04_wp, 4.0_wp,           & !parameter 12-17
-        0.01_wp, 0.001_wp, 0.78_wp,                                            & !parameter 18-20
-        0.005_wp, 0.01_wp, 0.31_wp, 0.43_wp, 2000000.0_wp,                     & !parameter 21-25
-        103000.0_wp, 900000.0_wp, 0.35_wp,                                     & !parameter 26-28                     
-        0.38_wp, 0.04_wp, 0.92_wp,                                             & !parameter 29-31       
-        27.0_wp, 0.22_wp, 0.003_wp, 0.006_wp, 0.012_wp, 0.018_wp,              & !parameter 32-37
-        1736000.0_wp, 1736000.0_wp, 1736000.0_wp,                              & !parameter 38-40
-        0.11_wp, 0.11_wp, 0.11_wp, 0.11_wp,                                    & !parameter 41-44
-        0.7_wp, 27.0_wp, 0.0_wp, 1.5_wp, 0.86_wp,                              & !parameter 45-49
-        5.0_wp, 0.001_wp, 0.0001_wp, 0.73_wp, 0.005_wp,                        & !parameter 50-54
-        0.01_wp, 0.31_wp, 0.43_wp, 2000000.0_wp,                               & !parameter 55-58
-        103000.0_wp, 900000.0_wp, 0.35_wp, 0.38_wp,                            & !parameter 59-62
-        0.04_wp, 0.92_wp, 27.0_wp, 0.27_wp,                                    & !parameter 63-66
-        0.003_wp, 0.006_wp, 0.012_wp, 0.018_wp,                                & !parameter 67-70
-        1736000.0_wp, 1736000.0_wp, 1736000.0_wp,                              & !parameter 71-73
-        0.11_wp, 0.11_wp, 0.11_wp, 0.87_wp, 0.7_wp,                            & !parameter 74-78
-        27.0_wp, 0.0_wp, 1.5_wp, 0.86_wp, 5.0_wp, 1.0_wp,                      & !parameter 79-84
-        0.005_wp, 0.01_wp, 0.5_wp, 0.79_wp, 2000000.0_wp, 103000.0_wp,         & !parameter 85-90
-        900000.0_wp, 0.35_wp, 0.38_wp, 0.04_wp, 0.93_wp, 27.0_wp, 0.0_wp,      & !parameter 91-97
-        0.003_wp, 0.006_wp, 0.012_wp, 0.018_wp, 1736000.0_wp,                  & !parameter 98-102 
-        1736000.0_wp, 1736000.0_wp, 0.11_wp, 0.11_wp, 0.11_wp,                 & !parameter 103-107
-        0.87_wp, 0.7_wp, 27.0_wp, 0.0_wp, 0.0_wp, 1.5_wp,                      & !parameter 108-113
-        0.86_wp, 5.0_wp, 0.0_wp,                                               & !parameter 114-116
-        299.15_wp, 293.15_wp, 0.8_wp, 0.6_wp, 3.0_wp,                          & !parameter 117-121
-        0.1_wp, 0.5_wp, 0.0_wp, 2.5_wp, 165000.0_wp, 4.5_wp,                   & !parameter 122-127
-        100000.0_wp, 0.0_wp, 4.0_wp, 8.0_wp, 3.0_wp, 0.2_wp,                   & !parameter 128-133- end of type 2
-        20000.0_wp, 20000.0_wp, 20000.0_wp, 23.0_wp, 23.0_wp, 10.0_wp,         & !parameter 0-5
-        1.0_wp, 0.005_wp, 0.01_wp, 0.41_wp, 0.7_wp, 2000000.0_wp,              & !parameter 6-11
-        103000.0_wp, 900000.0_wp, 0.35_wp, 0.14_wp, 0.035_wp, 4.0_wp,          & !parameter 12-17
-        0.01_wp, 0.001_wp, 0.75_wp,                                            & !parameter 18-20
-        0.005_wp, 0.01_wp, 0.41_wp, 0.7_wp, 2000000.0_wp,                      & !parameter 21-25
-        103000.0_wp, 900000.0_wp, 0.35_wp,                                     & !parameter 26-28                     
-        0.14_wp, 0.035_wp, 0.92_wp,                                            & !parameter 29-31       
-        27.0_wp, 0.25_wp, 0.003_wp, 0.006_wp, 0.012_wp, 0.018_wp,              & !parameter 32-37
-        1736000.0_wp, 1736000.0_wp, 1736000.0_wp,                              & !parameter 38-40
-        0.037_wp, 0.037_wp, 0.037_wp, 0.8_wp,                                  & !parameter 41-44
-        0.6_wp, 27.0_wp, 0.0_wp, 1.5_wp, 0.86_wp,                              & !parameter 45-49
-        5.0_wp, 0.001_wp, 0.0001_wp, 0.7_wp, 0.005_wp,                         & !parameter 50-54
-        0.01_wp, 0.41_wp, 0.7_wp, 2000000.0_wp,                                & !parameter 55-58
-        103000.0_wp, 900000.0_wp, 0.35_wp, 0.14_wp,                            & !parameter 59-62
-        0.035_wp, 0.92_wp, 27.0_wp, 0.3_wp,                                    & !parameter 63-66
-        0.003_wp, 0.006_wp, 0.012_wp, 0.018_wp,                                & !parameter 67-70
-        1736000.0_wp, 1736000.0_wp, 1736000.0_wp,                              & !parameter 71-73
-        0.037_wp, 0.037_wp, 0.037_wp, 0.8_wp, 0.6_wp,                          & !parameter 74-78
-        27.0_wp, 0.0_wp, 1.5_wp, 0.86_wp, 5.0_wp, 1.0_wp,                      & !parameter 79-84
-        0.005_wp, 0.01_wp, 0.41_wp, 0.7_wp, 2000000.0_wp, 103000.0_wp,         & !parameter 85-90
-        900000.0_wp, 0.35_wp, 0.14_wp, 0.035_wp, 0.93_wp, 27.0_wp, 0.0_wp,     & !parameter 91-97
-        0.003_wp, 0.006_wp, 0.012_wp, 0.018_wp, 1736000.0_wp,                  & !parameter 98-102 
-        1736000.0_wp, 1736000.0_wp, 0.037_wp, 0.037_wp, 0.037_wp,              & !parameter 103-107
-        0.8_wp, 0.6_wp, 27.0_wp, 0.0_wp, 0.0_wp, 1.5_wp,                       & !parameter 108-113
-        0.86_wp, 5.0_wp, 0.0_wp,                                               & !parameter 114-116
-        299.15_wp, 293.15_wp, 0.8_wp, 0.5_wp, 0.6_wp,                          & !parameter 117-121
-        0.1_wp, 0.5_wp, 0.8_wp, 2.5_wp, 80000.0_wp, 4.5_wp,                    & !parameter 122-127
-        100000.0_wp, 0.0_wp, 3.0_wp, 8.0_wp, 3.0_wp, 0.2_wp,                   & !parameter 128-133- end of type 3
-        20000.0_wp, 20000.0_wp, 20000.0_wp, 23.0_wp, 23.0_wp, 10.0_wp,         & !parameter 0-5
-        1.0_wp, 0.005_wp, 0.01_wp, 0.39_wp, 0.63_wp, 2200000.0_wp,             & !parameter 6-11
-        1400000.0_wp, 1300000.0_wp, 0.35_wp, 0.8_wp, 2.1_wp, 4.0_wp,           & !parameter 12-17
-        0.01_wp, 0.001_wp, 0.55_wp,                                            & !parameter 18-20
-        0.005_wp, 0.01_wp, 0.39_wp, 0.63_wp, 2200000.0_wp,                     & !parameter 21-25
-        1400000.0_wp, 1300000.0_wp, 0.35_wp,                                   & !parameter 26-28                     
-        0.8_wp, 2.1_wp, 0.93_wp,                                               & !parameter 29-31       
-        27.0_wp, 0.45_wp, 0.003_wp, 0.006_wp, 0.012_wp, 0.018_wp,              & !parameter 32-37
-        1736000.0_wp, 1736000.0_wp, 1736000.0_wp,                              & !parameter 38-40
-        0.57_wp, 0.57_wp, 0.57_wp, 0.91_wp,                                    & !parameter 41-44
-        0.75_wp, 27.0_wp, 0.0_wp, 1.5_wp, 0.86_wp,                             & !parameter 45-49
-        5.0_wp, 0.001_wp, 0.0001_wp, 0.5_wp, 0.005_wp,                         & !parameter 50-54
-        0.01_wp, 0.39_wp, 0.63_wp, 2200000.0_wp,                               & !parameter 55-58
-        1400000.0_wp, 1300000.0_wp, 0.35_wp, 0.8_wp,                           & !parameter 59-62
-        2.1_wp, 0.93_wp, 27.0_wp, 0.5_wp,                                      & !parameter 63-66
-        0.003_wp, 0.006_wp, 0.012_wp, 0.018_wp,                                & !parameter 67-70
-        1736000.0_wp, 1736000.0_wp, 1736000.0_wp,                              & !parameter 71-73
-        0.57_wp, 0.57_wp, 0.57_wp, 0.91_wp, 0.75_wp,                           & !parameter 74-78
-        27.0_wp, 0.0_wp, 1.5_wp, 0.86_wp, 5.0_wp, 1.0_wp,                      & !parameter 79-84
-        0.005_wp, 0.01_wp, 0.39_wp, 0.63_wp, 2200000.0_wp, 1400000.0_wp,       & !parameter 85-90
-        1300000.0_wp, 0.35_wp, 0.8_wp, 2.1_wp, 0.93_wp, 27.0_wp, 0.0_wp,       & !parameter 91-97
-        0.003_wp, 0.006_wp, 0.012_wp, 0.018_wp, 1736000.0_wp,                  & !parameter 98-102 
-        1736000.0_wp, 1736000.0_wp, 0.57_wp, 0.57_wp, 0.57_wp,                 & !parameter 103-107
-        0.91_wp, 0.75_wp, 27.0_wp, 0.0_wp, 0.0_wp, 1.5_wp,                     & !parameter 108-113
-        0.86_wp, 5.0_wp, 0.0_wp,                                               & !parameter 114-116
-        299.15_wp, 293.15_wp, 0.8_wp, 0.76_wp, 5.0_wp,                         & !parameter 117-121
-        0.1_wp, 1.5_wp, 0.0_wp, 3.5_wp, 370000.0_wp, 4.5_wp,                   & !parameter 122-127
-        100000.0_wp, 0.0_wp, 3.0_wp, 10.0_wp, 3.0_wp, 0.2_wp,                  & !parameter 128-133- end of type 4
-        20000.0_wp, 20000.0_wp, 20000.0_wp, 23.0_wp, 23.0_wp, 10.0_wp,         & !parameter 0-5
-        1.0_wp, 0.005_wp, 0.01_wp, 0.31_wp, 0.43_wp, 2000000.0_wp,             & !parameter 6-11
-        103000.0_wp, 900000.0_wp, 0.35_wp, 0.38_wp, 0.04_wp, 4.0_wp,           & !parameter 12-17
-        0.01_wp, 0.001_wp, 0.55_wp,                                            & !parameter 18-20
-        0.005_wp, 0.01_wp, 0.31_wp, 0.43_wp, 2000000.0_wp,                     & !parameter 21-25
-        103000.0_wp, 900000.0_wp, 0.35_wp,                                     & !parameter 26-28                     
-        0.38_wp, 0.04_wp, 0.92_wp,                                             & !parameter 29-31       
-        27.0_wp, 0.45_wp, 0.003_wp, 0.006_wp, 0.012_wp, 0.018_wp,              & !parameter 32-37
-        1736000.0_wp, 1736000.0_wp, 1736000.0_wp,                              & !parameter 38-40
-        0.11_wp, 0.11_wp, 0.11_wp, 0.87_wp,                                    & !parameter 41-44
-        0.7_wp, 27.0_wp, 0.0_wp, 1.5_wp, 0.86_wp,                              & !parameter 45-49
-        5.0_wp, 0.001_wp, 0.0001_wp, 0.5_wp, 0.005_wp,                         & !parameter 50-54
-        0.01_wp, 0.31_wp, 0.43_wp, 2000000.0_wp,                               & !parameter 55-58
-        103000.0_wp, 900000.0_wp, 0.35_wp, 0.38_wp,                            & !parameter 59-62
-        0.04_wp, 0.92_wp, 27.0_wp, 0.5_wp,                                     & !parameter 63-66
-        0.003_wp, 0.006_wp, 0.012_wp, 0.018_wp,                                & !parameter 67-70
-        1736000.0_wp, 1736000.0_wp, 1736000.0_wp,                              & !parameter 71-73
-        0.11_wp, 0.11_wp, 0.11_wp, 0.87_wp, 0.7_wp,                            & !parameter 74-78
-        27.0_wp, 0.0_wp, 1.5_wp, 0.86_wp, 5.0_wp, 1.0_wp,                      & !parameter 79-84
-        0.005_wp, 0.01_wp, 0.31_wp, 0.43_wp, 2000000.0_wp, 103000.0_wp,        & !parameter 85-90
-        900000.0_wp, 0.35_wp, 0.38_wp, 0.04_wp, 0.91_wp, 27.0_wp, 0.0_wp,      & !parameter 91-97
-        0.003_wp, 0.006_wp, 0.012_wp, 0.018_wp, 1736000.0_wp,                  & !parameter 98-102 
-        1736000.0_wp, 1736000.0_wp, 0.11_wp, 0.11_wp, 0.11_wp,                 & !parameter 103-107
-        0.87_wp, 0.7_wp, 27.0_wp, 0.0_wp, 0.0_wp, 1.5_wp,                      & !parameter 108-113
-        0.86_wp, 5.0_wp, 0.0_wp,                                               & !parameter 114-116
-        299.15_wp, 293.15_wp, 0.8_wp, 0.6_wp, 3.0_wp,                          & !parameter 117-121
-        0.1_wp, 1.5_wp, 0.65_wp, 2.5_wp, 165000.0_wp, 4.5_wp,                  & !parameter 122-127
-        100000.0_wp, 0.0_wp, 7.0_wp, 20.0_wp, 3.0_wp, 0.2_wp,                  & !parameter 128-133- end of type 5
-        20000.0_wp, 20000.0_wp, 20000.0_wp, 23.0_wp, 23.0_wp, 10.0_wp,         & !parameter 0-5
-        1.0_wp, 0.005_wp, 0.01_wp, 0.41_wp, 0.7_wp, 2000000.0_wp,              & !parameter 6-11
-        103000.0_wp, 900000.0_wp, 0.35_wp, 0.14_wp, 0.035_wp, 4.0_wp,          & !parameter 12-17
-        0.01_wp, 0.001_wp, 0.475_wp,                                           & !parameter 18-20
-        0.005_wp, 0.01_wp, 0.41_wp, 0.7_wp, 2000000.0_wp,                      & !parameter 21-25
-        103000.0_wp, 900000.0_wp, 0.35_wp,                                     & !parameter 26-28                     
-        0.14_wp, 0.035_wp, 0.92_wp,                                            & !parameter 29-31       
-        27.0_wp, 0.525_wp, 0.003_wp, 0.006_wp, 0.012_wp, 0.018_wp,             & !parameter 32-37
-        1736000.0_wp, 1736000.0_wp, 1736000.0_wp,                              & !parameter 38-40
-        0.037_wp, 0.037_wp, 0.037_wp, 0.8_wp,                                  & !parameter 41-44
-        0.6_wp, 27.0_wp, 0.0_wp, 1.5_wp, 0.86_wp,                              & !parameter 45-49
-        5.0_wp, 0.001_wp, 0.0001_wp, 0.425_wp, 0.005_wp,                       & !parameter 50-54
-        0.01_wp, 0.41_wp, 0.7_wp, 2000000.0_wp,                                & !parameter 55-58
-        103000.0_wp, 900000.0_wp, 0.35_wp, 0.14_wp,                            & !parameter 59-62
-        0.035_wp, 0.92_wp, 27.0_wp, 0.575_wp,                                  & !parameter 63-66
-        0.003_wp, 0.006_wp, 0.012_wp, 0.018_wp,                                & !parameter 67-70
-        1736000.0_wp, 1736000.0_wp, 1736000.0_wp,                              & !parameter 71-73
-        0.037_wp, 0.037_wp, 0.037_wp, 0.8_wp, 0.6_wp,                          & !parameter 74-78
-        27.0_wp, 0.0_wp, 1.5_wp, 0.86_wp, 5.0_wp, 1.0_wp,                      & !parameter 79-84
-        0.005_wp, 0.01_wp, 0.41_wp, 0.7_wp, 2000000.0_wp, 103000.0_wp,         & !parameter 85-90
-        900000.0_wp, 0.35_wp, 0.14_wp, 0.035_wp, 0.91_wp, 27.0_wp, 0.0_wp,     & !parameter 91-97
-        0.003_wp, 0.006_wp, 0.012_wp, 0.018_wp, 1736000.0_wp,                  & !parameter 98-102 
-        1736000.0_wp, 1736000.0_wp, 0.037_wp, 0.037_wp, 0.037_wp,              & !parameter 103-107
-        0.8_wp, 0.6_wp, 27.0_wp, 0.0_wp, 0.0_wp, 1.5_wp,                       & !parameter 108-113
-        0.86_wp, 5.0_wp, 0.0_wp,                                               & !parameter 114-116
-        299.15_wp, 293.15_wp, 0.8_wp, 0.5_wp, 0.6_wp,                          & !parameter 117-121
-        0.1_wp, 1.5_wp, 0.9_wp, 2.5_wp, 80000.0_wp, 4.5_wp,                    & !parameter 122-127
-        100000.0_wp, 0.0_wp, 5.0_wp, 15.0_wp, 3.0_wp, 0.2_wp,                  & !parameter 128-133- end of type 6    
-        20000.0_wp, 20000.0_wp, 20000.0_wp, 23.0_wp, 23.0_wp, 10.0_wp,         & !parameter 0-5
-        1.0_wp, 0.29_wp, 0.295_wp, 0.695_wp, 0.985_wp, 1950400.0_wp,           & !parameter 6-11
-        1848000.0_wp, 1848000.0_wp, 0.7_wp, 1.0_wp, 1.0_wp, 4.0_wp,            & !parameter 12-17
-        0.01_wp, 0.001_wp, 1.0_wp,                                             & !parameter 18-20
-        0.29_wp, 0.295_wp, 0.695_wp, 0.985_wp, 1950400.0_wp,                   & !parameter 21-25
-        1848000.0_wp, 1848000.0_wp, 0.7_wp,                                    & !parameter 26-28                     
-        1.0_wp, 1.0_wp, 0.9_wp,                                                & !parameter 29-31       
-        27.0_wp, 0.0_wp, 0.003_wp, 0.006_wp, 0.012_wp, 0.018_wp,               & !parameter 32-37
-        1736000.0_wp, 1736000.0_wp, 1736000.0_wp,                              & !parameter 38-40
-        0.57_wp, 0.57_wp, 0.57_wp, 0.8_wp,                                     & !parameter 41-44
-        0.6_wp, 27.0_wp, 0.0_wp, 1.5_wp, 0.86_wp,                              & !parameter 45-49
-        5.0_wp, 0.001_wp, 0.0001_wp, 1.0_wp, 0.29_wp,                          & !parameter 50-54
-        0.295_wp, 0.695_wp, 0.985_wp, 1950400.0_wp,                            & !parameter 55-58
-        1848000.0_wp, 1848000.0_wp, 0.7_wp, 1.0_wp,                            & !parameter 59-62
-        1.0_wp, 0.9_wp, 27.0_wp, 0.0_wp,                                       & !parameter 63-66
-        0.003_wp, 0.006_wp, 0.012_wp, 0.018_wp,                                & !parameter 67-70
-        1736000.0_wp, 1736000.0_wp, 1736000.0_wp,                              & !parameter 71-73
-        0.57_wp, 0.57_wp, 0.57_wp, 0.8_wp, 0.6_wp,                             & !parameter 74-78
-        27.0_wp, 0.0_wp, 1.5_wp, 0.86_wp, 5.0_wp, 1.0_wp,                      & !parameter 79-84
-        0.29_wp, 0.295_wp, 0.695_wp, 0.985_wp, 1950400.0_wp, 1848000.0_wp,     & !parameter 85-90
-        1848000.0_wp, 0.7_wp, 1.0_wp, 1.0_wp, 0.9_wp, 27.0_wp, 0.0_wp,         & !parameter 91-97
-        0.003_wp, 0.006_wp, 0.012_wp, 0.018_wp, 1736000.0_wp,                  & !parameter 98-102 
-        1736000.0_wp, 1736000.0_wp, 0.57_wp, 0.57_wp, 0.57_wp,                 & !parameter 103-107
-        0.8_wp, 0.6_wp, 27.0_wp, 0.0_wp, 0.0_wp, 1.5_wp,                       & !parameter 108-113
-        0.86_wp, 5.0_wp, 0.0_wp,                                               & !parameter 114-116
-        299.15_wp, 293.15_wp, 0.8_wp, 100.0_wp, 100.0_wp,                      & !parameter 117-121
-        20.0_wp, 20.0_wp, 0.0_wp, 1.0_wp, 1.0_wp, 4.5_wp,                      & !parameter 122-127
-        100000.0_wp, 0.0_wp, 0.0_wp, 0.0_wp, 3.0_wp, 0.2_wp                    & !parameter 128-133- end of type 7 (bridge)
-                                                                       /),     &
-                                                               (/134, 7/) )
+       0.7_wp,         &  !< parameter 0   - wall fraction above ground floor level
+       0.3_wp,         &  !< parameter 1   - window fraction above ground floor level
+       0.0_wp,         &  !< parameter 2   - green fraction above ground floor level
+       0.0_wp,         &  !< parameter 3   - green fraction roof above ground floor level
+       1.5_wp,         &  !< parameter 4   - LAI roof
+       1.5_wp,         &  !< parameter 5   - LAI on wall above ground floor level
+       2200000.0_wp,   &  !< parameter 6   - heat capacity 1st/2nd wall layer above ground floor level
+       1400000.0_wp,   &  !< parameter 7   - heat capacity 3rd wall layer above ground floor level
+       1300000.0_wp,   &  !< parameter 8   - heat capacity 4th wall layer above ground floor level
+       0.35_wp,        &  !< parameter 9   - thermal conductivity 1st/2nd wall layer above ground floor level
+       0.8_wp,         &  !< parameter 10  - thermal conductivity 3rd wall layer above ground floor level
+       2.1_wp,         &  !< parameter 11  - thermal conductivity 4th wall layer above ground floor level
+       299.15_wp,      &  !< parameter 12  - indoor target summer temperature
+       293.15_wp,      &  !< parameter 13  - indoor target winter temperature
+       0.93_wp,        &  !< parameter 14  - wall emissivity above ground floor level
+       0.86_wp,        &  !< parameter 15  - green emissivity above ground floor level
+       0.91_wp,        &  !< parameter 16  - window emissivity above ground floor level
+       0.75_wp,        &  !< parameter 17  - window transmissivity above ground floor level
+       0.001_wp,       &  !< parameter 18  - z0 roughness above ground floor level
+       0.0001_wp,      &  !< parameter 19  - z0h/z0g roughness heat/humidity above ground floor level
+       4.0_wp,         &  !< parameter 20  - ground floor level height
+       0.75_wp,        &  !< parameter 21  - wall fraction ground floor level
+       0.25_wp,        &  !< parameter 22  - window fraction ground floor level
+       0.0_wp,         &  !< parameter 23  - green fraction ground floor level
+       0.0_wp,         &  !< parameter 24  - green fraction roof ground floor level
+       1.5_wp,         &  !< parameter 25  - LAI on wall ground floor level
+       2200000.0_wp,   &  !< parameter 26  - heat capacity 1st/2nd wall layer ground floor level
+       1400000.0_wp,   &  !< parameter 27  - heat capacity 3rd wall layer ground floor level
+       1300000.0_wp,   &  !< parameter 28  - heat capacity 4th wall layer ground floor level
+       0.35_wp,        &  !< parameter 29  - thermal conductivity 1st/2nd wall layer ground floor level
+       0.8_wp,         &  !< parameter 30  - thermal conductivity 3rd wall layer ground floor level
+       2.1_wp,         &  !< parameter 31  - thermal conductivity 4th wall layer ground floor level
+       0.93_wp,        &  !< parameter 32  - wall emissivity ground floor level
+       0.91_wp,        &  !< parameter 33  - window emissivity ground floor level
+       0.86_wp,        &  !< parameter 34  - green emissivity ground floor level
+       0.75_wp,        &  !< parameter 35  - window transmissivity ground floor level
+       0.01_wp,        &  !< parameter 36  - z0 roughness ground floor level
+       0.001_wp,       &  !< parameter 37  - z0h/z0q roughness heat/humidity
+       27.0_wp,        &  !< parameter 38  - wall albedo above ground floor level
+       5.0_wp,         &  !< parameter 39  - green albedo above ground floor level
+       27.0_wp,        &  !< parameter 40  - window albedo above ground floor level
+       0.005_wp,       &  !< parameter 41  - 1st wall layer thickness above ground floor level
+       0.01_wp,        &  !< parameter 42  - 2nd wall layer thickness above ground floor level
+       0.39_wp,        &  !< parameter 43  - 3rd wall layer thickness above ground floor level
+       0.63_wp,        &  !< parameter 44  - 4th wall layer thickness above ground floor level
+       20000.0_wp,     &  !< parameter 45  - heat capacity wall surface
+       23.0_wp,        &  !< parameter 46  - thermal conductivity of wall surface
+       20000.0_wp,     &  !< parameter 47  - heat capacity of window surface
+       20000.0_wp,     &  !< parameter 48  - heat capacity of green surface
+       23.0_wp,        &  !< parameter 49  - thermal conductivity of window surface
+       10.0_wp,        &  !< parameter 50  - thermal conductivty of green surface
+       1.0_wp,         &  !< parameter 51  - wall fraction ground plate
+       0.005_wp,       &  !< parameter 52  - 1st wall layer thickness ground plate
+       0.01_wp,        &  !< parameter 53  - 2nd wall layer thickness ground plate
+       0.39_wp,        &  !< parameter 54  - 3rd wall layer thickness ground plate
+       0.63_wp,        &  !< parameter 55  - 4th wall layer thickness ground plate
+       2200000.0_wp,   &  !< parameter 56  - heat capacity 1st/2nd wall layer ground plate
+       1400000.0_wp,   &  !< parameter 57  - heat capacity 3rd wall layer ground plate
+       1300000.0_wp,   &  !< parameter 58  - heat capacity 4th wall layer ground plate
+       0.35_wp,        &  !< parameter 59  - thermal conductivity 1st/2nd wall layer ground plate
+       0.8_wp,         &  !< parameter 60  - thermal conductivity 3rd wall layer ground plate
+       2.1_wp,         &  !< parameter 61  - thermal conductivity 4th wall layer ground plate
+       0.005_wp,       &  !< parameter 62  - 1st wall layer thickness ground floor level
+       0.01_wp,        &  !< parameter 63  - 2nd wall layer thickness ground floor level
+       0.39_wp,        &  !< parameter 64  - 3rd wall layer thickness ground floor level
+       0.63_wp,        &  !< parameter 65  - 4th wall layer thickness ground floor level
+       27.0_wp,        &  !< parameter 66  - wall albedo ground floor level
+       0.003_wp,       &  !< parameter 67  - 1st window layer thickness ground floor level
+       0.006_wp,       &  !< parameter 68  - 2nd window layer thickness ground floor level
+       0.012_wp,       &  !< parameter 69  - 3rd window layer thickness ground floor level
+       0.018_wp,       &  !< parameter 70  - 4th window layer thickness ground floor level
+       1736000.0_wp,   &  !< parameter 71  - heat capacity 1st/2nd window layer ground floor level
+       1736000.0_wp,   &  !< parameter 72  - heat capacity 3rd window layer ground floor level
+       1736000.0_wp,   &  !< parameter 73  - heat capacity 4th window layer ground floor level
+       0.57_wp,        &  !< parameter 74  - thermal conductivity 1st/2nd window layer ground floor level
+       0.57_wp,        &  !< parameter 75  - thermal conductivity 3rd window layer ground floor level
+       0.57_wp,        &  !< parameter 76  - thermal conductivity 4th window layer ground floor level
+       27.0_wp,        &  !< parameter 77  - window albedo ground floor level
+       5.0_wp,         &  !< parameter 78  - green albedo ground floor level
+       0.003_wp,       &  !< parameter 79  - 1st window layer thickness above ground floor level
+       0.006_wp,       &  !< parameter 80  - 2nd thickness window layer above ground floor level
+       0.012_wp,       &  !< parameter 81  - 3rd window layer thickness above ground floor level
+       0.018_wp,       &  !< parameter 82  - 4th window layer thickness above ground floor level
+       1736000.0_wp,   &  !< parameter 83  - heat capacity 1st/2nd window layer above ground floor level
+       1736000.0_wp,   &  !< parameter 84  - heat capacity 3rd window layer above ground floor level
+       1736000.0_wp,   &  !< parameter 85  - heat capacity 4th window layer above ground floor level
+       0.57_wp,        &  !< parameter 86  - thermal conductivity 1st/2nd window layer above ground floor level
+       0.57_wp,        &  !< parameter 87  - thermal conductivity 3rd window layer above ground floor level
+       0.57_wp,        &  !< parameter 88  - thermal conductivity 4th window layer above ground floor level
+       1.0_wp,         &  !< parameter 89  - wall fraction roof
+       0.005_wp,       &  !< parameter 90  - 1st wall layer thickness roof
+       0.01_wp,        &  !< parameter 91  - 2nd wall layer thickness roof
+       0.31_wp,        &  !< parameter 92  - 3rd wall layer thickness roof
+       0.63_wp,        &  !< parameter 93  - 4th wall layer thickness roof
+       2200000.0_wp,   &  !< parameter 94  - heat capacity 1st/2nd wall layer roof
+       1400000.0_wp,   &  !< parameter 95  - heat capacity 3rd wall layer roof
+       1300000.0_wp,   &  !< parameter 96  - heat capacity 4th wall layer roof
+       0.35_wp,        &  !< parameter 97  - thermal conductivity 1st/2nd wall layer roof
+       0.8_wp,         &  !< parameter 98  - thermal conductivity 3rd wall layer roof
+       2.1_wp,         &  !< parameter 99  - thermal conductivity 4th wall layer roof
+       0.93_wp,        &  !< parameter 100 - wall emissivity roof
+       27.0_wp,        &  !< parameter 101 - wall albedo roof
+       0.0_wp,         &  !< parameter 102 - window fraction roof
+       0.003_wp,       &  !< parameter 103 - window 1st layer thickness roof
+       0.006_wp,       &  !< parameter 104 - window 2nd layer thickness roof
+       0.012_wp,       &  !< parameter 105 - window 3rd layer thickness roof
+       0.018_wp,       &  !< parameter 106 - window 4th layer thickness roof
+       1736000.0_wp,   &  !< parameter 107 - heat capacity 1st/2nd window layer roof
+       1736000.0_wp,   &  !< parameter 108 - heat capacity 3rd window layer roof
+       1736000.0_wp,   &  !< parameter 109 - heat capacity 4th window layer roof
+       0.57_wp,        &  !< parameter 110 - thermal conductivity 1st/2nd window layer roof
+       0.57_wp,        &  !< parameter 111 - thermal conductivity 3rd window layer roof
+       0.57_wp,        &  !< parameter 112 - thermal conductivity 4th window layer roof
+       0.91_wp,        &  !< parameter 113 - window emissivity roof
+       0.75_wp,        &  !< parameter 114 - window transmissivity roof
+       27.0_wp,        &  !< parameter 115 - window albedo roof
+       0.86_wp,        &  !< parameter 116 - green emissivity roof
+       5.0_wp,         &  !< parameter 117 - green albedo roof
+       0.0_wp,         &  !< parameter 118 - green type roof
+       0.8_wp,         &  !< parameter 119 - shading factor
+       0.76_wp,        &  !< parameter 120 - g-value windows
+       5.0_wp,         &  !< parameter 121 - u-value windows
+       0.1_wp,         &  !< parameter 122 - basical airflow without occupancy of the room
+       0.5_wp,         &  !< parameter 123 - additional airflow depend of occupancy of the room
+       0.0_wp,         &  !< parameter 124 - heat recovery efficiency
+       3.5_wp,         &  !< parameter 125 - dynamic parameter specific effective surface
+       370000.0_wp,    &  !< parameter 126 - dynamic parameter innner heatstorage
+       4.5_wp,         &  !< parameter 127 - ratio internal surface/floor area
+       100000.0_wp,    &  !< parameter 128 - maximal heating capacity
+       0.0_wp,         &  !< parameter 129 - maximal cooling capacity
+       3.0_wp,         &  !< parameter 130 - additional internal heat gains dependent on occupancy of the room
+       10.0_wp,        &  !< parameter 131 - basic internal heat gains without occupancy of the room
+       3.0_wp,         &  !< parameter 132 - storey height
+       0.2_wp,         &  !< parameter 133 - ceiling construction height
+       0.73_wp,        &  !< parameter 0   - wall fraction above ground floor level
+       0.27_wp,        &  !< parameter 1   - window fraction above ground floor level
+       0.0_wp,         &  !< parameter 2   - green fraction above ground floor level
+       0.0_wp,         &  !< parameter 3   - green fraction roof above ground floor level
+       1.5_wp,         &  !< parameter 4   - LAI roof
+       1.5_wp,         &  !< parameter 5   - LAI on wall above ground floor level
+       2000000.0_wp,   &  !< parameter 6   - heat capacity 1st/2nd wall layer above ground floor level
+       103000.0_wp,    &  !< parameter 7   - heat capacity 3rd wall layer above ground floor level
+       900000.0_wp,    &  !< parameter 8   - heat capacity 4th wall layer above ground floor level
+       0.35_wp,        &  !< parameter 9   - thermal conductivity 1st/2nd wall layer above ground floor level
+       0.38_wp,        &  !< parameter 10  - thermal conductivity 3rd wall layer above ground floor level
+       0.04_wp,        &  !< parameter 11  - thermal conductivity 4th wall layer above ground floor level
+       299.15_wp,      &  !< parameter 12  - indoor target summer temperature
+       293.15_wp,      &  !< parameter 13  - indoor target winter temperature
+       0.92_wp,        &  !< parameter 14  - wall emissivity above ground floor level
+       0.86_wp,        &  !< parameter 15  - green emissivity above ground floor level
+       0.87_wp,        &  !< parameter 16  - window emissivity above ground floor level
+       0.7_wp,         &  !< parameter 17  - window transmissivity above ground floor level
+       0.001_wp,       &  !< parameter 18  - z0 roughness above ground floor level
+       0.0001_wp,      &  !< parameter 19  - z0h/z0g roughness heat/humidity above ground floor level
+       4.0_wp,         &  !< parameter 20  - ground floor level height
+       0.78_wp,        &  !< parameter 21  - wall fraction ground floor level
+       0.22_wp,        &  !< parameter 22  - window fraction ground floor level
+       0.0_wp,         &  !< parameter 23  - green fraction ground floor level
+       0.0_wp,         &  !< parameter 24  - green fraction roof ground floor level
+       1.5_wp,         &  !< parameter 25  - LAI on wall ground floor level
+       2000000.0_wp,   &  !< parameter 26  - heat capacity 1st/2nd wall layer ground floor level
+       103000.0_wp,    &  !< parameter 27  - heat capacity 3rd wall layer ground floor level
+       900000.0_wp,    &  !< parameter 28  - heat capacity 4th wall layer ground floor level
+       0.35_wp,        &  !< parameter 29  - thermal conductivity 1st/2nd wall layer ground floor level
+       0.38_wp,        &  !< parameter 30  - thermal conductivity 3rd wall layer ground floor level
+       0.04_wp,        &  !< parameter 31  - thermal conductivity 4th wall layer ground floor level
+       0.92_wp,        &  !< parameter 32  - wall emissivity ground floor level
+       0.11_wp,        &  !< parameter 33  - window emissivity ground floor level
+       0.86_wp,        &  !< parameter 34  - green emissivity ground floor level
+       0.7_wp,         &  !< parameter 35  - window transmissivity ground floor level
+       0.01_wp,        &  !< parameter 36  - z0 roughness ground floor level
+       0.001_wp,       &  !< parameter 37  - z0h/z0q roughness heat/humidity
+       27.0_wp,        &  !< parameter 38  - wall albedo above ground floor level
+       5.0_wp,         &  !< parameter 39  - green albedo above ground floor level
+       27.0_wp,        &  !< parameter 40  - window albedo above ground floor level
+       0.005_wp,       &  !< parameter 41  - 1st wall layer thickness above ground floor level
+       0.01_wp,        &  !< parameter 42  - 2nd wall layer thickness above ground floor level
+       0.31_wp,        &  !< parameter 43  - 3rd wall layer thickness above ground floor level
+       0.43_wp,        &  !< parameter 44  - 4th wall layer thickness above ground floor level
+       20000.0_wp,     &  !< parameter 45  - heat capacity wall surface
+       23.0_wp,        &  !< parameter 46  - thermal conductivity of wall surface
+       20000.0_wp,     &  !< parameter 47  - heat capacity of window surface
+       20000.0_wp,     &  !< parameter 48  - heat capacity of green surface
+       23.0_wp,        &  !< parameter 49  - thermal conductivity of window surface
+       10.0_wp,        &  !< parameter 50  - thermal conductivty of green surface
+       1.0_wp,         &  !< parameter 51  - wall fraction ground plate
+       0.005_wp,       &  !< parameter 52  - 1st wall layer thickness ground plate
+       0.01_wp,        &  !< parameter 53  - 2nd wall layer thickness ground plate
+       0.31_wp,        &  !< parameter 54  - 3rd wall layer thickness ground plate
+       0.42_wp,        &  !< parameter 55  - 4th wall layer thickness ground plate
+       2000000.0_wp,   &  !< parameter 56  - heat capacity 1st/2nd wall layer ground plate
+       103000.0_wp,    &  !< parameter 57  - heat capacity 3rd wall layer ground plate
+       900000.0_wp,    &  !< parameter 58  - heat capacity 4th wall layer ground plate
+       0.35_wp,        &  !< parameter 59  - thermal conductivity 1st/2nd wall layer ground plate
+       0.38_wp,        &  !< parameter 60  - thermal conductivity 3rd wall layer ground plate
+       0.04_wp,        &  !< parameter 61  - thermal conductivity 4th wall layer ground plate
+       0.005_wp,       &  !< parameter 62  - 1st wall layer thickness ground floor level
+       0.01_wp,        &  !< parameter 63  - 2nd wall layer thickness ground floor level
+       0.31_wp,        &  !< parameter 64  - 3rd wall layer thickness ground floor level
+       0.43_wp,        &  !< parameter 65  - 4th wall layer thickness ground floor level
+       27.0_wp,        &  !< parameter 66  - wall albedo ground floor level
+       0.003_wp,       &  !< parameter 67  - 1st window layer thickness ground floor level
+       0.006_wp,       &  !< parameter 68  - 2nd window layer thickness ground floor level
+       0.012_wp,       &  !< parameter 69  - 3rd window layer thickness ground floor level
+       0.018_wp,       &  !< parameter 70  - 4th window layer thickness ground floor level
+       1736000.0_wp,   &  !< parameter 71  - heat capacity 1st/2nd window layer ground floor level
+       1736000.0_wp,   &  !< parameter 72  - heat capacity 3rd window layer ground floor level
+       1736000.0_wp,   &  !< parameter 73  - heat capacity 4th window layer ground floor level
+       0.11_wp,        &  !< parameter 74  - thermal conductivity 1st/2nd window layer ground floor level
+       0.11_wp,        &  !< parameter 75  - thermal conductivity 3rd window layer ground floor level
+       0.11_wp,        &  !< parameter 76  - thermal conductivity 4th window layer ground floor level
+       27.0_wp,        &  !< parameter 77  - window albedo ground floor level
+       5.0_wp,         &  !< parameter 78  - green albedo ground floor level
+       0.003_wp,       &  !< parameter 79  - 1st window layer thickness above ground floor level
+       0.006_wp,       &  !< parameter 80  - 2nd thickness window layer above ground floor level
+       0.012_wp,       &  !< parameter 81  - 3rd window layer thickness above ground floor level
+       0.018_wp,       &  !< parameter 82  - 4th window layer thickness above ground floor level
+       1736000.0_wp,   &  !< parameter 83  - heat capacity 1st/2nd window layer above ground floor level
+       1736000.0_wp,   &  !< parameter 84  - heat capacity 3rd window layer above ground floor level
+       1736000.0_wp,   &  !< parameter 85  - heat capacity 4th window layer above ground floor level
+       0.11_wp,        &  !< parameter 86  - thermal conductivity 1st/2nd window layer above ground floor level
+       0.11_wp,        &  !< parameter 87  - thermal conductivity 3rd window layer above ground floor level
+       0.11_wp,        &  !< parameter 88  - thermal conductivity 4th window layer above ground floor level
+       1.0_wp,         &  !< parameter 89  - wall fraction roof
+       0.005_wp,       &  !< parameter 90  - 1st wall layer thickness roof
+       0.01_wp,        &  !< parameter 91  - 2nd wall layer thickness roof
+       0.5_wp,         &  !< parameter 92  - 3rd wall layer thickness roof
+       0.79_wp,        &  !< parameter 93  - 4th wall layer thickness roof
+       2000000.0_wp,   &  !< parameter 94  - heat capacity 1st/2nd wall layer roof
+       103000.0_wp,    &  !< parameter 95  - heat capacity 3rd wall layer roof
+       900000.0_wp,    &  !< parameter 96  - heat capacity 4th wall layer roof
+       0.35_wp,        &  !< parameter 97  - thermal conductivity 1st/2nd wall layer roof
+       0.38_wp,        &  !< parameter 98  - thermal conductivity 3rd wall layer roof
+       0.04_wp,        &  !< parameter 99  - thermal conductivity 4th wall layer roof
+       0.93_wp,        &  !< parameter 100 - wall emissivity roof
+       27.0_wp,        &  !< parameter 101 - wall albedo roof
+       0.0_wp,         &  !< parameter 102 - window fraction roof
+       0.003_wp,       &  !< parameter 103 - window 1st layer thickness roof
+       0.006_wp,       &  !< parameter 104 - window 2nd layer thickness roof
+       0.012_wp,       &  !< parameter 105 - window 3rd layer thickness roof
+       0.018_wp,       &  !< parameter 106 - window 4th layer thickness roof
+       1736000.0_wp,   &  !< parameter 107 - heat capacity 1st/2nd window layer roof
+       1736000.0_wp,   &  !< parameter 108 - heat capacity 3rd window layer roof
+       1736000.0_wp,   &  !< parameter 109 - heat capacity 4th window layer roof
+       0.11_wp,        &  !< parameter 110 - thermal conductivity 1st/2nd window layer roof
+       0.11_wp,        &  !< parameter 111 - thermal conductivity 3rd window layer roof
+       0.11_wp,        &  !< parameter 112 - thermal conductivity 4th window layer roof
+       0.87_wp,        &  !< parameter 113 - window emissivity roof
+       0.7_wp,         &  !< parameter 114 - window transmissivity roof
+       27.0_wp,        &  !< parameter 115 - window albedo roof
+       0.86_wp,        &  !< parameter 116 - green emissivity roof
+       5.0_wp,         &  !< parameter 117 - green albedo roof
+       0.0_wp,         &  !< parameter 118 - green type roof
+       0.8_wp,         &  !< parameter 119 - shading factor
+       0.6_wp,         &  !< parameter 120 - g-value windows
+       3.0_wp,         &  !< parameter 121 - u-value windows
+       0.1_wp,         &  !< parameter 122 - basical airflow without occupancy of the room
+       0.5_wp,         &  !< parameter 123 - additional airflow depend of occupancy of the room
+       0.0_wp,         &  !< parameter 124 - heat recovery efficiency
+       2.5_wp,         &  !< parameter 125 - dynamic parameter specific effective surface
+       165000.0_wp,    &  !< parameter 126 - dynamic parameter innner heatstorage
+       4.5_wp,         &  !< parameter 127 - ratio internal surface/floor area
+       100000.0_wp,    &  !< parameter 128 - maximal heating capacity
+       0.0_wp,         &  !< parameter 129 - maximal cooling capacity
+       4.0_wp,         &  !< parameter 130 - additional internal heat gains dependent on occupancy of the room
+       8.0_wp,         &  !< parameter 131 - basic internal heat gains without occupancy of the room
+       3.0_wp,         &  !< parameter 132 - storey height
+       0.2_wp,         &  !< parameter 133 - ceiling construction height
+       0.7_wp,         &  !< parameter 0   - wall fraction above ground floor level
+       0.3_wp,         &  !< parameter 1   - window fraction above ground floor level
+       0.0_wp,         &  !< parameter 2   - green fraction above ground floor level
+       0.0_wp,         &  !< parameter 3   - green fraction roof above ground floor level
+       1.5_wp,         &  !< parameter 4   - LAI roof
+       1.5_wp,         &  !< parameter 5   - LAI on wall above ground floor level
+       2000000.0_wp,   &  !< parameter 6   - heat capacity 1st/2nd wall layer above ground floor level
+       103000.0_wp,    &  !< parameter 7   - heat capacity 3rd wall layer above ground floor level
+       900000.0_wp,    &  !< parameter 8   - heat capacity 4th wall layer above ground floor level
+       0.35_wp,        &  !< parameter 9   - thermal conductivity 1st/2nd wall layer above ground floor level
+       0.14_wp,        &  !< parameter 10  - thermal conductivity 3rd wall layer above ground floor level
+       0.035_wp,       &  !< parameter 11  - thermal conductivity 4th wall layer above ground floor level
+       299.15_wp,      &  !< parameter 12  - indoor target summer temperature
+       293.15_wp,      &  !< parameter 13  - indoor target winter temperature
+       0.92_wp,        &  !< parameter 14  - wall emissivity above ground floor level
+       0.86_wp,        &  !< parameter 15  - green emissivity above ground floor level
+       0.8_wp,         &  !< parameter 16  - window emissivity above ground floor level
+       0.6_wp,         &  !< parameter 17  - window transmissivity above ground floor level
+       0.001_wp,       &  !< parameter 18  - z0 roughness above ground floor level
+       0.0001_wp,      &  !< parameter 19  - z0h/z0g roughness heat/humidity above ground floor level
+       4.0_wp,         &  !< parameter 20  - ground floor level height
+       0.75_wp,        &  !< parameter 21  - wall fraction ground floor level
+       0.25_wp,        &  !< parameter 22  - window fraction ground floor level
+       0.0_wp,         &  !< parameter 23  - green fraction ground floor level
+       0.0_wp,         &  !< parameter 24  - green fraction roof ground floor level
+       1.5_wp,         &  !< parameter 25  - LAI on wall ground floor level
+       2000000.0_wp,   &  !< parameter 26  - heat capacity 1st/2nd wall layer ground floor level
+       103000.0_wp,    &  !< parameter 27  - heat capacity 3rd wall layer ground floor level
+       900000.0_wp,    &  !< parameter 28  - heat capacity 4th wall layer ground floor level
+       0.35_wp,        &  !< parameter 29  - thermal conductivity 1st/2nd wall layer ground floor level
+       0.14_wp,        &  !< parameter 30  - thermal conductivity 3rd wall layer ground floor level
+       0.035_wp,       &  !< parameter 31  - thermal conductivity 4th wall layer ground floor level
+       0.92_wp,        &  !< parameter 32  - wall emissivity ground floor level
+       0.8_wp,         &  !< parameter 33  - window emissivity ground floor level
+       0.86_wp,        &  !< parameter 34  - green emissivity ground floor level
+       0.6_wp,         &  !< parameter 35  - window transmissivity ground floor level
+       0.01_wp,        &  !< parameter 36  - z0 roughness ground floor level
+       0.001_wp,       &  !< parameter 37  - z0h/z0q roughness heat/humidity
+       27.0_wp,        &  !< parameter 38  - wall albedo above ground floor level
+       5.0_wp,         &  !< parameter 39  - green albedo above ground floor level
+       27.0_wp,        &  !< parameter 40  - window albedo above ground floor level
+       0.005_wp,       &  !< parameter 41  - 1st wall layer thickness above ground floor level
+       0.01_wp,        &  !< parameter 42  - 2nd wall layer thickness above ground floor level
+       0.41_wp,        &  !< parameter 43  - 3rd wall layer thickness above ground floor level
+       0.7_wp,         &  !< parameter 44  - 4th wall layer thickness above ground floor level
+       20000.0_wp,     &  !< parameter 45  - heat capacity wall surface
+       23.0_wp,        &  !< parameter 46  - thermal conductivity of wall surface
+       20000.0_wp,     &  !< parameter 47  - heat capacity of window surface
+       20000.0_wp,     &  !< parameter 48  - heat capacity of green surface
+       23.0_wp,        &  !< parameter 49  - thermal conductivity of window surface
+       10.0_wp,        &  !< parameter 50  - thermal conductivty of green surface
+       1.0_wp,         &  !< parameter 51  - wall fraction ground plate
+       0.005_wp,       &  !< parameter 52  - 1st wall layer thickness ground plate
+       0.01_wp,        &  !< parameter 53  - 2nd wall layer thickness ground plate
+       0.41_wp,        &  !< parameter 54  - 3rd wall layer thickness ground plate
+       0.7_wp,         &  !< parameter 55  - 4th wall layer thickness ground plate
+       2000000.0_wp,   &  !< parameter 56  - heat capacity 1st/2nd wall layer ground plate
+       103000.0_wp,    &  !< parameter 57  - heat capacity 3rd wall layer ground plate
+       900000.0_wp,    &  !< parameter 58  - heat capacity 4th wall layer ground plate
+       0.35_wp,        &  !< parameter 59  - thermal conductivity 1st/2nd wall layer ground plate
+       0.14_wp,        &  !< parameter 60  - thermal conductivity 3rd wall layer ground plate
+       0.035_wp,       &  !< parameter 61  - thermal conductivity 4th wall layer ground plate
+       0.005_wp,       &  !< parameter 62  - 1st wall layer thickness ground floor level
+       0.01_wp,        &  !< parameter 63  - 2nd wall layer thickness ground floor level
+       0.41_wp,        &  !< parameter 64  - 3rd wall layer thickness ground floor level
+       0.7_wp,         &  !< parameter 65  - 4th wall layer thickness ground floor level
+       27.0_wp,        &  !< parameter 66  - wall albedo ground floor level
+       0.003_wp,       &  !< parameter 67  - 1st window layer thickness ground floor level
+       0.006_wp,       &  !< parameter 68  - 2nd window layer thickness ground floor level
+       0.012_wp,       &  !< parameter 69  - 3rd window layer thickness ground floor level
+       0.018_wp,       &  !< parameter 70  - 4th window layer thickness ground floor level
+       1736000.0_wp,   &  !< parameter 71  - heat capacity 1st/2nd window layer ground floor level
+       1736000.0_wp,   &  !< parameter 72  - heat capacity 3rd window layer ground floor level
+       1736000.0_wp,   &  !< parameter 73  - heat capacity 4th window layer ground floor level
+       0.037_wp,       &  !< parameter 74  - thermal conductivity 1st/2nd window layer ground floor level
+       0.037_wp,       &  !< parameter 75  - thermal conductivity 3rd window layer ground floor level
+       0.037_wp,       &  !< parameter 76  - thermal conductivity 4th window layer ground floor level
+       27.0_wp,        &  !< parameter 77  - window albedo ground floor level
+       5.0_wp,         &  !< parameter 78  - green albedo ground floor level
+       0.003_wp,       &  !< parameter 79  - 1st window layer thickness above ground floor level
+       0.006_wp,       &  !< parameter 80  - 2nd thickness window layer above ground floor level
+       0.012_wp,       &  !< parameter 81  - 3rd window layer thickness above ground floor level
+       0.018_wp,       &  !< parameter 82  - 4th window layer thickness above ground floor level
+       1736000.0_wp,   &  !< parameter 83  - heat capacity 1st/2nd window layer above ground floor level
+       1736000.0_wp,   &  !< parameter 84  - heat capacity 3rd window layer above ground floor level
+       1736000.0_wp,   &  !< parameter 85  - heat capacity 4th window layer above ground floor level
+       0.037_wp,       &  !< parameter 86  - thermal conductivity 1st/2nd window layer above ground floor level
+       0.037_wp,       &  !< parameter 87  - thermal conductivity 3rd window layer above ground floor level
+       0.037_wp,       &  !< parameter 88  - thermal conductivity 4th window layer above ground floor level
+       1.0_wp,         &  !< parameter 89  - wall fraction roof
+       0.005_wp,       &  !< parameter 90  - 1st wall layer thickness roof
+       0.01_wp,        &  !< parameter 91  - 2nd wall layer thickness roof
+       0.41_wp,        &  !< parameter 92  - 3rd wall layer thickness roof
+       0.7_wp,         &  !< parameter 93  - 4th wall layer thickness roof
+       2000000.0_wp,   &  !< parameter 94  - heat capacity 1st/2nd wall layer roof
+       103000.0_wp,    &  !< parameter 95  - heat capacity 3rd wall layer roof
+       900000.0_wp,    &  !< parameter 96  - heat capacity 4th wall layer roof
+       0.35_wp,        &  !< parameter 97  - thermal conductivity 1st/2nd wall layer roof
+       0.14_wp,        &  !< parameter 98  - thermal conductivity 3rd wall layer roof
+       0.035_wp,       &  !< parameter 99  - thermal conductivity 4th wall layer roof
+       0.93_wp,        &  !< parameter 100 - wall emissivity roof
+       27.0_wp,        &  !< parameter 101 - wall albedo roof
+       0.0_wp,         &  !< parameter 102 - window fraction roof
+       0.003_wp,       &  !< parameter 103 - window 1st layer thickness roof
+       0.006_wp,       &  !< parameter 104 - window 2nd layer thickness roof
+       0.012_wp,       &  !< parameter 105 - window 3rd layer thickness roof
+       0.018_wp,       &  !< parameter 106 - window 4th layer thickness roof
+       1736000.0_wp,   &  !< parameter 107 - heat capacity 1st/2nd window layer roof
+       1736000.0_wp,   &  !< parameter 108 - heat capacity 3rd window layer roof
+       1736000.0_wp,   &  !< parameter 109 - heat capacity 4th window layer roof
+       0.037_wp,       &  !< parameter 110 - thermal conductivity 1st/2nd window layer roof
+       0.037_wp,       &  !< parameter 111 - thermal conductivity 3rd window layer roof
+       0.037_wp,       &  !< parameter 112 - thermal conductivity 4th window layer roof
+       0.8_wp,         &  !< parameter 113 - window emissivity roof
+       0.6_wp,         &  !< parameter 114 - window transmissivity roof
+       27.0_wp,        &  !< parameter 115 - window albedo roof
+       0.86_wp,        &  !< parameter 116 - green emissivity roof
+       5.0_wp,         &  !< parameter 117 - green albedo roof
+       0.0_wp,         &  !< parameter 118 - green type roof
+       0.8_wp,         &  !< parameter 119 - shading factor
+       0.5_wp,         &  !< parameter 120 - g-value windows
+       0.6_wp,         &  !< parameter 121 - u-value windows
+       0.1_wp,         &  !< parameter 122 - basical airflow without occupancy of the room
+       0.5_wp,         &  !< parameter 123 - additional airflow depend of occupancy of the room
+       0.8_wp,         &  !< parameter 124 - heat recovery efficiency
+       2.5_wp,         &  !< parameter 125 - dynamic parameter specific effective surface
+       80000.0_wp,     &  !< parameter 126 - dynamic parameter innner heatstorage
+       4.5_wp,         &  !< parameter 127 - ratio internal surface/floor area
+       100000.0_wp,    &  !< parameter 128 - maximal heating capacity
+       0.0_wp,         &  !< parameter 129 - maximal cooling capacity
+       3.0_wp,         &  !< parameter 130 - additional internal heat gains dependent on occupancy of the room
+       8.0_wp,         &  !< parameter 131 - basic internal heat gains without occupancy of the room
+       3.0_wp,         &  !< parameter 132 - storey height
+       0.2_wp,         &  !< parameter 133 - ceiling construction height
+       0.5_wp,         &  !< parameter 0   - wall fraction above ground floor level
+       0.5_wp,         &  !< parameter 1   - window fraction above ground floor level
+       0.0_wp,         &  !< parameter 2   - green fraction above ground floor level
+       0.0_wp,         &  !< parameter 3   - green fraction roof above ground floor level
+       1.5_wp,         &  !< parameter 4   - LAI roof
+       1.5_wp,         &  !< parameter 5   - LAI on wall above ground floor level
+       2200000.0_wp,   &  !< parameter 6   - heat capacity 1st/2nd wall layer above ground floor level
+       1400000.0_wp,   &  !< parameter 7   - heat capacity 3rd wall layer above ground floor level
+       1300000.0_wp,   &  !< parameter 8   - heat capacity 4th wall layer above ground floor level
+       0.35_wp,        &  !< parameter 9   - thermal conductivity 1st/2nd wall layer above ground floor level
+       0.8_wp,         &  !< parameter 10  - thermal conductivity 3rd wall layer above ground floor level
+       2.1_wp,         &  !< parameter 11  - thermal conductivity 4th wall layer above ground floor level
+       299.15_wp,      &  !< parameter 12  - indoor target summer temperature
+       293.15_wp,      &  !< parameter 13  - indoor target winter temperature
+       0.93_wp,        &  !< parameter 14  - wall emissivity above ground floor level
+       0.86_wp,        &  !< parameter 15  - green emissivity above ground floor level
+       0.91_wp,        &  !< parameter 16  - window emissivity above ground floor level
+       0.75_wp,        &  !< parameter 17  - window transmissivity above ground floor level
+       0.001_wp,       &  !< parameter 18  - z0 roughness above ground floor level
+       0.0001_wp,      &  !< parameter 19  - z0h/z0g roughness heat/humidity above ground floor level
+       4.0_wp,         &  !< parameter 20  - ground floor level height
+       0.55_wp,        &  !< parameter 21  - wall fraction ground floor level
+       0.45_wp,        &  !< parameter 22  - window fraction ground floor level
+       0.0_wp,         &  !< parameter 23  - green fraction ground floor level
+       0.0_wp,         &  !< parameter 24  - green fraction roof ground floor level
+       1.5_wp,         &  !< parameter 25  - LAI on wall ground floor level
+       2200000.0_wp,   &  !< parameter 26  - heat capacity 1st/2nd wall layer ground floor level
+       1400000.0_wp,   &  !< parameter 27  - heat capacity 3rd wall layer ground floor level
+       1300000.0_wp,   &  !< parameter 28  - heat capacity 4th wall layer ground floor level
+       0.35_wp,        &  !< parameter 29  - thermal conductivity 1st/2nd wall layer ground floor level
+       0.8_wp,         &  !< parameter 30  - thermal conductivity 3rd wall layer ground floor level
+       2.1_wp,         &  !< parameter 31  - thermal conductivity 4th wall layer ground floor level
+       0.93_wp,        &  !< parameter 32  - wall emissivity ground floor level
+       0.91_wp,        &  !< parameter 33  - window emissivity ground floor level
+       0.86_wp,        &  !< parameter 34  - green emissivity ground floor level
+       0.75_wp,        &  !< parameter 35  - window transmissivity ground floor level
+       0.01_wp,        &  !< parameter 36  - z0 roughness ground floor level
+       0.001_wp,       &  !< parameter 37  - z0h/z0q roughness heat/humidity
+       27.0_wp,        &  !< parameter 38  - wall albedo above ground floor level
+       5.0_wp,         &  !< parameter 39  - green albedo above ground floor level
+       27.0_wp,        &  !< parameter 40  - window albedo above ground floor level
+       0.005_wp,       &  !< parameter 41  - 1st wall layer thickness above ground floor level
+       0.01_wp,        &  !< parameter 42  - 2nd wall layer thickness above ground floor level
+       0.39_wp,        &  !< parameter 43  - 3rd wall layer thickness above ground floor level
+       0.63_wp,        &  !< parameter 44  - 4th wall layer thickness above ground floor level
+       20000.0_wp,     &  !< parameter 45  - heat capacity wall surface
+       23.0_wp,        &  !< parameter 46  - thermal conductivity of wall surface
+       20000.0_wp,     &  !< parameter 47  - heat capacity of window surface
+       20000.0_wp,     &  !< parameter 48  - heat capacity of green surface
+       23.0_wp,        &  !< parameter 49  - thermal conductivity of window surface
+       10.0_wp,        &  !< parameter 50  - thermal conductivty of green surface
+       1.0_wp,         &  !< parameter 51  - wall fraction ground plate
+       0.005_wp,       &  !< parameter 52  - 1st wall layer thickness ground plate
+       0.01_wp,        &  !< parameter 53  - 2nd wall layer thickness ground plate
+       0.39_wp,        &  !< parameter 54  - 3rd wall layer thickness ground plate
+       0.63_wp,        &  !< parameter 55  - 4th wall layer thickness ground plate
+       2200000.0_wp,   &  !< parameter 56  - heat capacity 1st/2nd wall layer ground plate
+       1400000.0_wp,   &  !< parameter 57  - heat capacity 3rd wall layer ground plate
+       1300000.0_wp,   &  !< parameter 58  - heat capacity 4th wall layer ground plate
+       0.35_wp,        &  !< parameter 59  - thermal conductivity 1st/2nd wall layer ground plate
+       0.8_wp,         &  !< parameter 60  - thermal conductivity 3rd wall layer ground plate
+       2.1_wp,         &  !< parameter 61  - thermal conductivity 4th wall layer ground plate
+       0.005_wp,       &  !< parameter 62  - 1st wall layer thickness ground floor level
+       0.01_wp,        &  !< parameter 63  - 2nd wall layer thickness ground floor level
+       0.39_wp,        &  !< parameter 64  - 3rd wall layer thickness ground floor level
+       0.63_wp,        &  !< parameter 65  - 4th wall layer thickness ground floor level
+       27.0_wp,        &  !< parameter 66  - wall albedo ground floor level
+       0.003_wp,       &  !< parameter 67  - 1st window layer thickness ground floor level
+       0.006_wp,       &  !< parameter 68  - 2nd window layer thickness ground floor level
+       0.012_wp,       &  !< parameter 69  - 3rd window layer thickness ground floor level
+       0.018_wp,       &  !< parameter 70  - 4th window layer thickness ground floor level
+       1736000.0_wp,   &  !< parameter 71  - heat capacity 1st/2nd window layer ground floor level
+       1736000.0_wp,   &  !< parameter 72  - heat capacity 3rd window layer ground floor level
+       1736000.0_wp,   &  !< parameter 73  - heat capacity 4th window layer ground floor level
+       0.57_wp,        &  !< parameter 74  - thermal conductivity 1st/2nd window layer ground floor level
+       0.57_wp,        &  !< parameter 75  - thermal conductivity 3rd window layer ground floor level
+       0.57_wp,        &  !< parameter 76  - thermal conductivity 4th window layer ground floor level
+       27.0_wp,        &  !< parameter 77  - window albedo ground floor level
+       5.0_wp,         &  !< parameter 78  - green albedo ground floor level
+       0.003_wp,       &  !< parameter 79  - 1st window layer thickness above ground floor level
+       0.006_wp,       &  !< parameter 80  - 2nd thickness window layer above ground floor level
+       0.012_wp,       &  !< parameter 81  - 3rd window layer thickness above ground floor level
+       0.018_wp,       &  !< parameter 82  - 4th window layer thickness above ground floor level
+       1736000.0_wp,   &  !< parameter 83  - heat capacity 1st/2nd window layer above ground floor level
+       1736000.0_wp,   &  !< parameter 84  - heat capacity 3rd window layer above ground floor level
+       1736000.0_wp,   &  !< parameter 85  - heat capacity 4th window layer above ground floor level
+       0.57_wp,        &  !< parameter 86  - thermal conductivity 1st/2nd window layer above ground floor level
+       0.57_wp,        &  !< parameter 87  - thermal conductivity 3rd window layer above ground floor level
+       0.57_wp,        &  !< parameter 88  - thermal conductivity 4th window layer above ground floor level
+       1.0_wp,         &  !< parameter 89  - wall fraction roof
+       0.005_wp,       &  !< parameter 90  - 1st wall layer thickness roof
+       0.01_wp,        &  !< parameter 91  - 2nd wall layer thickness roof
+       0.39_wp,        &  !< parameter 92  - 3rd wall layer thickness roof
+       0.63_wp,        &  !< parameter 93  - 4th wall layer thickness roof
+       2200000.0_wp,   &  !< parameter 94  - heat capacity 1st/2nd wall layer roof
+       1400000.0_wp,   &  !< parameter 95  - heat capacity 3rd wall layer roof
+       1300000.0_wp,   &  !< parameter 96  - heat capacity 4th wall layer roof
+       0.35_wp,        &  !< parameter 97  - thermal conductivity 1st/2nd wall layer roof
+       0.8_wp,         &  !< parameter 98  - thermal conductivity 3rd wall layer roof
+       2.1_wp,         &  !< parameter 99  - thermal conductivity 4th wall layer roof
+       0.93_wp,        &  !< parameter 100 - wall emissivity roof
+       27.0_wp,        &  !< parameter 101 - wall albedo roof
+       0.0_wp,         &  !< parameter 102 - window fraction roof
+       0.003_wp,       &  !< parameter 103 - window 1st layer thickness roof
+       0.006_wp,       &  !< parameter 104 - window 2nd layer thickness roof
+       0.012_wp,       &  !< parameter 105 - window 3rd layer thickness roof
+       0.018_wp,       &  !< parameter 106 - window 4th layer thickness roof
+       1736000.0_wp,   &  !< parameter 107 - heat capacity 1st/2nd window layer roof
+       1736000.0_wp,   &  !< parameter 108 - heat capacity 3rd window layer roof
+       1736000.0_wp,   &  !< parameter 109 - heat capacity 4th window layer roof
+       0.57_wp,        &  !< parameter 110 - thermal conductivity 1st/2nd window layer roof
+       0.57_wp,        &  !< parameter 111 - thermal conductivity 3rd window layer roof
+       0.57_wp,        &  !< parameter 112 - thermal conductivity 4th window layer roof
+       0.91_wp,        &  !< parameter 113 - window emissivity roof
+       0.75_wp,        &  !< parameter 114 - window transmissivity roof
+       27.0_wp,        &  !< parameter 115 - window albedo roof
+       0.86_wp,        &  !< parameter 116 - green emissivity roof
+       5.0_wp,         &  !< parameter 117 - green albedo roof
+       0.0_wp,         &  !< parameter 118 - green type roof
+       0.8_wp,         &  !< parameter 119 - shading factor
+       0.76_wp,        &  !< parameter 120 - g-value windows
+       5.0_wp,         &  !< parameter 121 - u-value windows
+       0.1_wp,         &  !< parameter 122 - basical airflow without occupancy of the room
+       1.5_wp,         &  !< parameter 123 - additional airflow depend of occupancy of the room
+       0.0_wp,         &  !< parameter 124 - heat recovery efficiency
+       3.5_wp,         &  !< parameter 125 - dynamic parameter specific effective surface
+       370000.0_wp,    &  !< parameter 126 - dynamic parameter innner heatstorage
+       4.5_wp,         &  !< parameter 127 - ratio internal surface/floor area
+       100000.0_wp,    &  !< parameter 128 - maximal heating capacity
+       0.0_wp,         &  !< parameter 129 - maximal cooling capacity
+       3.0_wp,         &  !< parameter 130 - additional internal heat gains dependent on occupancy of the room
+       10.0_wp,        &  !< parameter 131 - basic internal heat gains without occupancy of the room
+       3.0_wp,         &  !< parameter 132 - storey height
+       0.2_wp,         &  !< parameter 133 - ceiling construction height
+       0.5_wp,         &  !< parameter 0   - wall fraction above ground floor level
+       0.5_wp,         &  !< parameter 1   - window fraction above ground floor level
+       0.0_wp,         &  !< parameter 2   - green fraction above ground floor level
+       0.0_wp,         &  !< parameter 3   - green fraction roof above ground floor level
+       1.5_wp,         &  !< parameter 4   - LAI roof
+       1.5_wp,         &  !< parameter 5   - LAI on wall above ground floor level
+       2000000.0_wp,   &  !< parameter 6   - heat capacity 1st/2nd wall layer above ground floor level
+       103000.0_wp,    &  !< parameter 7   - heat capacity 3rd wall layer above ground floor level
+       900000.0_wp,    &  !< parameter 8   - heat capacity 4th wall layer above ground floor level
+       0.35_wp,        &  !< parameter 9   - thermal conductivity 1st/2nd wall layer above ground floor level
+       0.38_wp,        &  !< parameter 10  - thermal conductivity 3rd wall layer above ground floor level
+       0.04_wp,        &  !< parameter 11  - thermal conductivity 4th wall layer above ground floor level
+       299.15_wp,      &  !< parameter 12  - indoor target summer temperature
+       293.15_wp,      &  !< parameter 13  - indoor target winter temperature
+       0.92_wp,        &  !< parameter 14  - wall emissivity above ground floor level
+       0.86_wp,        &  !< parameter 15  - green emissivity above ground floor level
+       0.87_wp,        &  !< parameter 16  - window emissivity above ground floor level
+       0.7_wp,         &  !< parameter 17  - window transmissivity above ground floor level
+       0.001_wp,       &  !< parameter 18  - z0 roughness above ground floor level
+       0.0001_wp,      &  !< parameter 19  - z0h/z0g roughness heat/humidity above ground floor level
+       4.0_wp,         &  !< parameter 20  - ground floor level height
+       0.55_wp,        &  !< parameter 21  - wall fraction ground floor level
+       0.45_wp,        &  !< parameter 22  - window fraction ground floor level
+       0.0_wp,         &  !< parameter 23  - green fraction ground floor level
+       0.0_wp,         &  !< parameter 24  - green fraction roof ground floor level
+       1.5_wp,         &  !< parameter 25  - LAI on wall ground floor level
+       2000000.0_wp,   &  !< parameter 26  - heat capacity 1st/2nd wall layer ground floor level
+       103000.0_wp,    &  !< parameter 27  - heat capacity 3rd wall layer ground floor level
+       900000.0_wp,    &  !< parameter 28  - heat capacity 4th wall layer ground floor level
+       0.35_wp,        &  !< parameter 29  - thermal conductivity 1st/2nd wall layer ground floor level
+       0.38_wp,        &  !< parameter 30  - thermal conductivity 3rd wall layer ground floor level
+       0.04_wp,        &  !< parameter 31  - thermal conductivity 4th wall layer ground floor level
+       0.92_wp,        &  !< parameter 32  - wall emissivity ground floor level
+       0.87_wp,        &  !< parameter 33  - window emissivity ground floor level
+       0.86_wp,        &  !< parameter 34  - green emissivity ground floor level
+       0.7_wp,         &  !< parameter 35  - window transmissivity ground floor level
+       0.01_wp,        &  !< parameter 36  - z0 roughness ground floor level
+       0.001_wp,       &  !< parameter 37  - z0h/z0q roughness heat/humidity
+       27.0_wp,        &  !< parameter 38  - wall albedo above ground floor level
+       5.0_wp,         &  !< parameter 39  - green albedo above ground floor level
+       27.0_wp,        &  !< parameter 40  - window albedo above ground floor level
+       0.005_wp,       &  !< parameter 41  - 1st wall layer thickness above ground floor level
+       0.01_wp,        &  !< parameter 42  - 2nd wall layer thickness above ground floor level
+       0.31_wp,        &  !< parameter 43  - 3rd wall layer thickness above ground floor level
+       0.43_wp,        &  !< parameter 44  - 4th wall layer thickness above ground floor level
+       20000.0_wp,     &  !< parameter 45  - heat capacity wall surface
+       23.0_wp,        &  !< parameter 46  - thermal conductivity of wall surface
+       20000.0_wp,     &  !< parameter 47  - heat capacity of window surface
+       20000.0_wp,     &  !< parameter 48  - heat capacity of green surface
+       23.0_wp,        &  !< parameter 49  - thermal conductivity of window surface
+       10.0_wp,        &  !< parameter 50  - thermal conductivty of green surface
+       1.0_wp,         &  !< parameter 51  - wall fraction ground plate
+       0.005_wp,       &  !< parameter 52  - 1st wall layer thickness ground plate
+       0.01_wp,        &  !< parameter 53  - 2nd wall layer thickness ground plate
+       0.31_wp,        &  !< parameter 54  - 3rd wall layer thickness ground plate
+       0.43_wp,        &  !< parameter 55  - 4th wall layer thickness ground plate
+       2000000.0_wp,   &  !< parameter 56  - heat capacity 1st/2nd wall layer ground plate
+       103000.0_wp,    &  !< parameter 57  - heat capacity 3rd wall layer ground plate
+       900000.0_wp,    &  !< parameter 58  - heat capacity 4th wall layer ground plate
+       0.35_wp,        &  !< parameter 59  - thermal conductivity 1st/2nd wall layer ground plate
+       0.38_wp,        &  !< parameter 60  - thermal conductivity 3rd wall layer ground plate
+       0.04_wp,        &  !< parameter 61  - thermal conductivity 4th wall layer ground plate
+       0.005_wp,       &  !< parameter 62  - 1st wall layer thickness ground floor level
+       0.01_wp,        &  !< parameter 63  - 2nd wall layer thickness ground floor level
+       0.31_wp,        &  !< parameter 64  - 3rd wall layer thickness ground floor level
+       0.43_wp,        &  !< parameter 65  - 4th wall layer thickness ground floor level
+       27.0_wp,        &  !< parameter 66  - wall albedo ground floor level
+       0.003_wp,       &  !< parameter 67  - 1st window layer thickness ground floor level
+       0.006_wp,       &  !< parameter 68  - 2nd window layer thickness ground floor level
+       0.012_wp,       &  !< parameter 69  - 3rd window layer thickness ground floor level
+       0.018_wp,       &  !< parameter 70  - 4th window layer thickness ground floor level
+       1736000.0_wp,   &  !< parameter 71  - heat capacity 1st/2nd window layer ground floor level
+       1736000.0_wp,   &  !< parameter 72  - heat capacity 3rd window layer ground floor level
+       1736000.0_wp,   &  !< parameter 73  - heat capacity 4th window layer ground floor level
+       0.11_wp,        &  !< parameter 74  - thermal conductivity 1st/2nd window layer ground floor level
+       0.11_wp,        &  !< parameter 75  - thermal conductivity 3rd window layer ground floor level
+       0.11_wp,        &  !< parameter 76  - thermal conductivity 4th window layer ground floor level
+       27.0_wp,        &  !< parameter 77  - window albedo ground floor level
+       5.0_wp,         &  !< parameter 78  - green albedo ground floor level
+       0.003_wp,       &  !< parameter 79  - 1st window layer thickness above ground floor level
+       0.006_wp,       &  !< parameter 80  - 2nd thickness window layer above ground floor level
+       0.012_wp,       &  !< parameter 81  - 3rd window layer thickness above ground floor level
+       0.018_wp,       &  !< parameter 82  - 4th window layer thickness above ground floor level
+       1736000.0_wp,   &  !< parameter 83  - heat capacity 1st/2nd window layer above ground floor level
+       1736000.0_wp,   &  !< parameter 84  - heat capacity 3rd window layer above ground floor level
+       1736000.0_wp,   &  !< parameter 85  - heat capacity 4th window layer above ground floor level
+       0.11_wp,        &  !< parameter 86  - thermal conductivity 1st/2nd window layer above ground floor level
+       0.11_wp,        &  !< parameter 87  - thermal conductivity 3rd window layer above ground floor level
+       0.11_wp,        &  !< parameter 88  - thermal conductivity 4th window layer above ground floor level
+       1.0_wp,         &  !< parameter 89  - wall fraction roof
+       0.005_wp,       &  !< parameter 90  - 1st wall layer thickness roof
+       0.01_wp,        &  !< parameter 91  - 2nd wall layer thickness roof
+       0.31_wp,        &  !< parameter 92  - 3rd wall layer thickness roof
+       0.43_wp,        &  !< parameter 93  - 4th wall layer thickness roof
+       2000000.0_wp,   &  !< parameter 94  - heat capacity 1st/2nd wall layer roof
+       103000.0_wp,    &  !< parameter 95  - heat capacity 3rd wall layer roof
+       900000.0_wp,    &  !< parameter 96  - heat capacity 4th wall layer roof
+       0.35_wp,        &  !< parameter 97  - thermal conductivity 1st/2nd wall layer roof
+       0.38_wp,        &  !< parameter 98  - thermal conductivity 3rd wall layer roof
+       0.04_wp,        &  !< parameter 99  - thermal conductivity 4th wall layer roof
+       0.91_wp,        &  !< parameter 100 - wall emissivity roof
+       27.0_wp,        &  !< parameter 101 - wall albedo roof
+       0.0_wp,         &  !< parameter 102 - window fraction roof
+       0.003_wp,       &  !< parameter 103 - window 1st layer thickness roof
+       0.006_wp,       &  !< parameter 104 - window 2nd layer thickness roof
+       0.012_wp,       &  !< parameter 105 - window 3rd layer thickness roof
+       0.018_wp,       &  !< parameter 106 - window 4th layer thickness roof
+       1736000.0_wp,   &  !< parameter 107 - heat capacity 1st/2nd window layer roof
+       1736000.0_wp,   &  !< parameter 108 - heat capacity 3rd window layer roof
+       1736000.0_wp,   &  !< parameter 109 - heat capacity 4th window layer roof
+       0.11_wp,        &  !< parameter 110 - thermal conductivity 1st/2nd window layer roof
+       0.11_wp,        &  !< parameter 111 - thermal conductivity 3rd window layer roof
+       0.11_wp,        &  !< parameter 112 - thermal conductivity 4th window layer roof
+       0.87_wp,        &  !< parameter 113 - window emissivity roof
+       0.7_wp,         &  !< parameter 114 - window transmissivity roof
+       27.0_wp,        &  !< parameter 115 - window albedo roof
+       0.86_wp,        &  !< parameter 116 - green emissivity roof
+       5.0_wp,         &  !< parameter 117 - green albedo roof
+       0.0_wp,         &  !< parameter 118 - green type roof
+       0.8_wp,         &  !< parameter 119 - shading factor
+       0.6_wp,         &  !< parameter 120 - g-value windows
+       3.0_wp,         &  !< parameter 121 - u-value windows
+       0.1_wp,         &  !< parameter 122 - basical airflow without occupancy of the room
+       1.5_wp,         &  !< parameter 123 - additional airflow depend of occupancy of the room
+       0.65_wp,        &  !< parameter 124 - heat recovery efficiency
+       2.5_wp,         &  !< parameter 125 - dynamic parameter specific effective surface
+       165000.0_wp,    &  !< parameter 126 - dynamic parameter innner heatstorage
+       4.5_wp,         &  !< parameter 127 - ratio internal surface/floor area
+       100000.0_wp,    &  !< parameter 128 - maximal heating capacity
+       0.0_wp,         &  !< parameter 129 - maximal cooling capacity
+       7.0_wp,         &  !< parameter 130 - additional internal heat gains dependent on occupancy of the room
+       20.0_wp,        &  !< parameter 131 - basic internal heat gains without occupancy of the room
+       3.0_wp,         &  !< parameter 132 - storey height
+       0.2_wp,         &  !< parameter 133 - ceiling construction height
+       0.425_wp,       &  !< parameter 0   - wall fraction above ground floor level
+       0.575_wp,       &  !< parameter 1   - window fraction above ground floor level
+       0.0_wp,         &  !< parameter 2   - green fraction above ground floor level
+       0.0_wp,         &  !< parameter 3   - green fraction roof above ground floor level
+       1.5_wp,         &  !< parameter 4   - LAI roof
+       1.5_wp,         &  !< parameter 5   - LAI on wall above ground floor level
+       2000000.0_wp,   &  !< parameter 6   - heat capacity 1st/2nd wall layer above ground floor level
+       103000.0_wp,    &  !< parameter 7   - heat capacity 3rd wall layer above ground floor level
+       900000.0_wp,    &  !< parameter 8   - heat capacity 4th wall layer above ground floor level
+       0.35_wp,        &  !< parameter 9   - thermal conductivity 1st/2nd wall layer above ground floor level
+       0.14_wp,        &  !< parameter 10  - thermal conductivity 3rd wall layer above ground floor level
+       0.035_wp,       &  !< parameter 11  - thermal conductivity 4th wall layer above ground floor level
+       299.15_wp,      &  !< parameter 12  - indoor target summer temperature
+       293.15_wp,      &  !< parameter 13  - indoor target winter temperature
+       0.92_wp,        &  !< parameter 14  - wall emissivity above ground floor level
+       0.86_wp,        &  !< parameter 15  - green emissivity above ground floor level
+       0.8_wp,         &  !< parameter 16  - window emissivity above ground floor level
+       0.6_wp,         &  !< parameter 17  - window transmissivity above ground floor level
+       0.001_wp,       &  !< parameter 18  - z0 roughness above ground floor level
+       0.0001_wp,      &  !< parameter 19  - z0h/z0g roughness heat/humidity above ground floor level
+       4.0_wp,         &  !< parameter 20  - ground floor level height
+       0.475_wp,       &  !< parameter 21  - wall fraction ground floor level
+       0.525_wp,       &  !< parameter 22  - window fraction ground floor level
+       0.0_wp,         &  !< parameter 23  - green fraction ground floor level
+       0.0_wp,         &  !< parameter 24  - green fraction roof ground floor level
+       1.5_wp,         &  !< parameter 25  - LAI on wall ground floor level
+       2000000.0_wp,   &  !< parameter 26  - heat capacity 1st/2nd wall layer ground floor level
+       103000.0_wp,    &  !< parameter 27  - heat capacity 3rd wall layer ground floor level
+       900000.0_wp,    &  !< parameter 28  - heat capacity 4th wall layer ground floor level
+       0.35_wp,        &  !< parameter 29  - thermal conductivity 1st/2nd wall layer ground floor level
+       0.14_wp,        &  !< parameter 30  - thermal conductivity 3rd wall layer ground floor level
+       0.035_wp,       &  !< parameter 31  - thermal conductivity 4th wall layer ground floor level
+       0.92_wp,        &  !< parameter 32  - wall emissivity ground floor level
+       0.8_wp,         &  !< parameter 33  - window emissivity ground floor level
+       0.86_wp,        &  !< parameter 34  - green emissivity ground floor level
+       0.6_wp,         &  !< parameter 35  - window transmissivity ground floor level
+       0.01_wp,        &  !< parameter 36  - z0 roughness ground floor level
+       0.001_wp,       &  !< parameter 37  - z0h/z0q roughness heat/humidity
+       27.0_wp,        &  !< parameter 38  - wall albedo above ground floor level
+       5.0_wp,         &  !< parameter 39  - green albedo above ground floor level
+       27.0_wp,        &  !< parameter 40  - window albedo above ground floor level
+       0.005_wp,       &  !< parameter 41  - 1st wall layer thickness above ground floor level
+       0.01_wp,        &  !< parameter 42  - 2nd wall layer thickness above ground floor level
+       0.41_wp,        &  !< parameter 43  - 3rd wall layer thickness above ground floor level
+       0.7_wp,         &  !< parameter 44  - 4th wall layer thickness above ground floor level
+       20000.0_wp,     &  !< parameter 45  - heat capacity wall surface
+       23.0_wp,        &  !< parameter 46  - thermal conductivity of wall surface
+       20000.0_wp,     &  !< parameter 47  - heat capacity of window surface
+       20000.0_wp,     &  !< parameter 48  - heat capacity of green surface
+       23.0_wp,        &  !< parameter 49  - thermal conductivity of window surface
+       10.0_wp,        &  !< parameter 50  - thermal conductivty of green surface
+       1.0_wp,         &  !< parameter 51  - wall fraction ground plate
+       0.005_wp,       &  !< parameter 52  - 1st wall layer thickness ground plate
+       0.01_wp,        &  !< parameter 53  - 2nd wall layer thickness ground plate
+       0.41_wp,        &  !< parameter 54  - 3rd wall layer thickness ground plate
+       0.7_wp,         &  !< parameter 55  - 4th wall layer thickness ground plate
+       2000000.0_wp,   &  !< parameter 56  - heat capacity 1st/2nd wall layer ground plate
+       103000.0_wp,    &  !< parameter 57  - heat capacity 3rd wall layer ground plate
+       900000.0_wp,    &  !< parameter 58  - heat capacity 4th wall layer ground plate
+       0.35_wp,        &  !< parameter 59  - thermal conductivity 1st/2nd wall layer ground plate
+       0.14_wp,        &  !< parameter 60  - thermal conductivity 3rd wall layer ground plate
+       0.035_wp,       &  !< parameter 61  - thermal conductivity 4th wall layer ground plate
+       0.005_wp,       &  !< parameter 62  - 1st wall layer thickness ground floor level
+       0.01_wp,        &  !< parameter 63  - 2nd wall layer thickness ground floor level
+       0.41_wp,        &  !< parameter 64  - 3rd wall layer thickness ground floor level
+       0.7_wp,         &  !< parameter 65  - 4th wall layer thickness ground floor level
+       27.0_wp,        &  !< parameter 66  - wall albedo ground floor level
+       0.003_wp,       &  !< parameter 67  - 1st window layer thickness ground floor level
+       0.006_wp,       &  !< parameter 68  - 2nd window layer thickness ground floor level
+       0.012_wp,       &  !< parameter 69  - 3rd window layer thickness ground floor level
+       0.018_wp,       &  !< parameter 70  - 4th window layer thickness ground floor level
+       1736000.0_wp,   &  !< parameter 71  - heat capacity 1st/2nd window layer ground floor level
+       1736000.0_wp,   &  !< parameter 72  - heat capacity 3rd window layer ground floor level
+       1736000.0_wp,   &  !< parameter 73  - heat capacity 4th window layer ground floor level
+       0.037_wp,       &  !< parameter 74  - thermal conductivity 1st/2nd window layer ground floor level
+       0.037_wp,       &  !< parameter 75  - thermal conductivity 3rd window layer ground floor level
+       0.037_wp,       &  !< parameter 76  - thermal conductivity 4th window layer ground floor level
+       27.0_wp,        &  !< parameter 77  - window albedo ground floor level
+       5.0_wp,         &  !< parameter 78  - green albedo ground floor level
+       0.003_wp,       &  !< parameter 79  - 1st window layer thickness above ground floor level
+       0.006_wp,       &  !< parameter 80  - 2nd thickness window layer above ground floor level
+       0.012_wp,       &  !< parameter 81  - 3rd window layer thickness above ground floor level
+       0.018_wp,       &  !< parameter 82  - 4th window layer thickness above ground floor level
+       1736000.0_wp,   &  !< parameter 83  - heat capacity 1st/2nd window layer above ground floor level
+       1736000.0_wp,   &  !< parameter 84  - heat capacity 3rd window layer above ground floor level
+       1736000.0_wp,   &  !< parameter 85  - heat capacity 4th window layer above ground floor level
+       0.037_wp,       &  !< parameter 86  - thermal conductivity 1st/2nd window layer above ground floor level
+       0.037_wp,       &  !< parameter 87  - thermal conductivity 3rd window layer above ground floor level
+       0.037_wp,       &  !< parameter 88  - thermal conductivity 4th window layer above ground floor level
+       1.0_wp,         &  !< parameter 89  - wall fraction roof
+       0.005_wp,       &  !< parameter 90  - 1st wall layer thickness roof
+       0.01_wp,        &  !< parameter 91  - 2nd wall layer thickness roof
+       0.41_wp,        &  !< parameter 92  - 3rd wall layer thickness roof
+       0.7_wp,         &  !< parameter 93  - 4th wall layer thickness roof
+       2000000.0_wp,   &  !< parameter 94  - heat capacity 1st/2nd wall layer roof
+       103000.0_wp,    &  !< parameter 95  - heat capacity 3rd wall layer roof
+       900000.0_wp,    &  !< parameter 96  - heat capacity 4th wall layer roof
+       0.35_wp,        &  !< parameter 97  - thermal conductivity 1st/2nd wall layer roof
+       0.14_wp,        &  !< parameter 98  - thermal conductivity 3rd wall layer roof
+       0.035_wp,       &  !< parameter 99  - thermal conductivity 4th wall layer roof
+       0.91_wp,        &  !< parameter 100 - wall emissivity roof
+       27.0_wp,        &  !< parameter 101 - wall albedo roof
+       0.0_wp,         &  !< parameter 102 - window fraction roof
+       0.003_wp,       &  !< parameter 103 - window 1st layer thickness roof
+       0.006_wp,       &  !< parameter 104 - window 2nd layer thickness roof
+       0.012_wp,       &  !< parameter 105 - window 3rd layer thickness roof
+       0.018_wp,       &  !< parameter 106 - window 4th layer thickness roof
+       1736000.0_wp,   &  !< parameter 107 - heat capacity 1st/2nd window layer roof
+       1736000.0_wp,   &  !< parameter 108 - heat capacity 3rd window layer roof
+       1736000.0_wp,   &  !< parameter 109 - heat capacity 4th window layer roof
+       0.037_wp,       &  !< parameter 110 - thermal conductivity 1st/2nd window layer roof
+       0.037_wp,       &  !< parameter 111 - thermal conductivity 3rd window layer roof
+       0.037_wp,       &  !< parameter 112 - thermal conductivity 4th window layer roof
+       0.8_wp,         &  !< parameter 113 - window emissivity roof
+       0.6_wp,         &  !< parameter 114 - window transmissivity roof
+       27.0_wp,        &  !< parameter 115 - window albedo roof
+       0.86_wp,        &  !< parameter 116 - green emissivity roof
+       5.0_wp,         &  !< parameter 117 - green albedo roof
+       0.0_wp,         &  !< parameter 118 - green type roof
+       0.8_wp,         &  !< parameter 119 - shading factor
+       0.5_wp,         &  !< parameter 120 - g-value windows
+       0.6_wp,         &  !< parameter 121 - u-value windows
+       0.1_wp,         &  !< parameter 122 - basical airflow without occupancy of the room
+       1.5_wp,         &  !< parameter 123 - additional airflow depend of occupancy of the room
+       0.9_wp,         &  !< parameter 124 - heat recovery efficiency
+       2.5_wp,         &  !< parameter 125 - dynamic parameter specific effective surface
+       80000.0_wp,     &  !< parameter 126 - dynamic parameter innner heatstorage
+       4.5_wp,         &  !< parameter 127 - ratio internal surface/floor area
+       100000.0_wp,    &  !< parameter 128 - maximal heating capacity
+       0.0_wp,         &  !< parameter 129 - maximal cooling capacity
+       5.0_wp,         &  !< parameter 130 - additional internal heat gains dependent on occupancy of the room
+       15.0_wp,        &  !< parameter 131 - basic internal heat gains without occupancy of the room
+       3.0_wp,         &  !< parameter 132 - storey height
+       0.2_wp,         &  !< parameter 133 - ceiling construction height
+       1.0_wp,         &  !< parameter 0   - wall fraction above ground floor level
+       0.0_wp,         &  !< parameter 1   - window fraction above ground floor level
+       0.0_wp,         &  !< parameter 2   - green fraction above ground floor level
+       0.0_wp,         &  !< parameter 3   - green fraction roof above ground floor level
+       1.5_wp,         &  !< parameter 4   - LAI roof
+       1.5_wp,         &  !< parameter 5   - LAI on wall above ground floor level
+       1950400.0_wp,   &  !< parameter 6   - heat capacity 1st/2nd wall layer above ground floor level
+       1848000.0_wp,   &  !< parameter 7   - heat capacity 3rd wall layer above ground floor level
+       1848000.0_wp,   &  !< parameter 8   - heat capacity 4th wall layer above ground floor level
+       0.7_wp,         &  !< parameter 9   - thermal conductivity 1st/2nd wall layer above ground floor level
+       1.0_wp,         &  !< parameter 10  - thermal conductivity 3rd wall layer above ground floor level
+       1.0_wp,         &  !< parameter 11  - thermal conductivity 4th wall layer above ground floor level
+       299.15_wp,      &  !< parameter 12  - indoor target summer temperature
+       293.15_wp,      &  !< parameter 13  - indoor target winter temperature
+       0.9_wp,         &  !< parameter 14  - wall emissivity above ground floor level
+       0.86_wp,        &  !< parameter 15  - green emissivity above ground floor level
+       0.8_wp,         &  !< parameter 16  - window emissivity above ground floor level
+       0.6_wp,         &  !< parameter 17  - window transmissivity above ground floor level
+       0.001_wp,       &  !< parameter 18  - z0 roughness above ground floor level
+       0.0001_wp,      &  !< parameter 19  - z0h/z0g roughness heat/humidity above ground floor level
+       4.0_wp,         &  !< parameter 20  - ground floor level height
+       1.0_wp,         &  !< parameter 21  - wall fraction ground floor level
+       0.0_wp,         &  !< parameter 22  - window fraction ground floor level
+       0.0_wp,         &  !< parameter 23  - green fraction ground floor level
+       0.0_wp,         &  !< parameter 24  - green fraction roof ground floor level
+       1.5_wp,         &  !< parameter 25  - LAI on wall ground floor level
+       1950400.0_wp,   &  !< parameter 26  - heat capacity 1st/2nd wall layer ground floor level
+       1848000.0_wp,   &  !< parameter 27  - heat capacity 3rd wall layer ground floor level
+       1848000.0_wp,   &  !< parameter 28  - heat capacity 4th wall layer ground floor level
+       0.7_wp,         &  !< parameter 29  - thermal conductivity 1st/2nd wall layer ground floor level
+       1.0_wp,         &  !< parameter 30  - thermal conductivity 3rd wall layer ground floor level
+       1.0_wp,         &  !< parameter 31  - thermal conductivity 4th wall layer ground floor level
+       0.9_wp,         &  !< parameter 32  - wall emissivity ground floor level
+       0.8_wp,         &  !< parameter 33  - window emissivity ground floor level
+       0.86_wp,        &  !< parameter 34  - green emissivity ground floor level
+       0.6_wp,         &  !< parameter 35  - window transmissivity ground floor level
+       0.01_wp,        &  !< parameter 36  - z0 roughness ground floor level
+       0.001_wp,       &  !< parameter 37  - z0h/z0q roughness heat/humidity
+       27.0_wp,        &  !< parameter 38  - wall albedo above ground floor level
+       5.0_wp,         &  !< parameter 39  - green albedo above ground floor level
+       27.0_wp,        &  !< parameter 40  - window albedo above ground floor level
+       0.29_wp,        &  !< parameter 41  - 1st wall layer thickness above ground floor level
+       0.295_wp,       &  !< parameter 42  - 2nd wall layer thickness above ground floor level
+       0.695_wp,       &  !< parameter 43  - 3rd wall layer thickness above ground floor level
+       0.985_wp,       &  !< parameter 44  - 4th wall layer thickness above ground floor level
+       20000.0_wp,     &  !< parameter 45  - heat capacity wall surface
+       23.0_wp,        &  !< parameter 46  - thermal conductivity of wall surface
+       20000.0_wp,     &  !< parameter 47  - heat capacity of window surface
+       20000.0_wp,     &  !< parameter 48  - heat capacity of green surface
+       23.0_wp,        &  !< parameter 49  - thermal conductivity of window surface
+       10.0_wp,        &  !< parameter 50  - thermal conductivty of green surface
+       1.0_wp,         &  !< parameter 51  - wall fraction ground plate
+       0.29_wp,        &  !< parameter 52  - 1st wall layer thickness ground plate
+       0.295_wp,       &  !< parameter 53  - 2nd wall layer thickness ground plate
+       0.695_wp,       &  !< parameter 54  - 3rd wall layer thickness ground plate
+       0.985_wp,       &  !< parameter 55  - 4th wall layer thickness ground plate
+       1950400.0_wp,   &  !< parameter 56  - heat capacity 1st/2nd wall layer ground plate
+       1848000.0_wp,   &  !< parameter 57  - heat capacity 3rd wall layer ground plate
+       1848000.0_wp,   &  !< parameter 58  - heat capacity 4th wall layer ground plate
+       0.7_wp,         &  !< parameter 59  - thermal conductivity 1st/2nd wall layer ground plate
+       1.0_wp,         &  !< parameter 60  - thermal conductivity 3rd wall layer ground plate
+       1.0_wp,         &  !< parameter 61  - thermal conductivity 4th wall layer ground plate
+       0.29_wp,        &  !< parameter 62  - 1st wall layer thickness ground floor level
+       0.295_wp,       &  !< parameter 63  - 2nd wall layer thickness ground floor level
+       0.695_wp,       &  !< parameter 64  - 3rd wall layer thickness ground floor level
+       0.985_wp,       &  !< parameter 65  - 4th wall layer thickness ground floor level
+       27.0_wp,        &  !< parameter 66  - wall albedo ground floor level
+       0.003_wp,       &  !< parameter 67  - 1st window layer thickness ground floor level
+       0.006_wp,       &  !< parameter 68  - 2nd window layer thickness ground floor level
+       0.012_wp,       &  !< parameter 69  - 3rd window layer thickness ground floor level
+       0.018_wp,       &  !< parameter 70  - 4th window layer thickness ground floor level
+       1736000.0_wp,   &  !< parameter 71  - heat capacity 1st/2nd window layer ground floor level
+       1736000.0_wp,   &  !< parameter 72  - heat capacity 3rd window layer ground floor level
+       1736000.0_wp,   &  !< parameter 73  - heat capacity 4th window layer ground floor level
+       0.57_wp,        &  !< parameter 74  - thermal conductivity 1st/2nd window layer ground floor level
+       0.57_wp,        &  !< parameter 75  - thermal conductivity 3rd window layer ground floor level
+       0.57_wp,        &  !< parameter 76  - thermal conductivity 4th window layer ground floor level
+       27.0_wp,        &  !< parameter 77  - window albedo ground floor level
+       5.0_wp,         &  !< parameter 78  - green albedo ground floor level
+       0.003_wp,       &  !< parameter 79  - 1st window layer thickness above ground floor level
+       0.006_wp,       &  !< parameter 80  - 2nd thickness window layer above ground floor level
+       0.012_wp,       &  !< parameter 81  - 3rd window layer thickness above ground floor level
+       0.018_wp,       &  !< parameter 82  - 4th window layer thickness above ground floor level
+       1736000.0_wp,   &  !< parameter 83  - heat capacity 1st/2nd window layer above ground floor level
+       1736000.0_wp,   &  !< parameter 84  - heat capacity 3rd window layer above ground floor level
+       1736000.0_wp,   &  !< parameter 85  - heat capacity 4th window layer above ground floor level
+       0.57_wp,        &  !< parameter 86  - thermal conductivity 1st/2nd window layer above ground floor level
+       0.57_wp,        &  !< parameter 87  - thermal conductivity 3rd window layer above ground floor level
+       0.57_wp,        &  !< parameter 88  - thermal conductivity 4th window layer above ground floor level
+       1.0_wp,         &  !< parameter 89  - wall fraction roof
+       0.29_wp,        &  !< parameter 90  - 1st wall layer thickness roof
+       0.295_wp,       &  !< parameter 91  - 2nd wall layer thickness roof
+       0.695_wp,       &  !< parameter 92  - 3rd wall layer thickness roof
+       0.985_wp,       &  !< parameter 93  - 4th wall layer thickness roof
+       1950400.0_wp,   &  !< parameter 94  - heat capacity 1st/2nd wall layer roof
+       1848000.0_wp,   &  !< parameter 95  - heat capacity 3rd wall layer roof
+       1848000.0_wp,   &  !< parameter 96  - heat capacity 4th wall layer roof
+       0.7_wp,         &  !< parameter 97  - thermal conductivity 1st/2nd wall layer roof
+       1.0_wp,         &  !< parameter 98  - thermal conductivity 3rd wall layer roof
+       1.0_wp,         &  !< parameter 99  - thermal conductivity 4th wall layer roof
+       0.9_wp,         &  !< parameter 100 - wall emissivity roof
+       27.0_wp,        &  !< parameter 101 - wall albedo roof
+       0.0_wp,         &  !< parameter 102 - window fraction roof
+       0.003_wp,       &  !< parameter 103 - window 1st layer thickness roof
+       0.006_wp,       &  !< parameter 104 - window 2nd layer thickness roof
+       0.012_wp,       &  !< parameter 105 - window 3rd layer thickness roof
+       0.018_wp,       &  !< parameter 106 - window 4th layer thickness roof
+       1736000.0_wp,   &  !< parameter 107 - heat capacity 1st/2nd window layer roof
+       1736000.0_wp,   &  !< parameter 108 - heat capacity 3rd window layer roof
+       1736000.0_wp,   &  !< parameter 109 - heat capacity 4th window layer roof
+       0.57_wp,        &  !< parameter 110 - thermal conductivity 1st/2nd window layer roof
+       0.57_wp,        &  !< parameter 111 - thermal conductivity 3rd window layer roof
+       0.57_wp,        &  !< parameter 112 - thermal conductivity 4th window layer roof
+       0.8_wp,         &  !< parameter 113 - window emissivity roof
+       0.6_wp,         &  !< parameter 114 - window transmissivity roof
+       27.0_wp,        &  !< parameter 115 - window albedo roof
+       0.86_wp,        &  !< parameter 116 - green emissivity roof
+       5.0_wp,         &  !< parameter 117 - green albedo roof
+       0.0_wp,         &  !< parameter 118 - green type roof
+       0.8_wp,         &  !< parameter 119 - shading factor
+       100.0_wp,       &  !< parameter 120 - g-value windows
+       100.0_wp,       &  !< parameter 121 - u-value windows
+       20.0_wp,        &  !< parameter 122 - basical airflow without occupancy of the room
+       20.0_wp,        &  !< parameter 123 - additional airflow depend of occupancy of the room
+       0.0_wp,         &  !< parameter 124 - heat recovery efficiency
+       1.0_wp,         &  !< parameter 125 - dynamic parameter specific effective surface
+       1.0_wp,         &  !< parameter 126 - dynamic parameter innner heatstorage
+       4.5_wp,         &  !< parameter 127 - ratio internal surface/floor area
+       100000.0_wp,    &  !< parameter 128 - maximal heating capacity
+       0.0_wp,         &  !< parameter 129 - maximal cooling capacity
+       0.0_wp,         &  !< parameter 130 - additional internal heat gains dependent on occupancy of the room
+       0.0_wp,         &  !< parameter 131 - basic internal heat gains without occupancy of the room
+       3.0_wp,         &  !< parameter 132 - storey height
+       0.2_wp          &  !< parameter 133 - ceiling construction height
+       /),             &
+       (/134, 7/) )
 
 !
@@ -2653,5 +3336,5 @@
         CHARACTER(LEN=2)                              :: ls
         CHARACTER(LEN=varnamelength)                  :: var           !< TRIM(variable)
-        INTEGER(iwp), PARAMETER                       :: nl1 = 16      !< number of directional usm variables
+        INTEGER(iwp), PARAMETER                       :: nl1 = 14      !< number of directional usm variables
         CHARACTER(LEN=varnamelength), DIMENSION(nl1)  :: varlist1 = &  !< list of directional usm variables
                   (/'usm_wshf                      ', &
@@ -2664,6 +3347,4 @@
                     'usm_surfwintrans              ', &
                     'usm_surfcat                   ', &
-                    'usm_surfalb                   ', &
-                    'usm_surfemis                  ', &
                     'usm_t_surf_wall               ', &
                     'usm_t_surf_window             ', &
@@ -2734,6 +3415,5 @@
                   var(1:14) == 'usm_theta_10cm' )  THEN
             unit = 'K'
-        ELSE IF ( var(1:9) == 'usm_surfz'  .OR.  var(1:11) == 'usm_surfcat'  .OR.         &
-                  var(1:11) == 'usm_surfalb'  .OR.  var(1:12) == 'usm_surfemis'  )  THEN
+        ELSE IF ( var(1:9) == 'usm_surfz'  .OR.  var(1:11) == 'usm_surfcat' )  THEN
             unit = '1'
         ELSE
@@ -2851,5 +3531,5 @@
         LOGICAL, INTENT(OUT)           ::  found     !< 
         REAL(sp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do) ::  local_pf   !< sp - it has to correspond to module data_output_3d
-        REAL(wp), DIMENSION(nzb:nzt+1,nys:nyn,nxl:nxr)     ::  temp_pf    !< temp array for urban surface output procedure
+        REAL(sp), DIMENSION(nzb:nzt+1,nys:nyn,nxl:nxr)     ::  temp_pf    !< temp array for urban surface output procedure
         
         CHARACTER (len=varnamelength)                      :: var     !< trimmed variable name
@@ -2956,64 +3636,4 @@
               ENDIF
               
-          CASE ( 'usm_surfalb' )
-!
-!--           surface albedo, weighted average
-              IF ( idsint == iup_u )  THEN
-                 DO  m = 1, surf_usm_h%ns
-                    i = surf_usm_h%i(m)
-                    j = surf_usm_h%j(m)
-                    k = surf_usm_h%k(m)
-                    temp_pf(k,j,i) = surf_usm_h%frac(ind_veg_wall,m)     *     &
-                                     surf_usm_h%albedo(ind_veg_wall,m)  +      &
-                                     surf_usm_h%frac(ind_pav_green,m)    *     &
-                                     surf_usm_h%albedo(ind_pav_green,m) +      &
-                                     surf_usm_h%frac(ind_wat_win,m)      *     &
-                                     surf_usm_h%albedo(ind_wat_win,m)
-                 ENDDO
-              ELSE
-                 l = idsidx
-                 DO  m = 1, surf_usm_v(l)%ns
-                    i = surf_usm_v(l)%i(m)
-                    j = surf_usm_v(l)%j(m)
-                    k = surf_usm_v(l)%k(m)
-                    temp_pf(k,j,i) = surf_usm_v(l)%frac(ind_veg_wall,m)     *  &
-                                     surf_usm_v(l)%albedo(ind_veg_wall,m)  +   &
-                                     surf_usm_v(l)%frac(ind_pav_green,m)    *  &
-                                     surf_usm_v(l)%albedo(ind_pav_green,m) +   &
-                                     surf_usm_v(l)%frac(ind_wat_win,m)      *  &
-                                     surf_usm_v(l)%albedo(ind_wat_win,m)
-                 ENDDO
-              ENDIF
-              
-          CASE ( 'usm_surfemis' )
-!
-!--           surface emissivity, weighted average
-              IF ( idsint == iup_u )  THEN
-                 DO  m = 1, surf_usm_h%ns
-                    i = surf_usm_h%i(m)
-                    j = surf_usm_h%j(m)
-                    k = surf_usm_h%k(m)
-                    temp_pf(k,j,i) =  surf_usm_h%frac(ind_veg_wall,m)      *   &
-                                      surf_usm_h%emissivity(ind_veg_wall,m)  + &
-                                      surf_usm_h%frac(ind_pav_green,m)     *   &
-                                      surf_usm_h%emissivity(ind_pav_green,m) + &
-                                      surf_usm_h%frac(ind_wat_win,m)       *   &
-                                      surf_usm_h%emissivity(ind_wat_win,m)
-                 ENDDO
-              ELSE
-                 l = idsidx
-                 DO  m = 1, surf_usm_v(l)%ns
-                    i = surf_usm_v(l)%i(m)
-                    j = surf_usm_v(l)%j(m)
-                    k = surf_usm_v(l)%k(m)
-                    temp_pf(k,j,i) = surf_usm_v(l)%frac(ind_veg_wall,m)       *&
-                                     surf_usm_v(l)%emissivity(ind_veg_wall,m) +&
-                                     surf_usm_v(l)%frac(ind_pav_green,m)      *&
-                                     surf_usm_v(l)%emissivity(ind_pav_green,m)+&
-                                     surf_usm_v(l)%frac(ind_wat_win,m)        *&
-                                     surf_usm_v(l)%emissivity(ind_wat_win,m)
-                 ENDDO
-              ENDIF
-
           CASE ( 'usm_surfwintrans' )
 !
@@ -3759,5 +4379,4 @@
              var(1:15) == 'usm_theta_10cm' .OR.                                             &
              var(1:9) == 'usm_surfz'  .OR.  var(1:11) == 'usm_surfcat'  .OR.                &
-             var(1:11) == 'usm_surfalb'  .OR.  var(1:12) == 'usm_surfemis'  .OR.            &
              var(1:16) == 'usm_surfwintrans'  .OR. var(1:7) == 'usm_swc' ) THEN
 
@@ -4058,29 +4677,4 @@
 !--     ground-floor-level properties. 
         surf_usm_h%ground_level = .FALSE. 
-!         DO  m = 1, surf_usm_h%ns
-!            i = surf_usm_h%i(m)
-!            j = surf_usm_h%j(m)
-!            k = surf_usm_h%k(m)
-! !
-! !--        Get local ground level. If no ground level is given in input file,
-! !--        use default value.
-!            ground_floor_level_l = ground_floor_level
-!            IF ( building_pars_f%from_file )  THEN
-!               IF ( building_pars_f%pars_xy(ind_gflh,j,i) /=                    &
-!                    building_pars_f%fill )  &
-!                  ground_floor_level_l = building_pars_f%pars_xy(ind_gflh,j,i)          
-!            ENDIF
-! !
-! !--        Determine height of surface element above ground level
-!            IF (  terrain_height_f%from_file )  THEN
-!               z_agl = zw(k) - terrain_height_f%var(j,i)
-!            ELSE
-!               z_agl = zw(k)
-!            ENDIF
-! !
-! !--        Set flag for ground level
-!            IF ( z_agl <= ground_floor_level_l )                                &
-!               surf_usm_h%ground_level(m) = .TRUE.
-!         ENDDO
 
         DO  l = 0, 3
@@ -4091,10 +4685,17 @@
               k = surf_usm_v(l)%k(m)
 !
-!--           Get local ground level. If no ground level is given in input file,
-!--           use default value.
+!--           Determine local ground level. Level 1 - default value,
+!--           level 2 - initialization according to building type, 
+!--           level 3 - initialization from value read from file.
               ground_floor_level_l = ground_floor_level
+              
+              IF ( building_type_f%from_file )  THEN
+                  ground_floor_level_l =                                       &
+                              building_pars(ind_gflh,building_type_f%var(j,i))
+              ENDIF
+              
               IF ( building_pars_f%from_file )  THEN
                  IF ( building_pars_f%pars_xy(ind_gflh,j,i) /=                 &
-                      building_pars_f%fill ) &
+                      building_pars_f%fill )                                   &
                     ground_floor_level_l = building_pars_f%pars_xy(ind_gflh,j,i)
               ENDIF
@@ -4102,6 +4703,8 @@
 !--           Determine height of surface element above ground level. Please 
 !--           note, height of surface element is determined with respect to
-!--           its height of the adjoing atmospheric grid point. 
-              IF (  terrain_height_f%from_file )  THEN
+!--           its height above ground of the reference grid point in atmosphere,
+!--           Therefore, substract the offset values when assessing the terrain
+!--           height.
+              IF ( terrain_height_f%from_file )  THEN
                  z_agl = zw(k) - terrain_height_f%var(j-surf_usm_v(l)%joff,    &
                                                       i-surf_usm_v(l)%ioff)
@@ -4149,5 +4752,5 @@
 !
 !--     Initialize urban-type surface attribute. According to initialization in 
-!--     land-surface model, follow a 3-level approach.
+!--     land-surface model, follow a 3-level approach. 
 !--     Level 1 - initialization via default attributes 
         DO  m = 1, surf_usm_h%ns
@@ -4204,6 +4807,6 @@
            surf_usm_h%lambda_h_window(nzb_wall+3,m) = building_pars(ind_tc3_win_r,building_type)    
 
-           surf_usm_h%target_temp_summer(m)  = building_pars(117,building_type)    
-           surf_usm_h%target_temp_winter(m)  = building_pars(118,building_type)    
+           surf_usm_h%target_temp_summer(m)  = building_pars(ind_indoor_target_temp_summer,building_type)    
+           surf_usm_h%target_temp_winter(m)  = building_pars(ind_indoor_target_temp_winter,building_type)    
 !
 !--        emissivity of wall-, green- and window fraction 
@@ -4238,10 +4841,10 @@
            surf_usm_h%zw_window(nzb_wall+3,m)       = building_pars(ind_thick_4_win_r,building_type)
 
-           surf_usm_h%c_surface(m)           = building_pars(0,building_type)  
-           surf_usm_h%lambda_surf(m)         = building_pars(3,building_type)  
-           surf_usm_h%c_surface_green(m)     = building_pars(2,building_type)  
-           surf_usm_h%lambda_surf_green(m)   = building_pars(5,building_type)  
-           surf_usm_h%c_surface_window(m)    = building_pars(1,building_type)  
-           surf_usm_h%lambda_surf_window(m)  = building_pars(4,building_type)  
+           surf_usm_h%c_surface(m)           = building_pars(ind_c_surface,building_type)  
+           surf_usm_h%lambda_surf(m)         = building_pars(ind_lambda_surf,building_type)  
+           surf_usm_h%c_surface_green(m)     = building_pars(ind_c_surface_green,building_type)  
+           surf_usm_h%lambda_surf_green(m)   = building_pars(ind_lambda_surf_green,building_type)  
+           surf_usm_h%c_surface_window(m)    = building_pars(ind_c_surface_win,building_type)  
+           surf_usm_h%lambda_surf_window(m)  = building_pars(ind_lambda_surf_win,building_type)  
            
            surf_usm_h%green_type_roof(m)     = building_pars(ind_green_type_roof,building_type)
@@ -4367,6 +4970,6 @@
               surf_usm_v(l)%lambda_h_window(nzb_wall+3,m) = building_pars(ind_tc3_win,building_type)    
 
-              surf_usm_v(l)%target_temp_summer(m)  = building_pars(117,building_type)    
-              surf_usm_v(l)%target_temp_winter(m)  = building_pars(118,building_type)    
+              surf_usm_v(l)%target_temp_summer(m)  = building_pars(ind_indoor_target_temp_summer,building_type)    
+              surf_usm_v(l)%target_temp_winter(m)  = building_pars(ind_indoor_target_temp_winter,building_type)    
 !
 !--           emissivity of wall-, green- and window fraction 
@@ -4400,10 +5003,10 @@
               surf_usm_v(l)%zw_window(nzb_wall+3,m)       = building_pars(ind_thick_4_win,building_type)
 
-              surf_usm_v(l)%c_surface(m)           = building_pars(0,building_type)  
-              surf_usm_v(l)%lambda_surf(m)         = building_pars(3,building_type)
-              surf_usm_v(l)%c_surface_green(m)     = building_pars(2,building_type)  
-              surf_usm_v(l)%lambda_surf_green(m)   = building_pars(5,building_type)
-              surf_usm_v(l)%c_surface_window(m)    = building_pars(1,building_type)  
-              surf_usm_v(l)%lambda_surf_window(m)  = building_pars(4,building_type)
+              surf_usm_v(l)%c_surface(m)           = building_pars(ind_c_surface,building_type)  
+              surf_usm_v(l)%lambda_surf(m)         = building_pars(ind_lambda_surf,building_type)
+              surf_usm_v(l)%c_surface_green(m)     = building_pars(ind_c_surface_green,building_type)  
+              surf_usm_v(l)%lambda_surf_green(m)   = building_pars(ind_lambda_surf_green,building_type)
+              surf_usm_v(l)%c_surface_window(m)    = building_pars(ind_c_surface_win,building_type)  
+              surf_usm_v(l)%lambda_surf_window(m)  = building_pars(ind_lambda_surf_win,building_type)
 
            ENDDO
@@ -4470,6 +5073,6 @@
                  surf_usm_h%lambda_h_window(nzb_wall+3,m) = building_pars(ind_tc3_win_r,st)    
 
-                 surf_usm_h%target_temp_summer(m)  = building_pars(117,st)    
-                 surf_usm_h%target_temp_winter(m)  = building_pars(118,st)    
+                 surf_usm_h%target_temp_summer(m)  = building_pars(ind_indoor_target_temp_summer,st)    
+                 surf_usm_h%target_temp_winter(m)  = building_pars(ind_indoor_target_temp_winter,st)    
 !
 !--              emissivity of wall-, green- and window fraction 
@@ -4504,10 +5107,10 @@
                  surf_usm_h%zw_window(nzb_wall+3,m)       = building_pars(ind_thick_4_win_r,st)
 
-                 surf_usm_h%c_surface(m)           = building_pars(0,st)  
-                 surf_usm_h%lambda_surf(m)         = building_pars(3,st)
-                 surf_usm_h%c_surface_green(m)     = building_pars(2,st)  
-                 surf_usm_h%lambda_surf_green(m)   = building_pars(5,st)
-                 surf_usm_h%c_surface_window(m)    = building_pars(1,st)  
-                 surf_usm_h%lambda_surf_window(m)  = building_pars(4,st)
+                 surf_usm_h%c_surface(m)           = building_pars(ind_c_surface,st)  
+                 surf_usm_h%lambda_surf(m)         = building_pars(ind_lambda_surf,st)
+                 surf_usm_h%c_surface_green(m)     = building_pars(ind_c_surface_green,st)  
+                 surf_usm_h%lambda_surf_green(m)   = building_pars(ind_lambda_surf_green,st)
+                 surf_usm_h%c_surface_window(m)    = building_pars(ind_c_surface_win,st)  
+                 surf_usm_h%lambda_surf_window(m)  = building_pars(ind_lambda_surf_win,st)
                  
                  surf_usm_h%green_type_roof(m)     = building_pars(ind_green_type_roof,st)
@@ -4640,6 +5243,6 @@
                     surf_usm_v(l)%lambda_h_window(nzb_wall+3,m) = building_pars(ind_tc3_win,st) 
 
-                    surf_usm_v(l)%target_temp_summer(m)  = building_pars(117,st)    
-                    surf_usm_v(l)%target_temp_winter(m)  = building_pars(118,st)    
+                    surf_usm_v(l)%target_temp_summer(m)  = building_pars(ind_indoor_target_temp_summer,st)    
+                    surf_usm_v(l)%target_temp_winter(m)  = building_pars(ind_indoor_target_temp_winter,st)    
 !
 !--                 emissivity of wall-, green- and window fraction 
@@ -4673,10 +5276,10 @@
                     surf_usm_v(l)%zw_window(nzb_wall+3,m)       = building_pars(ind_thick_4_win,st)
 
-                    surf_usm_v(l)%c_surface(m)           = building_pars(0,st)  
-                    surf_usm_v(l)%lambda_surf(m)         = building_pars(3,st) 
-                    surf_usm_v(l)%c_surface_green(m)     = building_pars(2,st)  
-                    surf_usm_v(l)%lambda_surf_green(m)   = building_pars(5,st) 
-                    surf_usm_v(l)%c_surface_window(m)    = building_pars(1,st)  
-                    surf_usm_v(l)%lambda_surf_window(m)  = building_pars(4,st) 
+                    surf_usm_v(l)%c_surface(m)           = building_pars(ind_c_surface,st)  
+                    surf_usm_v(l)%lambda_surf(m)         = building_pars(ind_lambda_surf,st) 
+                    surf_usm_v(l)%c_surface_green(m)     = building_pars(ind_c_surface_green,st)  
+                    surf_usm_v(l)%lambda_surf_green(m)   = building_pars(ind_lambda_surf_green,st) 
+                    surf_usm_v(l)%c_surface_window(m)    = building_pars(ind_c_surface_win,st)  
+                    surf_usm_v(l)%lambda_surf_window(m)  = building_pars(ind_lambda_surf_win,st) 
 
 
@@ -4687,5 +5290,7 @@
         
 !
-!--     Level 3 - initialization via building_pars read from file
+!--     Level 3 - initialization via building_pars read from file. Note, only
+!--     variables that are also defined in the input-standard can be initialized
+!--     via file. Other variables will be initialized on level 1 or 2. 
         IF ( building_pars_f%from_file )  THEN
            DO  m = 1, surf_usm_h%ns
@@ -4696,145 +5301,283 @@
 !--           In order to distinguish between ground floor level and 
 !--           above-ground-floor level surfaces, set input indices.
-              ind_green_frac_r = MERGE( ind_green_frac_r_gfl, ind_green_frac_r_agfl, &
+              ind_wall_frac    = MERGE( ind_wall_frac_gfl,                     &
+                                        ind_wall_frac_agfl,                    &
                                         surf_usm_h%ground_level(m) )
-              ind_lai_r        = MERGE( ind_lai_r_gfl,        ind_lai_r_agfl,        &
+              ind_green_frac_r = MERGE( ind_green_frac_r_gfl,                  &
+                                        ind_green_frac_r_agfl,                 &
                                         surf_usm_h%ground_level(m) )
-              ind_z0           = MERGE( ind_z0_gfl,           ind_z0_agfl,           &
+              ind_win_frac     = MERGE( ind_win_frac_gfl,                      &
+                                        ind_win_frac_agfl,                     &
                                         surf_usm_h%ground_level(m) )
-              ind_z0qh         = MERGE( ind_z0qh_gfl,         ind_z0qh_agfl,         &
+              ind_lai_r        = MERGE( ind_lai_r_gfl,                         &
+                                        ind_lai_r_agfl,                        &
                                         surf_usm_h%ground_level(m) )
+              ind_z0           = MERGE( ind_z0_gfl,                            &
+                                        ind_z0_agfl,                           &
+                                        surf_usm_h%ground_level(m) )
+              ind_z0qh         = MERGE( ind_z0qh_gfl,                          &
+                                        ind_z0qh_agfl,                         &
+                                        surf_usm_h%ground_level(m) )
+              ind_hc1          = MERGE( ind_hc1_gfl,                           &
+                                        ind_hc1_agfl,                          &
+                                        surf_usm_h%ground_level(m) )
+              ind_hc2          = MERGE( ind_hc2_gfl,                           &
+                                        ind_hc2_agfl,                          &
+                                        surf_usm_h%ground_level(m) )
+              ind_hc3          = MERGE( ind_hc3_gfl,                           &
+                                        ind_hc3_agfl,                          &
+                                        surf_usm_h%ground_level(m) )
+              ind_tc1          = MERGE( ind_tc1_gfl,                           &
+                                        ind_tc1_agfl,                          &
+                                        surf_usm_h%ground_level(m) )
+              ind_tc2          = MERGE( ind_tc2_gfl,                           &
+                                        ind_tc2_agfl,                          &
+                                        surf_usm_h%ground_level(m) )
+              ind_tc3          = MERGE( ind_tc3_gfl,                           &
+                                        ind_tc3_agfl,                          &
+                                        surf_usm_h%ground_level(m) )
+              ind_emis_wall    = MERGE( ind_emis_wall_gfl,                     &
+                                        ind_emis_wall_agfl,                    &
+                                        surf_usm_h%ground_level(m) )
+              ind_emis_green   = MERGE( ind_emis_green_gfl,                    &
+                                        ind_emis_green_agfl,                   &
+                                        surf_usm_h%ground_level(m) )
+              ind_emis_win     = MERGE( ind_emis_win_gfl,                      &
+                                        ind_emis_win_agfl,                     &
+                                        surf_usm_h%ground_level(m) )
+              ind_trans        = MERGE( ind_trans_gfl,                         &
+                                        ind_trans_agfl,                        &
+                                        surf_usm_h%ground_level(m) )
 
 !
 !--           Initialize relatvie wall- (0), green- (1) and window (2) fractions
-              IF ( building_pars_f%pars_xy(ind_wall_frac_r,j,i) /= building_pars_f%fill )    &
-                 surf_usm_h%frac(ind_veg_wall,m)  = building_pars_f%pars_xy(ind_wall_frac_r,j,i)   
-              IF ( building_pars_f%pars_xy(ind_green_frac_r,j,i) /= building_pars_f%fill )   & 
-                 surf_usm_h%frac(ind_pav_green,m) = building_pars_f%pars_xy(ind_green_frac_r,j,i) 
-              IF ( building_pars_f%pars_xy(ind_win_frac_r,j,i) /= building_pars_f%fill )     & 
-                 surf_usm_h%frac(ind_wat_win,m)   = building_pars_f%pars_xy(ind_win_frac_r,j,i)
-
+              IF ( building_pars_f%pars_xy(ind_wall_frac,j,i) /=               &
+                   building_pars_f%fill )                                      &
+                 surf_usm_h%frac(ind_veg_wall,m)  =                            &
+                                    building_pars_f%pars_xy(ind_wall_frac,j,i)   
+                 
+              IF ( building_pars_f%pars_xy(ind_green_frac_r,j,i) /=            &          
+                   building_pars_f%fill )                                      & 
+                 surf_usm_h%frac(ind_pav_green,m) =                            &
+                                    building_pars_f%pars_xy(ind_green_frac_r,j,i) 
+                 
+              IF ( building_pars_f%pars_xy(ind_win_frac,j,i) /=                &
+                   building_pars_f%fill )                                      & 
+                 surf_usm_h%frac(ind_wat_win,m)   =                            &
+                                    building_pars_f%pars_xy(ind_win_frac,j,i)
  
-              IF ( building_pars_f%pars_xy(ind_lai_r,j,i) /= building_pars_f%fill )          &
-                 surf_usm_h%lai(m)             = building_pars_f%pars_xy(ind_lai_r,j,i)
-
-              IF ( building_pars_f%pars_xy(ind_hc1_wall_r,j,i) /= building_pars_f%fill )  THEN 
-                 surf_usm_h%rho_c_wall(nzb_wall,m)   = building_pars_f%pars_xy(ind_hc1_wall_r,j,i)  
-                 surf_usm_h%rho_c_wall(nzb_wall+1,m) = building_pars_f%pars_xy(ind_hc1_wall_r,j,i)
+              IF ( building_pars_f%pars_xy(ind_lai_r,j,i) /=                   &
+                   building_pars_f%fill )                                      &
+                 surf_usm_h%lai(m)  = building_pars_f%pars_xy(ind_lai_r,j,i)
+
+              IF ( building_pars_f%pars_xy(ind_hc1,j,i) /=                     &
+                   building_pars_f%fill )  THEN
+                 surf_usm_h%rho_c_wall(nzb_wall,m)   =                         &
+                                    building_pars_f%pars_xy(ind_hc1,j,i)  
+                 surf_usm_h%rho_c_wall(nzb_wall+1,m) =                         &
+                                    building_pars_f%pars_xy(ind_hc1,j,i)
               ENDIF
-              IF ( building_pars_f%pars_xy(ind_hc2_wall_r,j,i) /= building_pars_f%fill )    &
-                 surf_usm_h%rho_c_wall(nzb_wall+2,m) = building_pars_f%pars_xy(ind_hc2_wall_r,j,i)
-              IF ( building_pars_f%pars_xy(ind_hc3_wall_r,j,i) /= building_pars_f%fill )    & 
-                 surf_usm_h%rho_c_wall(nzb_wall+3,m) = building_pars_f%pars_xy(ind_hc3_wall_r,j,i)
-              IF ( building_pars_f%pars_xy(ind_hc1_wall_r,j,i) /= building_pars_f%fill )  THEN 
-                 surf_usm_h%rho_c_green(nzb_wall,m)   = rho_c_soil !building_pars_f%pars_xy(ind_hc1_wall_r,j,i)  
-                 surf_usm_h%rho_c_green(nzb_wall+1,m) = rho_c_soil !building_pars_f%pars_xy(ind_hc1_wall_r,j,i)
+              
+              
+              IF ( building_pars_f%pars_xy(ind_hc2,j,i) /=                     &
+                   building_pars_f%fill )                                      &
+                 surf_usm_h%rho_c_wall(nzb_wall+2,m) =                         &
+                                    building_pars_f%pars_xy(ind_hc2,j,i)
+                 
+              IF ( building_pars_f%pars_xy(ind_hc3,j,i) /=                     &
+                   building_pars_f%fill )                                      &
+                 surf_usm_h%rho_c_wall(nzb_wall+3,m) =                         &
+                                    building_pars_f%pars_xy(ind_hc3,j,i)
+                 
+              IF ( building_pars_f%pars_xy(ind_hc1,j,i) /=                     &
+                   building_pars_f%fill )  THEN 
+                 surf_usm_h%rho_c_green(nzb_wall,m)   =                        &
+                                    building_pars_f%pars_xy(ind_hc1,j,i)  
+                 surf_usm_h%rho_c_green(nzb_wall+1,m) =                        &
+                                    building_pars_f%pars_xy(ind_hc1,j,i)
               ENDIF
-              IF ( building_pars_f%pars_xy(ind_hc2_wall_r,j,i) /= building_pars_f%fill )    &
-                 surf_usm_h%rho_c_green(nzb_wall+2,m) = rho_c_soil !building_pars_f%pars_xy(ind_hc2_wall_r,j,i)
-              IF ( building_pars_f%pars_xy(ind_hc3_wall_r,j,i) /= building_pars_f%fill )    & 
-                 surf_usm_h%rho_c_green(nzb_wall+3,m) = rho_c_soil !building_pars_f%pars_xy(ind_hc3_wall_r,j,i)
-              IF ( building_pars_f%pars_xy(ind_hc1_win_r,j,i) /= building_pars_f%fill )  THEN 
-                 surf_usm_h%rho_c_window(nzb_wall,m)   = building_pars_f%pars_xy(ind_hc1_win_r,j,i)  
-                 surf_usm_h%rho_c_window(nzb_wall+1,m) = building_pars_f%pars_xy(ind_hc1_win_r,j,i)
+              IF ( building_pars_f%pars_xy(ind_hc2,j,i) /=                     &
+                   building_pars_f%fill )                                      &
+                 surf_usm_h%rho_c_green(nzb_wall+2,m) =                        &
+                                    building_pars_f%pars_xy(ind_hc2,j,i)
+                 
+              IF ( building_pars_f%pars_xy(ind_hc3,j,i) /=                     &
+                   building_pars_f%fill )                                      &
+                 surf_usm_h%rho_c_green(nzb_wall+3,m) =                        &
+                                    building_pars_f%pars_xy(ind_hc3,j,i)
+                 
+              IF ( building_pars_f%pars_xy(ind_hc1,j,i) /=                     &
+                   building_pars_f%fill )  THEN 
+                 surf_usm_h%rho_c_window(nzb_wall,m)   =                       &
+                                    building_pars_f%pars_xy(ind_hc1,j,i)  
+                 surf_usm_h%rho_c_window(nzb_wall+1,m) =                       &
+                                    building_pars_f%pars_xy(ind_hc1,j,i)
               ENDIF
-              IF ( building_pars_f%pars_xy(ind_hc2_win_r,j,i) /= building_pars_f%fill )    &
-                 surf_usm_h%rho_c_window(nzb_wall+2,m) = building_pars_f%pars_xy(ind_hc2_win_r,j,i)
-              IF ( building_pars_f%pars_xy(ind_hc3_win_r,j,i) /= building_pars_f%fill )    & 
-                 surf_usm_h%rho_c_window(nzb_wall+3,m) = building_pars_f%pars_xy(ind_hc3_win_r,j,i)
-
-              IF ( building_pars_f%pars_xy(ind_tc1_wall_r,j,i) /= building_pars_f%fill )  THEN
-                 surf_usm_h%lambda_h(nzb_wall,m)   = building_pars_f%pars_xy(ind_tc1_wall_r,j,i)         
-                 surf_usm_h%lambda_h(nzb_wall+1,m) = building_pars_f%pars_xy(ind_tc1_wall_r,j,i)        
+              IF ( building_pars_f%pars_xy(ind_hc2,j,i) /=                     &
+                   building_pars_f%fill )                                      &
+                 surf_usm_h%rho_c_window(nzb_wall+2,m) =                       &
+                                    building_pars_f%pars_xy(ind_hc2,j,i)
+                 
+              IF ( building_pars_f%pars_xy(ind_hc3,j,i) /=                     &
+                   building_pars_f%fill )                                      &
+                 surf_usm_h%rho_c_window(nzb_wall+3,m) =                       &
+                                    building_pars_f%pars_xy(ind_hc3,j,i)
+
+              IF ( building_pars_f%pars_xy(ind_tc1,j,i) /=                     &
+                   building_pars_f%fill )  THEN
+                 surf_usm_h%lambda_h(nzb_wall,m)   =                           &
+                                    building_pars_f%pars_xy(ind_tc1,j,i)         
+                 surf_usm_h%lambda_h(nzb_wall+1,m) =                           &
+                                    building_pars_f%pars_xy(ind_tc1,j,i)        
               ENDIF
-              IF ( building_pars_f%pars_xy(ind_tc2_wall_r,j,i) /= building_pars_f%fill )    &
-                 surf_usm_h%lambda_h(nzb_wall+2,m) = building_pars_f%pars_xy(ind_tc2_wall_r,j,i)
-              IF ( building_pars_f%pars_xy(ind_tc3_wall_r,j,i) /= building_pars_f%fill )    & 
-                 surf_usm_h%lambda_h(nzb_wall+3,m) = building_pars_f%pars_xy(ind_tc3_wall_r,j,i)    
-              IF ( building_pars_f%pars_xy(ind_tc1_wall_r,j,i) /= building_pars_f%fill )  THEN
-                 surf_usm_h%lambda_h_green(nzb_wall,m)   = lambda_h_green_sm !building_pars_f%pars_xy(ind_tc1_wall_r,j,i)         
-                 surf_usm_h%lambda_h_green(nzb_wall+1,m) = lambda_h_green_sm !building_pars_f%pars_xy(ind_tc1_wall_r,j,i)        
+              IF ( building_pars_f%pars_xy(ind_tc2,j,i) /=                     &  
+                   building_pars_f%fill )                                      &
+                 surf_usm_h%lambda_h(nzb_wall+2,m) =                           &
+                                    building_pars_f%pars_xy(ind_tc2,j,i)
+                 
+              IF ( building_pars_f%pars_xy(ind_tc3,j,i) /=                     &
+                   building_pars_f%fill )                                      & 
+                 surf_usm_h%lambda_h(nzb_wall+3,m) =                           &
+                                    building_pars_f%pars_xy(ind_tc3,j,i)    
+                 
+              IF ( building_pars_f%pars_xy(ind_tc1,j,i) /=                     &
+                   building_pars_f%fill )  THEN
+                 surf_usm_h%lambda_h_green(nzb_wall,m)   =                     &
+                                     building_pars_f%pars_xy(ind_tc1,j,i)         
+                 surf_usm_h%lambda_h_green(nzb_wall+1,m) =                     &
+                                     building_pars_f%pars_xy(ind_tc1,j,i)        
               ENDIF
-              IF ( building_pars_f%pars_xy(ind_tc2_wall_r,j,i) /= building_pars_f%fill )    &
-                 surf_usm_h%lambda_h_green(nzb_wall+2,m) = lambda_h_green_sm !building_pars_f%pars_xy(ind_tc2_wall_r,j,i)
-              IF ( building_pars_f%pars_xy(ind_tc3_wall_r,j,i) /= building_pars_f%fill )    & 
-                 surf_usm_h%lambda_h_green(nzb_wall+3,m) = lambda_h_green_sm !building_pars_f%pars_xy(ind_tc3_wall_r,j,i)    
-              IF ( building_pars_f%pars_xy(ind_tc1_win_r,j,i) /= building_pars_f%fill )  THEN
-                 surf_usm_h%lambda_h_window(nzb_wall,m)   = building_pars_f%pars_xy(ind_tc1_win_r,j,i)         
-                 surf_usm_h%lambda_h_window(nzb_wall+1,m) = building_pars_f%pars_xy(ind_tc1_win_r,j,i)        
+              IF ( building_pars_f%pars_xy(ind_tc2,j,i) /=                     & 
+                   building_pars_f%fill )                                      &
+                 surf_usm_h%lambda_h_green(nzb_wall+2,m) =                     &
+                                    building_pars_f%pars_xy(ind_tc2,j,i)
+                 
+              IF ( building_pars_f%pars_xy(ind_tc3,j,i) /=                     &       
+                   building_pars_f%fill )                                      &
+                 surf_usm_h%lambda_h_green(nzb_wall+3,m) =                     &
+                                    building_pars_f%pars_xy(ind_tc3,j,i)    
+                 
+              IF ( building_pars_f%pars_xy(ind_tc1_win_r,j,i) /=               &
+                   building_pars_f%fill )  THEN
+                 surf_usm_h%lambda_h_window(nzb_wall,m)   =                    &
+                                     building_pars_f%pars_xy(ind_tc1_win_r,j,i)         
+                 surf_usm_h%lambda_h_window(nzb_wall+1,m) =                    &
+                                     building_pars_f%pars_xy(ind_tc1_win_r,j,i)        
               ENDIF
-              IF ( building_pars_f%pars_xy(ind_tc2_win_r,j,i) /= building_pars_f%fill )    &
-                 surf_usm_h%lambda_h_window(nzb_wall+2,m) = building_pars_f%pars_xy(ind_tc2_win_r,j,i)
-              IF ( building_pars_f%pars_xy(ind_tc3_win_r,j,i) /= building_pars_f%fill )    & 
-                 surf_usm_h%lambda_h_window(nzb_wall+3,m) = building_pars_f%pars_xy(ind_tc3_win_r,j,i)    
-
-              IF ( building_pars_f%pars_xy(117,j,i) /= building_pars_f%fill )         & 
-                 surf_usm_h%target_temp_summer(m)  = building_pars_f%pars_xy(117,j,i)   
-              IF ( building_pars_f%pars_xy(118,j,i) /= building_pars_f%fill )         & 
-                 surf_usm_h%target_temp_winter(m)  = building_pars_f%pars_xy(118,j,i)   
-
-              IF ( building_pars_f%pars_xy(ind_emis_wall_r,j,i) /= building_pars_f%fill ) & 
-                 surf_usm_h%emissivity(ind_veg_wall,m)  = building_pars_f%pars_xy(ind_emis_wall_r,j,i)
-              IF ( building_pars_f%pars_xy(ind_emis_green_r,j,i) /= building_pars_f%fill )& 
-                 surf_usm_h%emissivity(ind_pav_green,m) = building_pars_f%pars_xy(ind_emis_green_r,j,i)
-              IF ( building_pars_f%pars_xy(ind_emis_win_r,j,i) /= building_pars_f%fill )  & 
-                 surf_usm_h%emissivity(ind_wat_win,m)   = building_pars_f%pars_xy(ind_emis_win_r,j,i)
-
-              IF ( building_pars_f%pars_xy(ind_trans_r,j,i) /= building_pars_f%fill )    & 
-                 surf_usm_h%transmissivity(m)      = building_pars_f%pars_xy(ind_trans_r,j,i)
-
-              IF ( building_pars_f%pars_xy(ind_z0,j,i) /= building_pars_f%fill )    & 
-                 surf_usm_h%z0(m)                  = building_pars_f%pars_xy(ind_z0,j,i)
-              IF ( building_pars_f%pars_xy(ind_z0qh,j,i) /= building_pars_f%fill )  & 
-                 surf_usm_h%z0h(m)                 = building_pars_f%pars_xy(ind_z0qh,j,i)
-              IF ( building_pars_f%pars_xy(ind_z0qh,j,i) /= building_pars_f%fill )  & 
-                 surf_usm_h%z0q(m)                 = building_pars_f%pars_xy(ind_z0qh,j,i)
-
-              IF ( building_pars_f%pars_xy(ind_alb_wall_r,j,i) /= building_pars_f%fill )    & 
-                 surf_usm_h%albedo_type(ind_veg_wall,m)  = building_pars_f%pars_xy(ind_alb_wall_r,j,i)
-              IF ( building_pars_f%pars_xy(ind_alb_green_r,j,i) /= building_pars_f%fill )   & 
-                 surf_usm_h%albedo_type(ind_pav_green,m) = building_pars_f%pars_xy(ind_alb_green_r,j,i)
-              IF ( building_pars_f%pars_xy(ind_alb_win_r,j,i) /= building_pars_f%fill )     & 
-                 surf_usm_h%albedo_type(ind_wat_win,m)   = building_pars_f%pars_xy(ind_alb_win_r,j,i)
-
-              IF ( building_pars_f%pars_xy(ind_thick_1_wall_r,j,i) /= building_pars_f%fill )    & 
-                 surf_usm_h%zw(nzb_wall,m) = building_pars_f%pars_xy(ind_thick_1_wall_r,j,i)
-              IF ( building_pars_f%pars_xy(ind_thick_2_wall_r,j,i) /= building_pars_f%fill )    & 
-                 surf_usm_h%zw(nzb_wall+1,m) = building_pars_f%pars_xy(ind_thick_2_wall_r,j,i)
-              IF ( building_pars_f%pars_xy(ind_thick_3_wall_r,j,i) /= building_pars_f%fill )    & 
-                 surf_usm_h%zw(nzb_wall+2,m) = building_pars_f%pars_xy(ind_thick_3_wall_r,j,i)
-              IF ( building_pars_f%pars_xy(ind_thick_4_wall_r,j,i) /= building_pars_f%fill )    & 
-                 surf_usm_h%zw(nzb_wall+3,m) = building_pars_f%pars_xy(ind_thick_4_wall_r,j,i)
-              IF ( building_pars_f%pars_xy(ind_thick_1_wall_r,j,i) /= building_pars_f%fill )    & 
-                 surf_usm_h%zw_green(nzb_wall,m) = building_pars_f%pars_xy(ind_thick_1_wall_r,j,i)
-              IF ( building_pars_f%pars_xy(ind_thick_2_wall_r,j,i) /= building_pars_f%fill )    & 
-                 surf_usm_h%zw_green(nzb_wall+1,m) = building_pars_f%pars_xy(ind_thick_2_wall_r,j,i)
-              IF ( building_pars_f%pars_xy(ind_thick_3_wall_r,j,i) /= building_pars_f%fill )    & 
-                 surf_usm_h%zw_green(nzb_wall+2,m) = building_pars_f%pars_xy(ind_thick_3_wall_r,j,i)
-              IF ( building_pars_f%pars_xy(ind_thick_4_wall_r,j,i) /= building_pars_f%fill )    & 
-                 surf_usm_h%zw_green(nzb_wall+3,m) = building_pars_f%pars_xy(ind_thick_4_wall_r,j,i)
-              IF ( building_pars_f%pars_xy(ind_thick_1_win_r,j,i) /= building_pars_f%fill )     & 
-                 surf_usm_h%zw_window(nzb_wall,m) = building_pars_f%pars_xy(ind_thick_1_win_r,j,i)
-              IF ( building_pars_f%pars_xy(ind_thick_2_win_r,j,i) /= building_pars_f%fill )     & 
-                 surf_usm_h%zw_window(nzb_wall+1,m) = building_pars_f%pars_xy(ind_thick_2_win_r,j,i)
-              IF ( building_pars_f%pars_xy(ind_thick_3_win_r,j,i) /= building_pars_f%fill )     & 
-                 surf_usm_h%zw_window(nzb_wall+2,m) = building_pars_f%pars_xy(ind_thick_3_win_r,j,i)
-              IF ( building_pars_f%pars_xy(ind_thick_4_win_r,j,i) /= building_pars_f%fill )     & 
-                 surf_usm_h%zw_window(nzb_wall+3,m) = building_pars_f%pars_xy(ind_thick_4_win_r,j,i)
-
-              IF ( building_pars_f%pars_xy(0,j,i) /= building_pars_f%fill )    & 
-                 surf_usm_h%c_surface(m)           = building_pars_f%pars_xy(0,j,i)
-              IF ( building_pars_f%pars_xy(3,j,i) /= building_pars_f%fill )    & 
-                 surf_usm_h%lambda_surf(m)         = building_pars_f%pars_xy(3,j,i)
-              IF ( building_pars_f%pars_xy(2,j,i) /= building_pars_f%fill )    & 
-                 surf_usm_h%c_surface_green(m)           = building_pars_f%pars_xy(2,j,i)
-              IF ( building_pars_f%pars_xy(5,j,i) /= building_pars_f%fill )    & 
-                 surf_usm_h%lambda_surf_green(m)         = building_pars_f%pars_xy(5,j,i)
-              IF ( building_pars_f%pars_xy(1,j,i) /= building_pars_f%fill )    & 
-                 surf_usm_h%c_surface_window(m)           = building_pars_f%pars_xy(1,j,i)
-              IF ( building_pars_f%pars_xy(4,j,i) /= building_pars_f%fill )    & 
-                 surf_usm_h%lambda_surf_window(m)         = building_pars_f%pars_xy(4,j,i)
-
-              IF ( building_pars_f%pars_xy(ind_green_type_roof,j,i) /= building_pars_f%fill )    & 
-                 surf_usm_h%green_type_roof(m)            = building_pars_f%pars_xy(ind_green_type_roof,j,i)
+              IF ( building_pars_f%pars_xy(ind_tc2,j,i) /=                     &     
+                   building_pars_f%fill )                                      &
+                 surf_usm_h%lambda_h_window(nzb_wall+2,m) =                    &
+                                     building_pars_f%pars_xy(ind_tc2,j,i)
+                 
+              IF ( building_pars_f%pars_xy(ind_tc3,j,i) /=                     &   
+                   building_pars_f%fill )                                      &
+                 surf_usm_h%lambda_h_window(nzb_wall+3,m) =                    &
+                                    building_pars_f%pars_xy(ind_tc3,j,i)    
+
+              IF ( building_pars_f%pars_xy(ind_indoor_target_temp_summer,j,i) /=&            
+                   building_pars_f%fill )                                      & 
+                 surf_usm_h%target_temp_summer(m)  =                           &
+                      building_pars_f%pars_xy(ind_indoor_target_temp_summer,j,i)   
+              IF ( building_pars_f%pars_xy(ind_indoor_target_temp_winter,j,i) /=&            
+                   building_pars_f%fill )                                      & 
+                 surf_usm_h%target_temp_winter(m)  =                           &
+                      building_pars_f%pars_xy(ind_indoor_target_temp_winter,j,i)   
+
+              IF ( building_pars_f%pars_xy(ind_emis_wall,j,i) /=               &   
+                   building_pars_f%fill )                                      &
+                 surf_usm_h%emissivity(ind_veg_wall,m)  =                      &
+                                    building_pars_f%pars_xy(ind_emis_wall,j,i)
+                 
+              IF ( building_pars_f%pars_xy(ind_emis_green,j,i) /=              &            
+                   building_pars_f%fill )                                      &
+                 surf_usm_h%emissivity(ind_pav_green,m) =                      &
+                                     building_pars_f%pars_xy(ind_emis_green,j,i)
+                 
+              IF ( building_pars_f%pars_xy(ind_emis_win,j,i) /=                & 
+                   building_pars_f%fill )                                      &
+                 surf_usm_h%emissivity(ind_wat_win,m)   =                      &
+                                     building_pars_f%pars_xy(ind_emis_win,j,i)
+                 
+              IF ( building_pars_f%pars_xy(ind_trans,j,i) /=                   &    
+                   building_pars_f%fill )                                      &
+                 surf_usm_h%transmissivity(m) =                                &
+                                    building_pars_f%pars_xy(ind_trans,j,i)
+
+              IF ( building_pars_f%pars_xy(ind_z0,j,i) /=                      &          
+                   building_pars_f%fill )                                      &
+                 surf_usm_h%z0(m) = building_pars_f%pars_xy(ind_z0,j,i)
+                 
+              IF ( building_pars_f%pars_xy(ind_z0qh,j,i) /=                    &            
+                   building_pars_f%fill )                                      &
+                 surf_usm_h%z0h(m) = building_pars_f%pars_xy(ind_z0qh,j,i)
+              IF ( building_pars_f%pars_xy(ind_z0qh,j,i) /=                    &           
+                   building_pars_f%fill )                                      &
+                 surf_usm_h%z0q(m) = building_pars_f%pars_xy(ind_z0qh,j,i)
+
+              IF ( building_pars_f%pars_xy(ind_alb_wall_agfl,j,i) /=           &          
+                   building_pars_f%fill )                                      & 
+                 surf_usm_h%albedo_type(ind_veg_wall,m)  =                     &
+                                 building_pars_f%pars_xy(ind_alb_wall_agfl,j,i)
+                 
+              IF ( building_pars_f%pars_xy(ind_alb_green_agfl,j,i) /=          &           
+                   building_pars_f%fill )                                      &
+                 surf_usm_h%albedo_type(ind_pav_green,m) =                     &
+                                building_pars_f%pars_xy(ind_alb_green_agfl,j,i)
+              IF ( building_pars_f%pars_xy(ind_alb_win_agfl,j,i) /=            &          
+                   building_pars_f%fill )                                      &
+                 surf_usm_h%albedo_type(ind_wat_win,m)   =                     &
+                                   building_pars_f%pars_xy(ind_alb_win_agfl,j,i)
+
+              IF ( building_pars_f%pars_xy(ind_thick_1_agfl,j,i) /=            &          
+                   building_pars_f%fill )                                      & 
+                 surf_usm_h%zw(nzb_wall,m) =                                   &
+                                  building_pars_f%pars_xy(ind_thick_1_agfl,j,i)
+                 
+              IF ( building_pars_f%pars_xy(ind_thick_2_agfl,j,i) /=            &          
+                   building_pars_f%fill )                                      &
+                 surf_usm_h%zw(nzb_wall+1,m) =                                 &
+                                  building_pars_f%pars_xy(ind_thick_2_agfl,j,i)
+                 
+              IF ( building_pars_f%pars_xy(ind_thick_3_agfl,j,i) /=            &         
+                   building_pars_f%fill )                                      &
+                 surf_usm_h%zw(nzb_wall+2,m) =                                 &
+                                  building_pars_f%pars_xy(ind_thick_3_agfl,j,i)
+                 
+                 
+              IF ( building_pars_f%pars_xy(ind_thick_4_agfl,j,i) /=            &         
+                   building_pars_f%fill )                                      & 
+                 surf_usm_h%zw(nzb_wall+3,m) =                                 &
+                                  building_pars_f%pars_xy(ind_thick_4_agfl,j,i)
+                 
+              IF ( building_pars_f%pars_xy(ind_thick_1_agfl,j,i) /=            &           
+                   building_pars_f%fill )                                      &
+                 surf_usm_h%zw_green(nzb_wall,m) =                             &
+                                  building_pars_f%pars_xy(ind_thick_1_agfl,j,i)
+                 
+              IF ( building_pars_f%pars_xy(ind_thick_2_agfl,j,i) /=            &          
+                   building_pars_f%fill )                                      &
+                 surf_usm_h%zw_green(nzb_wall+1,m) =                           &
+                                   building_pars_f%pars_xy(ind_thick_2_agfl,j,i)
+                 
+              IF ( building_pars_f%pars_xy(ind_thick_3_agfl,j,i) /=            &          
+                   building_pars_f%fill )                                      & 
+                 surf_usm_h%zw_green(nzb_wall+2,m) =                           &
+                                   building_pars_f%pars_xy(ind_thick_3_agfl,j,i)
+                 
+              IF ( building_pars_f%pars_xy(ind_thick_4_agfl,j,i) /=            &          
+                   building_pars_f%fill )                                      &
+                 surf_usm_h%zw_green(nzb_wall+3,m) =                           &
+                                   building_pars_f%pars_xy(ind_thick_4_agfl,j,i)
+
+              IF ( building_pars_f%pars_xy(ind_c_surface,j,i) /=               &       
+                   building_pars_f%fill )                                      & 
+                 surf_usm_h%c_surface(m) =                                     &
+                                    building_pars_f%pars_xy(ind_c_surface,j,i)
+                 
+              IF ( building_pars_f%pars_xy(ind_lambda_surf,j,i) /=             &        
+                   building_pars_f%fill )                                      &
+                 surf_usm_h%lambda_surf(m) =                                   &
+                                    building_pars_f%pars_xy(ind_lambda_surf,j,i)
+              
            ENDDO
 
@@ -4847,219 +5590,288 @@
                 
 !
-!--              In order to distinguish between ground floor level and 
-!--              above-ground-floor level surfaces, set input indices.
-                 ind_alb_green    = MERGE( ind_alb_green_gfl,    ind_alb_green_agfl,    &
-                                           surf_usm_v(l)%ground_level(m) )
-                 ind_alb_wall     = MERGE( ind_alb_wall_gfl,     ind_alb_wall_agfl,     &
-                                           surf_usm_v(l)%ground_level(m) )
-                 ind_alb_win      = MERGE( ind_alb_win_gfl,      ind_alb_win_agfl,      &
-                                           surf_usm_v(l)%ground_level(m) )
-                 ind_wall_frac    = MERGE( ind_wall_frac_gfl,    ind_wall_frac_agfl,    &
-                                           surf_usm_v(l)%ground_level(m) )
-                 ind_win_frac     = MERGE( ind_win_frac_gfl,     ind_win_frac_agfl,     &
-                                           surf_usm_v(l)%ground_level(m) )
-                 ind_green_frac_w = MERGE( ind_green_frac_w_gfl, ind_green_frac_w_agfl, &
-                                           surf_usm_v(l)%ground_level(m) )
-                 ind_green_frac_r = MERGE( ind_green_frac_r_gfl, ind_green_frac_r_agfl, &
-                                           surf_usm_v(l)%ground_level(m) )
-                 ind_lai_r        = MERGE( ind_lai_r_gfl,        ind_lai_r_agfl,        &
-                                           surf_usm_v(l)%ground_level(m) )
-                 ind_lai_w        = MERGE( ind_lai_w_gfl,        ind_lai_w_agfl,        &
-                                           surf_usm_v(l)%ground_level(m) )
-                 ind_hc1          = MERGE( ind_hc1_gfl,          ind_hc1_agfl,          &
-                                           surf_usm_v(l)%ground_level(m) )
-                 ind_hc1_win      = MERGE( ind_hc1_win_gfl,      ind_hc1_win_agfl,      &
-                                           surf_usm_v(l)%ground_level(m) )
-                 ind_hc2          = MERGE( ind_hc2_gfl,          ind_hc2_agfl,          &
-                                           surf_usm_v(l)%ground_level(m) )
-                 ind_hc2_win      = MERGE( ind_hc2_win_gfl,      ind_hc2_win_agfl,      &
-                                           surf_usm_v(l)%ground_level(m) )
-                 ind_hc3          = MERGE( ind_hc3_gfl,          ind_hc3_agfl,          &
-                                           surf_usm_v(l)%ground_level(m) )
-                 ind_hc3_win      = MERGE( ind_hc3_win_gfl,      ind_hc3_win_agfl,      &
-                                           surf_usm_v(l)%ground_level(m) )
-                 ind_tc1          = MERGE( ind_tc1_gfl,          ind_tc1_agfl,          &
-                                           surf_usm_v(l)%ground_level(m) )
-                 ind_tc1_win      = MERGE( ind_tc1_win_gfl,      ind_tc1_win_agfl,      &
-                                           surf_usm_v(l)%ground_level(m) )
-                 ind_tc2          = MERGE( ind_tc2_gfl,          ind_tc2_agfl,          &
-                                           surf_usm_v(l)%ground_level(m) )
-                 ind_tc2_win      = MERGE( ind_tc2_win_gfl,      ind_tc2_win_agfl,      &
-                                           surf_usm_v(l)%ground_level(m) )
-                 ind_tc3          = MERGE( ind_tc3_gfl,          ind_tc3_agfl,          &
-                                           surf_usm_v(l)%ground_level(m) )
-                 ind_tc3_win      = MERGE( ind_tc3_win_gfl,      ind_tc3_win_agfl,      &
-                                           surf_usm_v(l)%ground_level(m) )
-                 ind_thick_1      = MERGE( ind_thick_1_gfl,      ind_thick_1_agfl,      &
-                                           surf_usm_v(l)%ground_level(m) )
-                 ind_thick_1_win  = MERGE( ind_thick_1_win_gfl,  ind_thick_1_win_agfl,  &
-                                           surf_usm_v(l)%ground_level(m) )
-                 ind_thick_2      = MERGE( ind_thick_2_gfl,      ind_thick_2_agfl,      &
-                                           surf_usm_v(l)%ground_level(m) )
-                 ind_thick_2_win  = MERGE( ind_thick_2_win_gfl,  ind_thick_2_win_agfl,  &
-                                           surf_usm_v(l)%ground_level(m) )
-                 ind_thick_3      = MERGE( ind_thick_3_gfl,      ind_thick_3_agfl,      &
-                                           surf_usm_v(l)%ground_level(m) )
-                 ind_thick_3_win  = MERGE( ind_thick_3_win_gfl,  ind_thick_3_win_agfl,  &
-                                           surf_usm_v(l)%ground_level(m) )
-                 ind_thick_4      = MERGE( ind_thick_4_gfl,      ind_thick_4_agfl,      &
-                                           surf_usm_v(l)%ground_level(m) )
-                 ind_thick_4_win  = MERGE( ind_thick_4_win_gfl,  ind_thick_4_win_agfl,  &
-                                           surf_usm_v(l)%ground_level(m) )
-                 ind_emis_wall    = MERGE( ind_emis_wall_gfl,    ind_emis_wall_agfl,    &
-                                           surf_usm_v(l)%ground_level(m) )
-                 ind_emis_green   = MERGE( ind_emis_green_gfl,   ind_emis_green_agfl,   &
-                                           surf_usm_v(l)%ground_level(m) )
-                 ind_emis_win     = MERGE( ind_emis_win_gfl,     ind_emis_win_agfl,     &
-                                           surf_usm_v(l)%ground_level(m) )
-                 ind_trans        = MERGE( ind_trans_gfl,       ind_trans_agfl,         &
-                                           surf_usm_v(l)%ground_level(m) )
-                 ind_z0           = MERGE( ind_z0_gfl,           ind_z0_agfl,           &
-                                           surf_usm_v(l)%ground_level(m) )
-                 ind_z0qh         = MERGE( ind_z0qh_gfl,         ind_z0qh_agfl,         &
+!--                 In order to distinguish between ground floor level and 
+!--                 above-ground-floor level surfaces, set input indices.
+                    ind_wall_frac    = MERGE( ind_wall_frac_gfl,               &
+                                              ind_wall_frac_agfl,              &
                                               surf_usm_v(l)%ground_level(m) )
-
-!
-!--              Initialize relatvie wall- (0), green- (1) and window (2) fractions
-                 IF ( building_pars_f%pars_xy(ind_wall_frac,j,i) /=                     &
-                      building_pars_f%fill )                                            &
-                    surf_usm_v(l)%frac(ind_veg_wall,m)  =                               &
-                                      building_pars_f%pars_xy(ind_wall_frac,j,i)   
-                 IF ( building_pars_f%pars_xy(ind_green_frac_w,j,i) /=                  &
-                      building_pars_f%fill )                                            & 
-                    surf_usm_v(l)%frac(ind_pav_green,m) =                               &
-                                      building_pars_f%pars_xy(ind_green_frac_w,j,i) 
-                 IF ( building_pars_f%pars_xy(ind_win_frac,j,i) /=                      &
-                      building_pars_f%fill )                                            & 
-                    surf_usm_v(l)%frac(ind_wat_win,m)   =                               &
-                                      building_pars_f%pars_xy(ind_win_frac,j,i)
-  
-                 IF ( building_pars_f%pars_xy(ind_lai_w,j,i) /= building_pars_f%fill )  & 
-                    surf_usm_v(l)%lai(m) = building_pars_f%pars_xy(ind_lai_w,j,i)
-
-                 IF ( building_pars_f%pars_xy(ind_hc1,j,i) /= building_pars_f%fill )    &
-                 THEN 
-                    surf_usm_v(l)%rho_c_wall(nzb_wall,m)   =                            &
-                                                    building_pars_f%pars_xy(ind_hc1,j,i)  
-                    surf_usm_v(l)%rho_c_wall(nzb_wall+1,m) =                            &
-                                                    building_pars_f%pars_xy(ind_hc1,j,i)
-                 ENDIF
-                 IF ( building_pars_f%pars_xy(ind_hc2,j,i) /= building_pars_f%fill )    &
-                    surf_usm_v(l)%rho_c_wall(nzb_wall+2,m) =                            &                            
-                                                    building_pars_f%pars_xy(ind_hc2,j,i)
-                 IF ( building_pars_f%pars_xy(ind_hc3,j,i) /= building_pars_f%fill )    & 
-                    surf_usm_v(l)%rho_c_wall(nzb_wall+3,m) =                            &
-                                                    building_pars_f%pars_xy(ind_hc3,j,i)
-                 IF ( building_pars_f%pars_xy(ind_hc1,j,i) /= building_pars_f%fill )  THEN 
-                    surf_usm_v(l)%rho_c_green(nzb_wall,m)   =                           &
-                                                    rho_c_soil !building_pars_f%pars_xy(ind_hc1,j,i)  
-                    surf_usm_v(l)%rho_c_green(nzb_wall+1,m) =                           &
-                                                    rho_c_soil !building_pars_f%pars_xy(ind_hc1,j,i)
-                 ENDIF
-                 IF ( building_pars_f%pars_xy(ind_hc2,j,i) /= building_pars_f%fill )    &
-                    surf_usm_v(l)%rho_c_green(nzb_wall+2,m) = rho_c_soil !building_pars_f%pars_xy(ind_hc2,j,i)
-                 IF ( building_pars_f%pars_xy(ind_hc3,j,i) /= building_pars_f%fill )    & 
-                    surf_usm_v(l)%rho_c_green(nzb_wall+3,m) = rho_c_soil !building_pars_f%pars_xy(ind_hc3,j,i)    
-                 IF ( building_pars_f%pars_xy(ind_hc1_win,j,i) /= building_pars_f%fill )  THEN 
-                    surf_usm_v(l)%rho_c_window(nzb_wall,m)   = building_pars_f%pars_xy(ind_hc1_win,j,i)  
-                    surf_usm_v(l)%rho_c_window(nzb_wall+1,m) = building_pars_f%pars_xy(ind_hc1_win,j,i)
-                 ENDIF
-                 IF ( building_pars_f%pars_xy(ind_hc2_win,j,i) /= building_pars_f%fill )    &
-                    surf_usm_v(l)%rho_c_window(nzb_wall+2,m) = building_pars_f%pars_xy(ind_hc2_win,j,i)
-                 IF ( building_pars_f%pars_xy(ind_hc3_win,j,i) /= building_pars_f%fill )    &
-                    surf_usm_v(l)%rho_c_window(nzb_wall+3,m) = building_pars_f%pars_xy(ind_hc3_win,j,i)
-
-                 IF ( building_pars_f%pars_xy(ind_tc1,j,i) /= building_pars_f%fill )  THEN
-                    surf_usm_v(l)%lambda_h(nzb_wall,m)   = building_pars_f%pars_xy(ind_tc1,j,i)         
-                    surf_usm_v(l)%lambda_h(nzb_wall+1,m) = building_pars_f%pars_xy(ind_tc1,j,i)        
-                 ENDIF
-                 IF ( building_pars_f%pars_xy(ind_tc2,j,i) /= building_pars_f%fill )    &
-                    surf_usm_v(l)%lambda_h(nzb_wall+2,m) = building_pars_f%pars_xy(ind_tc2,j,i)
-                 IF ( building_pars_f%pars_xy(ind_tc3,j,i) /= building_pars_f%fill )    & 
-                    surf_usm_v(l)%lambda_h(nzb_wall+3,m) = building_pars_f%pars_xy(ind_tc3,j,i)    
-                 IF ( building_pars_f%pars_xy(ind_tc1,j,i) /= building_pars_f%fill )  THEN
-                    surf_usm_v(l)%lambda_h_green(nzb_wall,m)   = lambda_h_green_sm !building_pars_f%pars_xy(ind_tc1,j,i)         
-                    surf_usm_v(l)%lambda_h_green(nzb_wall+1,m) = lambda_h_green_sm !building_pars_f%pars_xy(ind_tc1,j,i)        
-                 ENDIF
-                 IF ( building_pars_f%pars_xy(ind_tc2,j,i) /= building_pars_f%fill )    &
-                    surf_usm_v(l)%lambda_h_green(nzb_wall+2,m) = lambda_h_green_sm !building_pars_f%pars_xy(ind_tc2,j,i)
-                 IF ( building_pars_f%pars_xy(ind_tc3,j,i) /= building_pars_f%fill )    & 
-                    surf_usm_v(l)%lambda_h_green(nzb_wall+3,m) = lambda_h_green_sm !building_pars_f%pars_xy(ind_tc3,j,i)    
-                 IF ( building_pars_f%pars_xy(ind_tc1_win,j,i) /= building_pars_f%fill )  THEN
-                    surf_usm_v(l)%lambda_h_window(nzb_wall,m)   = building_pars_f%pars_xy(ind_tc1_win,j,i)         
-                    surf_usm_v(l)%lambda_h_window(nzb_wall+1,m) = building_pars_f%pars_xy(ind_tc1_win,j,i)        
-                 ENDIF
-                 IF ( building_pars_f%pars_xy(ind_tc2_win,j,i) /= building_pars_f%fill )    &
-                    surf_usm_v(l)%lambda_h_window(nzb_wall+2,m) = building_pars_f%pars_xy(ind_tc2_win,j,i)
-                 IF ( building_pars_f%pars_xy(ind_tc3_win,j,i) /= building_pars_f%fill )    & 
-                    surf_usm_v(l)%lambda_h_window(nzb_wall+3,m) = building_pars_f%pars_xy(ind_tc3_win,j,i)
-
-                 IF ( building_pars_f%pars_xy(117,j,i) /= building_pars_f%fill )         & 
-                    surf_usm_v(l)%target_temp_summer(m)  = building_pars_f%pars_xy(117,j,i)   
-                 IF ( building_pars_f%pars_xy(118,j,i) /= building_pars_f%fill )         & 
-                    surf_usm_v(l)%target_temp_winter(m)  = building_pars_f%pars_xy(118,j,i)   
-
-                 IF ( building_pars_f%pars_xy(ind_emis_wall,j,i) /= building_pars_f%fill ) & 
-                    surf_usm_v(l)%emissivity(ind_veg_wall,m)  = building_pars_f%pars_xy(ind_emis_wall,j,i)
-                 IF ( building_pars_f%pars_xy(ind_emis_green,j,i) /= building_pars_f%fill )& 
-                    surf_usm_v(l)%emissivity(ind_pav_green,m) = building_pars_f%pars_xy(ind_emis_green,j,i)
-                 IF ( building_pars_f%pars_xy(ind_emis_win,j,i) /= building_pars_f%fill )  & 
-                    surf_usm_v(l)%emissivity(ind_wat_win,m)   = building_pars_f%pars_xy(ind_emis_win,j,i)
-
-                 IF ( building_pars_f%pars_xy(ind_trans,j,i) /= building_pars_f%fill )    & 
-                    surf_usm_v(l)%transmissivity(m)      = building_pars_f%pars_xy(ind_trans,j,i)
-
-                 IF ( building_pars_f%pars_xy(ind_z0,j,i) /= building_pars_f%fill )    & 
-                    surf_usm_v(l)%z0(m)                  = building_pars_f%pars_xy(ind_z0,j,i)
-                 IF ( building_pars_f%pars_xy(ind_z0qh,j,i) /= building_pars_f%fill )  & 
-                    surf_usm_v(l)%z0h(m)                 = building_pars_f%pars_xy(ind_z0qh,j,i)
-                 IF ( building_pars_f%pars_xy(ind_z0qh,j,i) /= building_pars_f%fill )  & 
-                    surf_usm_v(l)%z0q(m)                 = building_pars_f%pars_xy(ind_z0qh,j,i)
-
-                 IF ( building_pars_f%pars_xy(ind_alb_wall,j,i) /= building_pars_f%fill )    & 
-                    surf_usm_v(l)%albedo_type(ind_veg_wall,m)  = building_pars_f%pars_xy(ind_alb_wall,j,i)
-                 IF ( building_pars_f%pars_xy(ind_alb_green,j,i) /= building_pars_f%fill )   & 
-                    surf_usm_v(l)%albedo_type(ind_pav_green,m) = building_pars_f%pars_xy(ind_alb_green,j,i)
-                 IF ( building_pars_f%pars_xy(ind_alb_win,j,i) /= building_pars_f%fill )     & 
-                    surf_usm_v(l)%albedo_type(ind_wat_win,m)   = building_pars_f%pars_xy(ind_alb_win,j,i)
-
-                 IF ( building_pars_f%pars_xy(ind_thick_1,j,i) /= building_pars_f%fill )    & 
-                    surf_usm_v(l)%zw(nzb_wall,m) = building_pars_f%pars_xy(ind_thick_1,j,i)
-                 IF ( building_pars_f%pars_xy(ind_thick_2,j,i) /= building_pars_f%fill )    & 
-                    surf_usm_v(l)%zw(nzb_wall+1,m) = building_pars_f%pars_xy(ind_thick_2,j,i)
-                 IF ( building_pars_f%pars_xy(ind_thick_3,j,i) /= building_pars_f%fill )    & 
-                    surf_usm_v(l)%zw(nzb_wall+2,m) = building_pars_f%pars_xy(ind_thick_3,j,i)
-                 IF ( building_pars_f%pars_xy(ind_thick_4,j,i) /= building_pars_f%fill )    & 
-                    surf_usm_v(l)%zw(nzb_wall+3,m) = building_pars_f%pars_xy(ind_thick_4,j,i)
-                 IF ( building_pars_f%pars_xy(ind_thick_1,j,i) /= building_pars_f%fill )    & 
-                    surf_usm_v(l)%zw_green(nzb_wall,m) = building_pars_f%pars_xy(ind_thick_1,j,i)
-                 IF ( building_pars_f%pars_xy(ind_thick_2,j,i) /= building_pars_f%fill )    & 
-                    surf_usm_v(l)%zw_green(nzb_wall+1,m) = building_pars_f%pars_xy(ind_thick_2,j,i)
-                 IF ( building_pars_f%pars_xy(ind_thick_3,j,i) /= building_pars_f%fill )    & 
-                    surf_usm_v(l)%zw_green(nzb_wall+2,m) = building_pars_f%pars_xy(ind_thick_3,j,i)
-                 IF ( building_pars_f%pars_xy(ind_thick_4,j,i) /= building_pars_f%fill )    & 
-                    surf_usm_v(l)%zw_green(nzb_wall+3,m) = building_pars_f%pars_xy(ind_thick_4,j,i)
-                 IF ( building_pars_f%pars_xy(ind_thick_1_win,j,i) /= building_pars_f%fill )    & 
-                    surf_usm_v(l)%zw_window(nzb_wall,m) = building_pars_f%pars_xy(ind_thick_1_win,j,i)
-                 IF ( building_pars_f%pars_xy(ind_thick_2_win,j,i) /= building_pars_f%fill )    & 
-                    surf_usm_v(l)%zw_window(nzb_wall+1,m) = building_pars_f%pars_xy(ind_thick_2_win,j,i)
-                 IF ( building_pars_f%pars_xy(ind_thick_3_win,j,i) /= building_pars_f%fill )    & 
-                    surf_usm_v(l)%zw_window(nzb_wall+2,m) = building_pars_f%pars_xy(ind_thick_3_win,j,i)
-                 IF ( building_pars_f%pars_xy(ind_thick_4_win,j,i) /= building_pars_f%fill )    & 
-                    surf_usm_v(l)%zw_window(nzb_wall+3,m) = building_pars_f%pars_xy(ind_thick_4_win,j,i)
-
-                 IF ( building_pars_f%pars_xy(0,j,i) /= building_pars_f%fill )    & 
-                    surf_usm_v(l)%c_surface(m)           = building_pars_f%pars_xy(0,j,i)
-                 IF ( building_pars_f%pars_xy(3,j,i) /= building_pars_f%fill )    & 
-                    surf_usm_v(l)%lambda_surf(m)         = building_pars_f%pars_xy(3,j,i)
-                 IF ( building_pars_f%pars_xy(2,j,i) /= building_pars_f%fill )    & 
-                    surf_usm_v(l)%c_surface_green(m)     = building_pars_f%pars_xy(2,j,i)
-                 IF ( building_pars_f%pars_xy(5,j,i) /= building_pars_f%fill )    & 
-                    surf_usm_v(l)%lambda_surf_green(m)   = building_pars_f%pars_xy(5,j,i)
-                 IF ( building_pars_f%pars_xy(1,j,i) /= building_pars_f%fill )    & 
-                    surf_usm_v(l)%c_surface_window(m)    = building_pars_f%pars_xy(1,j,i)
-                 IF ( building_pars_f%pars_xy(4,j,i) /= building_pars_f%fill )    & 
-                    surf_usm_v(l)%lambda_surf_window(m)  = building_pars_f%pars_xy(4,j,i)
-
+                    ind_green_frac_w = MERGE( ind_green_frac_w_gfl,            &
+                                              ind_green_frac_w_agfl,           &
+                                              surf_usm_v(l)%ground_level(m) )
+                    ind_win_frac     = MERGE( ind_win_frac_gfl,                &
+                                              ind_win_frac_agfl,               &
+                                              surf_usm_v(l)%ground_level(m) )
+                    ind_lai_w        = MERGE( ind_lai_w_gfl,                   &
+                                              ind_lai_w_agfl,                  &
+                                              surf_usm_v(l)%ground_level(m) )
+                    ind_z0           = MERGE( ind_z0_gfl,                      &
+                                              ind_z0_agfl,                     &
+                                              surf_usm_v(l)%ground_level(m) )
+                    ind_z0qh         = MERGE( ind_z0qh_gfl,                    &
+                                              ind_z0qh_agfl,                   &
+                                              surf_usm_v(l)%ground_level(m) )
+                    ind_hc1          = MERGE( ind_hc1_gfl,                     &
+                                              ind_hc1_agfl,                    &
+                                              surf_usm_v(l)%ground_level(m) )
+                    ind_hc2          = MERGE( ind_hc2_gfl,                     &
+                                              ind_hc2_agfl,                    &
+                                              surf_usm_v(l)%ground_level(m) )
+                    ind_hc3          = MERGE( ind_hc3_gfl,                     &
+                                              ind_hc3_agfl,                    &
+                                              surf_usm_v(l)%ground_level(m) )
+                    ind_tc1          = MERGE( ind_tc1_gfl,                     &
+                                              ind_tc1_agfl,                    &
+                                              surf_usm_v(l)%ground_level(m) )
+                    ind_tc2          = MERGE( ind_tc2_gfl,                     &
+                                              ind_tc2_agfl,                    &
+                                              surf_usm_v(l)%ground_level(m) )
+                    ind_tc3          = MERGE( ind_tc3_gfl,                     &
+                                              ind_tc3_agfl,                    &
+                                              surf_usm_v(l)%ground_level(m) )
+                    ind_emis_wall    = MERGE( ind_emis_wall_gfl,               &
+                                              ind_emis_wall_agfl,              &
+                                              surf_usm_v(l)%ground_level(m) )
+                    ind_emis_green   = MERGE( ind_emis_green_gfl,              &
+                                              ind_emis_green_agfl,             &
+                                              surf_usm_v(l)%ground_level(m) )
+                    ind_emis_win     = MERGE( ind_emis_win_gfl,                &
+                                              ind_emis_win_agfl,               &
+                                              surf_usm_v(l)%ground_level(m) )
+                    ind_trans        = MERGE( ind_trans_gfl,                   &
+                                              ind_trans_agfl,                  &
+                                              surf_usm_v(l)%ground_level(m) )
+                    
+!                   
+!--                 Initialize relatvie wall- (0), green- (1) and window (2) fractions
+                    IF ( building_pars_f%pars_xy(ind_wall_frac,j,i) /=         &
+                         building_pars_f%fill )                                &
+                       surf_usm_v(l)%frac(ind_veg_wall,m)  =                   &
+                                          building_pars_f%pars_xy(ind_wall_frac,j,i)   
+                       
+                    IF ( building_pars_f%pars_xy(ind_green_frac_w,j,i) /=      &          
+                         building_pars_f%fill )                                & 
+                       surf_usm_v(l)%frac(ind_pav_green,m) =                   &
+                                  building_pars_f%pars_xy(ind_green_frac_w,j,i) 
+                       
+                    IF ( building_pars_f%pars_xy(ind_win_frac,j,i) /=          &
+                         building_pars_f%fill )                                & 
+                       surf_usm_v(l)%frac(ind_wat_win,m)   =                   &
+                                       building_pars_f%pars_xy(ind_win_frac,j,i)
+                    
+                    IF ( building_pars_f%pars_xy(ind_lai_w,j,i) /=             &
+                         building_pars_f%fill )                                &
+                       surf_usm_v(l)%lai(m)  =                                 &
+                                       building_pars_f%pars_xy(ind_lai_w,j,i)
+                    
+                    IF ( building_pars_f%pars_xy(ind_hc1,j,i) /=               &
+                         building_pars_f%fill )  THEN
+                       surf_usm_v(l)%rho_c_wall(nzb_wall,m)   =                &
+                                          building_pars_f%pars_xy(ind_hc1,j,i) 
+                       surf_usm_v(l)%rho_c_wall(nzb_wall+1,m) =                &
+                                          building_pars_f%pars_xy(ind_hc1,j,i)
+                    ENDIF
+                    
+                    
+                    IF ( building_pars_f%pars_xy(ind_hc2,j,i) /=               &
+                         building_pars_f%fill )                                &
+                       surf_usm_v(l)%rho_c_wall(nzb_wall+2,m) =                &
+                                          building_pars_f%pars_xy(ind_hc2,j,i)
+                       
+                    IF ( building_pars_f%pars_xy(ind_hc3,j,i) /=               &         
+                         building_pars_f%fill )                                &
+                       surf_usm_v(l)%rho_c_wall(nzb_wall+3,m) =                &
+                                          building_pars_f%pars_xy(ind_hc3,j,i)
+                       
+                    IF ( building_pars_f%pars_xy(ind_hc1,j,i) /=               &
+                         building_pars_f%fill )  THEN 
+                       surf_usm_v(l)%rho_c_green(nzb_wall,m)   =               &
+                                          building_pars_f%pars_xy(ind_hc1,j,i) 
+                       surf_usm_v(l)%rho_c_green(nzb_wall+1,m) =               &
+                                          building_pars_f%pars_xy(ind_hc1,j,i)
+                    ENDIF
+                    IF ( building_pars_f%pars_xy(ind_hc2,j,i) /=               &
+                         building_pars_f%fill )                                &
+                       surf_usm_v(l)%rho_c_green(nzb_wall+2,m) =               &
+                                          building_pars_f%pars_xy(ind_hc2,j,i)
+                       
+                    IF ( building_pars_f%pars_xy(ind_hc3,j,i) /=               &
+                         building_pars_f%fill )                                &
+                       surf_usm_v(l)%rho_c_green(nzb_wall+3,m) =               &
+                                          building_pars_f%pars_xy(ind_hc3,j,i)
+                       
+                    IF ( building_pars_f%pars_xy(ind_hc1,j,i) /=               &
+                         building_pars_f%fill )  THEN 
+                       surf_usm_v(l)%rho_c_window(nzb_wall,m)   =              &
+                                          building_pars_f%pars_xy(ind_hc1,j,i) 
+                       surf_usm_v(l)%rho_c_window(nzb_wall+1,m) =              &
+                                          building_pars_f%pars_xy(ind_hc1,j,i)
+                    ENDIF
+                    IF ( building_pars_f%pars_xy(ind_hc2,j,i) /=               &
+                         building_pars_f%fill )                                &
+                       surf_usm_v(l)%rho_c_window(nzb_wall+2,m) =              &
+                                          building_pars_f%pars_xy(ind_hc2,j,i)
+                       
+                    IF ( building_pars_f%pars_xy(ind_hc3,j,i) /=               &
+                         building_pars_f%fill )                                &
+                       surf_usm_v(l)%rho_c_window(nzb_wall+3,m) =              &
+                                          building_pars_f%pars_xy(ind_hc3,j,i)
+                    
+                    IF ( building_pars_f%pars_xy(ind_tc1,j,i) /=               &
+                         building_pars_f%fill )  THEN
+                       surf_usm_v(l)%lambda_h(nzb_wall,m)   =                  &
+                                          building_pars_f%pars_xy(ind_tc1,j,i)   
+                       surf_usm_v(l)%lambda_h(nzb_wall+1,m) =                  &
+                                          building_pars_f%pars_xy(ind_tc1,j,i)  
+                    ENDIF
+                    IF ( building_pars_f%pars_xy(ind_tc2,j,i) /=               &  
+                         building_pars_f%fill )                                &
+                       surf_usm_v(l)%lambda_h(nzb_wall+2,m) =                  &
+                                          building_pars_f%pars_xy(ind_tc2,j,i)
+                       
+                    IF ( building_pars_f%pars_xy(ind_tc3,j,i) /=               &
+                         building_pars_f%fill )                                & 
+                       surf_usm_v(l)%lambda_h(nzb_wall+3,m) =                  &
+                                          building_pars_f%pars_xy(ind_tc3,j,i) 
+                       
+                    IF ( building_pars_f%pars_xy(ind_tc1,j,i) /=               &
+                         building_pars_f%fill )  THEN
+                       surf_usm_v(l)%lambda_h_green(nzb_wall,m)   =            &
+                                           building_pars_f%pars_xy(ind_tc1,j,i)   
+                       surf_usm_v(l)%lambda_h_green(nzb_wall+1,m) =            &
+                                           building_pars_f%pars_xy(ind_tc1,j,i)  
+                    ENDIF
+                    IF ( building_pars_f%pars_xy(ind_tc2,j,i) /=               & 
+                         building_pars_f%fill )                                &
+                       surf_usm_v(l)%lambda_h_green(nzb_wall+2,m) =            &
+                                          building_pars_f%pars_xy(ind_tc2,j,i)
+                       
+                    IF ( building_pars_f%pars_xy(ind_tc3,j,i) /=               &       
+                         building_pars_f%fill )                                &
+                       surf_usm_v(l)%lambda_h_green(nzb_wall+3,m) =            &
+                                          building_pars_f%pars_xy(ind_tc3,j,i) 
+                       
+                    IF ( building_pars_f%pars_xy(ind_tc1_win_r,j,i) /=         &
+                         building_pars_f%fill )  THEN
+                       surf_usm_v(l)%lambda_h_window(nzb_wall,m)   =           &
+                                     building_pars_f%pars_xy(ind_tc1_win_r,j,i)         
+                       surf_usm_v(l)%lambda_h_window(nzb_wall+1,m) =           &
+                                     building_pars_f%pars_xy(ind_tc1_win_r,j,i)        
+                    ENDIF
+                    IF ( building_pars_f%pars_xy(ind_tc2,j,i) /=               &     
+                         building_pars_f%fill )                                &
+                       surf_usm_v(l)%lambda_h_window(nzb_wall+2,m) =           &
+                                           building_pars_f%pars_xy(ind_tc2,j,i)
+                       
+                    IF ( building_pars_f%pars_xy(ind_tc3,j,i) /=               &   
+                         building_pars_f%fill )                                &
+                       surf_usm_v(l)%lambda_h_window(nzb_wall+3,m) =           &
+                                          building_pars_f%pars_xy(ind_tc3,j,i)    
+                    
+                    IF ( building_pars_f%pars_xy(ind_indoor_target_temp_summer,j,i) /=&            
+                         building_pars_f%fill )                                & 
+                       surf_usm_v(l)%target_temp_summer(m)  =                  &
+                            building_pars_f%pars_xy(ind_indoor_target_temp_summer,j,i)   
+                    IF ( building_pars_f%pars_xy(ind_indoor_target_temp_winter,j,i) /=&            
+                         building_pars_f%fill )                                & 
+                       surf_usm_v(l)%target_temp_winter(m)  =                  &
+                            building_pars_f%pars_xy(ind_indoor_target_temp_winter,j,i)   
+                    
+                    IF ( building_pars_f%pars_xy(ind_emis_wall,j,i) /=         &   
+                         building_pars_f%fill )                                &
+                       surf_usm_v(l)%emissivity(ind_veg_wall,m)  =             &
+                                      building_pars_f%pars_xy(ind_emis_wall,j,i)
+                       
+                    IF ( building_pars_f%pars_xy(ind_emis_green,j,i) /=        &            
+                         building_pars_f%fill )                                &
+                       surf_usm_v(l)%emissivity(ind_pav_green,m) =             &
+                                      building_pars_f%pars_xy(ind_emis_green,j,i)
+                       
+                    IF ( building_pars_f%pars_xy(ind_emis_win,j,i) /=          & 
+                         building_pars_f%fill )                                &
+                       surf_usm_v(l)%emissivity(ind_wat_win,m)   =             &
+                                      building_pars_f%pars_xy(ind_emis_win,j,i)
+                       
+                    IF ( building_pars_f%pars_xy(ind_trans,j,i) /=             &    
+                         building_pars_f%fill )                                &
+                       surf_usm_v(l)%transmissivity(m) =                       &
+                                          building_pars_f%pars_xy(ind_trans,j,i)
+                    
+                    IF ( building_pars_f%pars_xy(ind_z0,j,i) /=                &          
+                         building_pars_f%fill )                                &
+                       surf_usm_v(l)%z0(m) = building_pars_f%pars_xy(ind_z0,j,i)
+                       
+                    IF ( building_pars_f%pars_xy(ind_z0qh,j,i) /=              &            
+                         building_pars_f%fill )                                &
+                       surf_usm_v(l)%z0h(m) =                                  &
+                                       building_pars_f%pars_xy(ind_z0qh,j,i)
+                    IF ( building_pars_f%pars_xy(ind_z0qh,j,i) /=              &           
+                         building_pars_f%fill )                                &
+                       surf_usm_v(l)%z0q(m) =                                  &
+                                       building_pars_f%pars_xy(ind_z0qh,j,i)
+                    
+                    IF ( building_pars_f%pars_xy(ind_alb_wall_agfl,j,i) /=     &          
+                         building_pars_f%fill )                                & 
+                       surf_usm_v(l)%albedo_type(ind_veg_wall,m)  =            &
+                                 building_pars_f%pars_xy(ind_alb_wall_agfl,j,i)
+                       
+                    IF ( building_pars_f%pars_xy(ind_alb_green_agfl,j,i) /=    &           
+                         building_pars_f%fill )                                &
+                       surf_usm_v(l)%albedo_type(ind_pav_green,m) =            &
+                                 building_pars_f%pars_xy(ind_alb_green_agfl,j,i)
+                    IF ( building_pars_f%pars_xy(ind_alb_win_agfl,j,i) /=      &          
+                         building_pars_f%fill )                                &
+                       surf_usm_v(l)%albedo_type(ind_wat_win,m)   =            &
+                                   building_pars_f%pars_xy(ind_alb_win_agfl,j,i)
+                    
+                    IF ( building_pars_f%pars_xy(ind_thick_1_agfl,j,i) /=      &          
+                         building_pars_f%fill )                                & 
+                       surf_usm_v(l)%zw(nzb_wall,m) =                          &
+                                   building_pars_f%pars_xy(ind_thick_1_agfl,j,i)
+                       
+                    IF ( building_pars_f%pars_xy(ind_thick_2_agfl,j,i) /=      &          
+                         building_pars_f%fill )                                &
+                       surf_usm_v(l)%zw(nzb_wall+1,m) =                        &
+                                   building_pars_f%pars_xy(ind_thick_2_agfl,j,i)
+                       
+                    IF ( building_pars_f%pars_xy(ind_thick_3_agfl,j,i) /=      &         
+                         building_pars_f%fill )                                &
+                       surf_usm_v(l)%zw(nzb_wall+2,m) =                        &
+                                   building_pars_f%pars_xy(ind_thick_3_agfl,j,i)
+                       
+                       
+                    IF ( building_pars_f%pars_xy(ind_thick_4_agfl,j,i) /=      &         
+                         building_pars_f%fill )                                & 
+                       surf_usm_v(l)%zw(nzb_wall+3,m) =                        &
+                                   building_pars_f%pars_xy(ind_thick_4_agfl,j,i)
+                       
+                    IF ( building_pars_f%pars_xy(ind_thick_1_agfl,j,i) /=      &           
+                         building_pars_f%fill )                                &
+                       surf_usm_v(l)%zw_green(nzb_wall,m) =                    &
+                                   building_pars_f%pars_xy(ind_thick_1_agfl,j,i)
+                       
+                    IF ( building_pars_f%pars_xy(ind_thick_2_agfl,j,i) /=      &          
+                         building_pars_f%fill )                                &
+                       surf_usm_v(l)%zw_green(nzb_wall+1,m) =                  &
+                                   building_pars_f%pars_xy(ind_thick_2_agfl,j,i)
+                       
+                    IF ( building_pars_f%pars_xy(ind_thick_3_agfl,j,i) /=      &          
+                         building_pars_f%fill )                                & 
+                       surf_usm_v(l)%zw_green(nzb_wall+2,m) =                  &
+                                   building_pars_f%pars_xy(ind_thick_3_agfl,j,i)
+                       
+                    IF ( building_pars_f%pars_xy(ind_thick_4_agfl,j,i) /=      &          
+                         building_pars_f%fill )                                &
+                       surf_usm_v(l)%zw_green(nzb_wall+3,m) =                  &
+                                   building_pars_f%pars_xy(ind_thick_4_agfl,j,i)
+                    
+                    IF ( building_pars_f%pars_xy(ind_c_surface,j,i) /=         &       
+                         building_pars_f%fill )                                & 
+                       surf_usm_v(l)%c_surface(m) =                            &
+                                     building_pars_f%pars_xy(ind_c_surface,j,i)
+                       
+                    IF ( building_pars_f%pars_xy(ind_lambda_surf,j,i) /=       &        
+                         building_pars_f%fill )                                &
+                       surf_usm_v(l)%lambda_surf(m) =                          &
+                                    building_pars_f%pars_xy(ind_lambda_surf,j,i)
+                    
               ENDDO
            ENDDO
@@ -5773,4 +6585,5 @@
    
               ENDDO
+              lambda_green_temp(nzt_wall+1) = lambda_green_temp(nzt_wall)
    
    
