Changeset 3881


Ignore:
Timestamp:
Apr 10, 2019 9:31:22 AM (2 years ago)
Author:
suehring
Message:

Make level 3 initialization of urban-surfaces consistent to input data standard; revise flagging of ground-floor level surfaces at buidlings; bugfixes in level 3 initialization of albedo; bugfix in OpenMP directive; check for zero output timestep in surface output; assure that Obukhov lenght does not become zero

Location:
palm/trunk/SOURCE
Files:
6 edited

Legend:

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

    r3868 r3881  
    2525! -----------------
    2626! $Id$
     27! Bugfix in level 3 initialization of pavement albedo type and pavement
     28! emissivity
     29!
     30! 3868 2019-04-08 11:52:36Z suehring
    2731! More strict limitation of roughness length when it is in the order of the
    2832! vertical grid spacing
     
    41864190                     pavement_pars_f%fill )                                    &
    41874191                   surf_lsm_h%albedo_type(ind_pav_green,m) =                   &
    4188                                               INT( pavement_pars(ind_p_at,st) )
     4192                                   INT( pavement_pars_f%pars_xy(ind_p_at,j,i) )
    41894193                IF ( pavement_pars_f%pars_xy(ind_p_emis,j,i) /=                &
    41904194                     pavement_pars_f%fill )                                    &
    41914195                   surf_lsm_h%emissivity(ind_pav_green,m)  =                   &
    4192                                               pavement_pars(ind_p_emis,st)
     4196                                   pavement_pars_f%pars_xy(ind_p_emis,j,i)
    41934197             ENDIF
    41944198
     
    42464250                        pavement_pars_f%fill )                                 &
    42474251                      surf_lsm_v(l)%albedo_type(ind_pav_green,m) =             &
    4248                                             INT( pavement_pars(ind_p_at,st) )
     4252                                   INT( pavement_pars_f%pars_xy(ind_p_at,j,i) )
    42494253
    42504254                   IF ( pavement_pars_f%pars_xy(ind_p_emis,j,i) /=             &
    42514255                        pavement_pars_f%fill )                                 &
    42524256                      surf_lsm_v(l)%emissivity(ind_pav_green,m)  =             &
    4253                                             pavement_pars(ind_p_emis,st
     4257                                   pavement_pars_f%pars_xy(ind_p_emis,j,i
    42544258                ENDIF
    42554259             ENDDO
  • palm/trunk/SOURCE/prognostic_equations.f90

    r3880 r3881  
    2525! -----------------
    2626! $Id$
     27! Bugfix in OpenMP directive
     28!
     29! 3880 2019-04-08 21:43:02Z knoop
    2730! Moved wtm_tendencies to module_interface_actions
    2831!
     
    589592!
    590593!-- Loop over all prognostic equations
    591 !-- b, c ang g added for SALSA
    592     !$OMP PARALLEL PRIVATE (i,i_omp_start,j,k,loop_start,tn,b,c,g)
     594    !$OMP PARALLEL PRIVATE (i,i_omp_start,j,k,loop_start,tn)
    593595
    594596    !$  tn = omp_get_thread_num()
  • palm/trunk/SOURCE/radiation_model_mod.f90

    r3861 r3881  
    2828! -----------------
    2929! $Id$
     30! Output of albedo and emissivity moved from USM, bugfixes in initialization
     31! of albedo
     32!
     33! 3861 2019-04-04 06:27:41Z maronga
    3034! Bugfix: factor of 4.0 required instead of 3.0 in calculation of rad_lw_out_change_0
    3135!
     
    13971401          unit = 'W/m2'
    13981402       ELSE IF ( var(1:7) == 'rtm_svf'  .OR.  var(1:7) == 'rtm_dif'  .OR.                &
    1399                  var(1:9) == 'rtm_skyvf' .OR. var(1:9) == 'rtm_skyvft' )  THEN
     1403                 var(1:9) == 'rtm_skyvf' .OR. var(1:9) == 'rtm_skyvft'  .OR.             &
     1404                 var(1:12) == 'rtm_surfalb_'  .OR.  var(1:13) == 'rtm_surfemis_'  ) THEN
    14001405          IF ( .NOT.  radiation ) THEN
    14011406                message_string = 'output of "' // TRIM( var ) // '" require'&
     
    24432448                DO  ind_type = 0, 2
    24442449                   IF ( surf_lsm_h%albedo_type(ind_type,m) == 0 )  THEN
    2445                       IF ( albedo_pars_f%pars_xy(1,j,i) /= albedo_pars_f%fill )&
     2450                      IF ( albedo_pars_f%pars_xy(0,j,i) /= albedo_pars_f%fill )&
    24462451                         surf_lsm_h%albedo(ind_type,m) =                       &
    2447                                                 albedo_pars_f%pars_xy(1,j,i)
     2452                                                albedo_pars_f%pars_xy(0,j,i)
    24482453                      IF ( albedo_pars_f%pars_xy(1,j,i) /= albedo_pars_f%fill )&
    24492454                         surf_lsm_h%aldir(ind_type,m) =                        &
    24502455                                                albedo_pars_f%pars_xy(1,j,i)
     2456                      IF ( albedo_pars_f%pars_xy(1,j,i) /= albedo_pars_f%fill )&
     2457                         surf_lsm_h%aldif(ind_type,m) =                        &
     2458                                                albedo_pars_f%pars_xy(1,j,i)
    24512459                      IF ( albedo_pars_f%pars_xy(2,j,i) /= albedo_pars_f%fill )&
    2452                          surf_lsm_h%aldif(ind_type,m) =                        &
     2460                         surf_lsm_h%asdir(ind_type,m) =                        &
    24532461                                                albedo_pars_f%pars_xy(2,j,i)
    2454                       IF ( albedo_pars_f%pars_xy(3,j,i) /= albedo_pars_f%fill )&
    2455                          surf_lsm_h%asdir(ind_type,m) =                        &
    2456                                                 albedo_pars_f%pars_xy(3,j,i)
    2457                       IF ( albedo_pars_f%pars_xy(4,j,i) /= albedo_pars_f%fill )&
     2462                      IF ( albedo_pars_f%pars_xy(2,j,i) /= albedo_pars_f%fill )&
    24582463                         surf_lsm_h%asdif(ind_type,m) =                        &
    2459                                                 albedo_pars_f%pars_xy(4,j,i)
     2464                                                albedo_pars_f%pars_xy(2,j,i)
    24602465                   ENDIF
    24612466                ENDDO
     
    24692474                   j = surf_usm_h%j(m)
    24702475!
    2471 !--                Spectral albedos for wall/green/window surfaces
     2476!--                Broadband albedos for wall/green/window surfaces
    24722477                   DO  ind_type = 0, 2
    24732478                      IF ( surf_usm_h%albedo_type(ind_type,m) == 0 )  THEN
    2474                          IF ( albedo_pars_f%pars_xy(1,j,i) /= albedo_pars_f%fill )&
     2479                         IF ( albedo_pars_f%pars_xy(0,j,i) /= albedo_pars_f%fill )&
    24752480                            surf_usm_h%albedo(ind_type,m) =                       &
    2476                                                 albedo_pars_f%pars_xy(1,j,i)
    2477                          IF ( albedo_pars_f%pars_xy(1,j,i) /= albedo_pars_f%fill )&
    2478                             surf_usm_h%aldir(ind_type,m) =                        &
    2479                                                 albedo_pars_f%pars_xy(1,j,i)
    2480                          IF ( albedo_pars_f%pars_xy(2,j,i) /= albedo_pars_f%fill )&
    2481                             surf_usm_h%aldif(ind_type,m) =                        &
    2482                                                 albedo_pars_f%pars_xy(2,j,i)
    2483                          IF ( albedo_pars_f%pars_xy(3,j,i) /= albedo_pars_f%fill )&
    2484                             surf_usm_h%asdir(ind_type,m) =                        &
    2485                                                 albedo_pars_f%pars_xy(3,j,i)
    2486                          IF ( albedo_pars_f%pars_xy(4,j,i) /= albedo_pars_f%fill )&
    2487                             surf_usm_h%asdif(ind_type,m) =                        &
    2488                                                 albedo_pars_f%pars_xy(4,j,i)
     2481                                                albedo_pars_f%pars_xy(0,j,i)
    24892482                      ENDIF
    24902483                   ENDDO
     2484!
     2485!--                Spectral albedos especially for building wall surfaces
     2486                   IF ( albedo_pars_f%pars_xy(1,j,i) /= albedo_pars_f%fill )  THEN
     2487                      surf_usm_h%aldir(ind_veg_wall,m) =                       &
     2488                                                albedo_pars_f%pars_xy(1,j,i)
     2489                      surf_usm_h%aldif(ind_veg_wall,m) =                       &
     2490                                                albedo_pars_f%pars_xy(1,j,i)
     2491                   ENDIF
     2492                   IF ( albedo_pars_f%pars_xy(2,j,i) /= albedo_pars_f%fill )  THEN
     2493                      surf_usm_h%asdir(ind_veg_wall,m) =                       &
     2494                                                albedo_pars_f%pars_xy(2,j,i)
     2495                      surf_usm_h%asdif(ind_veg_wall,m) =                       &
     2496                                                albedo_pars_f%pars_xy(2,j,i)
     2497                   ENDIF
     2498!
     2499!--                Spectral albedos especially for building green surfaces
     2500                   IF ( albedo_pars_f%pars_xy(3,j,i) /= albedo_pars_f%fill )  THEN
     2501                      surf_usm_h%aldir(ind_pav_green,m) =                      &
     2502                                                albedo_pars_f%pars_xy(3,j,i)
     2503                      surf_usm_h%aldif(ind_pav_green,m) =                      &
     2504                                                albedo_pars_f%pars_xy(3,j,i)
     2505                   ENDIF
     2506                   IF ( albedo_pars_f%pars_xy(4,j,i) /= albedo_pars_f%fill )  THEN
     2507                      surf_usm_h%asdir(ind_pav_green,m) =                      &
     2508                                                albedo_pars_f%pars_xy(4,j,i)
     2509                      surf_usm_h%asdif(ind_pav_green,m) =                      &
     2510                                                albedo_pars_f%pars_xy(4,j,i)
     2511                   ENDIF
     2512!
     2513!--                Spectral albedos especially for building window surfaces
     2514                   IF ( albedo_pars_f%pars_xy(5,j,i) /= albedo_pars_f%fill )  THEN
     2515                      surf_usm_h%aldir(ind_wat_win,m) =                        &
     2516                                                albedo_pars_f%pars_xy(5,j,i)
     2517                      surf_usm_h%aldif(ind_wat_win,m) =                        &
     2518                                                albedo_pars_f%pars_xy(5,j,i)
     2519                   ENDIF
     2520                   IF ( albedo_pars_f%pars_xy(6,j,i) /= albedo_pars_f%fill )  THEN
     2521                      surf_usm_h%asdir(ind_wat_win,m) =                        &
     2522                                                albedo_pars_f%pars_xy(6,j,i)
     2523                      surf_usm_h%asdif(ind_wat_win,m) =                        &
     2524                                                albedo_pars_f%pars_xy(6,j,i)
     2525                   ENDIF
    24912526
    24922527                ENDDO
     
    25052540                   DO  ind_type = 0, 2
    25062541                      IF ( surf_lsm_v(l)%albedo_type(ind_type,m) == 0 )  THEN
    2507                          IF ( albedo_pars_f%pars_xy(1,j+joff,i+ioff) /=        &
     2542                         IF ( albedo_pars_f%pars_xy(0,j+joff,i+ioff) /=        &
    25082543                              albedo_pars_f%fill )                             &
    25092544                            surf_lsm_v(l)%albedo(ind_type,m) =                 &
    2510                                           albedo_pars_f%pars_xy(1,j+joff,i+ioff)
     2545                                          albedo_pars_f%pars_xy(0,j+joff,i+ioff)
    25112546                         IF ( albedo_pars_f%pars_xy(1,j+joff,i+ioff) /=        &
    25122547                              albedo_pars_f%fill )                             &
    25132548                            surf_lsm_v(l)%aldir(ind_type,m) =                  &
    25142549                                          albedo_pars_f%pars_xy(1,j+joff,i+ioff)
     2550                         IF ( albedo_pars_f%pars_xy(1,j+joff,i+ioff) /=        &
     2551                              albedo_pars_f%fill )                             &
     2552                            surf_lsm_v(l)%aldif(ind_type,m) =                  &
     2553                                          albedo_pars_f%pars_xy(1,j+joff,i+ioff)
    25152554                         IF ( albedo_pars_f%pars_xy(2,j+joff,i+ioff) /=        &
    25162555                              albedo_pars_f%fill )                             &
    2517                             surf_lsm_v(l)%aldif(ind_type,m) =                  &
     2556                            surf_lsm_v(l)%asdir(ind_type,m) =                  &
    25182557                                          albedo_pars_f%pars_xy(2,j+joff,i+ioff)
    2519                          IF ( albedo_pars_f%pars_xy(3,j+joff,i+ioff) /=        &
    2520                               albedo_pars_f%fill )                             &
    2521                             surf_lsm_v(l)%asdir(ind_type,m) =                  &
    2522                                           albedo_pars_f%pars_xy(3,j+joff,i+ioff)
    2523                          IF ( albedo_pars_f%pars_xy(4,j+joff,i+ioff) /=        &
     2558                         IF ( albedo_pars_f%pars_xy(2,j+joff,i+ioff) /=        &
    25242559                              albedo_pars_f%fill )                             &
    25252560                            surf_lsm_v(l)%asdif(ind_type,m) =                  &
    2526                                           albedo_pars_f%pars_xy(4,j+joff,i+ioff)
     2561                                          albedo_pars_f%pars_xy(2,j+joff,i+ioff)
    25272562                      ENDIF
    25282563                   ENDDO
     
    25392574                      j = surf_usm_v(l)%j(m)
    25402575!
    2541 !--                   Spectral albedos for wall/green/window surfaces
     2576!--                   Broadband albedos for wall/green/window surfaces
    25422577                      DO  ind_type = 0, 2
    25432578                         IF ( surf_usm_v(l)%albedo_type(ind_type,m) == 0 )  THEN
    2544                             IF ( albedo_pars_f%pars_xy(1,j+joff,i+ioff) /=        &
    2545                                  albedo_pars_f%fill )                             &
    2546                                surf_usm_v(l)%albedo(ind_type,m) =                 &
    2547                                              albedo_pars_f%pars_xy(1,j+joff,i+ioff)
    2548                             IF ( albedo_pars_f%pars_xy(1,j+joff,i+ioff) /=        &
    2549                                  albedo_pars_f%fill )                             &
    2550                                surf_usm_v(l)%aldir(ind_type,m) =                  &
    2551                                              albedo_pars_f%pars_xy(1,j+joff,i+ioff)
    2552                             IF ( albedo_pars_f%pars_xy(2,j+joff,i+ioff) /=        &
    2553                                  albedo_pars_f%fill )                             &
    2554                                surf_usm_v(l)%aldif(ind_type,m) =                  &
    2555                                              albedo_pars_f%pars_xy(2,j+joff,i+ioff)
    2556                             IF ( albedo_pars_f%pars_xy(3,j+joff,i+ioff) /=        &
    2557                                  albedo_pars_f%fill )                             &
    2558                                surf_usm_v(l)%asdir(ind_type,m) =                  &
    2559                                              albedo_pars_f%pars_xy(3,j+joff,i+ioff)
    2560                             IF ( albedo_pars_f%pars_xy(4,j+joff,i+ioff) /=        &
    2561                                  albedo_pars_f%fill )                             &
    2562                                surf_usm_v(l)%asdif(ind_type,m) =                  &
    2563                                              albedo_pars_f%pars_xy(4,j+joff,i+ioff)
     2579                            IF ( albedo_pars_f%pars_xy(0,j+joff,i+ioff) /=     &
     2580                                 albedo_pars_f%fill )                          &
     2581                               surf_usm_v(l)%albedo(ind_type,m) =              &
     2582                                             albedo_pars_f%pars_xy(0,j+joff,i+ioff)
    25642583                         ENDIF
    25652584                      ENDDO
    2566 
     2585!
     2586!--                   Spectral albedos especially for building wall surfaces
     2587                      IF ( albedo_pars_f%pars_xy(1,j+joff,i+ioff) /=           &
     2588                           albedo_pars_f%fill )  THEN
     2589                         surf_usm_v(l)%aldir(ind_veg_wall,m) =                 &
     2590                                         albedo_pars_f%pars_xy(1,j+joff,i+ioff)
     2591                         surf_usm_v(l)%aldif(ind_veg_wall,m) =                 &
     2592                                         albedo_pars_f%pars_xy(1,j+joff,i+ioff)
     2593                      ENDIF
     2594                      IF ( albedo_pars_f%pars_xy(2,j+joff,i+ioff) /=           &
     2595                           albedo_pars_f%fill )  THEN
     2596                         surf_usm_v(l)%asdir(ind_veg_wall,m) =                 &
     2597                                         albedo_pars_f%pars_xy(2,j+joff,i+ioff)
     2598                         surf_usm_v(l)%asdif(ind_veg_wall,m) =                 &
     2599                                         albedo_pars_f%pars_xy(2,j+joff,i+ioff)
     2600                      ENDIF
     2601!                     
     2602!--                   Spectral albedos especially for building green surfaces
     2603                      IF ( albedo_pars_f%pars_xy(3,j+joff,i+ioff) /=           &
     2604                           albedo_pars_f%fill )  THEN
     2605                         surf_usm_v(l)%aldir(ind_pav_green,m) =                &
     2606                                         albedo_pars_f%pars_xy(3,j+joff,i+ioff)
     2607                         surf_usm_v(l)%aldif(ind_pav_green,m) =                &
     2608                                         albedo_pars_f%pars_xy(3,j+joff,i+ioff)
     2609                      ENDIF
     2610                      IF ( albedo_pars_f%pars_xy(4,j+joff,i+ioff) /=           &
     2611                           albedo_pars_f%fill )  THEN
     2612                         surf_usm_v(l)%asdir(ind_pav_green,m) =                &
     2613                                         albedo_pars_f%pars_xy(4,j+joff,i+ioff)
     2614                         surf_usm_v(l)%asdif(ind_pav_green,m) =                &
     2615                                         albedo_pars_f%pars_xy(4,j+joff,i+ioff)
     2616                      ENDIF
     2617!                     
     2618!--                   Spectral albedos especially for building window surfaces
     2619                      IF ( albedo_pars_f%pars_xy(5,j+joff,i+ioff) /=           &
     2620                           albedo_pars_f%fill )  THEN
     2621                         surf_usm_v(l)%aldir(ind_wat_win,m) =                  &
     2622                                         albedo_pars_f%pars_xy(5,j+joff,i+ioff)
     2623                         surf_usm_v(l)%aldif(ind_wat_win,m) =                  &
     2624                                         albedo_pars_f%pars_xy(5,j+joff,i+ioff)
     2625                      ENDIF
     2626                      IF ( albedo_pars_f%pars_xy(6,j+joff,i+ioff) /=           &
     2627                           albedo_pars_f%fill )  THEN
     2628                         surf_usm_v(l)%asdir(ind_wat_win,m) =                  &
     2629                                         albedo_pars_f%pars_xy(6,j+joff,i+ioff)
     2630                         surf_usm_v(l)%asdif(ind_wat_win,m) =                  &
     2631                                         albedo_pars_f%pars_xy(6,j+joff,i+ioff)
     2632                      ENDIF
    25672633                   ENDDO
    25682634                ENDIF
     
    98099875         var(1:7) == 'rtm_svf'  .OR.  var(1:7) == 'rtm_dif'  .OR.                       &
    98109876         var(1:9) == 'rtm_skyvf' .OR. var(1:10) == 'rtm_skyvft'  .OR.                   &
     9877         var(1:12) == 'rtm_surfalb_'  .OR.  var(1:13) == 'rtm_surfemis_'  .OR.          &
    98119878         var == 'rtm_mrt'  .OR.  var ==  'rtm_mrt_sw'  .OR.  var == 'rtm_mrt_lw' )  THEN
    98129879
     
    1051110578         ENDIF
    1051210579
    10513 !--   block of RTM output variables
    10514 !--   variables are intended mainly for debugging and detailed analyse purposes
    10515       CASE ( 'rtm_skyvf' )
    10516 !--        sky view factor
    10517          DO isurf = dirstart(ids), dirend(ids)
    10518             IF ( surfl(id,isurf) == idsint_u  .OR.  surfl(id,isurf) == idsint_l )  THEN
    10519                local_pf(surfl(ix,isurf),surfl(iy,isurf),surfl(iz,isurf)) = skyvf(isurf)
    10520             ENDIF
    10521          ENDDO
    10522 
    10523       CASE ( 'rtm_skyvft' )
    10524 !--      sky view factor
    10525          DO isurf = dirstart(ids), dirend(ids)
    10526             IF ( surfl(id,isurf) == ids )  THEN
    10527                local_pf(surfl(ix,isurf),surfl(iy,isurf),surfl(iz,isurf)) = skyvft(isurf)
    10528             ENDIF
    10529          ENDDO
    10530 
    10531       CASE ( 'rtm_svf', 'rtm_dif' )
    10532 !--      shape view factors or iradiance factors to selected surface
    10533          IF ( TRIM(var)=='rtm_svf' )  THEN
    10534              k = 1
    10535          ELSE
    10536              k = 2
    10537          ENDIF
    10538          DO isvf = 1, nsvfl
    10539             isurflt = svfsurf(1, isvf)
    10540             isurfs = svfsurf(2, isvf)
    10541 
    10542             IF ( surf(ix,isurfs) == is  .AND.  surf(iy,isurfs) == js  .AND. surf(iz,isurfs) == ks  .AND. &
    10543                  (surf(id,isurfs) == idsint_u .OR. surfl(id,isurfs) == idsint_l ) ) THEN
    10544 !--            correct source surface
    10545                local_pf(surfl(ix,isurflt),surfl(iy,isurflt),surfl(iz,isurflt)) = svf(k,isvf)
    10546             ENDIF
    10547          ENDDO
    10548 
    1054910580      CASE ( 'rtm_rad_net' )
    1055010581!--     array of complete radiation balance
     
    1078510816            ENDIF
    1078610817         ENDIF
    10787 
    10788        CASE DEFAULT
    10789           found = .FALSE.
     10818!         
     10819!--   block of RTM output variables
     10820!--   variables are intended mainly for debugging and detailed analyse purposes
     10821      CASE ( 'rtm_skyvf' )
     10822!     
     10823!--      sky view factor
     10824         DO isurf = dirstart(ids), dirend(ids)
     10825            IF ( surfl(id,isurf) == idsint_u  .OR.  surfl(id,isurf) == idsint_l )  THEN
     10826               local_pf(surfl(ix,isurf),surfl(iy,isurf),surfl(iz,isurf)) = skyvf(isurf)
     10827            ENDIF
     10828         ENDDO
     10829
     10830      CASE ( 'rtm_skyvft' )
     10831!
     10832!--      sky view factor
     10833         DO isurf = dirstart(ids), dirend(ids)
     10834            IF ( surfl(id,isurf) == ids )  THEN
     10835               local_pf(surfl(ix,isurf),surfl(iy,isurf),surfl(iz,isurf)) = skyvft(isurf)
     10836            ENDIF
     10837         ENDDO
     10838
     10839      CASE ( 'rtm_svf', 'rtm_dif' )
     10840!
     10841!--      shape view factors or iradiance factors to selected surface
     10842         IF ( TRIM(var)=='rtm_svf' )  THEN
     10843             k = 1
     10844         ELSE
     10845             k = 2
     10846         ENDIF
     10847         DO isvf = 1, nsvfl
     10848            isurflt = svfsurf(1, isvf)
     10849            isurfs = svfsurf(2, isvf)
     10850
     10851            IF ( surf(ix,isurfs) == is  .AND.  surf(iy,isurfs) == js  .AND. surf(iz,isurfs) == ks  .AND. &
     10852                 (surf(id,isurfs) == idsint_u .OR. surfl(id,isurfs) == idsint_l ) ) THEN
     10853!
     10854!--            correct source surface
     10855               local_pf(surfl(ix,isurflt),surfl(iy,isurflt),surfl(iz,isurflt)) = svf(k,isvf)
     10856            ENDIF
     10857         ENDDO
     10858
     10859      CASE ( 'rtm_surfalb' )
     10860!
     10861!--      surface albedo
     10862         DO isurf = dirstart(ids), dirend(ids)
     10863            IF ( surfl(id,isurf) == idsint_u  .OR.  surfl(id,isurf) == idsint_l )  THEN
     10864               local_pf(surfl(ix,isurf),surfl(iy,isurf),surfl(iz,isurf)) = albedo_surf(isurf)
     10865            ENDIF
     10866         ENDDO
     10867
     10868      CASE ( 'rtm_surfemis' )
     10869!
     10870!--      surface emissivity, weighted average
     10871         DO isurf = dirstart(ids), dirend(ids)
     10872            IF ( surfl(id,isurf) == idsint_u  .OR.  surfl(id,isurf) == idsint_l )  THEN
     10873               local_pf(surfl(ix,isurf),surfl(iy,isurf),surfl(iz,isurf)) = emiss_surf(isurf)
     10874            ENDIF
     10875         ENDDO
     10876
     10877      CASE DEFAULT
     10878         found = .FALSE.
    1079010879
    1079110880    END SELECT
  • palm/trunk/SOURCE/surface_data_output_mod.f90

    r3817 r3881  
    2525! -----------------
    2626! $Id$
     27! Check for zero output timestep (not allowed in parallel NetCDF output mode)
     28!
     29! 3817 2019-03-26 13:53:57Z suehring
    2730! Correct output coordinates of vertical surface elements
    2831!
     
    125128   USE netcdf_interface,                                                       &
    126129       ONLY:  netcdf_create_att, netcdf_create_dim, netcdf_create_file,        &
    127               netcdf_create_global_atts, netcdf_create_var, netcdf_handle_error
     130              netcdf_create_global_atts, netcdf_create_var, netcdf_data_format,&
     131              netcdf_handle_error
    128132
    129133   USE pegrid
     
    43164320#endif
    43174321!
     4322!--    In case of parallel NetCDF output the output timestep must not be zero.
     4323!--    This is because the number of requiered output timesteps is
     4324!--    pre-calculated, which is not possible with zero output timestep.
     4325       IF ( netcdf_data_format > 4 )  THEN
     4326          IF ( dt_dosurf == 0.0_wp )  THEN
     4327             message_string = 'dt_dosurf = 0.0 while using a ' //             &
     4328                              'variable timestep and parallel netCDF4 ' //    &
     4329                              'is not allowed.'
     4330             CALL message( 'surface_data_output_check_parameters', 'PA0081',  &
     4331                           1, 2, 0, 6, 0 )
     4332          ENDIF
     4333         
     4334          IF ( dt_dosurf_av == 0.0_wp )  THEN
     4335             message_string = 'dt_dosurf_av = 0.0 while using a ' //          &
     4336                              'variable timestep and parallel netCDF4 ' //    &
     4337                              'is not allowed.'
     4338             CALL message( 'surface_data_output_check_parameters', 'PA0081',  &
     4339                           1, 2, 0, 6, 0 )
     4340          ENDIF
     4341       ENDIF
     4342
     4343!
    43184344!--   In case of restart runs, check it the number of cores has been changed.
    43194345!--   With surface output this is not allowed.
  • palm/trunk/SOURCE/surface_layer_fluxes_mod.f90

    r3834 r3881  
    2626! -----------------
    2727! $Id$
     28! Assure that Obukhov length does not become zero
     29!
     30! 3834 2019-03-28 15:40:15Z forkel
    2831! added USE chem_gasphase_mod
    2932!
     
    11881191                CYCLE
    11891192             ENDIF
     1193!
     1194!--          Assure that Obukhov length does not become zero
     1195             IF ( ABS( surf%ol(m) ) < 1E-6_wp )                                &
     1196                surf%ol(m) = SIGN( 10E-6_wp, surf%ol(m) )             
    11901197
    11911198          ENDDO
  • palm/trunk/SOURCE/urban_surface_mod.f90

    r3832 r3881  
    2828! -----------------
    2929! $Id$
     30! Revise determination of local ground-floor level height.
     31! Make level 3 initalization conform with Palm-input-data standard
     32! Move output of albedo and emissivity to radiation module
     33!
     34! 3832 2019-03-28 13:16:58Z raasch
    3035! instrumented with openmp directives
    3136!
     
    564569    REAL(wp)     ::  roughness_concrete = 0.001_wp   !< roughness length of average concrete surface
    565570!
    566 !-- Indices of input attributes for (above) ground floor level
    567     INTEGER(iwp) ::  ind_alb_wall_agfl     = 65   !< index in input list for albedo_type of wall above ground floor level
    568     INTEGER(iwp) ::  ind_alb_wall_gfl      = 32   !< index in input list for albedo_type of wall ground floor level
    569     INTEGER(iwp) ::  ind_alb_wall_r        = 96   !< index in input list for albedo_type of wall roof
    570     INTEGER(iwp) ::  ind_alb_green_agfl    = 83   !< index in input list for albedo_type of green above ground floor level
    571     INTEGER(iwp) ::  ind_alb_green_gfl     = 50   !< index in input list for albedo_type of green ground floor level
    572     INTEGER(iwp) ::  ind_alb_green_r       = 115  !< index in input list for albedo_type of green roof
    573     INTEGER(iwp) ::  ind_alb_win_agfl      = 79   !< index in input list for albedo_type of window fraction
     571!-- Indices of input attributes in building_pars for (above) ground floor level
     572    INTEGER(iwp) ::  ind_alb_wall_agfl     = 38   !< index in input list for albedo_type of wall above ground floor level
     573    INTEGER(iwp) ::  ind_alb_wall_gfl      = 66   !< index in input list for albedo_type of wall ground floor level
     574    INTEGER(iwp) ::  ind_alb_wall_r        = 101  !< index in input list for albedo_type of wall roof
     575    INTEGER(iwp) ::  ind_alb_green_agfl    = 39   !< index in input list for albedo_type of green above ground floor level
     576    INTEGER(iwp) ::  ind_alb_green_gfl     = 78   !< index in input list for albedo_type of green ground floor level
     577    INTEGER(iwp) ::  ind_alb_green_r       = 117  !< index in input list for albedo_type of green roof
     578    INTEGER(iwp) ::  ind_alb_win_agfl      = 40   !< index in input list for albedo_type of window fraction above ground floor level
     579    INTEGER(iwp) ::  ind_alb_win_gfl       = 77   !< index in input list for albedo_type of window fraction ground floor level
     580    INTEGER(iwp) ::  ind_alb_win_r         = 115  !< index in input list for albedo_type of window fraction roof
     581    INTEGER(iwp) ::  ind_c_surface         = 45   !< index in input list for heat capacity wall surface
     582    INTEGER(iwp) ::  ind_c_surface_green   = 48   !< index in input list for heat capacity green surface
     583    INTEGER(iwp) ::  ind_c_surface_win     = 47   !< index in input list for heat capacity window surface
     584    INTEGER(iwp) ::  ind_emis_wall_agfl    = 14   !< index in input list for wall emissivity, above ground floor level
     585    INTEGER(iwp) ::  ind_emis_wall_gfl     = 32   !< index in input list for wall emissivity, ground floor level
     586    INTEGER(iwp) ::  ind_emis_wall_r       = 100  !< index in input list for wall emissivity, roof
     587    INTEGER(iwp) ::  ind_emis_green_agfl   = 15   !< index in input list for green emissivity, above ground floor level
     588    INTEGER(iwp) ::  ind_emis_green_gfl    = 34   !< index in input list for green emissivity, ground floor level
     589    INTEGER(iwp) ::  ind_emis_green_r      = 116  !< index in input list for green emissivity, roof
     590    INTEGER(iwp) ::  ind_emis_win_agfl     = 16   !< index in input list for window emissivity, above ground floor level
     591    INTEGER(iwp) ::  ind_emis_win_gfl      = 33   !< index in input list for window emissivity, ground floor level
     592    INTEGER(iwp) ::  ind_emis_win_r        = 113  !< index in input list for window emissivity, roof
     593    INTEGER(iwp) ::  ind_gflh              = 20   !< index in input list for ground floor level height
     594    INTEGER(iwp) ::  ind_green_frac_w_agfl = 2    !< index in input list for green fraction on wall, above ground floor level
     595    INTEGER(iwp) ::  ind_green_frac_w_gfl  = 23   !< index in input list for green fraction on wall, ground floor level
     596    INTEGER(iwp) ::  ind_green_frac_r_agfl = 3    !< index in input list for green fraction on roof, above ground floor level
     597    INTEGER(iwp) ::  ind_green_frac_r_gfl  = 24   !< index in input list for green fraction on roof, ground floor level
     598    INTEGER(iwp) ::  ind_hc1_agfl          = 6    !< index in input list for heat capacity at first wall layer,
    574599                                                  !< above ground floor level
    575     INTEGER(iwp) ::  ind_alb_win_gfl       = 46   !< index in input list for albedo_type of window fraction ground floor level
    576     INTEGER(iwp) ::  ind_alb_win_r         = 110  !< index in input list for albedo_type of window fraction roof
    577     INTEGER(iwp) ::  ind_emis_wall_agfl    = 64   !< index in input list for wall emissivity, above ground floor level
    578     INTEGER(iwp) ::  ind_emis_wall_gfl     = 31   !< index in input list for wall emissivity, ground floor level
    579     INTEGER(iwp) ::  ind_emis_wall_r       = 95   !< index in input list for wall emissivity, roof
    580     INTEGER(iwp) ::  ind_emis_green_agfl   = 82   !< index in input list for green emissivity, above ground floor level
    581     INTEGER(iwp) ::  ind_emis_green_gfl    = 49   !< index in input list for green emissivity, ground floor level
    582     INTEGER(iwp) ::  ind_emis_green_r      = 114  !< index in input list for green emissivity, roof
    583     INTEGER(iwp) ::  ind_emis_win_agfl     = 77   !< index in input list for window emissivity, above ground floor level
    584     INTEGER(iwp) ::  ind_emis_win_gfl      = 44   !< index in input list for window emissivity, ground floor level
    585     INTEGER(iwp) ::  ind_emis_win_r        = 108  !< index in input list for window emissivity, roof
    586     INTEGER(iwp) ::  ind_green_frac_w_agfl = 80   !< index in input list for green fraction on wall, above ground floor level
    587     INTEGER(iwp) ::  ind_green_frac_w_gfl  = 47   !< index in input list for green fraction on wall, ground floor level
    588     INTEGER(iwp) ::  ind_green_frac_r_agfl = 112  !< index in input list for green fraction on roof, above ground floor level
    589     INTEGER(iwp) ::  ind_green_frac_r_gfl  = 111  !< index in input list for green fraction on roof, ground floor level
    590     INTEGER(iwp) ::  ind_hc1_agfl          = 58   !< index in input list for heat capacity at first wall layer,
     600    INTEGER(iwp) ::  ind_hc1_gfl           = 26   !< index in input list for heat capacity at first wall layer, ground floor level
     601    INTEGER(iwp) ::  ind_hc1_wall_r        = 94   !< index in input list for heat capacity at first wall layer, roof
     602    INTEGER(iwp) ::  ind_hc1_win_agfl      = 83   !< index in input list for heat capacity at first window layer,
    591603                                                  !< above ground floor level
    592     INTEGER(iwp) ::  ind_hc1_gfl           = 25   !< index in input list for heat capacity at first wall layer, ground floor level
    593     INTEGER(iwp) ::  ind_hc1_wall_r        = 89   !< index in input list for heat capacity at first wall layer, roof
    594     INTEGER(iwp) ::  ind_hc1_win_agfl      = 71   !< index in input list for heat capacity at first window layer,
     604    INTEGER(iwp) ::  ind_hc1_win_gfl       = 71   !< index in input list for heat capacity at first window layer,
     605                                                  !< ground floor level
     606    INTEGER(iwp) ::  ind_hc1_win_r         = 107  !< index in input list for heat capacity at first window layer, roof
     607    INTEGER(iwp) ::  ind_hc2_agfl          = 7    !< index in input list for heat capacity at second wall layer,
    595608                                                  !< above ground floor level
    596     INTEGER(iwp) ::  ind_hc1_win_gfl       = 38   !< index in input list for heat capacity at first window layer,
     609    INTEGER(iwp) ::  ind_hc2_gfl           = 27   !< index in input list for heat capacity at second wall layer, ground floor level
     610    INTEGER(iwp) ::  ind_hc2_wall_r        = 95   !< index in input list for heat capacity at second wall layer, roof
     611    INTEGER(iwp) ::  ind_hc2_win_agfl      = 84   !< index in input list for heat capacity at second window layer,
     612                                                  !< above ground floor level
     613    INTEGER(iwp) ::  ind_hc2_win_gfl       = 72   !< index in input list for heat capacity at second window layer,
    597614                                                  !< ground floor level
    598     INTEGER(iwp) ::  ind_hc1_win_r         = 102  !< index in input list for heat capacity at first window layer, roof
    599     INTEGER(iwp) ::  ind_hc2_agfl          = 59   !< index in input list for heat capacity at second wall layer,
     615    INTEGER(iwp) ::  ind_hc2_win_r         = 108  !< index in input list for heat capacity at second window layer, roof
     616    INTEGER(iwp) ::  ind_hc3_agfl          = 8    !< index in input list for heat capacity at third wall layer,
    600617                                                  !< above ground floor level
    601     INTEGER(iwp) ::  ind_hc2_gfl           = 26   !< index in input list for heat capacity at second wall layer, ground floor level
    602     INTEGER(iwp) ::  ind_hc2_wall_r        = 90   !< index in input list for heat capacity at second wall layer, roof
    603     INTEGER(iwp) ::  ind_hc2_win_agfl      = 72   !< index in input list for heat capacity at second window layer,
     618    INTEGER(iwp) ::  ind_hc3_gfl           = 28   !< index in input list for heat capacity at third wall layer, ground floor level
     619    INTEGER(iwp) ::  ind_hc3_wall_r        = 96   !< index in input list for heat capacity at third wall layer, roof
     620    INTEGER(iwp) ::  ind_hc3_win_agfl      = 85   !< index in input list for heat capacity at third window layer,
    604621                                                  !< above ground floor level
    605     INTEGER(iwp) ::  ind_hc2_win_gfl       = 39   !< index in input list for heat capacity at second window layer,
     622    INTEGER(iwp) ::  ind_hc3_win_gfl       = 73   !< index in input list for heat capacity at third window layer,
    606623                                                  !< ground floor level
    607     INTEGER(iwp) ::  ind_hc2_win_r         = 103  !< index in input list for heat capacity at second window layer, roof
    608     INTEGER(iwp) ::  ind_hc3_agfl          = 60   !< index in input list for heat capacity at third wall layer,
     624    INTEGER(iwp) ::  ind_hc3_win_r         = 109  !< index in input list for heat capacity at third window layer, roof
     625    INTEGER(iwp) ::  ind_indoor_target_temp_summer = 12
     626    INTEGER(iwp) ::  ind_indoor_target_temp_winter = 13
     627    INTEGER(iwp) ::  ind_lai_r_agfl        = 4    !< index in input list for LAI on roof, above ground floor level
     628    INTEGER(iwp) ::  ind_lai_r_gfl         = 4  !< index in input list for LAI on roof, ground floor level
     629    INTEGER(iwp) ::  ind_lai_w_agfl        = 5    !< index in input list for LAI on wall, above ground floor level
     630    INTEGER(iwp) ::  ind_lai_w_gfl         = 25   !< index in input list for LAI on wall, ground floor level
     631    INTEGER(iwp) ::  ind_lambda_surf       = 46   !< index in input list for thermal conductivity of wall surface
     632    INTEGER(iwp) ::  ind_lambda_surf_green = 50   !< index in input list for thermal conductivity of green surface
     633    INTEGER(iwp) ::  ind_lambda_surf_win   = 49   !< index in input list for thermal conductivity of window surface
     634    INTEGER(iwp) ::  ind_tc1_agfl          = 9    !< index in input list for thermal conductivity at first wall layer,
    609635                                                  !< above ground floor level
    610     INTEGER(iwp) ::  ind_hc3_gfl           = 27   !< index in input list for heat capacity at third wall layer, ground floor level
    611     INTEGER(iwp) ::  ind_hc3_wall_r        = 91   !< index in input list for heat capacity at third wall layer, roof
    612     INTEGER(iwp) ::  ind_hc3_win_agfl      = 73   !< index in input list for heat capacity at third window layer,
     636    INTEGER(iwp) ::  ind_tc1_gfl           = 29   !< index in input list for thermal conductivity at first wall layer,
     637                                                  !< ground floor level
     638    INTEGER(iwp) ::  ind_tc1_wall_r        = 97   !< index in input list for thermal conductivity at first wall layer, roof
     639    INTEGER(iwp) ::  ind_tc1_win_agfl      = 86   !< index in input list for thermal conductivity at first window layer,
    613640                                                  !< above ground floor level
    614     INTEGER(iwp) ::  ind_hc3_win_gfl       = 40   !< index in input list for heat capacity at third window layer,
     641    INTEGER(iwp) ::  ind_tc1_win_gfl       = 74   !< index in input list for thermal conductivity at first window layer,
    615642                                                  !< ground floor level
    616     INTEGER(iwp) ::  ind_hc3_win_r         = 104  !< index in input list for heat capacity at third window layer, roof
    617     INTEGER(iwp) ::  ind_gflh              = 17   !< index in input list for ground floor level height
    618     INTEGER(iwp) ::  ind_lai_r_agfl        = 113  !< index in input list for LAI on roof, above ground floor level
    619     INTEGER(iwp) ::  ind_lai_r_gfl         = 113  !< index in input list for LAI on roof, ground floor level
    620     INTEGER(iwp) ::  ind_lai_w_agfl        = 81   !< index in input list for LAI on wall, above ground floor level
    621     INTEGER(iwp) ::  ind_lai_w_gfl         = 48   !< index in input list for LAI on wall, ground floor level
    622     INTEGER(iwp) ::  ind_tc1_agfl          = 61   !< index in input list for thermal conductivity at first wall layer,
     643    INTEGER(iwp) ::  ind_tc1_win_r         = 110  !< index in input list for thermal conductivity at first window layer, roof
     644    INTEGER(iwp) ::  ind_tc2_agfl          = 10   !< index in input list for thermal conductivity at second wall layer,
    623645                                                  !< above ground floor level
    624     INTEGER(iwp) ::  ind_tc1_gfl           = 28   !< index in input list for thermal conductivity at first wall layer,
     646    INTEGER(iwp) ::  ind_tc2_gfl           = 30   !< index in input list for thermal conductivity at second wall layer,
    625647                                                  !< ground floor level
    626     INTEGER(iwp) ::  ind_tc1_wall_r        = 92   !< index in input list for thermal conductivity at first wall layer, roof
    627     INTEGER(iwp) ::  ind_tc1_win_agfl      = 74   !< index in input list for thermal conductivity at first window layer,
     648    INTEGER(iwp) ::  ind_tc2_wall_r        = 98   !< index in input list for thermal conductivity at second wall layer, roof
     649    INTEGER(iwp) ::  ind_tc2_win_agfl      = 87   !< index in input list for thermal conductivity at second window layer,
    628650                                                  !< above ground floor level
    629     INTEGER(iwp) ::  ind_tc1_win_gfl       = 41   !< index in input list for thermal conductivity at first window layer,
     651    INTEGER(iwp) ::  ind_tc2_win_gfl       = 75   !< index in input list for thermal conductivity at second window layer,
    630652                                                  !< ground floor level
    631     INTEGER(iwp) ::  ind_tc1_win_r         = 105  !< index in input list for thermal conductivity at first window layer, roof
    632     INTEGER(iwp) ::  ind_tc2_agfl          = 62   !< index in input list for thermal conductivity at second wall layer,
     653    INTEGER(iwp) ::  ind_tc2_win_r         = 111  !< index in input list for thermal conductivity at second window layer,
     654                                                  !< ground floor level
     655    INTEGER(iwp) ::  ind_tc3_agfl          = 11   !< index in input list for thermal conductivity at third wall layer,
    633656                                                  !< above ground floor level
    634     INTEGER(iwp) ::  ind_tc2_gfl           = 29   !< index in input list for thermal conductivity at second wall layer,
     657    INTEGER(iwp) ::  ind_tc3_gfl           = 31   !< index in input list for thermal conductivity at third wall layer,
    635658                                                  !< ground floor level
    636     INTEGER(iwp) ::  ind_tc2_wall_r        = 93   !< index in input list for thermal conductivity at second wall layer, roof
    637     INTEGER(iwp) ::  ind_tc2_win_agfl      = 75   !< index in input list for thermal conductivity at second window layer,
     659    INTEGER(iwp) ::  ind_tc3_wall_r        = 99   !< index in input list for thermal conductivity at third wall layer, roof
     660    INTEGER(iwp) ::  ind_tc3_win_agfl      = 88   !< index in input list for thermal conductivity at third window layer,
    638661                                                  !< above ground floor level
    639     INTEGER(iwp) ::  ind_tc2_win_gfl       = 42   !< index in input list for thermal conductivity at second window layer,
     662    INTEGER(iwp) ::  ind_tc3_win_gfl       = 76   !< index in input list for thermal conductivity at third window layer,
    640663                                                  !< ground floor level
    641     INTEGER(iwp) ::  ind_tc2_win_r         = 106  !< index in input list for thermal conductivity at second window layer,
    642                                                   !< ground floor level
    643     INTEGER(iwp) ::  ind_tc3_agfl          = 63   !< index in input list for thermal conductivity at third wall layer,
    644                                                   !< above ground floor level
    645     INTEGER(iwp) ::  ind_tc3_gfl           = 30   !< index in input list for thermal conductivity at third wall layer,
    646                                                   !< ground floor level
    647     INTEGER(iwp) ::  ind_tc3_wall_r        = 94   !< index in input list for thermal conductivity at third wall layer, roof
    648     INTEGER(iwp) ::  ind_tc3_win_agfl      = 76   !< index in input list for thermal conductivity at third window layer,
    649                                                   !< above ground floor level
    650     INTEGER(iwp) ::  ind_tc3_win_gfl       = 43   !< index in input list for thermal conductivity at third window layer,
    651                                                   !< ground floor level
    652     INTEGER(iwp) ::  ind_tc3_win_r         = 107  !< index in input list for thermal conductivity at third window layer, roof
    653     INTEGER(iwp) ::  ind_thick_1_agfl      = 54   !< index for wall layer thickness - 1st layer above ground floor level
    654     INTEGER(iwp) ::  ind_thick_1_gfl       = 21   !< index for wall layer thickness - 1st layer ground floor level
    655     INTEGER(iwp) ::  ind_thick_1_wall_r    = 85   !< index for wall layer thickness - 1st layer roof
    656     INTEGER(iwp) ::  ind_thick_1_win_agfl  = 67   !< index for window layer thickness - 1st layer above ground floor level
    657     INTEGER(iwp) ::  ind_thick_1_win_gfl   = 34   !< index for window layer thickness - 1st layer ground floor level
    658     INTEGER(iwp) ::  ind_thick_1_win_r     = 98   !< index for window layer thickness - 1st layer roof
    659     INTEGER(iwp) ::  ind_thick_2_agfl      = 55   !< index for wall layer thickness - 2nd layer above ground floor level
    660     INTEGER(iwp) ::  ind_thick_2_gfl       = 22   !< index for wall layer thickness - 2nd layer ground floor level
    661     INTEGER(iwp) ::  ind_thick_2_wall_r    = 86   !< index for wall layer thickness - 2nd layer roof
    662     INTEGER(iwp) ::  ind_thick_2_win_agfl  = 68   !< index for window layer thickness - 2nd layer above ground floor level
    663     INTEGER(iwp) ::  ind_thick_2_win_gfl   = 35   !< index for window layer thickness - 2nd layer ground floor level
    664     INTEGER(iwp) ::  ind_thick_2_win_r     = 99   !< index for window layer thickness - 2nd layer roof
    665     INTEGER(iwp) ::  ind_thick_3_agfl      = 56   !< index for wall layer thickness - 3rd layer above ground floor level
    666     INTEGER(iwp) ::  ind_thick_3_gfl       = 23   !< index for wall layer thickness - 3rd layer ground floor level
    667     INTEGER(iwp) ::  ind_thick_3_wall_r    = 87   !< index for wall layer thickness - 3rd layer roof
    668     INTEGER(iwp) ::  ind_thick_3_win_agfl  = 69   !< index for window layer thickness - 3rd layer above ground floor level
    669     INTEGER(iwp) ::  ind_thick_3_win_gfl   = 36   !< index for window layer thickness - 3rd layer ground floor level 
    670     INTEGER(iwp) ::  ind_thick_3_win_r     = 100  !< index for window layer thickness - 3rd layer roof
    671     INTEGER(iwp) ::  ind_thick_4_agfl      = 57   !< index for wall layer thickness - 4th layer above ground floor level
    672     INTEGER(iwp) ::  ind_thick_4_gfl       = 24   !< index for wall layer thickness - 4th layer ground floor level
    673     INTEGER(iwp) ::  ind_thick_4_wall_r    = 88   !< index for wall layer thickness - 4st layer roof
    674     INTEGER(iwp) ::  ind_thick_4_win_agfl  = 70   !< index for window layer thickness - 4th layer above ground floor level
    675     INTEGER(iwp) ::  ind_thick_4_win_gfl   = 37   !< index for window layer thickness - 4th layer ground floor level
    676     INTEGER(iwp) ::  ind_thick_4_win_r     = 101  !< index for window layer thickness - 4th layer roof
    677     INTEGER(iwp) ::  ind_trans_agfl        = 78   !< index in input list for window transmissivity, above ground floor level
    678     INTEGER(iwp) ::  ind_trans_gfl         = 45   !< index in input list for window transmissivity, ground floor level
    679     INTEGER(iwp) ::  ind_trans_r           = 109  !< index in input list for window transmissivity, roof
    680     INTEGER(iwp) ::  ind_wall_frac_agfl    = 53   !< index in input list for wall fraction, above ground floor level
    681     INTEGER(iwp) ::  ind_wall_frac_gfl     = 20   !< index in input list for wall fraction, ground floor level
    682     INTEGER(iwp) ::  ind_wall_frac_r       = 84   !< index in input list for wall fraction, roof
    683     INTEGER(iwp) ::  ind_win_frac_agfl     = 66   !< index in input list for window fraction, above ground floor level
    684     INTEGER(iwp) ::  ind_win_frac_gfl      = 33   !< index in input list for window fraction, ground floor level
    685     INTEGER(iwp) ::  ind_win_frac_r        = 97   !< index in input list for window fraction, roof
    686     INTEGER(iwp) ::  ind_z0_agfl           = 51   !< index in input list for z0, above ground floor level
    687     INTEGER(iwp) ::  ind_z0_gfl            = 18   !< index in input list for z0, ground floor level
    688     INTEGER(iwp) ::  ind_z0qh_agfl         = 52   !< index in input list for z0h / z0q, above ground floor level
    689     INTEGER(iwp) ::  ind_z0qh_gfl          = 19   !< index in input list for z0h / z0q, ground floor level
    690     INTEGER(iwp) ::  ind_green_type_roof   = 116  !< index in input list for type of green roof
     664    INTEGER(iwp) ::  ind_tc3_win_r         = 112  !< index in input list for thermal conductivity at third window layer, roof
     665    INTEGER(iwp) ::  ind_thick_1_agfl      = 41   !< index for wall layer thickness - 1st layer above ground floor level
     666    INTEGER(iwp) ::  ind_thick_1_gfl       = 62   !< index for wall layer thickness - 1st layer ground floor level
     667    INTEGER(iwp) ::  ind_thick_1_wall_r    = 90   !< index for wall layer thickness - 1st layer roof
     668    INTEGER(iwp) ::  ind_thick_1_win_agfl  = 79   !< index for window layer thickness - 1st layer above ground floor level
     669    INTEGER(iwp) ::  ind_thick_1_win_gfl   = 67   !< index for window layer thickness - 1st layer ground floor level
     670    INTEGER(iwp) ::  ind_thick_1_win_r     = 103  !< index for window layer thickness - 1st layer roof
     671    INTEGER(iwp) ::  ind_thick_2_agfl      = 42   !< index for wall layer thickness - 2nd layer above ground floor level
     672    INTEGER(iwp) ::  ind_thick_2_gfl       = 63   !< index for wall layer thickness - 2nd layer ground floor level
     673    INTEGER(iwp) ::  ind_thick_2_wall_r    = 91   !< index for wall layer thickness - 2nd layer roof
     674    INTEGER(iwp) ::  ind_thick_2_win_agfl  = 80   !< index for window layer thickness - 2nd layer above ground floor level
     675    INTEGER(iwp) ::  ind_thick_2_win_gfl   = 68   !< index for window layer thickness - 2nd layer ground floor level
     676    INTEGER(iwp) ::  ind_thick_2_win_r     = 104  !< index for window layer thickness - 2nd layer roof
     677    INTEGER(iwp) ::  ind_thick_3_agfl      = 43   !< index for wall layer thickness - 3rd layer above ground floor level
     678    INTEGER(iwp) ::  ind_thick_3_gfl       = 64   !< index for wall layer thickness - 3rd layer ground floor level
     679    INTEGER(iwp) ::  ind_thick_3_wall_r    = 92   !< index for wall layer thickness - 3rd layer roof
     680    INTEGER(iwp) ::  ind_thick_3_win_agfl  = 81   !< index for window layer thickness - 3rd layer above ground floor level
     681    INTEGER(iwp) ::  ind_thick_3_win_gfl   = 69   !< index for window layer thickness - 3rd layer ground floor level 
     682    INTEGER(iwp) ::  ind_thick_3_win_r     = 105  !< index for window layer thickness - 3rd layer roof
     683    INTEGER(iwp) ::  ind_thick_4_agfl      = 44   !< index for wall layer thickness - 4th layer above ground floor level
     684    INTEGER(iwp) ::  ind_thick_4_gfl       = 65   !< index for wall layer thickness - 4th layer ground floor level
     685    INTEGER(iwp) ::  ind_thick_4_wall_r    = 93   !< index for wall layer thickness - 4st layer roof
     686    INTEGER(iwp) ::  ind_thick_4_win_agfl  = 82   !< index for window layer thickness - 4th layer above ground floor level
     687    INTEGER(iwp) ::  ind_thick_4_win_gfl   = 70   !< index for window layer thickness - 4th layer ground floor level
     688    INTEGER(iwp) ::  ind_thick_4_win_r     = 106  !< index for window layer thickness - 4th layer roof
     689    INTEGER(iwp) ::  ind_trans_agfl        = 17   !< index in input list for window transmissivity, above ground floor level
     690    INTEGER(iwp) ::  ind_trans_gfl         = 35   !< index in input list for window transmissivity, ground floor level
     691    INTEGER(iwp) ::  ind_trans_r           = 114  !< index in input list for window transmissivity, roof
     692    INTEGER(iwp) ::  ind_wall_frac_agfl    = 0    !< index in input list for wall fraction, above ground floor level
     693    INTEGER(iwp) ::  ind_wall_frac_gfl     = 21   !< index in input list for wall fraction, ground floor level
     694    INTEGER(iwp) ::  ind_wall_frac_r       = 89   !< index in input list for wall fraction, roof
     695    INTEGER(iwp) ::  ind_win_frac_agfl     = 1    !< index in input list for window fraction, above ground floor level
     696    INTEGER(iwp) ::  ind_win_frac_gfl      = 22   !< index in input list for window fraction, ground floor level
     697    INTEGER(iwp) ::  ind_win_frac_r        = 102  !< index in input list for window fraction, roof
     698    INTEGER(iwp) ::  ind_z0_agfl           = 18   !< index in input list for z0, above ground floor level
     699    INTEGER(iwp) ::  ind_z0_gfl            = 36   !< index in input list for z0, ground floor level
     700    INTEGER(iwp) ::  ind_z0qh_agfl         = 19   !< index in input list for z0h / z0q, above ground floor level
     701    INTEGER(iwp) ::  ind_z0qh_gfl          = 37   !< index in input list for z0h / z0q, ground floor level
     702    INTEGER(iwp) ::  ind_green_type_roof   = 118  !< index in input list for type of green roof
    691703
    692704
     
    705717                                   'bridges                              '  &  !< type  7
    706718                                                                     /)
    707 !
    708 !-- building parameters, 6 different types
    709 !-- Parameter for urban surface model
    710 !-- 0 - heat capacity wall surface, 1 - heat capacity of window surface, 2 - heat capacity of green surface
    711 !-- 3 - thermal conductivity of wall surface, 4 - thermal conductivity of window surface,
    712 !-- 5 - thermal conductivty of green surface, 6 - wall fraction ground plate,
    713 !-- 7 - 1st wall layer thickness ground plate, 8 - 2nd wall layer thickness ground plate
    714 !-- 9 - 3rd wall layer thickness ground plate, 10 - 4th wall layer thickness ground plate,
    715 !-- 11 - heat capacity 1st/2nd wall layer ground plate, 12 - heat capacity 3rd wall layer ground plate
    716 !-- 13 - heat capacity 4th wall layer ground plate, 14 - thermal conductivity 1st/2nd wall layer ground plate,
    717 !-- 15 - thermal conductivity 3rd wall layer ground plate, 16 - thermal conductivity 4th wall layer ground plate
    718 !-- 17 - ground floor level height, 18 - z0 roughness ground floor level, 19 - z0h/z0g roughness heaat/humidity,
    719 !-- 20 - wall fraction ground floor level, 21 - 1st wall layer thickness ground floor level,
    720 !-- 22 - 2nd wall layer thickness ground floor level, 23 - 3rd wall layer thickness ground floor level,
    721 !-- 24 - 4th wall layer thickness ground floor level, 25 - heat capacity 1st/2nd wall layer ground floor level,
    722 !-- 26 - heat capacity 3rd wall layer ground floor level, 27 - heat capacity 4th wall layer ground floor level,
    723 !-- 28 - thermal conductivity 1st/2nd wall layer ground floor level,
    724 !-- 29 - thermal conductivity 3rd wall layer ground floor level, 30 - thermal conductivity 4th wall layer ground floor level
    725 !-- 31 - wall emissivity ground floor level, 32 - wall albedo ground floor level, 33 - window fraction ground floor level,
    726 !-- 34 - 1st window layer thickness ground floor level, 35 - 2nd window layer thickness ground floor level,
    727 !-- 36 - 3rd window layer thickness ground floor level, 37 - 4th window layer thickness ground floor level,
    728 !-- 38 - heat capacity 1st/2nd window layer ground floor level, 39 - heat capacity 3rd window layer ground floor level,
    729 !-- 40 - heat capacity 4th window layer ground floor level,
    730 !-- 41 - thermal conductivity 1st/2nd window layer ground floor level,
    731 !-- 42 - thermal conductivity 3rd window layer ground floor level,
    732 !-- 43 - thermal conductivity 4th window layer ground floor level, 44 - window emissivity ground floor level,
    733 !-- 45 - window transmissivity ground floor level, 46 - window albedo ground floor level,
    734 !-- 47 - green fraction ground floor level, 48 - LAI on wall ground floor level, 49 - green emissivity ground floor level,
    735 !-- 50 - green albedo ground floor level, 51 - z0 roughness above ground floor level,
    736 !-- 52 - z0h/z0g roughness heat/humidity above ground floor level, 53 - wall fraction above ground floor level
    737 !-- 54 - 1st wall layer thickness above ground floor level, 55 - 2nd wall layer thickness above ground floor level
    738 !-- 56 - 3rd wall layer thickness above ground floor level, 57 - 4th wall layer thickness above ground floor level
    739 !-- 58 - heat capacity 1st/2nd wall layer above ground floor level,
    740 !-- 59 - heat capacity 3rd wall layer above ground floor level,
    741 !-- 60 - heat capacity 4th wall layer above ground floor level,
    742 !-- 61 - thermal conductivity 1st/2nd wall layer above ground floor level,
    743 !-- 62 - thermal conductivity 3rd wall layer above ground floor level,
    744 !-- 63 - thermal conductivity 4th wall layer above ground floor level,
    745 !-- 64 - wall emissivity above ground floor level, 65 - wall albedo above ground floor level,
    746 !-- 66 - window fraction above ground floor level, 67 - 1st window layer thickness above ground floor level,
    747 !-- 68 - 2nd thickness window layer above ground floor level, 69 - 3rd window layer thickness above ground floor level,
    748 !-- 70 - 4th window layer thickness above ground floor level,
    749 !-- 71 - heat capacity 1st/2nd window layer above ground floor level,
    750 !-- 72 - heat capacity 3rd window layer above ground floor level,
    751 !-- 73 - heat capacity 4th window layer above ground floor level,
    752 !-- 74 - conductivity 1st/2nd window layer above ground floor level,
    753 !-- 75 - thermal conductivity 3rd window layer above ground floor level,
    754 !-- 76 - thermal conductivity 4th window layer above ground floor level, 77 - window emissivity above ground floor level,
    755 !-- 78 - window transmissivity above ground floor level, 79 - window albedo above ground floor level,
    756 !-- 80 - green fraction above ground floor level, 81 - LAI on wall above ground floor level,
    757 !-- 82 - green emissivity above ground floor level, 83 - green albedo above ground floor level,
    758 !-- 84 - wall fraction roof, 85 - 1st wall layer thickness roof, 86 - 2nd wall layer thickness roof,
    759 !-- 87 - 3rd wall layer thickness roof, 88 - 4th wall layer thickness roof,
    760 !-- 89 - heat capacity 1st/2nd wall layer roof, 90 - heat capacity 3rd wall layer roof,
    761 !-- 91 - heat capacity 4th wall layer roof, 92 - thermal conductivity 1st/2nd wall layer roof,
    762 !-- 93 - thermal conductivity 3rd wall layer roof, 94 - thermal conductivity 4th wall layer roof,
    763 !-- 95 - wall emissivity roof, 96 - wall albedo roof, 97 - window fraction roof,
    764 !-- 98 - window 1st layer thickness roof, 99 - window 2nd layer thickness roof, 100 - window 3rd layer thickness roof,
    765 !-- 101 - window 4th layer thickness, 102 - heat capacity 1st/2nd window layer roof,
    766 !-- 103 - heat capacity 3rd window layer roof, 104 - heat capacity 4th window layer roof,
    767 !-- 105 - thermal conductivity 1st/2nd window layer roof, 106 - thermal conductivity 3rd window layer roof,
    768 !-- 107 - thermal conductivity 4th window layer roof, 108 - window emissivity roof, 109 - window transmissivity roof,
    769 !-- 110 - window albedo roof, 111 - green fraction roof ground floor level,
    770 !-- 112 - green fraction roof above ground floor level, 113 - LAI roof, 114 - green emissivity roof,
    771 !-- 115 - green albedo roof, 116 - green type roof,
    772 !-- Parameter for indoor model
    773 !-- 117 - indoor target summer temperature, 118 - indoor target winter temperature,
    774 !-- 119 - shading factor, 120 - g-value windows, 121 - u-value windows, 122 - basical airflow without occupancy of the room,
    775 !-- 123 - additional airflow depend of occupancy of the room, 124 - heat recovery efficiency,
    776 !-- 125 - dynamic parameter specific effective surface, 126 - dynamic parameter innner heatstorage,
    777 !-- 127 - ratio internal surface/floor area, 128 - maximal heating capacity, 129 - maximal cooling capacity,
    778 !-- 130 - additional internal heat gains dependent on occupancy of the room,
    779 !-- 131 - basic internal heat gains without occupancy of the room, 132 - storey height, 133 - ceiling construction height
    780 
    781 
     719
     720
     721!
     722!-- Building facade/wall/green/window properties (partly according to PIDS)
    782723    REAL(wp), DIMENSION(0:133,1:7), PARAMETER :: building_pars = RESHAPE( (/   &
    783         20000.0_wp, 20000.0_wp, 20000.0_wp, 23.0_wp, 23.0_wp, 10.0_wp,         & !parameter 0-5
    784         1.0_wp, 0.005_wp, 0.01_wp, 0.39_wp, 0.63_wp, 2200000.0_wp,             & !parameter 6-11
    785         1400000.0_wp, 1300000.0_wp, 0.35_wp, 0.8_wp, 2.1_wp, 4.0_wp,           & !parameter 12-17
    786         0.01_wp, 0.001_wp, 0.75_wp,                                            & !parameter 18-20
    787         0.005_wp, 0.01_wp, 0.39_wp, 0.63_wp, 2200000.0_wp,                     & !parameter 21-25
    788         1400000.0_wp, 1300000.0_wp, 0.35_wp,                                   & !parameter 26-28                     
    789         0.8_wp, 2.1_wp, 0.93_wp,                                               & !parameter 29-31       
    790         27.0_wp, 0.25_wp, 0.003_wp, 0.006_wp, 0.012_wp, 0.018_wp,              & !parameter 32-37
    791         1736000.0_wp, 1736000.0_wp, 1736000.0_wp,                              & !parameter 38-40
    792         0.57_wp, 0.57_wp, 0.57_wp, 0.91_wp,                                    & !parameter 41-44
    793         0.75_wp, 27.0_wp, 0.0_wp, 1.5_wp, 0.86_wp,                             & !parameter 45-49
    794         5.0_wp, 0.001_wp, 0.0001_wp, 0.7_wp, 0.005_wp,                         & !parameter 50-54
    795         0.01_wp, 0.39_wp, 0.63_wp, 2200000.0_wp,                               & !parameter 55-58
    796         1400000.0_wp, 1300000.0_wp, 0.35_wp, 0.8_wp,                           & !parameter 59-62
    797         2.1_wp, 0.93_wp, 27.0_wp, 0.3_wp,                                      & !parameter 63-66
    798         0.003_wp, 0.006_wp, 0.012_wp, 0.018_wp,                                & !parameter 67-70
    799         1736000.0_wp, 1736000.0_wp, 1736000.0_wp,                              & !parameter 71-73
    800         0.57_wp, 0.57_wp, 0.57_wp, 0.91_wp, 0.75_wp,                           & !parameter 74-78
    801         27.0_wp, 0.0_wp, 1.5_wp, 0.86_wp, 5.0_wp, 1.0_wp,                      & !parameter 79-84
    802         0.005_wp, 0.01_wp, 0.31_wp, 0.63_wp, 2200000.0_wp, 1400000.0_wp,       & !parameter 85-90
    803         1300000.0_wp, 0.35_wp, 0.8_wp, 2.1_wp, 0.93_wp, 27.0_wp, 0.0_wp,       & !parameter 91-97
    804         0.003_wp, 0.006_wp, 0.012_wp, 0.018_wp, 1736000.0_wp,                  & !parameter 98-102
    805         1736000.0_wp, 1736000.0_wp, 0.57_wp, 0.57_wp, 0.57_wp,                 & !parameter 103-107
    806         0.91_wp, 0.75_wp, 27.0_wp, 0.0_wp, 0.0_wp, 1.5_wp,                     & !parameter 108-113
    807         0.86_wp, 5.0_wp, 0.0_wp,                                               & !parameter 114-116
    808         299.15_wp, 293.15_wp, 0.8_wp, 0.76_wp, 5.0_wp,                         & !parameter 117-121
    809         0.1_wp, 0.5_wp, 0.0_wp, 3.5_wp, 370000.0_wp, 4.5_wp,                   & !parameter 122-127
    810         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
    811         20000.0_wp, 20000.0_wp, 20000.0_wp, 23.0_wp, 23.0_wp, 10.0_wp,         & !parameter 0-5
    812         1.0_wp, 0.005_wp, 0.01_wp, 0.31_wp, 0.42_wp, 2000000.0_wp,             & !parameter 6-11
    813         103000.0_wp, 900000.0_wp, 0.35_wp, 0.38_wp, 0.04_wp, 4.0_wp,           & !parameter 12-17
    814         0.01_wp, 0.001_wp, 0.78_wp,                                            & !parameter 18-20
    815         0.005_wp, 0.01_wp, 0.31_wp, 0.43_wp, 2000000.0_wp,                     & !parameter 21-25
    816         103000.0_wp, 900000.0_wp, 0.35_wp,                                     & !parameter 26-28                     
    817         0.38_wp, 0.04_wp, 0.92_wp,                                             & !parameter 29-31       
    818         27.0_wp, 0.22_wp, 0.003_wp, 0.006_wp, 0.012_wp, 0.018_wp,              & !parameter 32-37
    819         1736000.0_wp, 1736000.0_wp, 1736000.0_wp,                              & !parameter 38-40
    820         0.11_wp, 0.11_wp, 0.11_wp, 0.11_wp,                                    & !parameter 41-44
    821         0.7_wp, 27.0_wp, 0.0_wp, 1.5_wp, 0.86_wp,                              & !parameter 45-49
    822         5.0_wp, 0.001_wp, 0.0001_wp, 0.73_wp, 0.005_wp,                        & !parameter 50-54
    823         0.01_wp, 0.31_wp, 0.43_wp, 2000000.0_wp,                               & !parameter 55-58
    824         103000.0_wp, 900000.0_wp, 0.35_wp, 0.38_wp,                            & !parameter 59-62
    825         0.04_wp, 0.92_wp, 27.0_wp, 0.27_wp,                                    & !parameter 63-66
    826         0.003_wp, 0.006_wp, 0.012_wp, 0.018_wp,                                & !parameter 67-70
    827         1736000.0_wp, 1736000.0_wp, 1736000.0_wp,                              & !parameter 71-73
    828         0.11_wp, 0.11_wp, 0.11_wp, 0.87_wp, 0.7_wp,                            & !parameter 74-78
    829         27.0_wp, 0.0_wp, 1.5_wp, 0.86_wp, 5.0_wp, 1.0_wp,                      & !parameter 79-84
    830         0.005_wp, 0.01_wp, 0.5_wp, 0.79_wp, 2000000.0_wp, 103000.0_wp,         & !parameter 85-90
    831         900000.0_wp, 0.35_wp, 0.38_wp, 0.04_wp, 0.93_wp, 27.0_wp, 0.0_wp,      & !parameter 91-97
    832         0.003_wp, 0.006_wp, 0.012_wp, 0.018_wp, 1736000.0_wp,                  & !parameter 98-102
    833         1736000.0_wp, 1736000.0_wp, 0.11_wp, 0.11_wp, 0.11_wp,                 & !parameter 103-107
    834         0.87_wp, 0.7_wp, 27.0_wp, 0.0_wp, 0.0_wp, 1.5_wp,                      & !parameter 108-113
    835         0.86_wp, 5.0_wp, 0.0_wp,                                               & !parameter 114-116
    836         299.15_wp, 293.15_wp, 0.8_wp, 0.6_wp, 3.0_wp,                          & !parameter 117-121
    837         0.1_wp, 0.5_wp, 0.0_wp, 2.5_wp, 165000.0_wp, 4.5_wp,                   & !parameter 122-127
    838         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
    839         20000.0_wp, 20000.0_wp, 20000.0_wp, 23.0_wp, 23.0_wp, 10.0_wp,         & !parameter 0-5
    840         1.0_wp, 0.005_wp, 0.01_wp, 0.41_wp, 0.7_wp, 2000000.0_wp,              & !parameter 6-11
    841         103000.0_wp, 900000.0_wp, 0.35_wp, 0.14_wp, 0.035_wp, 4.0_wp,          & !parameter 12-17
    842         0.01_wp, 0.001_wp, 0.75_wp,                                            & !parameter 18-20
    843         0.005_wp, 0.01_wp, 0.41_wp, 0.7_wp, 2000000.0_wp,                      & !parameter 21-25
    844         103000.0_wp, 900000.0_wp, 0.35_wp,                                     & !parameter 26-28                     
    845         0.14_wp, 0.035_wp, 0.92_wp,                                            & !parameter 29-31       
    846         27.0_wp, 0.25_wp, 0.003_wp, 0.006_wp, 0.012_wp, 0.018_wp,              & !parameter 32-37
    847         1736000.0_wp, 1736000.0_wp, 1736000.0_wp,                              & !parameter 38-40
    848         0.037_wp, 0.037_wp, 0.037_wp, 0.8_wp,                                  & !parameter 41-44
    849         0.6_wp, 27.0_wp, 0.0_wp, 1.5_wp, 0.86_wp,                              & !parameter 45-49
    850         5.0_wp, 0.001_wp, 0.0001_wp, 0.7_wp, 0.005_wp,                         & !parameter 50-54
    851         0.01_wp, 0.41_wp, 0.7_wp, 2000000.0_wp,                                & !parameter 55-58
    852         103000.0_wp, 900000.0_wp, 0.35_wp, 0.14_wp,                            & !parameter 59-62
    853         0.035_wp, 0.92_wp, 27.0_wp, 0.3_wp,                                    & !parameter 63-66
    854         0.003_wp, 0.006_wp, 0.012_wp, 0.018_wp,                                & !parameter 67-70
    855         1736000.0_wp, 1736000.0_wp, 1736000.0_wp,                              & !parameter 71-73
    856         0.037_wp, 0.037_wp, 0.037_wp, 0.8_wp, 0.6_wp,                          & !parameter 74-78
    857         27.0_wp, 0.0_wp, 1.5_wp, 0.86_wp, 5.0_wp, 1.0_wp,                      & !parameter 79-84
    858         0.005_wp, 0.01_wp, 0.41_wp, 0.7_wp, 2000000.0_wp, 103000.0_wp,         & !parameter 85-90
    859         900000.0_wp, 0.35_wp, 0.14_wp, 0.035_wp, 0.93_wp, 27.0_wp, 0.0_wp,     & !parameter 91-97
    860         0.003_wp, 0.006_wp, 0.012_wp, 0.018_wp, 1736000.0_wp,                  & !parameter 98-102
    861         1736000.0_wp, 1736000.0_wp, 0.037_wp, 0.037_wp, 0.037_wp,              & !parameter 103-107
    862         0.8_wp, 0.6_wp, 27.0_wp, 0.0_wp, 0.0_wp, 1.5_wp,                       & !parameter 108-113
    863         0.86_wp, 5.0_wp, 0.0_wp,                                               & !parameter 114-116
    864         299.15_wp, 293.15_wp, 0.8_wp, 0.5_wp, 0.6_wp,                          & !parameter 117-121
    865         0.1_wp, 0.5_wp, 0.8_wp, 2.5_wp, 80000.0_wp, 4.5_wp,                    & !parameter 122-127
    866         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
    867         20000.0_wp, 20000.0_wp, 20000.0_wp, 23.0_wp, 23.0_wp, 10.0_wp,         & !parameter 0-5
    868         1.0_wp, 0.005_wp, 0.01_wp, 0.39_wp, 0.63_wp, 2200000.0_wp,             & !parameter 6-11
    869         1400000.0_wp, 1300000.0_wp, 0.35_wp, 0.8_wp, 2.1_wp, 4.0_wp,           & !parameter 12-17
    870         0.01_wp, 0.001_wp, 0.55_wp,                                            & !parameter 18-20
    871         0.005_wp, 0.01_wp, 0.39_wp, 0.63_wp, 2200000.0_wp,                     & !parameter 21-25
    872         1400000.0_wp, 1300000.0_wp, 0.35_wp,                                   & !parameter 26-28                     
    873         0.8_wp, 2.1_wp, 0.93_wp,                                               & !parameter 29-31       
    874         27.0_wp, 0.45_wp, 0.003_wp, 0.006_wp, 0.012_wp, 0.018_wp,              & !parameter 32-37
    875         1736000.0_wp, 1736000.0_wp, 1736000.0_wp,                              & !parameter 38-40
    876         0.57_wp, 0.57_wp, 0.57_wp, 0.91_wp,                                    & !parameter 41-44
    877         0.75_wp, 27.0_wp, 0.0_wp, 1.5_wp, 0.86_wp,                             & !parameter 45-49
    878         5.0_wp, 0.001_wp, 0.0001_wp, 0.5_wp, 0.005_wp,                         & !parameter 50-54
    879         0.01_wp, 0.39_wp, 0.63_wp, 2200000.0_wp,                               & !parameter 55-58
    880         1400000.0_wp, 1300000.0_wp, 0.35_wp, 0.8_wp,                           & !parameter 59-62
    881         2.1_wp, 0.93_wp, 27.0_wp, 0.5_wp,                                      & !parameter 63-66
    882         0.003_wp, 0.006_wp, 0.012_wp, 0.018_wp,                                & !parameter 67-70
    883         1736000.0_wp, 1736000.0_wp, 1736000.0_wp,                              & !parameter 71-73
    884         0.57_wp, 0.57_wp, 0.57_wp, 0.91_wp, 0.75_wp,                           & !parameter 74-78
    885         27.0_wp, 0.0_wp, 1.5_wp, 0.86_wp, 5.0_wp, 1.0_wp,                      & !parameter 79-84
    886         0.005_wp, 0.01_wp, 0.39_wp, 0.63_wp, 2200000.0_wp, 1400000.0_wp,       & !parameter 85-90
    887         1300000.0_wp, 0.35_wp, 0.8_wp, 2.1_wp, 0.93_wp, 27.0_wp, 0.0_wp,       & !parameter 91-97
    888         0.003_wp, 0.006_wp, 0.012_wp, 0.018_wp, 1736000.0_wp,                  & !parameter 98-102
    889         1736000.0_wp, 1736000.0_wp, 0.57_wp, 0.57_wp, 0.57_wp,                 & !parameter 103-107
    890         0.91_wp, 0.75_wp, 27.0_wp, 0.0_wp, 0.0_wp, 1.5_wp,                     & !parameter 108-113
    891         0.86_wp, 5.0_wp, 0.0_wp,                                               & !parameter 114-116
    892         299.15_wp, 293.15_wp, 0.8_wp, 0.76_wp, 5.0_wp,                         & !parameter 117-121
    893         0.1_wp, 1.5_wp, 0.0_wp, 3.5_wp, 370000.0_wp, 4.5_wp,                   & !parameter 122-127
    894         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
    895         20000.0_wp, 20000.0_wp, 20000.0_wp, 23.0_wp, 23.0_wp, 10.0_wp,         & !parameter 0-5
    896         1.0_wp, 0.005_wp, 0.01_wp, 0.31_wp, 0.43_wp, 2000000.0_wp,             & !parameter 6-11
    897         103000.0_wp, 900000.0_wp, 0.35_wp, 0.38_wp, 0.04_wp, 4.0_wp,           & !parameter 12-17
    898         0.01_wp, 0.001_wp, 0.55_wp,                                            & !parameter 18-20
    899         0.005_wp, 0.01_wp, 0.31_wp, 0.43_wp, 2000000.0_wp,                     & !parameter 21-25
    900         103000.0_wp, 900000.0_wp, 0.35_wp,                                     & !parameter 26-28                     
    901         0.38_wp, 0.04_wp, 0.92_wp,                                             & !parameter 29-31       
    902         27.0_wp, 0.45_wp, 0.003_wp, 0.006_wp, 0.012_wp, 0.018_wp,              & !parameter 32-37
    903         1736000.0_wp, 1736000.0_wp, 1736000.0_wp,                              & !parameter 38-40
    904         0.11_wp, 0.11_wp, 0.11_wp, 0.87_wp,                                    & !parameter 41-44
    905         0.7_wp, 27.0_wp, 0.0_wp, 1.5_wp, 0.86_wp,                              & !parameter 45-49
    906         5.0_wp, 0.001_wp, 0.0001_wp, 0.5_wp, 0.005_wp,                         & !parameter 50-54
    907         0.01_wp, 0.31_wp, 0.43_wp, 2000000.0_wp,                               & !parameter 55-58
    908         103000.0_wp, 900000.0_wp, 0.35_wp, 0.38_wp,                            & !parameter 59-62
    909         0.04_wp, 0.92_wp, 27.0_wp, 0.5_wp,                                     & !parameter 63-66
    910         0.003_wp, 0.006_wp, 0.012_wp, 0.018_wp,                                & !parameter 67-70
    911         1736000.0_wp, 1736000.0_wp, 1736000.0_wp,                              & !parameter 71-73
    912         0.11_wp, 0.11_wp, 0.11_wp, 0.87_wp, 0.7_wp,                            & !parameter 74-78
    913         27.0_wp, 0.0_wp, 1.5_wp, 0.86_wp, 5.0_wp, 1.0_wp,                      & !parameter 79-84
    914         0.005_wp, 0.01_wp, 0.31_wp, 0.43_wp, 2000000.0_wp, 103000.0_wp,        & !parameter 85-90
    915         900000.0_wp, 0.35_wp, 0.38_wp, 0.04_wp, 0.91_wp, 27.0_wp, 0.0_wp,      & !parameter 91-97
    916         0.003_wp, 0.006_wp, 0.012_wp, 0.018_wp, 1736000.0_wp,                  & !parameter 98-102
    917         1736000.0_wp, 1736000.0_wp, 0.11_wp, 0.11_wp, 0.11_wp,                 & !parameter 103-107
    918         0.87_wp, 0.7_wp, 27.0_wp, 0.0_wp, 0.0_wp, 1.5_wp,                      & !parameter 108-113
    919         0.86_wp, 5.0_wp, 0.0_wp,                                               & !parameter 114-116
    920         299.15_wp, 293.15_wp, 0.8_wp, 0.6_wp, 3.0_wp,                          & !parameter 117-121
    921         0.1_wp, 1.5_wp, 0.65_wp, 2.5_wp, 165000.0_wp, 4.5_wp,                  & !parameter 122-127
    922         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
    923         20000.0_wp, 20000.0_wp, 20000.0_wp, 23.0_wp, 23.0_wp, 10.0_wp,         & !parameter 0-5
    924         1.0_wp, 0.005_wp, 0.01_wp, 0.41_wp, 0.7_wp, 2000000.0_wp,              & !parameter 6-11
    925         103000.0_wp, 900000.0_wp, 0.35_wp, 0.14_wp, 0.035_wp, 4.0_wp,          & !parameter 12-17
    926         0.01_wp, 0.001_wp, 0.475_wp,                                           & !parameter 18-20
    927         0.005_wp, 0.01_wp, 0.41_wp, 0.7_wp, 2000000.0_wp,                      & !parameter 21-25
    928         103000.0_wp, 900000.0_wp, 0.35_wp,                                     & !parameter 26-28                     
    929         0.14_wp, 0.035_wp, 0.92_wp,                                            & !parameter 29-31       
    930         27.0_wp, 0.525_wp, 0.003_wp, 0.006_wp, 0.012_wp, 0.018_wp,             & !parameter 32-37
    931         1736000.0_wp, 1736000.0_wp, 1736000.0_wp,                              & !parameter 38-40
    932         0.037_wp, 0.037_wp, 0.037_wp, 0.8_wp,                                  & !parameter 41-44
    933         0.6_wp, 27.0_wp, 0.0_wp, 1.5_wp, 0.86_wp,                              & !parameter 45-49
    934         5.0_wp, 0.001_wp, 0.0001_wp, 0.425_wp, 0.005_wp,                       & !parameter 50-54
    935         0.01_wp, 0.41_wp, 0.7_wp, 2000000.0_wp,                                & !parameter 55-58
    936         103000.0_wp, 900000.0_wp, 0.35_wp, 0.14_wp,                            & !parameter 59-62
    937         0.035_wp, 0.92_wp, 27.0_wp, 0.575_wp,                                  & !parameter 63-66
    938         0.003_wp, 0.006_wp, 0.012_wp, 0.018_wp,                                & !parameter 67-70
    939         1736000.0_wp, 1736000.0_wp, 1736000.0_wp,                              & !parameter 71-73
    940         0.037_wp, 0.037_wp, 0.037_wp, 0.8_wp, 0.6_wp,                          & !parameter 74-78
    941         27.0_wp, 0.0_wp, 1.5_wp, 0.86_wp, 5.0_wp, 1.0_wp,                      & !parameter 79-84
    942         0.005_wp, 0.01_wp, 0.41_wp, 0.7_wp, 2000000.0_wp, 103000.0_wp,         & !parameter 85-90
    943         900000.0_wp, 0.35_wp, 0.14_wp, 0.035_wp, 0.91_wp, 27.0_wp, 0.0_wp,     & !parameter 91-97
    944         0.003_wp, 0.006_wp, 0.012_wp, 0.018_wp, 1736000.0_wp,                  & !parameter 98-102
    945         1736000.0_wp, 1736000.0_wp, 0.037_wp, 0.037_wp, 0.037_wp,              & !parameter 103-107
    946         0.8_wp, 0.6_wp, 27.0_wp, 0.0_wp, 0.0_wp, 1.5_wp,                       & !parameter 108-113
    947         0.86_wp, 5.0_wp, 0.0_wp,                                               & !parameter 114-116
    948         299.15_wp, 293.15_wp, 0.8_wp, 0.5_wp, 0.6_wp,                          & !parameter 117-121
    949         0.1_wp, 1.5_wp, 0.9_wp, 2.5_wp, 80000.0_wp, 4.5_wp,                    & !parameter 122-127
    950         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   
    951         20000.0_wp, 20000.0_wp, 20000.0_wp, 23.0_wp, 23.0_wp, 10.0_wp,         & !parameter 0-5
    952         1.0_wp, 0.29_wp, 0.295_wp, 0.695_wp, 0.985_wp, 1950400.0_wp,           & !parameter 6-11
    953         1848000.0_wp, 1848000.0_wp, 0.7_wp, 1.0_wp, 1.0_wp, 4.0_wp,            & !parameter 12-17
    954         0.01_wp, 0.001_wp, 1.0_wp,                                             & !parameter 18-20
    955         0.29_wp, 0.295_wp, 0.695_wp, 0.985_wp, 1950400.0_wp,                   & !parameter 21-25
    956         1848000.0_wp, 1848000.0_wp, 0.7_wp,                                    & !parameter 26-28                     
    957         1.0_wp, 1.0_wp, 0.9_wp,                                                & !parameter 29-31       
    958         27.0_wp, 0.0_wp, 0.003_wp, 0.006_wp, 0.012_wp, 0.018_wp,               & !parameter 32-37
    959         1736000.0_wp, 1736000.0_wp, 1736000.0_wp,                              & !parameter 38-40
    960         0.57_wp, 0.57_wp, 0.57_wp, 0.8_wp,                                     & !parameter 41-44
    961         0.6_wp, 27.0_wp, 0.0_wp, 1.5_wp, 0.86_wp,                              & !parameter 45-49
    962         5.0_wp, 0.001_wp, 0.0001_wp, 1.0_wp, 0.29_wp,                          & !parameter 50-54
    963         0.295_wp, 0.695_wp, 0.985_wp, 1950400.0_wp,                            & !parameter 55-58
    964         1848000.0_wp, 1848000.0_wp, 0.7_wp, 1.0_wp,                            & !parameter 59-62
    965         1.0_wp, 0.9_wp, 27.0_wp, 0.0_wp,                                       & !parameter 63-66
    966         0.003_wp, 0.006_wp, 0.012_wp, 0.018_wp,                                & !parameter 67-70
    967         1736000.0_wp, 1736000.0_wp, 1736000.0_wp,                              & !parameter 71-73
    968         0.57_wp, 0.57_wp, 0.57_wp, 0.8_wp, 0.6_wp,                             & !parameter 74-78
    969         27.0_wp, 0.0_wp, 1.5_wp, 0.86_wp, 5.0_wp, 1.0_wp,                      & !parameter 79-84
    970         0.29_wp, 0.295_wp, 0.695_wp, 0.985_wp, 1950400.0_wp, 1848000.0_wp,     & !parameter 85-90
    971         1848000.0_wp, 0.7_wp, 1.0_wp, 1.0_wp, 0.9_wp, 27.0_wp, 0.0_wp,         & !parameter 91-97
    972         0.003_wp, 0.006_wp, 0.012_wp, 0.018_wp, 1736000.0_wp,                  & !parameter 98-102
    973         1736000.0_wp, 1736000.0_wp, 0.57_wp, 0.57_wp, 0.57_wp,                 & !parameter 103-107
    974         0.8_wp, 0.6_wp, 27.0_wp, 0.0_wp, 0.0_wp, 1.5_wp,                       & !parameter 108-113
    975         0.86_wp, 5.0_wp, 0.0_wp,                                               & !parameter 114-116
    976         299.15_wp, 293.15_wp, 0.8_wp, 100.0_wp, 100.0_wp,                      & !parameter 117-121
    977         20.0_wp, 20.0_wp, 0.0_wp, 1.0_wp, 1.0_wp, 4.5_wp,                      & !parameter 122-127
    978         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)
    979                                                                        /),     &
    980                                                                (/134, 7/) )
     724       0.7_wp,         &  !< parameter 0   - wall fraction above ground floor level
     725       0.3_wp,         &  !< parameter 1   - window fraction above ground floor level
     726       0.0_wp,         &  !< parameter 2   - green fraction above ground floor level
     727       0.0_wp,         &  !< parameter 3   - green fraction roof above ground floor level
     728       1.5_wp,         &  !< parameter 4   - LAI roof
     729       1.5_wp,         &  !< parameter 5   - LAI on wall above ground floor level
     730       2200000.0_wp,   &  !< parameter 6   - heat capacity 1st/2nd wall layer above ground floor level
     731       1400000.0_wp,   &  !< parameter 7   - heat capacity 3rd wall layer above ground floor level
     732       1300000.0_wp,   &  !< parameter 8   - heat capacity 4th wall layer above ground floor level
     733       0.35_wp,        &  !< parameter 9   - thermal conductivity 1st/2nd wall layer above ground floor level
     734       0.8_wp,         &  !< parameter 10  - thermal conductivity 3rd wall layer above ground floor level
     735       2.1_wp,         &  !< parameter 11  - thermal conductivity 4th wall layer above ground floor level
     736       299.15_wp,      &  !< parameter 12  - indoor target summer temperature
     737       293.15_wp,      &  !< parameter 13  - indoor target winter temperature
     738       0.93_wp,        &  !< parameter 14  - wall emissivity above ground floor level
     739       0.86_wp,        &  !< parameter 15  - green emissivity above ground floor level
     740       0.91_wp,        &  !< parameter 16  - window emissivity above ground floor level
     741       0.75_wp,        &  !< parameter 17  - window transmissivity above ground floor level
     742       0.001_wp,       &  !< parameter 18  - z0 roughness above ground floor level
     743       0.0001_wp,      &  !< parameter 19  - z0h/z0g roughness heat/humidity above ground floor level
     744       4.0_wp,         &  !< parameter 20  - ground floor level height
     745       0.75_wp,        &  !< parameter 21  - wall fraction ground floor level
     746       0.25_wp,        &  !< parameter 22  - window fraction ground floor level
     747       0.0_wp,         &  !< parameter 23  - green fraction ground floor level
     748       0.0_wp,         &  !< parameter 24  - green fraction roof ground floor level
     749       1.5_wp,         &  !< parameter 25  - LAI on wall ground floor level
     750       2200000.0_wp,   &  !< parameter 26  - heat capacity 1st/2nd wall layer ground floor level
     751       1400000.0_wp,   &  !< parameter 27  - heat capacity 3rd wall layer ground floor level
     752       1300000.0_wp,   &  !< parameter 28  - heat capacity 4th wall layer ground floor level
     753       0.35_wp,        &  !< parameter 29  - thermal conductivity 1st/2nd wall layer ground floor level
     754       0.8_wp,         &  !< parameter 30  - thermal conductivity 3rd wall layer ground floor level
     755       2.1_wp,         &  !< parameter 31  - thermal conductivity 4th wall layer ground floor level
     756       0.93_wp,        &  !< parameter 32  - wall emissivity ground floor level
     757       0.91_wp,        &  !< parameter 33  - window emissivity ground floor level
     758       0.86_wp,        &  !< parameter 34  - green emissivity ground floor level
     759       0.75_wp,        &  !< parameter 35  - window transmissivity ground floor level
     760       0.01_wp,        &  !< parameter 36  - z0 roughness ground floor level
     761       0.001_wp,       &  !< parameter 37  - z0h/z0q roughness heat/humidity
     762       27.0_wp,        &  !< parameter 38  - wall albedo above ground floor level
     763       5.0_wp,         &  !< parameter 39  - green albedo above ground floor level
     764       27.0_wp,        &  !< parameter 40  - window albedo above ground floor level
     765       0.005_wp,       &  !< parameter 41  - 1st wall layer thickness above ground floor level
     766       0.01_wp,        &  !< parameter 42  - 2nd wall layer thickness above ground floor level
     767       0.39_wp,        &  !< parameter 43  - 3rd wall layer thickness above ground floor level
     768       0.63_wp,        &  !< parameter 44  - 4th wall layer thickness above ground floor level
     769       20000.0_wp,     &  !< parameter 45  - heat capacity wall surface
     770       23.0_wp,        &  !< parameter 46  - thermal conductivity of wall surface
     771       20000.0_wp,     &  !< parameter 47  - heat capacity of window surface
     772       20000.0_wp,     &  !< parameter 48  - heat capacity of green surface
     773       23.0_wp,        &  !< parameter 49  - thermal conductivity of window surface
     774       10.0_wp,        &  !< parameter 50  - thermal conductivty of green surface
     775       1.0_wp,         &  !< parameter 51  - wall fraction ground plate
     776       0.005_wp,       &  !< parameter 52  - 1st wall layer thickness ground plate
     777       0.01_wp,        &  !< parameter 53  - 2nd wall layer thickness ground plate
     778       0.39_wp,        &  !< parameter 54  - 3rd wall layer thickness ground plate
     779       0.63_wp,        &  !< parameter 55  - 4th wall layer thickness ground plate
     780       2200000.0_wp,   &  !< parameter 56  - heat capacity 1st/2nd wall layer ground plate
     781       1400000.0_wp,   &  !< parameter 57  - heat capacity 3rd wall layer ground plate
     782       1300000.0_wp,   &  !< parameter 58  - heat capacity 4th wall layer ground plate
     783       0.35_wp,        &  !< parameter 59  - thermal conductivity 1st/2nd wall layer ground plate
     784       0.8_wp,         &  !< parameter 60  - thermal conductivity 3rd wall layer ground plate
     785       2.1_wp,         &  !< parameter 61  - thermal conductivity 4th wall layer ground plate
     786       0.005_wp,       &  !< parameter 62  - 1st wall layer thickness ground floor level
     787       0.01_wp,        &  !< parameter 63  - 2nd wall layer thickness ground floor level
     788       0.39_wp,        &  !< parameter 64  - 3rd wall layer thickness ground floor level
     789       0.63_wp,        &  !< parameter 65  - 4th wall layer thickness ground floor level
     790       27.0_wp,        &  !< parameter 66  - wall albedo ground floor level
     791       0.003_wp,       &  !< parameter 67  - 1st window layer thickness ground floor level
     792       0.006_wp,       &  !< parameter 68  - 2nd window layer thickness ground floor level
     793       0.012_wp,       &  !< parameter 69  - 3rd window layer thickness ground floor level
     794       0.018_wp,       &  !< parameter 70  - 4th window layer thickness ground floor level
     795       1736000.0_wp,   &  !< parameter 71  - heat capacity 1st/2nd window layer ground floor level
     796       1736000.0_wp,   &  !< parameter 72  - heat capacity 3rd window layer ground floor level
     797       1736000.0_wp,   &  !< parameter 73  - heat capacity 4th window layer ground floor level
     798       0.57_wp,        &  !< parameter 74  - thermal conductivity 1st/2nd window layer ground floor level
     799       0.57_wp,        &  !< parameter 75  - thermal conductivity 3rd window layer ground floor level
     800       0.57_wp,        &  !< parameter 76  - thermal conductivity 4th window layer ground floor level
     801       27.0_wp,        &  !< parameter 77  - window albedo ground floor level
     802       5.0_wp,         &  !< parameter 78  - green albedo ground floor level
     803       0.003_wp,       &  !< parameter 79  - 1st window layer thickness above ground floor level
     804       0.006_wp,       &  !< parameter 80  - 2nd thickness window layer above ground floor level
     805       0.012_wp,       &  !< parameter 81  - 3rd window layer thickness above ground floor level
     806       0.018_wp,       &  !< parameter 82  - 4th window layer thickness above ground floor level
     807       1736000.0_wp,   &  !< parameter 83  - heat capacity 1st/2nd window layer above ground floor level
     808       1736000.0_wp,   &  !< parameter 84  - heat capacity 3rd window layer above ground floor level
     809       1736000.0_wp,   &  !< parameter 85  - heat capacity 4th window layer above ground floor level
     810       0.57_wp,        &  !< parameter 86  - thermal conductivity 1st/2nd window layer above ground floor level
     811       0.57_wp,        &  !< parameter 87  - thermal conductivity 3rd window layer above ground floor level
     812       0.57_wp,        &  !< parameter 88  - thermal conductivity 4th window layer above ground floor level
     813       1.0_wp,         &  !< parameter 89  - wall fraction roof
     814       0.005_wp,       &  !< parameter 90  - 1st wall layer thickness roof
     815       0.01_wp,        &  !< parameter 91  - 2nd wall layer thickness roof
     816       0.31_wp,        &  !< parameter 92  - 3rd wall layer thickness roof
     817       0.63_wp,        &  !< parameter 93  - 4th wall layer thickness roof
     818       2200000.0_wp,   &  !< parameter 94  - heat capacity 1st/2nd wall layer roof
     819       1400000.0_wp,   &  !< parameter 95  - heat capacity 3rd wall layer roof
     820       1300000.0_wp,   &  !< parameter 96  - heat capacity 4th wall layer roof
     821       0.35_wp,        &  !< parameter 97  - thermal conductivity 1st/2nd wall layer roof
     822       0.8_wp,         &  !< parameter 98  - thermal conductivity 3rd wall layer roof
     823       2.1_wp,         &  !< parameter 99  - thermal conductivity 4th wall layer roof
     824       0.93_wp,        &  !< parameter 100 - wall emissivity roof
     825       27.0_wp,        &  !< parameter 101 - wall albedo roof
     826       0.0_wp,         &  !< parameter 102 - window fraction roof
     827       0.003_wp,       &  !< parameter 103 - window 1st layer thickness roof
     828       0.006_wp,       &  !< parameter 104 - window 2nd layer thickness roof
     829       0.012_wp,       &  !< parameter 105 - window 3rd layer thickness roof
     830       0.018_wp,       &  !< parameter 106 - window 4th layer thickness roof
     831       1736000.0_wp,   &  !< parameter 107 - heat capacity 1st/2nd window layer roof
     832       1736000.0_wp,   &  !< parameter 108 - heat capacity 3rd window layer roof
     833       1736000.0_wp,   &  !< parameter 109 - heat capacity 4th window layer roof
     834       0.57_wp,        &  !< parameter 110 - thermal conductivity 1st/2nd window layer roof
     835       0.57_wp,        &  !< parameter 111 - thermal conductivity 3rd window layer roof
     836       0.57_wp,        &  !< parameter 112 - thermal conductivity 4th window layer roof
     837       0.91_wp,        &  !< parameter 113 - window emissivity roof
     838       0.75_wp,        &  !< parameter 114 - window transmissivity roof
     839       27.0_wp,        &  !< parameter 115 - window albedo roof
     840       0.86_wp,        &  !< parameter 116 - green emissivity roof
     841       5.0_wp,         &  !< parameter 117 - green albedo roof
     842       0.0_wp,         &  !< parameter 118 - green type roof
     843       0.8_wp,         &  !< parameter 119 - shading factor
     844       0.76_wp,        &  !< parameter 120 - g-value windows
     845       5.0_wp,         &  !< parameter 121 - u-value windows
     846       0.1_wp,         &  !< parameter 122 - basical airflow without occupancy of the room
     847       0.5_wp,         &  !< parameter 123 - additional airflow depend of occupancy of the room
     848       0.0_wp,         &  !< parameter 124 - heat recovery efficiency
     849       3.5_wp,         &  !< parameter 125 - dynamic parameter specific effective surface
     850       370000.0_wp,    &  !< parameter 126 - dynamic parameter innner heatstorage
     851       4.5_wp,         &  !< parameter 127 - ratio internal surface/floor area
     852       100000.0_wp,    &  !< parameter 128 - maximal heating capacity
     853       0.0_wp,         &  !< parameter 129 - maximal cooling capacity
     854       3.0_wp,         &  !< parameter 130 - additional internal heat gains dependent on occupancy of the room
     855       10.0_wp,        &  !< parameter 131 - basic internal heat gains without occupancy of the room
     856       3.0_wp,         &  !< parameter 132 - storey height
     857       0.2_wp,         &  !< parameter 133 - ceiling construction height
     858       0.73_wp,        &  !< parameter 0   - wall fraction above ground floor level
     859       0.27_wp,        &  !< parameter 1   - window fraction above ground floor level
     860       0.0_wp,         &  !< parameter 2   - green fraction above ground floor level
     861       0.0_wp,         &  !< parameter 3   - green fraction roof above ground floor level
     862       1.5_wp,         &  !< parameter 4   - LAI roof
     863       1.5_wp,         &  !< parameter 5   - LAI on wall above ground floor level
     864       2000000.0_wp,   &  !< parameter 6   - heat capacity 1st/2nd wall layer above ground floor level
     865       103000.0_wp,    &  !< parameter 7   - heat capacity 3rd wall layer above ground floor level
     866       900000.0_wp,    &  !< parameter 8   - heat capacity 4th wall layer above ground floor level
     867       0.35_wp,        &  !< parameter 9   - thermal conductivity 1st/2nd wall layer above ground floor level
     868       0.38_wp,        &  !< parameter 10  - thermal conductivity 3rd wall layer above ground floor level
     869       0.04_wp,        &  !< parameter 11  - thermal conductivity 4th wall layer above ground floor level
     870       299.15_wp,      &  !< parameter 12  - indoor target summer temperature
     871       293.15_wp,      &  !< parameter 13  - indoor target winter temperature
     872       0.92_wp,        &  !< parameter 14  - wall emissivity above ground floor level
     873       0.86_wp,        &  !< parameter 15  - green emissivity above ground floor level
     874       0.87_wp,        &  !< parameter 16  - window emissivity above ground floor level
     875       0.7_wp,         &  !< parameter 17  - window transmissivity above ground floor level
     876       0.001_wp,       &  !< parameter 18  - z0 roughness above ground floor level
     877       0.0001_wp,      &  !< parameter 19  - z0h/z0g roughness heat/humidity above ground floor level
     878       4.0_wp,         &  !< parameter 20  - ground floor level height
     879       0.78_wp,        &  !< parameter 21  - wall fraction ground floor level
     880       0.22_wp,        &  !< parameter 22  - window fraction ground floor level
     881       0.0_wp,         &  !< parameter 23  - green fraction ground floor level
     882       0.0_wp,         &  !< parameter 24  - green fraction roof ground floor level
     883       1.5_wp,         &  !< parameter 25  - LAI on wall ground floor level
     884       2000000.0_wp,   &  !< parameter 26  - heat capacity 1st/2nd wall layer ground floor level
     885       103000.0_wp,    &  !< parameter 27  - heat capacity 3rd wall layer ground floor level
     886       900000.0_wp,    &  !< parameter 28  - heat capacity 4th wall layer ground floor level
     887       0.35_wp,        &  !< parameter 29  - thermal conductivity 1st/2nd wall layer ground floor level
     888       0.38_wp,        &  !< parameter 30  - thermal conductivity 3rd wall layer ground floor level
     889       0.04_wp,        &  !< parameter 31  - thermal conductivity 4th wall layer ground floor level
     890       0.92_wp,        &  !< parameter 32  - wall emissivity ground floor level
     891       0.11_wp,        &  !< parameter 33  - window emissivity ground floor level
     892       0.86_wp,        &  !< parameter 34  - green emissivity ground floor level
     893       0.7_wp,         &  !< parameter 35  - window transmissivity ground floor level
     894       0.01_wp,        &  !< parameter 36  - z0 roughness ground floor level
     895       0.001_wp,       &  !< parameter 37  - z0h/z0q roughness heat/humidity
     896       27.0_wp,        &  !< parameter 38  - wall albedo above ground floor level
     897       5.0_wp,         &  !< parameter 39  - green albedo above ground floor level
     898       27.0_wp,        &  !< parameter 40  - window albedo above ground floor level
     899       0.005_wp,       &  !< parameter 41  - 1st wall layer thickness above ground floor level
     900       0.01_wp,        &  !< parameter 42  - 2nd wall layer thickness above ground floor level
     901       0.31_wp,        &  !< parameter 43  - 3rd wall layer thickness above ground floor level
     902       0.43_wp,        &  !< parameter 44  - 4th wall layer thickness above ground floor level
     903       20000.0_wp,     &  !< parameter 45  - heat capacity wall surface
     904       23.0_wp,        &  !< parameter 46  - thermal conductivity of wall surface
     905       20000.0_wp,     &  !< parameter 47  - heat capacity of window surface
     906       20000.0_wp,     &  !< parameter 48  - heat capacity of green surface
     907       23.0_wp,        &  !< parameter 49  - thermal conductivity of window surface
     908       10.0_wp,        &  !< parameter 50  - thermal conductivty of green surface
     909       1.0_wp,         &  !< parameter 51  - wall fraction ground plate
     910       0.005_wp,       &  !< parameter 52  - 1st wall layer thickness ground plate
     911       0.01_wp,        &  !< parameter 53  - 2nd wall layer thickness ground plate
     912       0.31_wp,        &  !< parameter 54  - 3rd wall layer thickness ground plate
     913       0.42_wp,        &  !< parameter 55  - 4th wall layer thickness ground plate
     914       2000000.0_wp,   &  !< parameter 56  - heat capacity 1st/2nd wall layer ground plate
     915       103000.0_wp,    &  !< parameter 57  - heat capacity 3rd wall layer ground plate
     916       900000.0_wp,    &  !< parameter 58  - heat capacity 4th wall layer ground plate
     917       0.35_wp,        &  !< parameter 59  - thermal conductivity 1st/2nd wall layer ground plate
     918       0.38_wp,        &  !< parameter 60  - thermal conductivity 3rd wall layer ground plate
     919       0.04_wp,        &  !< parameter 61  - thermal conductivity 4th wall layer ground plate
     920       0.005_wp,       &  !< parameter 62  - 1st wall layer thickness ground floor level
     921       0.01_wp,        &  !< parameter 63  - 2nd wall layer thickness ground floor level
     922       0.31_wp,        &  !< parameter 64  - 3rd wall layer thickness ground floor level
     923       0.43_wp,        &  !< parameter 65  - 4th wall layer thickness ground floor level
     924       27.0_wp,        &  !< parameter 66  - wall albedo ground floor level
     925       0.003_wp,       &  !< parameter 67  - 1st window layer thickness ground floor level
     926       0.006_wp,       &  !< parameter 68  - 2nd window layer thickness ground floor level
     927       0.012_wp,       &  !< parameter 69  - 3rd window layer thickness ground floor level
     928       0.018_wp,       &  !< parameter 70  - 4th window layer thickness ground floor level
     929       1736000.0_wp,   &  !< parameter 71  - heat capacity 1st/2nd window layer ground floor level
     930       1736000.0_wp,   &  !< parameter 72  - heat capacity 3rd window layer ground floor level
     931       1736000.0_wp,   &  !< parameter 73  - heat capacity 4th window layer ground floor level
     932       0.11_wp,        &  !< parameter 74  - thermal conductivity 1st/2nd window layer ground floor level
     933       0.11_wp,        &  !< parameter 75  - thermal conductivity 3rd window layer ground floor level
     934       0.11_wp,        &  !< parameter 76  - thermal conductivity 4th window layer ground floor level
     935       27.0_wp,        &  !< parameter 77  - window albedo ground floor level
     936       5.0_wp,         &  !< parameter 78  - green albedo ground floor level
     937       0.003_wp,       &  !< parameter 79  - 1st window layer thickness above ground floor level
     938       0.006_wp,       &  !< parameter 80  - 2nd thickness window layer above ground floor level
     939       0.012_wp,       &  !< parameter 81  - 3rd window layer thickness above ground floor level
     940       0.018_wp,       &  !< parameter 82  - 4th window layer thickness above ground floor level
     941       1736000.0_wp,   &  !< parameter 83  - heat capacity 1st/2nd window layer above ground floor level
     942       1736000.0_wp,   &  !< parameter 84  - heat capacity 3rd window layer above ground floor level
     943       1736000.0_wp,   &  !< parameter 85  - heat capacity 4th window layer above ground floor level
     944       0.11_wp,        &  !< parameter 86  - thermal conductivity 1st/2nd window layer above ground floor level
     945       0.11_wp,        &  !< parameter 87  - thermal conductivity 3rd window layer above ground floor level
     946       0.11_wp,        &  !< parameter 88  - thermal conductivity 4th window layer above ground floor level
     947       1.0_wp,         &  !< parameter 89  - wall fraction roof
     948       0.005_wp,       &  !< parameter 90  - 1st wall layer thickness roof
     949       0.01_wp,        &  !< parameter 91  - 2nd wall layer thickness roof
     950       0.5_wp,         &  !< parameter 92  - 3rd wall layer thickness roof
     951       0.79_wp,        &  !< parameter 93  - 4th wall layer thickness roof
     952       2000000.0_wp,   &  !< parameter 94  - heat capacity 1st/2nd wall layer roof
     953       103000.0_wp,    &  !< parameter 95  - heat capacity 3rd wall layer roof
     954       900000.0_wp,    &  !< parameter 96  - heat capacity 4th wall layer roof
     955       0.35_wp,        &  !< parameter 97  - thermal conductivity 1st/2nd wall layer roof
     956       0.38_wp,        &  !< parameter 98  - thermal conductivity 3rd wall layer roof
     957       0.04_wp,        &  !< parameter 99  - thermal conductivity 4th wall layer roof
     958       0.93_wp,        &  !< parameter 100 - wall emissivity roof
     959       27.0_wp,        &  !< parameter 101 - wall albedo roof
     960       0.0_wp,         &  !< parameter 102 - window fraction roof
     961       0.003_wp,       &  !< parameter 103 - window 1st layer thickness roof
     962       0.006_wp,       &  !< parameter 104 - window 2nd layer thickness roof
     963       0.012_wp,       &  !< parameter 105 - window 3rd layer thickness roof
     964       0.018_wp,       &  !< parameter 106 - window 4th layer thickness roof
     965       1736000.0_wp,   &  !< parameter 107 - heat capacity 1st/2nd window layer roof
     966       1736000.0_wp,   &  !< parameter 108 - heat capacity 3rd window layer roof
     967       1736000.0_wp,   &  !< parameter 109 - heat capacity 4th window layer roof
     968       0.11_wp,        &  !< parameter 110 - thermal conductivity 1st/2nd window layer roof
     969       0.11_wp,        &  !< parameter 111 - thermal conductivity 3rd window layer roof
     970       0.11_wp,        &  !< parameter 112 - thermal conductivity 4th window layer roof
     971       0.87_wp,        &  !< parameter 113 - window emissivity roof
     972       0.7_wp,         &  !< parameter 114 - window transmissivity roof
     973       27.0_wp,        &  !< parameter 115 - window albedo roof
     974       0.86_wp,        &  !< parameter 116 - green emissivity roof
     975       5.0_wp,         &  !< parameter 117 - green albedo roof
     976       0.0_wp,         &  !< parameter 118 - green type roof
     977       0.8_wp,         &  !< parameter 119 - shading factor
     978       0.6_wp,         &  !< parameter 120 - g-value windows
     979       3.0_wp,         &  !< parameter 121 - u-value windows
     980       0.1_wp,         &  !< parameter 122 - basical airflow without occupancy of the room
     981       0.5_wp,         &  !< parameter 123 - additional airflow depend of occupancy of the room
     982       0.0_wp,         &  !< parameter 124 - heat recovery efficiency
     983       2.5_wp,         &  !< parameter 125 - dynamic parameter specific effective surface
     984       165000.0_wp,    &  !< parameter 126 - dynamic parameter innner heatstorage
     985       4.5_wp,         &  !< parameter 127 - ratio internal surface/floor area
     986       100000.0_wp,    &  !< parameter 128 - maximal heating capacity
     987       0.0_wp,         &  !< parameter 129 - maximal cooling capacity
     988       4.0_wp,         &  !< parameter 130 - additional internal heat gains dependent on occupancy of the room
     989       8.0_wp,         &  !< parameter 131 - basic internal heat gains without occupancy of the room
     990       3.0_wp,         &  !< parameter 132 - storey height
     991       0.2_wp,         &  !< parameter 133 - ceiling construction height
     992       0.7_wp,         &  !< parameter 0   - wall fraction above ground floor level
     993       0.3_wp,         &  !< parameter 1   - window fraction above ground floor level
     994       0.0_wp,         &  !< parameter 2   - green fraction above ground floor level
     995       0.0_wp,         &  !< parameter 3   - green fraction roof above ground floor level
     996       1.5_wp,         &  !< parameter 4   - LAI roof
     997       1.5_wp,         &  !< parameter 5   - LAI on wall above ground floor level
     998       2000000.0_wp,   &  !< parameter 6   - heat capacity 1st/2nd wall layer above ground floor level
     999       103000.0_wp,    &  !< parameter 7   - heat capacity 3rd wall layer above ground floor level
     1000       900000.0_wp,    &  !< parameter 8   - heat capacity 4th wall layer above ground floor level
     1001       0.35_wp,        &  !< parameter 9   - thermal conductivity 1st/2nd wall layer above ground floor level
     1002       0.14_wp,        &  !< parameter 10  - thermal conductivity 3rd wall layer above ground floor level
     1003       0.035_wp,       &  !< parameter 11  - thermal conductivity 4th wall layer above ground floor level
     1004       299.15_wp,      &  !< parameter 12  - indoor target summer temperature
     1005       293.15_wp,      &  !< parameter 13  - indoor target winter temperature
     1006       0.92_wp,        &  !< parameter 14  - wall emissivity above ground floor level
     1007       0.86_wp,        &  !< parameter 15  - green emissivity above ground floor level
     1008       0.8_wp,         &  !< parameter 16  - window emissivity above ground floor level
     1009       0.6_wp,         &  !< parameter 17  - window transmissivity above ground floor level
     1010       0.001_wp,       &  !< parameter 18  - z0 roughness above ground floor level
     1011       0.0001_wp,      &  !< parameter 19  - z0h/z0g roughness heat/humidity above ground floor level
     1012       4.0_wp,         &  !< parameter 20  - ground floor level height
     1013       0.75_wp,        &  !< parameter 21  - wall fraction ground floor level
     1014       0.25_wp,        &  !< parameter 22  - window fraction ground floor level
     1015       0.0_wp,         &  !< parameter 23  - green fraction ground floor level
     1016       0.0_wp,         &  !< parameter 24  - green fraction roof ground floor level
     1017       1.5_wp,         &  !< parameter 25  - LAI on wall ground floor level
     1018       2000000.0_wp,   &  !< parameter 26  - heat capacity 1st/2nd wall layer ground floor level
     1019       103000.0_wp,    &  !< parameter 27  - heat capacity 3rd wall layer ground floor level
     1020       900000.0_wp,    &  !< parameter 28  - heat capacity 4th wall layer ground floor level
     1021       0.35_wp,        &  !< parameter 29  - thermal conductivity 1st/2nd wall layer ground floor level
     1022       0.14_wp,        &  !< parameter 30  - thermal conductivity 3rd wall layer ground floor level
     1023       0.035_wp,       &  !< parameter 31  - thermal conductivity 4th wall layer ground floor level
     1024       0.92_wp,        &  !< parameter 32  - wall emissivity ground floor level
     1025       0.8_wp,         &  !< parameter 33  - window emissivity ground floor level
     1026       0.86_wp,        &  !< parameter 34  - green emissivity ground floor level
     1027       0.6_wp,         &  !< parameter 35  - window transmissivity ground floor level
     1028       0.01_wp,        &  !< parameter 36  - z0 roughness ground floor level
     1029       0.001_wp,       &  !< parameter 37  - z0h/z0q roughness heat/humidity
     1030       27.0_wp,        &  !< parameter 38  - wall albedo above ground floor level
     1031       5.0_wp,         &  !< parameter 39  - green albedo above ground floor level
     1032       27.0_wp,        &  !< parameter 40  - window albedo above ground floor level
     1033       0.005_wp,       &  !< parameter 41  - 1st wall layer thickness above ground floor level
     1034       0.01_wp,        &  !< parameter 42  - 2nd wall layer thickness above ground floor level
     1035       0.41_wp,        &  !< parameter 43  - 3rd wall layer thickness above ground floor level
     1036       0.7_wp,         &  !< parameter 44  - 4th wall layer thickness above ground floor level
     1037       20000.0_wp,     &  !< parameter 45  - heat capacity wall surface
     1038       23.0_wp,        &  !< parameter 46  - thermal conductivity of wall surface
     1039       20000.0_wp,     &  !< parameter 47  - heat capacity of window surface
     1040       20000.0_wp,     &  !< parameter 48  - heat capacity of green surface
     1041       23.0_wp,        &  !< parameter 49  - thermal conductivity of window surface
     1042       10.0_wp,        &  !< parameter 50  - thermal conductivty of green surface
     1043       1.0_wp,         &  !< parameter 51  - wall fraction ground plate
     1044       0.005_wp,       &  !< parameter 52  - 1st wall layer thickness ground plate
     1045       0.01_wp,        &  !< parameter 53  - 2nd wall layer thickness ground plate
     1046       0.41_wp,        &  !< parameter 54  - 3rd wall layer thickness ground plate
     1047       0.7_wp,         &  !< parameter 55  - 4th wall layer thickness ground plate
     1048       2000000.0_wp,   &  !< parameter 56  - heat capacity 1st/2nd wall layer ground plate
     1049       103000.0_wp,    &  !< parameter 57  - heat capacity 3rd wall layer ground plate
     1050       900000.0_wp,    &  !< parameter 58  - heat capacity 4th wall layer ground plate
     1051       0.35_wp,        &  !< parameter 59  - thermal conductivity 1st/2nd wall layer ground plate
     1052       0.14_wp,        &  !< parameter 60  - thermal conductivity 3rd wall layer ground plate
     1053       0.035_wp,       &  !< parameter 61  - thermal conductivity 4th wall layer ground plate
     1054       0.005_wp,       &  !< parameter 62  - 1st wall layer thickness ground floor level
     1055       0.01_wp,        &  !< parameter 63  - 2nd wall layer thickness ground floor level
     1056       0.41_wp,        &  !< parameter 64  - 3rd wall layer thickness ground floor level
     1057       0.7_wp,         &  !< parameter 65  - 4th wall layer thickness ground floor level
     1058       27.0_wp,        &  !< parameter 66  - wall albedo ground floor level
     1059       0.003_wp,       &  !< parameter 67  - 1st window layer thickness ground floor level
     1060       0.006_wp,       &  !< parameter 68  - 2nd window layer thickness ground floor level
     1061       0.012_wp,       &  !< parameter 69  - 3rd window layer thickness ground floor level
     1062       0.018_wp,       &  !< parameter 70  - 4th window layer thickness ground floor level
     1063       1736000.0_wp,   &  !< parameter 71  - heat capacity 1st/2nd window layer ground floor level
     1064       1736000.0_wp,   &  !< parameter 72  - heat capacity 3rd window layer ground floor level
     1065       1736000.0_wp,   &  !< parameter 73  - heat capacity 4th window layer ground floor level
     1066       0.037_wp,       &  !< parameter 74  - thermal conductivity 1st/2nd window layer ground floor level
     1067       0.037_wp,       &  !< parameter 75  - thermal conductivity 3rd window layer ground floor level
     1068       0.037_wp,       &  !< parameter 76  - thermal conductivity 4th window layer ground floor level
     1069       27.0_wp,        &  !< parameter 77  - window albedo ground floor level
     1070       5.0_wp,         &  !< parameter 78  - green albedo ground floor level
     1071       0.003_wp,       &  !< parameter 79  - 1st window layer thickness above ground floor level
     1072       0.006_wp,       &  !< parameter 80  - 2nd thickness window layer above ground floor level
     1073       0.012_wp,       &  !< parameter 81  - 3rd window layer thickness above ground floor level
     1074       0.018_wp,       &  !< parameter 82  - 4th window layer thickness above ground floor level
     1075       1736000.0_wp,   &  !< parameter 83  - heat capacity 1st/2nd window layer above ground floor level
     1076       1736000.0_wp,   &  !< parameter 84  - heat capacity 3rd window layer above ground floor level
     1077       1736000.0_wp,   &  !< parameter 85  - heat capacity 4th window layer above ground floor level
     1078       0.037_wp,       &  !< parameter 86  - thermal conductivity 1st/2nd window layer above ground floor level
     1079       0.037_wp,       &  !< parameter 87  - thermal conductivity 3rd window layer above ground floor level
     1080       0.037_wp,       &  !< parameter 88  - thermal conductivity 4th window layer above ground floor level
     1081       1.0_wp,         &  !< parameter 89  - wall fraction roof
     1082       0.005_wp,       &  !< parameter 90  - 1st wall layer thickness roof
     1083       0.01_wp,        &  !< parameter 91  - 2nd wall layer thickness roof
     1084       0.41_wp,        &  !< parameter 92  - 3rd wall layer thickness roof
     1085       0.7_wp,         &  !< parameter 93  - 4th wall layer thickness roof
     1086       2000000.0_wp,   &  !< parameter 94  - heat capacity 1st/2nd wall layer roof
     1087       103000.0_wp,    &  !< parameter 95  - heat capacity 3rd wall layer roof
     1088       900000.0_wp,    &  !< parameter 96  - heat capacity 4th wall layer roof
     1089       0.35_wp,        &  !< parameter 97  - thermal conductivity 1st/2nd wall layer roof
     1090       0.14_wp,        &  !< parameter 98  - thermal conductivity 3rd wall layer roof
     1091       0.035_wp,       &  !< parameter 99  - thermal conductivity 4th wall layer roof
     1092       0.93_wp,        &  !< parameter 100 - wall emissivity roof
     1093       27.0_wp,        &  !< parameter 101 - wall albedo roof
     1094       0.0_wp,         &  !< parameter 102 - window fraction roof
     1095       0.003_wp,       &  !< parameter 103 - window 1st layer thickness roof
     1096       0.006_wp,       &  !< parameter 104 - window 2nd layer thickness roof
     1097       0.012_wp,       &  !< parameter 105 - window 3rd layer thickness roof
     1098       0.018_wp,       &  !< parameter 106 - window 4th layer thickness roof
     1099       1736000.0_wp,   &  !< parameter 107 - heat capacity 1st/2nd window layer roof
     1100       1736000.0_wp,   &  !< parameter 108 - heat capacity 3rd window layer roof
     1101       1736000.0_wp,   &  !< parameter 109 - heat capacity 4th window layer roof
     1102       0.037_wp,       &  !< parameter 110 - thermal conductivity 1st/2nd window layer roof
     1103       0.037_wp,       &  !< parameter 111 - thermal conductivity 3rd window layer roof
     1104       0.037_wp,       &  !< parameter 112 - thermal conductivity 4th window layer roof
     1105       0.8_wp,         &  !< parameter 113 - window emissivity roof
     1106       0.6_wp,         &  !< parameter 114 - window transmissivity roof
     1107       27.0_wp,        &  !< parameter 115 - window albedo roof
     1108       0.86_wp,        &  !< parameter 116 - green emissivity roof
     1109       5.0_wp,         &  !< parameter 117 - green albedo roof
     1110       0.0_wp,         &  !< parameter 118 - green type roof
     1111       0.8_wp,         &  !< parameter 119 - shading factor
     1112       0.5_wp,         &  !< parameter 120 - g-value windows
     1113       0.6_wp,         &  !< parameter 121 - u-value windows
     1114       0.1_wp,         &  !< parameter 122 - basical airflow without occupancy of the room
     1115       0.5_wp,         &  !< parameter 123 - additional airflow depend of occupancy of the room
     1116       0.8_wp,         &  !< parameter 124 - heat recovery efficiency
     1117       2.5_wp,         &  !< parameter 125 - dynamic parameter specific effective surface
     1118       80000.0_wp,     &  !< parameter 126 - dynamic parameter innner heatstorage
     1119       4.5_wp,         &  !< parameter 127 - ratio internal surface/floor area
     1120       100000.0_wp,    &  !< parameter 128 - maximal heating capacity
     1121       0.0_wp,         &  !< parameter 129 - maximal cooling capacity
     1122       3.0_wp,         &  !< parameter 130 - additional internal heat gains dependent on occupancy of the room
     1123       8.0_wp,         &  !< parameter 131 - basic internal heat gains without occupancy of the room
     1124       3.0_wp,         &  !< parameter 132 - storey height
     1125       0.2_wp,         &  !< parameter 133 - ceiling construction height
     1126       0.5_wp,         &  !< parameter 0   - wall fraction above ground floor level
     1127       0.5_wp,         &  !< parameter 1   - window fraction above ground floor level
     1128       0.0_wp,         &  !< parameter 2   - green fraction above ground floor level
     1129       0.0_wp,         &  !< parameter 3   - green fraction roof above ground floor level
     1130       1.5_wp,         &  !< parameter 4   - LAI roof
     1131       1.5_wp,         &  !< parameter 5   - LAI on wall above ground floor level
     1132       2200000.0_wp,   &  !< parameter 6   - heat capacity 1st/2nd wall layer above ground floor level
     1133       1400000.0_wp,   &  !< parameter 7   - heat capacity 3rd wall layer above ground floor level
     1134       1300000.0_wp,   &  !< parameter 8   - heat capacity 4th wall layer above ground floor level
     1135       0.35_wp,        &  !< parameter 9   - thermal conductivity 1st/2nd wall layer above ground floor level
     1136       0.8_wp,         &  !< parameter 10  - thermal conductivity 3rd wall layer above ground floor level
     1137       2.1_wp,         &  !< parameter 11  - thermal conductivity 4th wall layer above ground floor level
     1138       299.15_wp,      &  !< parameter 12  - indoor target summer temperature
     1139       293.15_wp,      &  !< parameter 13  - indoor target winter temperature
     1140       0.93_wp,        &  !< parameter 14  - wall emissivity above ground floor level
     1141       0.86_wp,        &  !< parameter 15  - green emissivity above ground floor level
     1142       0.91_wp,        &  !< parameter 16  - window emissivity above ground floor level
     1143       0.75_wp,        &  !< parameter 17  - window transmissivity above ground floor level
     1144       0.001_wp,       &  !< parameter 18  - z0 roughness above ground floor level
     1145       0.0001_wp,      &  !< parameter 19  - z0h/z0g roughness heat/humidity above ground floor level
     1146       4.0_wp,         &  !< parameter 20  - ground floor level height
     1147       0.55_wp,        &  !< parameter 21  - wall fraction ground floor level
     1148       0.45_wp,        &  !< parameter 22  - window fraction ground floor level
     1149       0.0_wp,         &  !< parameter 23  - green fraction ground floor level
     1150       0.0_wp,         &  !< parameter 24  - green fraction roof ground floor level
     1151       1.5_wp,         &  !< parameter 25  - LAI on wall ground floor level
     1152       2200000.0_wp,   &  !< parameter 26  - heat capacity 1st/2nd wall layer ground floor level
     1153       1400000.0_wp,   &  !< parameter 27  - heat capacity 3rd wall layer ground floor level
     1154       1300000.0_wp,   &  !< parameter 28  - heat capacity 4th wall layer ground floor level
     1155       0.35_wp,        &  !< parameter 29  - thermal conductivity 1st/2nd wall layer ground floor level
     1156       0.8_wp,         &  !< parameter 30  - thermal conductivity 3rd wall layer ground floor level
     1157       2.1_wp,         &  !< parameter 31  - thermal conductivity 4th wall layer ground floor level
     1158       0.93_wp,        &  !< parameter 32  - wall emissivity ground floor level
     1159       0.91_wp,        &  !< parameter 33  - window emissivity ground floor level
     1160       0.86_wp,        &  !< parameter 34  - green emissivity ground floor level
     1161       0.75_wp,        &  !< parameter 35  - window transmissivity ground floor level
     1162       0.01_wp,        &  !< parameter 36  - z0 roughness ground floor level
     1163       0.001_wp,       &  !< parameter 37  - z0h/z0q roughness heat/humidity
     1164       27.0_wp,        &  !< parameter 38  - wall albedo above ground floor level
     1165       5.0_wp,         &  !< parameter 39  - green albedo above ground floor level
     1166       27.0_wp,        &  !< parameter 40  - window albedo above ground floor level
     1167       0.005_wp,       &  !< parameter 41  - 1st wall layer thickness above ground floor level
     1168       0.01_wp,        &  !< parameter 42  - 2nd wall layer thickness above ground floor level
     1169       0.39_wp,        &  !< parameter 43  - 3rd wall layer thickness above ground floor level
     1170       0.63_wp,        &  !< parameter 44  - 4th wall layer thickness above ground floor level
     1171       20000.0_wp,     &  !< parameter 45  - heat capacity wall surface
     1172       23.0_wp,        &  !< parameter 46  - thermal conductivity of wall surface
     1173       20000.0_wp,     &  !< parameter 47  - heat capacity of window surface
     1174       20000.0_wp,     &  !< parameter 48  - heat capacity of green surface
     1175       23.0_wp,        &  !< parameter 49  - thermal conductivity of window surface
     1176       10.0_wp,        &  !< parameter 50  - thermal conductivty of green surface
     1177       1.0_wp,         &  !< parameter 51  - wall fraction ground plate
     1178       0.005_wp,       &  !< parameter 52  - 1st wall layer thickness ground plate
     1179       0.01_wp,        &  !< parameter 53  - 2nd wall layer thickness ground plate
     1180       0.39_wp,        &  !< parameter 54  - 3rd wall layer thickness ground plate
     1181       0.63_wp,        &  !< parameter 55  - 4th wall layer thickness ground plate
     1182       2200000.0_wp,   &  !< parameter 56  - heat capacity 1st/2nd wall layer ground plate
     1183       1400000.0_wp,   &  !< parameter 57  - heat capacity 3rd wall layer ground plate
     1184       1300000.0_wp,   &  !< parameter 58  - heat capacity 4th wall layer ground plate
     1185       0.35_wp,        &  !< parameter 59  - thermal conductivity 1st/2nd wall layer ground plate
     1186       0.8_wp,         &  !< parameter 60  - thermal conductivity 3rd wall layer ground plate
     1187       2.1_wp,         &  !< parameter 61  - thermal conductivity 4th wall layer ground plate
     1188       0.005_wp,       &  !< parameter 62  - 1st wall layer thickness ground floor level
     1189       0.01_wp,        &  !< parameter 63  - 2nd wall layer thickness ground floor level
     1190       0.39_wp,        &  !< parameter 64  - 3rd wall layer thickness ground floor level
     1191       0.63_wp,        &  !< parameter 65  - 4th wall layer thickness ground floor level
     1192       27.0_wp,        &  !< parameter 66  - wall albedo ground floor level
     1193       0.003_wp,       &  !< parameter 67  - 1st window layer thickness ground floor level
     1194       0.006_wp,       &  !< parameter 68  - 2nd window layer thickness ground floor level
     1195       0.012_wp,       &  !< parameter 69  - 3rd window layer thickness ground floor level
     1196       0.018_wp,       &  !< parameter 70  - 4th window layer thickness ground floor level
     1197       1736000.0_wp,   &  !< parameter 71  - heat capacity 1st/2nd window layer ground floor level
     1198       1736000.0_wp,   &  !< parameter 72  - heat capacity 3rd window layer ground floor level
     1199       1736000.0_wp,   &  !< parameter 73  - heat capacity 4th window layer ground floor level
     1200       0.57_wp,        &  !< parameter 74  - thermal conductivity 1st/2nd window layer ground floor level
     1201       0.57_wp,        &  !< parameter 75  - thermal conductivity 3rd window layer ground floor level
     1202       0.57_wp,        &  !< parameter 76  - thermal conductivity 4th window layer ground floor level
     1203       27.0_wp,        &  !< parameter 77  - window albedo ground floor level
     1204       5.0_wp,         &  !< parameter 78  - green albedo ground floor level
     1205       0.003_wp,       &  !< parameter 79  - 1st window layer thickness above ground floor level
     1206       0.006_wp,       &  !< parameter 80  - 2nd thickness window layer above ground floor level
     1207       0.012_wp,       &  !< parameter 81  - 3rd window layer thickness above ground floor level
     1208       0.018_wp,       &  !< parameter 82  - 4th window layer thickness above ground floor level
     1209       1736000.0_wp,   &  !< parameter 83  - heat capacity 1st/2nd window layer above ground floor level
     1210       1736000.0_wp,   &  !< parameter 84  - heat capacity 3rd window layer above ground floor level
     1211       1736000.0_wp,   &  !< parameter 85  - heat capacity 4th window layer above ground floor level
     1212       0.57_wp,        &  !< parameter 86  - thermal conductivity 1st/2nd window layer above ground floor level
     1213       0.57_wp,        &  !< parameter 87  - thermal conductivity 3rd window layer above ground floor level
     1214       0.57_wp,        &  !< parameter 88  - thermal conductivity 4th window layer above ground floor level
     1215       1.0_wp,         &  !< parameter 89  - wall fraction roof
     1216       0.005_wp,       &  !< parameter 90  - 1st wall layer thickness roof
     1217       0.01_wp,        &  !< parameter 91  - 2nd wall layer thickness roof
     1218       0.39_wp,        &  !< parameter 92  - 3rd wall layer thickness roof
     1219       0.63_wp,        &  !< parameter 93  - 4th wall layer thickness roof
     1220       2200000.0_wp,   &  !< parameter 94  - heat capacity 1st/2nd wall layer roof
     1221       1400000.0_wp,   &  !< parameter 95  - heat capacity 3rd wall layer roof
     1222       1300000.0_wp,   &  !< parameter 96  - heat capacity 4th wall layer roof
     1223       0.35_wp,        &  !< parameter 97  - thermal conductivity 1st/2nd wall layer roof
     1224       0.8_wp,         &  !< parameter 98  - thermal conductivity 3rd wall layer roof
     1225       2.1_wp,         &  !< parameter 99  - thermal conductivity 4th wall layer roof
     1226       0.93_wp,        &  !< parameter 100 - wall emissivity roof
     1227       27.0_wp,        &  !< parameter 101 - wall albedo roof
     1228       0.0_wp,         &  !< parameter 102 - window fraction roof
     1229       0.003_wp,       &  !< parameter 103 - window 1st layer thickness roof
     1230       0.006_wp,       &  !< parameter 104 - window 2nd layer thickness roof
     1231       0.012_wp,       &  !< parameter 105 - window 3rd layer thickness roof
     1232       0.018_wp,       &  !< parameter 106 - window 4th layer thickness roof
     1233       1736000.0_wp,   &  !< parameter 107 - heat capacity 1st/2nd window layer roof
     1234       1736000.0_wp,   &  !< parameter 108 - heat capacity 3rd window layer roof
     1235       1736000.0_wp,   &  !< parameter 109 - heat capacity 4th window layer roof
     1236       0.57_wp,        &  !< parameter 110 - thermal conductivity 1st/2nd window layer roof
     1237       0.57_wp,        &  !< parameter 111 - thermal conductivity 3rd window layer roof
     1238       0.57_wp,        &  !< parameter 112 - thermal conductivity 4th window layer roof
     1239       0.91_wp,        &  !< parameter 113 - window emissivity roof
     1240       0.75_wp,        &  !< parameter 114 - window transmissivity roof
     1241       27.0_wp,        &  !< parameter 115 - window albedo roof
     1242       0.86_wp,        &  !< parameter 116 - green emissivity roof
     1243       5.0_wp,         &  !< parameter 117 - green albedo roof
     1244       0.0_wp,         &  !< parameter 118 - green type roof
     1245       0.8_wp,         &  !< parameter 119 - shading factor
     1246       0.76_wp,        &  !< parameter 120 - g-value windows
     1247       5.0_wp,         &  !< parameter 121 - u-value windows
     1248       0.1_wp,         &  !< parameter 122 - basical airflow without occupancy of the room
     1249       1.5_wp,         &  !< parameter 123 - additional airflow depend of occupancy of the room
     1250       0.0_wp,         &  !< parameter 124 - heat recovery efficiency
     1251       3.5_wp,         &  !< parameter 125 - dynamic parameter specific effective surface
     1252       370000.0_wp,    &  !< parameter 126 - dynamic parameter innner heatstorage
     1253       4.5_wp,         &  !< parameter 127 - ratio internal surface/floor area
     1254       100000.0_wp,    &  !< parameter 128 - maximal heating capacity
     1255       0.0_wp,         &  !< parameter 129 - maximal cooling capacity
     1256       3.0_wp,         &  !< parameter 130 - additional internal heat gains dependent on occupancy of the room
     1257       10.0_wp,        &  !< parameter 131 - basic internal heat gains without occupancy of the room
     1258       3.0_wp,         &  !< parameter 132 - storey height
     1259       0.2_wp,         &  !< parameter 133 - ceiling construction height
     1260       0.5_wp,         &  !< parameter 0   - wall fraction above ground floor level
     1261       0.5_wp,         &  !< parameter 1   - window fraction above ground floor level
     1262       0.0_wp,         &  !< parameter 2   - green fraction above ground floor level
     1263       0.0_wp,         &  !< parameter 3   - green fraction roof above ground floor level
     1264       1.5_wp,         &  !< parameter 4   - LAI roof
     1265       1.5_wp,         &  !< parameter 5   - LAI on wall above ground floor level
     1266       2000000.0_wp,   &  !< parameter 6   - heat capacity 1st/2nd wall layer above ground floor level
     1267       103000.0_wp,    &  !< parameter 7   - heat capacity 3rd wall layer above ground floor level
     1268       900000.0_wp,    &  !< parameter 8   - heat capacity 4th wall layer above ground floor level
     1269       0.35_wp,        &  !< parameter 9   - thermal conductivity 1st/2nd wall layer above ground floor level
     1270       0.38_wp,        &  !< parameter 10  - thermal conductivity 3rd wall layer above ground floor level
     1271       0.04_wp,        &  !< parameter 11  - thermal conductivity 4th wall layer above ground floor level
     1272       299.15_wp,      &  !< parameter 12  - indoor target summer temperature
     1273       293.15_wp,      &  !< parameter 13  - indoor target winter temperature
     1274       0.92_wp,        &  !< parameter 14  - wall emissivity above ground floor level
     1275       0.86_wp,        &  !< parameter 15  - green emissivity above ground floor level
     1276       0.87_wp,        &  !< parameter 16  - window emissivity above ground floor level
     1277       0.7_wp,         &  !< parameter 17  - window transmissivity above ground floor level
     1278       0.001_wp,       &  !< parameter 18  - z0 roughness above ground floor level
     1279       0.0001_wp,      &  !< parameter 19  - z0h/z0g roughness heat/humidity above ground floor level
     1280       4.0_wp,         &  !< parameter 20  - ground floor level height
     1281       0.55_wp,        &  !< parameter 21  - wall fraction ground floor level
     1282       0.45_wp,        &  !< parameter 22  - window fraction ground floor level
     1283       0.0_wp,         &  !< parameter 23  - green fraction ground floor level
     1284       0.0_wp,         &  !< parameter 24  - green fraction roof ground floor level
     1285       1.5_wp,         &  !< parameter 25  - LAI on wall ground floor level
     1286       2000000.0_wp,   &  !< parameter 26  - heat capacity 1st/2nd wall layer ground floor level
     1287       103000.0_wp,    &  !< parameter 27  - heat capacity 3rd wall layer ground floor level
     1288       900000.0_wp,    &  !< parameter 28  - heat capacity 4th wall layer ground floor level
     1289       0.35_wp,        &  !< parameter 29  - thermal conductivity 1st/2nd wall layer ground floor level
     1290       0.38_wp,        &  !< parameter 30  - thermal conductivity 3rd wall layer ground floor level
     1291       0.04_wp,        &  !< parameter 31  - thermal conductivity 4th wall layer ground floor level
     1292       0.92_wp,        &  !< parameter 32  - wall emissivity ground floor level
     1293       0.87_wp,        &  !< parameter 33  - window emissivity ground floor level
     1294       0.86_wp,        &  !< parameter 34  - green emissivity ground floor level
     1295       0.7_wp,         &  !< parameter 35  - window transmissivity ground floor level
     1296       0.01_wp,        &  !< parameter 36  - z0 roughness ground floor level
     1297       0.001_wp,       &  !< parameter 37  - z0h/z0q roughness heat/humidity
     1298       27.0_wp,        &  !< parameter 38  - wall albedo above ground floor level
     1299       5.0_wp,         &  !< parameter 39  - green albedo above ground floor level
     1300       27.0_wp,        &  !< parameter 40  - window albedo above ground floor level
     1301       0.005_wp,       &  !< parameter 41  - 1st wall layer thickness above ground floor level
     1302       0.01_wp,        &  !< parameter 42  - 2nd wall layer thickness above ground floor level
     1303       0.31_wp,        &  !< parameter 43  - 3rd wall layer thickness above ground floor level
     1304       0.43_wp,        &  !< parameter 44  - 4th wall layer thickness above ground floor level
     1305       20000.0_wp,     &  !< parameter 45  - heat capacity wall surface
     1306       23.0_wp,        &  !< parameter 46  - thermal conductivity of wall surface
     1307       20000.0_wp,     &  !< parameter 47  - heat capacity of window surface
     1308       20000.0_wp,     &  !< parameter 48  - heat capacity of green surface
     1309       23.0_wp,        &  !< parameter 49  - thermal conductivity of window surface
     1310       10.0_wp,        &  !< parameter 50  - thermal conductivty of green surface
     1311       1.0_wp,         &  !< parameter 51  - wall fraction ground plate
     1312       0.005_wp,       &  !< parameter 52  - 1st wall layer thickness ground plate
     1313       0.01_wp,        &  !< parameter 53  - 2nd wall layer thickness ground plate
     1314       0.31_wp,        &  !< parameter 54  - 3rd wall layer thickness ground plate
     1315       0.43_wp,        &  !< parameter 55  - 4th wall layer thickness ground plate
     1316       2000000.0_wp,   &  !< parameter 56  - heat capacity 1st/2nd wall layer ground plate
     1317       103000.0_wp,    &  !< parameter 57  - heat capacity 3rd wall layer ground plate
     1318       900000.0_wp,    &  !< parameter 58  - heat capacity 4th wall layer ground plate
     1319       0.35_wp,        &  !< parameter 59  - thermal conductivity 1st/2nd wall layer ground plate
     1320       0.38_wp,        &  !< parameter 60  - thermal conductivity 3rd wall layer ground plate
     1321       0.04_wp,        &  !< parameter 61  - thermal conductivity 4th wall layer ground plate
     1322       0.005_wp,       &  !< parameter 62  - 1st wall layer thickness ground floor level
     1323       0.01_wp,        &  !< parameter 63  - 2nd wall layer thickness ground floor level
     1324       0.31_wp,        &  !< parameter 64  - 3rd wall layer thickness ground floor level
     1325       0.43_wp,        &  !< parameter 65  - 4th wall layer thickness ground floor level
     1326       27.0_wp,        &  !< parameter 66  - wall albedo ground floor level
     1327       0.003_wp,       &  !< parameter 67  - 1st window layer thickness ground floor level
     1328       0.006_wp,       &  !< parameter 68  - 2nd window layer thickness ground floor level
     1329       0.012_wp,       &  !< parameter 69  - 3rd window layer thickness ground floor level
     1330       0.018_wp,       &  !< parameter 70  - 4th window layer thickness ground floor level
     1331       1736000.0_wp,   &  !< parameter 71  - heat capacity 1st/2nd window layer ground floor level
     1332       1736000.0_wp,   &  !< parameter 72  - heat capacity 3rd window layer ground floor level
     1333       1736000.0_wp,   &  !< parameter 73  - heat capacity 4th window layer ground floor level
     1334       0.11_wp,        &  !< parameter 74  - thermal conductivity 1st/2nd window layer ground floor level
     1335       0.11_wp,        &  !< parameter 75  - thermal conductivity 3rd window layer ground floor level
     1336       0.11_wp,        &  !< parameter 76  - thermal conductivity 4th window layer ground floor level
     1337       27.0_wp,        &  !< parameter 77  - window albedo ground floor level
     1338       5.0_wp,         &  !< parameter 78  - green albedo ground floor level
     1339       0.003_wp,       &  !< parameter 79  - 1st window layer thickness above ground floor level
     1340       0.006_wp,       &  !< parameter 80  - 2nd thickness window layer above ground floor level
     1341       0.012_wp,       &  !< parameter 81  - 3rd window layer thickness above ground floor level
     1342       0.018_wp,       &  !< parameter 82  - 4th window layer thickness above ground floor level
     1343       1736000.0_wp,   &  !< parameter 83  - heat capacity 1st/2nd window layer above ground floor level
     1344       1736000.0_wp,   &  !< parameter 84  - heat capacity 3rd window layer above ground floor level
     1345       1736000.0_wp,   &  !< parameter 85  - heat capacity 4th window layer above ground floor level
     1346       0.11_wp,        &  !< parameter 86  - thermal conductivity 1st/2nd window layer above ground floor level
     1347       0.11_wp,        &  !< parameter 87  - thermal conductivity 3rd window layer above ground floor level
     1348       0.11_wp,        &  !< parameter 88  - thermal conductivity 4th window layer above ground floor level
     1349       1.0_wp,         &  !< parameter 89  - wall fraction roof
     1350       0.005_wp,       &  !< parameter 90  - 1st wall layer thickness roof
     1351       0.01_wp,        &  !< parameter 91  - 2nd wall layer thickness roof
     1352       0.31_wp,        &  !< parameter 92  - 3rd wall layer thickness roof
     1353       0.43_wp,        &  !< parameter 93  - 4th wall layer thickness roof
     1354       2000000.0_wp,   &  !< parameter 94  - heat capacity 1st/2nd wall layer roof
     1355       103000.0_wp,    &  !< parameter 95  - heat capacity 3rd wall layer roof
     1356       900000.0_wp,    &  !< parameter 96  - heat capacity 4th wall layer roof
     1357       0.35_wp,        &  !< parameter 97  - thermal conductivity 1st/2nd wall layer roof
     1358       0.38_wp,        &  !< parameter 98  - thermal conductivity 3rd wall layer roof
     1359       0.04_wp,        &  !< parameter 99  - thermal conductivity 4th wall layer roof
     1360       0.91_wp,        &  !< parameter 100 - wall emissivity roof
     1361       27.0_wp,        &  !< parameter 101 - wall albedo roof
     1362       0.0_wp,         &  !< parameter 102 - window fraction roof
     1363       0.003_wp,       &  !< parameter 103 - window 1st layer thickness roof
     1364       0.006_wp,       &  !< parameter 104 - window 2nd layer thickness roof
     1365       0.012_wp,       &  !< parameter 105 - window 3rd layer thickness roof
     1366       0.018_wp,       &  !< parameter 106 - window 4th layer thickness roof
     1367       1736000.0_wp,   &  !< parameter 107 - heat capacity 1st/2nd window layer roof
     1368       1736000.0_wp,   &  !< parameter 108 - heat capacity 3rd window layer roof
     1369       1736000.0_wp,   &  !< parameter 109 - heat capacity 4th window layer roof
     1370       0.11_wp,        &  !< parameter 110 - thermal conductivity 1st/2nd window layer roof
     1371       0.11_wp,        &  !< parameter 111 - thermal conductivity 3rd window layer roof
     1372       0.11_wp,        &  !< parameter 112 - thermal conductivity 4th window layer roof
     1373       0.87_wp,        &  !< parameter 113 - window emissivity roof
     1374       0.7_wp,         &  !< parameter 114 - window transmissivity roof
     1375       27.0_wp,        &  !< parameter 115 - window albedo roof
     1376       0.86_wp,        &  !< parameter 116 - green emissivity roof
     1377       5.0_wp,         &  !< parameter 117 - green albedo roof
     1378       0.0_wp,         &  !< parameter 118 - green type roof
     1379       0.8_wp,         &  !< parameter 119 - shading factor
     1380       0.6_wp,         &  !< parameter 120 - g-value windows
     1381       3.0_wp,         &  !< parameter 121 - u-value windows
     1382       0.1_wp,         &  !< parameter 122 - basical airflow without occupancy of the room
     1383       1.5_wp,         &  !< parameter 123 - additional airflow depend of occupancy of the room
     1384       0.65_wp,        &  !< parameter 124 - heat recovery efficiency
     1385       2.5_wp,         &  !< parameter 125 - dynamic parameter specific effective surface
     1386       165000.0_wp,    &  !< parameter 126 - dynamic parameter innner heatstorage
     1387       4.5_wp,         &  !< parameter 127 - ratio internal surface/floor area
     1388       100000.0_wp,    &  !< parameter 128 - maximal heating capacity
     1389       0.0_wp,         &  !< parameter 129 - maximal cooling capacity
     1390       7.0_wp,         &  !< parameter 130 - additional internal heat gains dependent on occupancy of the room
     1391       20.0_wp,        &  !< parameter 131 - basic internal heat gains without occupancy of the room
     1392       3.0_wp,         &  !< parameter 132 - storey height
     1393       0.2_wp,         &  !< parameter 133 - ceiling construction height
     1394       0.425_wp,       &  !< parameter 0   - wall fraction above ground floor level
     1395       0.575_wp,       &  !< parameter 1   - window fraction above ground floor level
     1396       0.0_wp,         &  !< parameter 2   - green fraction above ground floor level
     1397       0.0_wp,         &  !< parameter 3   - green fraction roof above ground floor level
     1398       1.5_wp,         &  !< parameter 4   - LAI roof
     1399       1.5_wp,         &  !< parameter 5   - LAI on wall above ground floor level
     1400       2000000.0_wp,   &  !< parameter 6   - heat capacity 1st/2nd wall layer above ground floor level
     1401       103000.0_wp,    &  !< parameter 7   - heat capacity 3rd wall layer above ground floor level
     1402       900000.0_wp,    &  !< parameter 8   - heat capacity 4th wall layer above ground floor level
     1403       0.35_wp,        &  !< parameter 9   - thermal conductivity 1st/2nd wall layer above ground floor level
     1404       0.14_wp,        &  !< parameter 10  - thermal conductivity 3rd wall layer above ground floor level
     1405       0.035_wp,       &  !< parameter 11  - thermal conductivity 4th wall layer above ground floor level
     1406       299.15_wp,      &  !< parameter 12  - indoor target summer temperature
     1407       293.15_wp,      &  !< parameter 13  - indoor target winter temperature
     1408       0.92_wp,        &  !< parameter 14  - wall emissivity above ground floor level
     1409       0.86_wp,        &  !< parameter 15  - green emissivity above ground floor level
     1410       0.8_wp,         &  !< parameter 16  - window emissivity above ground floor level
     1411       0.6_wp,         &  !< parameter 17  - window transmissivity above ground floor level
     1412       0.001_wp,       &  !< parameter 18  - z0 roughness above ground floor level
     1413       0.0001_wp,      &  !< parameter 19  - z0h/z0g roughness heat/humidity above ground floor level
     1414       4.0_wp,         &  !< parameter 20  - ground floor level height
     1415       0.475_wp,       &  !< parameter 21  - wall fraction ground floor level
     1416       0.525_wp,       &  !< parameter 22  - window fraction ground floor level
     1417       0.0_wp,         &  !< parameter 23  - green fraction ground floor level
     1418       0.0_wp,         &  !< parameter 24  - green fraction roof ground floor level
     1419       1.5_wp,         &  !< parameter 25  - LAI on wall ground floor level
     1420       2000000.0_wp,   &  !< parameter 26  - heat capacity 1st/2nd wall layer ground floor level
     1421       103000.0_wp,    &  !< parameter 27  - heat capacity 3rd wall layer ground floor level
     1422       900000.0_wp,    &  !< parameter 28  - heat capacity 4th wall layer ground floor level
     1423       0.35_wp,        &  !< parameter 29  - thermal conductivity 1st/2nd wall layer ground floor level
     1424       0.14_wp,        &  !< parameter 30  - thermal conductivity 3rd wall layer ground floor level
     1425       0.035_wp,       &  !< parameter 31  - thermal conductivity 4th wall layer ground floor level
     1426       0.92_wp,        &  !< parameter 32  - wall emissivity ground floor level
     1427       0.8_wp,         &  !< parameter 33  - window emissivity ground floor level
     1428       0.86_wp,        &  !< parameter 34  - green emissivity ground floor level
     1429       0.6_wp,         &  !< parameter 35  - window transmissivity ground floor level
     1430       0.01_wp,        &  !< parameter 36  - z0 roughness ground floor level
     1431       0.001_wp,       &  !< parameter 37  - z0h/z0q roughness heat/humidity
     1432       27.0_wp,        &  !< parameter 38  - wall albedo above ground floor level
     1433       5.0_wp,         &  !< parameter 39  - green albedo above ground floor level
     1434       27.0_wp,        &  !< parameter 40  - window albedo above ground floor level
     1435       0.005_wp,       &  !< parameter 41  - 1st wall layer thickness above ground floor level
     1436       0.01_wp,        &  !< parameter 42  - 2nd wall layer thickness above ground floor level
     1437       0.41_wp,        &  !< parameter 43  - 3rd wall layer thickness above ground floor level
     1438       0.7_wp,         &  !< parameter 44  - 4th wall layer thickness above ground floor level
     1439       20000.0_wp,     &  !< parameter 45  - heat capacity wall surface
     1440       23.0_wp,        &  !< parameter 46  - thermal conductivity of wall surface
     1441       20000.0_wp,     &  !< parameter 47  - heat capacity of window surface
     1442       20000.0_wp,     &  !< parameter 48  - heat capacity of green surface
     1443       23.0_wp,        &  !< parameter 49  - thermal conductivity of window surface
     1444       10.0_wp,        &  !< parameter 50  - thermal conductivty of green surface
     1445       1.0_wp,         &  !< parameter 51  - wall fraction ground plate
     1446       0.005_wp,       &  !< parameter 52  - 1st wall layer thickness ground plate
     1447       0.01_wp,        &  !< parameter 53  - 2nd wall layer thickness ground plate
     1448       0.41_wp,        &  !< parameter 54  - 3rd wall layer thickness ground plate
     1449       0.7_wp,         &  !< parameter 55  - 4th wall layer thickness ground plate
     1450       2000000.0_wp,   &  !< parameter 56  - heat capacity 1st/2nd wall layer ground plate
     1451       103000.0_wp,    &  !< parameter 57  - heat capacity 3rd wall layer ground plate
     1452       900000.0_wp,    &  !< parameter 58  - heat capacity 4th wall layer ground plate
     1453       0.35_wp,        &  !< parameter 59  - thermal conductivity 1st/2nd wall layer ground plate
     1454       0.14_wp,        &  !< parameter 60  - thermal conductivity 3rd wall layer ground plate
     1455       0.035_wp,       &  !< parameter 61  - thermal conductivity 4th wall layer ground plate
     1456       0.005_wp,       &  !< parameter 62  - 1st wall layer thickness ground floor level
     1457       0.01_wp,        &  !< parameter 63  - 2nd wall layer thickness ground floor level
     1458       0.41_wp,        &  !< parameter 64  - 3rd wall layer thickness ground floor level
     1459       0.7_wp,         &  !< parameter 65  - 4th wall layer thickness ground floor level
     1460       27.0_wp,        &  !< parameter 66  - wall albedo ground floor level
     1461       0.003_wp,       &  !< parameter 67  - 1st window layer thickness ground floor level
     1462       0.006_wp,       &  !< parameter 68  - 2nd window layer thickness ground floor level
     1463       0.012_wp,       &  !< parameter 69  - 3rd window layer thickness ground floor level
     1464       0.018_wp,       &  !< parameter 70  - 4th window layer thickness ground floor level
     1465       1736000.0_wp,   &  !< parameter 71  - heat capacity 1st/2nd window layer ground floor level
     1466       1736000.0_wp,   &  !< parameter 72  - heat capacity 3rd window layer ground floor level
     1467       1736000.0_wp,   &  !< parameter 73  - heat capacity 4th window layer ground floor level
     1468       0.037_wp,       &  !< parameter 74  - thermal conductivity 1st/2nd window layer ground floor level
     1469       0.037_wp,       &  !< parameter 75  - thermal conductivity 3rd window layer ground floor level
     1470       0.037_wp,       &  !< parameter 76  - thermal conductivity 4th window layer ground floor level
     1471       27.0_wp,        &  !< parameter 77  - window albedo ground floor level
     1472       5.0_wp,         &  !< parameter 78  - green albedo ground floor level
     1473       0.003_wp,       &  !< parameter 79  - 1st window layer thickness above ground floor level
     1474       0.006_wp,       &  !< parameter 80  - 2nd thickness window layer above ground floor level
     1475       0.012_wp,       &  !< parameter 81  - 3rd window layer thickness above ground floor level
     1476       0.018_wp,       &  !< parameter 82  - 4th window layer thickness above ground floor level
     1477       1736000.0_wp,   &  !< parameter 83  - heat capacity 1st/2nd window layer above ground floor level
     1478       1736000.0_wp,   &  !< parameter 84  - heat capacity 3rd window layer above ground floor level
     1479       1736000.0_wp,   &  !< parameter 85  - heat capacity 4th window layer above ground floor level
     1480       0.037_wp,       &  !< parameter 86  - thermal conductivity 1st/2nd window layer above ground floor level
     1481       0.037_wp,       &  !< parameter 87  - thermal conductivity 3rd window layer above ground floor level
     1482       0.037_wp,       &  !< parameter 88  - thermal conductivity 4th window layer above ground floor level
     1483       1.0_wp,         &  !< parameter 89  - wall fraction roof
     1484       0.005_wp,       &  !< parameter 90  - 1st wall layer thickness roof
     1485       0.01_wp,        &  !< parameter 91  - 2nd wall layer thickness roof
     1486       0.41_wp,        &  !< parameter 92  - 3rd wall layer thickness roof
     1487       0.7_wp,         &  !< parameter 93  - 4th wall layer thickness roof
     1488       2000000.0_wp,   &  !< parameter 94  - heat capacity 1st/2nd wall layer roof
     1489       103000.0_wp,    &  !< parameter 95  - heat capacity 3rd wall layer roof
     1490       900000.0_wp,    &  !< parameter 96  - heat capacity 4th wall layer roof
     1491       0.35_wp,        &  !< parameter 97  - thermal conductivity 1st/2nd wall layer roof
     1492       0.14_wp,        &  !< parameter 98  - thermal conductivity 3rd wall layer roof
     1493       0.035_wp,       &  !< parameter 99  - thermal conductivity 4th wall layer roof
     1494       0.91_wp,        &  !< parameter 100 - wall emissivity roof
     1495       27.0_wp,        &  !< parameter 101 - wall albedo roof
     1496       0.0_wp,         &  !< parameter 102 - window fraction roof
     1497       0.003_wp,       &  !< parameter 103 - window 1st layer thickness roof
     1498       0.006_wp,       &  !< parameter 104 - window 2nd layer thickness roof
     1499       0.012_wp,       &  !< parameter 105 - window 3rd layer thickness roof
     1500       0.018_wp,       &  !< parameter 106 - window 4th layer thickness roof
     1501       1736000.0_wp,   &  !< parameter 107 - heat capacity 1st/2nd window layer roof
     1502       1736000.0_wp,   &  !< parameter 108 - heat capacity 3rd window layer roof
     1503       1736000.0_wp,   &  !< parameter 109 - heat capacity 4th window layer roof
     1504       0.037_wp,       &  !< parameter 110 - thermal conductivity 1st/2nd window layer roof
     1505       0.037_wp,       &  !< parameter 111 - thermal conductivity 3rd window layer roof
     1506       0.037_wp,       &  !< parameter 112 - thermal conductivity 4th window layer roof
     1507       0.8_wp,         &  !< parameter 113 - window emissivity roof
     1508       0.6_wp,         &  !< parameter 114 - window transmissivity roof
     1509       27.0_wp,        &  !< parameter 115 - window albedo roof
     1510       0.86_wp,        &  !< parameter 116 - green emissivity roof
     1511       5.0_wp,         &  !< parameter 117 - green albedo roof
     1512       0.0_wp,         &  !< parameter 118 - green type roof
     1513       0.8_wp,         &  !< parameter 119 - shading factor
     1514       0.5_wp,         &  !< parameter 120 - g-value windows
     1515       0.6_wp,         &  !< parameter 121 - u-value windows
     1516       0.1_wp,         &  !< parameter 122 - basical airflow without occupancy of the room
     1517       1.5_wp,         &  !< parameter 123 - additional airflow depend of occupancy of the room
     1518       0.9_wp,         &  !< parameter 124 - heat recovery efficiency
     1519       2.5_wp,         &  !< parameter 125 - dynamic parameter specific effective surface
     1520       80000.0_wp,     &  !< parameter 126 - dynamic parameter innner heatstorage
     1521       4.5_wp,         &  !< parameter 127 - ratio internal surface/floor area
     1522       100000.0_wp,    &  !< parameter 128 - maximal heating capacity
     1523       0.0_wp,         &  !< parameter 129 - maximal cooling capacity
     1524       5.0_wp,         &  !< parameter 130 - additional internal heat gains dependent on occupancy of the room
     1525       15.0_wp,        &  !< parameter 131 - basic internal heat gains without occupancy of the room
     1526       3.0_wp,         &  !< parameter 132 - storey height
     1527       0.2_wp,         &  !< parameter 133 - ceiling construction height
     1528       1.0_wp,         &  !< parameter 0   - wall fraction above ground floor level
     1529       0.0_wp,         &  !< parameter 1   - window fraction above ground floor level
     1530       0.0_wp,         &  !< parameter 2   - green fraction above ground floor level
     1531       0.0_wp,         &  !< parameter 3   - green fraction roof above ground floor level
     1532       1.5_wp,         &  !< parameter 4   - LAI roof
     1533       1.5_wp,         &  !< parameter 5   - LAI on wall above ground floor level
     1534       1950400.0_wp,   &  !< parameter 6   - heat capacity 1st/2nd wall layer above ground floor level
     1535       1848000.0_wp,   &  !< parameter 7   - heat capacity 3rd wall layer above ground floor level
     1536       1848000.0_wp,   &  !< parameter 8   - heat capacity 4th wall layer above ground floor level
     1537       0.7_wp,         &  !< parameter 9   - thermal conductivity 1st/2nd wall layer above ground floor level
     1538       1.0_wp,         &  !< parameter 10  - thermal conductivity 3rd wall layer above ground floor level
     1539       1.0_wp,         &  !< parameter 11  - thermal conductivity 4th wall layer above ground floor level
     1540       299.15_wp,      &  !< parameter 12  - indoor target summer temperature
     1541       293.15_wp,      &  !< parameter 13  - indoor target winter temperature
     1542       0.9_wp,         &  !< parameter 14  - wall emissivity above ground floor level
     1543       0.86_wp,        &  !< parameter 15  - green emissivity above ground floor level
     1544       0.8_wp,         &  !< parameter 16  - window emissivity above ground floor level
     1545       0.6_wp,         &  !< parameter 17  - window transmissivity above ground floor level
     1546       0.001_wp,       &  !< parameter 18  - z0 roughness above ground floor level
     1547       0.0001_wp,      &  !< parameter 19  - z0h/z0g roughness heat/humidity above ground floor level
     1548       4.0_wp,         &  !< parameter 20  - ground floor level height
     1549       1.0_wp,         &  !< parameter 21  - wall fraction ground floor level
     1550       0.0_wp,         &  !< parameter 22  - window fraction ground floor level
     1551       0.0_wp,         &  !< parameter 23  - green fraction ground floor level
     1552       0.0_wp,         &  !< parameter 24  - green fraction roof ground floor level
     1553       1.5_wp,         &  !< parameter 25  - LAI on wall ground floor level
     1554       1950400.0_wp,   &  !< parameter 26  - heat capacity 1st/2nd wall layer ground floor level
     1555       1848000.0_wp,   &  !< parameter 27  - heat capacity 3rd wall layer ground floor level
     1556       1848000.0_wp,   &  !< parameter 28  - heat capacity 4th wall layer ground floor level
     1557       0.7_wp,         &  !< parameter 29  - thermal conductivity 1st/2nd wall layer ground floor level
     1558       1.0_wp,         &  !< parameter 30  - thermal conductivity 3rd wall layer ground floor level
     1559       1.0_wp,         &  !< parameter 31  - thermal conductivity 4th wall layer ground floor level
     1560       0.9_wp,         &  !< parameter 32  - wall emissivity ground floor level
     1561       0.8_wp,         &  !< parameter 33  - window emissivity ground floor level
     1562       0.86_wp,        &  !< parameter 34  - green emissivity ground floor level
     1563       0.6_wp,         &  !< parameter 35  - window transmissivity ground floor level
     1564       0.01_wp,        &  !< parameter 36  - z0 roughness ground floor level
     1565       0.001_wp,       &  !< parameter 37  - z0h/z0q roughness heat/humidity
     1566       27.0_wp,        &  !< parameter 38  - wall albedo above ground floor level
     1567       5.0_wp,         &  !< parameter 39  - green albedo above ground floor level
     1568       27.0_wp,        &  !< parameter 40  - window albedo above ground floor level
     1569       0.29_wp,        &  !< parameter 41  - 1st wall layer thickness above ground floor level
     1570       0.295_wp,       &  !< parameter 42  - 2nd wall layer thickness above ground floor level
     1571       0.695_wp,       &  !< parameter 43  - 3rd wall layer thickness above ground floor level
     1572       0.985_wp,       &  !< parameter 44  - 4th wall layer thickness above ground floor level
     1573       20000.0_wp,     &  !< parameter 45  - heat capacity wall surface
     1574       23.0_wp,        &  !< parameter 46  - thermal conductivity of wall surface
     1575       20000.0_wp,     &  !< parameter 47  - heat capacity of window surface
     1576       20000.0_wp,     &  !< parameter 48  - heat capacity of green surface
     1577       23.0_wp,        &  !< parameter 49  - thermal conductivity of window surface
     1578       10.0_wp,        &  !< parameter 50  - thermal conductivty of green surface
     1579       1.0_wp,         &  !< parameter 51  - wall fraction ground plate
     1580       0.29_wp,        &  !< parameter 52  - 1st wall layer thickness ground plate
     1581       0.295_wp,       &  !< parameter 53  - 2nd wall layer thickness ground plate
     1582       0.695_wp,       &  !< parameter 54  - 3rd wall layer thickness ground plate
     1583       0.985_wp,       &  !< parameter 55  - 4th wall layer thickness ground plate
     1584       1950400.0_wp,   &  !< parameter 56  - heat capacity 1st/2nd wall layer ground plate
     1585       1848000.0_wp,   &  !< parameter 57  - heat capacity 3rd wall layer ground plate
     1586       1848000.0_wp,   &  !< parameter 58  - heat capacity 4th wall layer ground plate
     1587       0.7_wp,         &  !< parameter 59  - thermal conductivity 1st/2nd wall layer ground plate
     1588       1.0_wp,         &  !< parameter 60  - thermal conductivity 3rd wall layer ground plate
     1589       1.0_wp,         &  !< parameter 61  - thermal conductivity 4th wall layer ground plate
     1590       0.29_wp,        &  !< parameter 62  - 1st wall layer thickness ground floor level
     1591       0.295_wp,       &  !< parameter 63  - 2nd wall layer thickness ground floor level
     1592       0.695_wp,       &  !< parameter 64  - 3rd wall layer thickness ground floor level
     1593       0.985_wp,       &  !< parameter 65  - 4th wall layer thickness ground floor level
     1594       27.0_wp,        &  !< parameter 66  - wall albedo ground floor level
     1595       0.003_wp,       &  !< parameter 67  - 1st window layer thickness ground floor level
     1596       0.006_wp,       &  !< parameter 68  - 2nd window layer thickness ground floor level
     1597       0.012_wp,       &  !< parameter 69  - 3rd window layer thickness ground floor level
     1598       0.018_wp,       &  !< parameter 70  - 4th window layer thickness ground floor level
     1599       1736000.0_wp,   &  !< parameter 71  - heat capacity 1st/2nd window layer ground floor level
     1600       1736000.0_wp,   &  !< parameter 72  - heat capacity 3rd window layer ground floor level
     1601       1736000.0_wp,   &  !< parameter 73  - heat capacity 4th window layer ground floor level
     1602       0.57_wp,        &  !< parameter 74  - thermal conductivity 1st/2nd window layer ground floor level
     1603       0.57_wp,        &  !< parameter 75  - thermal conductivity 3rd window layer ground floor level
     1604       0.57_wp,        &  !< parameter 76  - thermal conductivity 4th window layer ground floor level
     1605       27.0_wp,        &  !< parameter 77  - window albedo ground floor level
     1606       5.0_wp,         &  !< parameter 78  - green albedo ground floor level
     1607       0.003_wp,       &  !< parameter 79  - 1st window layer thickness above ground floor level
     1608       0.006_wp,       &  !< parameter 80  - 2nd thickness window layer above ground floor level
     1609       0.012_wp,       &  !< parameter 81  - 3rd window layer thickness above ground floor level
     1610       0.018_wp,       &  !< parameter 82  - 4th window layer thickness above ground floor level
     1611       1736000.0_wp,   &  !< parameter 83  - heat capacity 1st/2nd window layer above ground floor level
     1612       1736000.0_wp,   &  !< parameter 84  - heat capacity 3rd window layer above ground floor level
     1613       1736000.0_wp,   &  !< parameter 85  - heat capacity 4th window layer above ground floor level
     1614       0.57_wp,        &  !< parameter 86  - thermal conductivity 1st/2nd window layer above ground floor level
     1615       0.57_wp,        &  !< parameter 87  - thermal conductivity 3rd window layer above ground floor level
     1616       0.57_wp,        &  !< parameter 88  - thermal conductivity 4th window layer above ground floor level
     1617       1.0_wp,         &  !< parameter 89  - wall fraction roof
     1618       0.29_wp,        &  !< parameter 90  - 1st wall layer thickness roof
     1619       0.295_wp,       &  !< parameter 91  - 2nd wall layer thickness roof
     1620       0.695_wp,       &  !< parameter 92  - 3rd wall layer thickness roof
     1621       0.985_wp,       &  !< parameter 93  - 4th wall layer thickness roof
     1622       1950400.0_wp,   &  !< parameter 94  - heat capacity 1st/2nd wall layer roof
     1623       1848000.0_wp,   &  !< parameter 95  - heat capacity 3rd wall layer roof
     1624       1848000.0_wp,   &  !< parameter 96  - heat capacity 4th wall layer roof
     1625       0.7_wp,         &  !< parameter 97  - thermal conductivity 1st/2nd wall layer roof
     1626       1.0_wp,         &  !< parameter 98  - thermal conductivity 3rd wall layer roof
     1627       1.0_wp,         &  !< parameter 99  - thermal conductivity 4th wall layer roof
     1628       0.9_wp,         &  !< parameter 100 - wall emissivity roof
     1629       27.0_wp,        &  !< parameter 101 - wall albedo roof
     1630       0.0_wp,         &  !< parameter 102 - window fraction roof
     1631       0.003_wp,       &  !< parameter 103 - window 1st layer thickness roof
     1632       0.006_wp,       &  !< parameter 104 - window 2nd layer thickness roof
     1633       0.012_wp,       &  !< parameter 105 - window 3rd layer thickness roof
     1634       0.018_wp,       &  !< parameter 106 - window 4th layer thickness roof
     1635       1736000.0_wp,   &  !< parameter 107 - heat capacity 1st/2nd window layer roof
     1636       1736000.0_wp,   &  !< parameter 108 - heat capacity 3rd window layer roof
     1637       1736000.0_wp,   &  !< parameter 109 - heat capacity 4th window layer roof
     1638       0.57_wp,        &  !< parameter 110 - thermal conductivity 1st/2nd window layer roof
     1639       0.57_wp,        &  !< parameter 111 - thermal conductivity 3rd window layer roof
     1640       0.57_wp,        &  !< parameter 112 - thermal conductivity 4th window layer roof
     1641       0.8_wp,         &  !< parameter 113 - window emissivity roof
     1642       0.6_wp,         &  !< parameter 114 - window transmissivity roof
     1643       27.0_wp,        &  !< parameter 115 - window albedo roof
     1644       0.86_wp,        &  !< parameter 116 - green emissivity roof
     1645       5.0_wp,         &  !< parameter 117 - green albedo roof
     1646       0.0_wp,         &  !< parameter 118 - green type roof
     1647       0.8_wp,         &  !< parameter 119 - shading factor
     1648       100.0_wp,       &  !< parameter 120 - g-value windows
     1649       100.0_wp,       &  !< parameter 121 - u-value windows
     1650       20.0_wp,        &  !< parameter 122 - basical airflow without occupancy of the room
     1651       20.0_wp,        &  !< parameter 123 - additional airflow depend of occupancy of the room
     1652       0.0_wp,         &  !< parameter 124 - heat recovery efficiency
     1653       1.0_wp,         &  !< parameter 125 - dynamic parameter specific effective surface
     1654       1.0_wp,         &  !< parameter 126 - dynamic parameter innner heatstorage
     1655       4.5_wp,         &  !< parameter 127 - ratio internal surface/floor area
     1656       100000.0_wp,    &  !< parameter 128 - maximal heating capacity
     1657       0.0_wp,         &  !< parameter 129 - maximal cooling capacity
     1658       0.0_wp,         &  !< parameter 130 - additional internal heat gains dependent on occupancy of the room
     1659       0.0_wp,         &  !< parameter 131 - basic internal heat gains without occupancy of the room
     1660       3.0_wp,         &  !< parameter 132 - storey height
     1661       0.2_wp          &  !< parameter 133 - ceiling construction height
     1662       /),             &
     1663       (/134, 7/) )
    9811664
    9821665!
     
    26533336        CHARACTER(LEN=2)                              :: ls
    26543337        CHARACTER(LEN=varnamelength)                  :: var           !< TRIM(variable)
    2655         INTEGER(iwp), PARAMETER                       :: nl1 = 16      !< number of directional usm variables
     3338        INTEGER(iwp), PARAMETER                       :: nl1 = 14      !< number of directional usm variables
    26563339        CHARACTER(LEN=varnamelength), DIMENSION(nl1)  :: varlist1 = &  !< list of directional usm variables
    26573340                  (/'usm_wshf                      ', &
     
    26643347                    'usm_surfwintrans              ', &
    26653348                    'usm_surfcat                   ', &
    2666                     'usm_surfalb                   ', &
    2667                     'usm_surfemis                  ', &
    26683349                    'usm_t_surf_wall               ', &
    26693350                    'usm_t_surf_window             ', &
     
    27343415                  var(1:14) == 'usm_theta_10cm' )  THEN
    27353416            unit = 'K'
    2736         ELSE IF ( var(1:9) == 'usm_surfz'  .OR.  var(1:11) == 'usm_surfcat'  .OR.         &
    2737                   var(1:11) == 'usm_surfalb'  .OR.  var(1:12) == 'usm_surfemis'  )  THEN
     3417        ELSE IF ( var(1:9) == 'usm_surfz'  .OR.  var(1:11) == 'usm_surfcat' )  THEN
    27383418            unit = '1'
    27393419        ELSE
     
    28513531        LOGICAL, INTENT(OUT)           ::  found     !<
    28523532        REAL(sp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do) ::  local_pf   !< sp - it has to correspond to module data_output_3d
    2853         REAL(wp), DIMENSION(nzb:nzt+1,nys:nyn,nxl:nxr)     ::  temp_pf    !< temp array for urban surface output procedure
     3533        REAL(sp), DIMENSION(nzb:nzt+1,nys:nyn,nxl:nxr)     ::  temp_pf    !< temp array for urban surface output procedure
    28543534       
    28553535        CHARACTER (len=varnamelength)                      :: var     !< trimmed variable name
     
    29563636              ENDIF
    29573637             
    2958           CASE ( 'usm_surfalb' )
    2959 !
    2960 !--           surface albedo, weighted average
    2961               IF ( idsint == iup_u )  THEN
    2962                  DO  m = 1, surf_usm_h%ns
    2963                     i = surf_usm_h%i(m)
    2964                     j = surf_usm_h%j(m)
    2965                     k = surf_usm_h%k(m)
    2966                     temp_pf(k,j,i) = surf_usm_h%frac(ind_veg_wall,m)     *     &
    2967                                      surf_usm_h%albedo(ind_veg_wall,m)  +      &
    2968                                      surf_usm_h%frac(ind_pav_green,m)    *     &
    2969                                      surf_usm_h%albedo(ind_pav_green,m) +      &
    2970                                      surf_usm_h%frac(ind_wat_win,m)      *     &
    2971                                      surf_usm_h%albedo(ind_wat_win,m)
    2972                  ENDDO
    2973               ELSE
    2974                  l = idsidx
    2975                  DO  m = 1, surf_usm_v(l)%ns
    2976                     i = surf_usm_v(l)%i(m)
    2977                     j = surf_usm_v(l)%j(m)
    2978                     k = surf_usm_v(l)%k(m)
    2979                     temp_pf(k,j,i) = surf_usm_v(l)%frac(ind_veg_wall,m)     *  &
    2980                                      surf_usm_v(l)%albedo(ind_veg_wall,m)  +   &
    2981                                      surf_usm_v(l)%frac(ind_pav_green,m)    *  &
    2982                                      surf_usm_v(l)%albedo(ind_pav_green,m) +   &
    2983                                      surf_usm_v(l)%frac(ind_wat_win,m)      *  &
    2984                                      surf_usm_v(l)%albedo(ind_wat_win,m)
    2985                  ENDDO
    2986               ENDIF
    2987              
    2988           CASE ( 'usm_surfemis' )
    2989 !
    2990 !--           surface emissivity, weighted average
    2991               IF ( idsint == iup_u )  THEN
    2992                  DO  m = 1, surf_usm_h%ns
    2993                     i = surf_usm_h%i(m)
    2994                     j = surf_usm_h%j(m)
    2995                     k = surf_usm_h%k(m)
    2996                     temp_pf(k,j,i) =  surf_usm_h%frac(ind_veg_wall,m)      *   &
    2997                                       surf_usm_h%emissivity(ind_veg_wall,m)  + &
    2998                                       surf_usm_h%frac(ind_pav_green,m)     *   &
    2999                                       surf_usm_h%emissivity(ind_pav_green,m) + &
    3000                                       surf_usm_h%frac(ind_wat_win,m)       *   &
    3001                                       surf_usm_h%emissivity(ind_wat_win,m)
    3002                  ENDDO
    3003               ELSE
    3004                  l = idsidx
    3005                  DO  m = 1, surf_usm_v(l)%ns
    3006                     i = surf_usm_v(l)%i(m)
    3007                     j = surf_usm_v(l)%j(m)
    3008                     k = surf_usm_v(l)%k(m)
    3009                     temp_pf(k,j,i) = surf_usm_v(l)%frac(ind_veg_wall,m)       *&
    3010                                      surf_usm_v(l)%emissivity(ind_veg_wall,m) +&
    3011                                      surf_usm_v(l)%frac(ind_pav_green,m)      *&
    3012                                      surf_usm_v(l)%emissivity(ind_pav_green,m)+&
    3013                                      surf_usm_v(l)%frac(ind_wat_win,m)        *&
    3014                                      surf_usm_v(l)%emissivity(ind_wat_win,m)
    3015                  ENDDO
    3016               ENDIF
    3017 
    30183638          CASE ( 'usm_surfwintrans' )
    30193639!
     
    37594379             var(1:15) == 'usm_theta_10cm' .OR.                                             &
    37604380             var(1:9) == 'usm_surfz'  .OR.  var(1:11) == 'usm_surfcat'  .OR.                &
    3761              var(1:11) == 'usm_surfalb'  .OR.  var(1:12) == 'usm_surfemis'  .OR.            &
    37624381             var(1:16) == 'usm_surfwintrans'  .OR. var(1:7) == 'usm_swc' ) THEN
    37634382
     
    40584677!--     ground-floor-level properties.
    40594678        surf_usm_h%ground_level = .FALSE.
    4060 !         DO  m = 1, surf_usm_h%ns
    4061 !            i = surf_usm_h%i(m)
    4062 !            j = surf_usm_h%j(m)
    4063 !            k = surf_usm_h%k(m)
    4064 ! !
    4065 ! !--        Get local ground level. If no ground level is given in input file,
    4066 ! !--        use default value.
    4067 !            ground_floor_level_l = ground_floor_level
    4068 !            IF ( building_pars_f%from_file )  THEN
    4069 !               IF ( building_pars_f%pars_xy(ind_gflh,j,i) /=                    &
    4070 !                    building_pars_f%fill )  &
    4071 !                  ground_floor_level_l = building_pars_f%pars_xy(ind_gflh,j,i)         
    4072 !            ENDIF
    4073 ! !
    4074 ! !--        Determine height of surface element above ground level
    4075 !            IF (  terrain_height_f%from_file )  THEN
    4076 !               z_agl = zw(k) - terrain_height_f%var(j,i)
    4077 !            ELSE
    4078 !               z_agl = zw(k)
    4079 !            ENDIF
    4080 ! !
    4081 ! !--        Set flag for ground level
    4082 !            IF ( z_agl <= ground_floor_level_l )                                &
    4083 !               surf_usm_h%ground_level(m) = .TRUE.
    4084 !         ENDDO
    40854679
    40864680        DO  l = 0, 3
     
    40914685              k = surf_usm_v(l)%k(m)
    40924686!
    4093 !--           Get local ground level. If no ground level is given in input file,
    4094 !--           use default value.
     4687!--           Determine local ground level. Level 1 - default value,
     4688!--           level 2 - initialization according to building type,
     4689!--           level 3 - initialization from value read from file.
    40954690              ground_floor_level_l = ground_floor_level
     4691             
     4692              IF ( building_type_f%from_file )  THEN
     4693                  ground_floor_level_l =                                       &
     4694                              building_pars(ind_gflh,building_type_f%var(j,i))
     4695              ENDIF
     4696             
    40964697              IF ( building_pars_f%from_file )  THEN
    40974698                 IF ( building_pars_f%pars_xy(ind_gflh,j,i) /=                 &
    4098                       building_pars_f%fill ) &
     4699                      building_pars_f%fill )                                   &
    40994700                    ground_floor_level_l = building_pars_f%pars_xy(ind_gflh,j,i)
    41004701              ENDIF
     
    41024703!--           Determine height of surface element above ground level. Please
    41034704!--           note, height of surface element is determined with respect to
    4104 !--           its height of the adjoing atmospheric grid point.
    4105               IF (  terrain_height_f%from_file )  THEN
     4705!--           its height above ground of the reference grid point in atmosphere,
     4706!--           Therefore, substract the offset values when assessing the terrain
     4707!--           height.
     4708              IF ( terrain_height_f%from_file )  THEN
    41064709                 z_agl = zw(k) - terrain_height_f%var(j-surf_usm_v(l)%joff,    &
    41074710                                                      i-surf_usm_v(l)%ioff)
     
    41494752!
    41504753!--     Initialize urban-type surface attribute. According to initialization in
    4151 !--     land-surface model, follow a 3-level approach.
     4754!--     land-surface model, follow a 3-level approach. 
    41524755!--     Level 1 - initialization via default attributes
    41534756        DO  m = 1, surf_usm_h%ns
     
    42044807           surf_usm_h%lambda_h_window(nzb_wall+3,m) = building_pars(ind_tc3_win_r,building_type)   
    42054808
    4206            surf_usm_h%target_temp_summer(m)  = building_pars(117,building_type)   
    4207            surf_usm_h%target_temp_winter(m)  = building_pars(118,building_type)   
     4809           surf_usm_h%target_temp_summer(m)  = building_pars(ind_indoor_target_temp_summer,building_type)   
     4810           surf_usm_h%target_temp_winter(m)  = building_pars(ind_indoor_target_temp_winter,building_type)   
    42084811!
    42094812!--        emissivity of wall-, green- and window fraction
     
    42384841           surf_usm_h%zw_window(nzb_wall+3,m)       = building_pars(ind_thick_4_win_r,building_type)
    42394842
    4240            surf_usm_h%c_surface(m)           = building_pars(0,building_type) 
    4241            surf_usm_h%lambda_surf(m)         = building_pars(3,building_type) 
    4242            surf_usm_h%c_surface_green(m)     = building_pars(2,building_type) 
    4243            surf_usm_h%lambda_surf_green(m)   = building_pars(5,building_type) 
    4244            surf_usm_h%c_surface_window(m)    = building_pars(1,building_type) 
    4245            surf_usm_h%lambda_surf_window(m)  = building_pars(4,building_type) 
     4843           surf_usm_h%c_surface(m)           = building_pars(ind_c_surface,building_type) 
     4844           surf_usm_h%lambda_surf(m)         = building_pars(ind_lambda_surf,building_type) 
     4845           surf_usm_h%c_surface_green(m)     = building_pars(ind_c_surface_green,building_type) 
     4846           surf_usm_h%lambda_surf_green(m)   = building_pars(ind_lambda_surf_green,building_type) 
     4847           surf_usm_h%c_surface_window(m)    = building_pars(ind_c_surface_win,building_type) 
     4848           surf_usm_h%lambda_surf_window(m)  = building_pars(ind_lambda_surf_win,building_type) 
    42464849           
    42474850           surf_usm_h%green_type_roof(m)     = building_pars(ind_green_type_roof,building_type)
     
    43674970              surf_usm_v(l)%lambda_h_window(nzb_wall+3,m) = building_pars(ind_tc3_win,building_type)   
    43684971
    4369               surf_usm_v(l)%target_temp_summer(m)  = building_pars(117,building_type)   
    4370               surf_usm_v(l)%target_temp_winter(m)  = building_pars(118,building_type)   
     4972              surf_usm_v(l)%target_temp_summer(m)  = building_pars(ind_indoor_target_temp_summer,building_type)   
     4973              surf_usm_v(l)%target_temp_winter(m)  = building_pars(ind_indoor_target_temp_winter,building_type)   
    43714974!
    43724975!--           emissivity of wall-, green- and window fraction
     
    44005003              surf_usm_v(l)%zw_window(nzb_wall+3,m)       = building_pars(ind_thick_4_win,building_type)
    44015004
    4402               surf_usm_v(l)%c_surface(m)           = building_pars(0,building_type) 
    4403               surf_usm_v(l)%lambda_surf(m)         = building_pars(3,building_type)
    4404               surf_usm_v(l)%c_surface_green(m)     = building_pars(2,building_type) 
    4405               surf_usm_v(l)%lambda_surf_green(m)   = building_pars(5,building_type)
    4406               surf_usm_v(l)%c_surface_window(m)    = building_pars(1,building_type) 
    4407               surf_usm_v(l)%lambda_surf_window(m)  = building_pars(4,building_type)
     5005              surf_usm_v(l)%c_surface(m)           = building_pars(ind_c_surface,building_type) 
     5006              surf_usm_v(l)%lambda_surf(m)         = building_pars(ind_lambda_surf,building_type)
     5007              surf_usm_v(l)%c_surface_green(m)     = building_pars(ind_c_surface_green,building_type) 
     5008              surf_usm_v(l)%lambda_surf_green(m)   = building_pars(ind_lambda_surf_green,building_type)
     5009              surf_usm_v(l)%c_surface_window(m)    = building_pars(ind_c_surface_win,building_type) 
     5010              surf_usm_v(l)%lambda_surf_window(m)  = building_pars(ind_lambda_surf_win,building_type)
    44085011
    44095012           ENDDO
     
    44705073                 surf_usm_h%lambda_h_window(nzb_wall+3,m) = building_pars(ind_tc3_win_r,st)   
    44715074
    4472                  surf_usm_h%target_temp_summer(m)  = building_pars(117,st)   
    4473                  surf_usm_h%target_temp_winter(m)  = building_pars(118,st)   
     5075                 surf_usm_h%target_temp_summer(m)  = building_pars(ind_indoor_target_temp_summer,st)   
     5076                 surf_usm_h%target_temp_winter(m)  = building_pars(ind_indoor_target_temp_winter,st)   
    44745077!
    44755078!--              emissivity of wall-, green- and window fraction
     
    45045107                 surf_usm_h%zw_window(nzb_wall+3,m)       = building_pars(ind_thick_4_win_r,st)
    45055108
    4506                  surf_usm_h%c_surface(m)           = building_pars(0,st) 
    4507                  surf_usm_h%lambda_surf(m)         = building_pars(3,st)
    4508                  surf_usm_h%c_surface_green(m)     = building_pars(2,st) 
    4509                  surf_usm_h%lambda_surf_green(m)   = building_pars(5,st)
    4510                  surf_usm_h%c_surface_window(m)    = building_pars(1,st) 
    4511                  surf_usm_h%lambda_surf_window(m)  = building_pars(4,st)
     5109                 surf_usm_h%c_surface(m)           = building_pars(ind_c_surface,st) 
     5110                 surf_usm_h%lambda_surf(m)         = building_pars(ind_lambda_surf,st)
     5111                 surf_usm_h%c_surface_green(m)     = building_pars(ind_c_surface_green,st) 
     5112                 surf_usm_h%lambda_surf_green(m)   = building_pars(ind_lambda_surf_green,st)
     5113                 surf_usm_h%c_surface_window(m)    = building_pars(ind_c_surface_win,st) 
     5114                 surf_usm_h%lambda_surf_window(m)  = building_pars(ind_lambda_surf_win,st)
    45125115                 
    45135116                 surf_usm_h%green_type_roof(m)     = building_pars(ind_green_type_roof,st)
     
    46405243                    surf_usm_v(l)%lambda_h_window(nzb_wall+3,m) = building_pars(ind_tc3_win,st)
    46415244
    4642                     surf_usm_v(l)%target_temp_summer(m)  = building_pars(117,st)   
    4643                     surf_usm_v(l)%target_temp_winter(m)  = building_pars(118,st)   
     5245                    surf_usm_v(l)%target_temp_summer(m)  = building_pars(ind_indoor_target_temp_summer,st)   
     5246                    surf_usm_v(l)%target_temp_winter(m)  = building_pars(ind_indoor_target_temp_winter,st)   
    46445247!
    46455248!--                 emissivity of wall-, green- and window fraction
     
    46735276                    surf_usm_v(l)%zw_window(nzb_wall+3,m)       = building_pars(ind_thick_4_win,st)
    46745277
    4675                     surf_usm_v(l)%c_surface(m)           = building_pars(0,st) 
    4676                     surf_usm_v(l)%lambda_surf(m)         = building_pars(3,st)
    4677                     surf_usm_v(l)%c_surface_green(m)     = building_pars(2,st) 
    4678                     surf_usm_v(l)%lambda_surf_green(m)   = building_pars(5,st)
    4679                     surf_usm_v(l)%c_surface_window(m)    = building_pars(1,st) 
    4680                     surf_usm_v(l)%lambda_surf_window(m)  = building_pars(4,st)
     5278                    surf_usm_v(l)%c_surface(m)           = building_pars(ind_c_surface,st) 
     5279                    surf_usm_v(l)%lambda_surf(m)         = building_pars(ind_lambda_surf,st)
     5280                    surf_usm_v(l)%c_surface_green(m)     = building_pars(ind_c_surface_green,st) 
     5281                    surf_usm_v(l)%lambda_surf_green(m)   = building_pars(ind_lambda_surf_green,st)
     5282                    surf_usm_v(l)%c_surface_window(m)    = building_pars(ind_c_surface_win,st) 
     5283                    surf_usm_v(l)%lambda_surf_window(m)  = building_pars(ind_lambda_surf_win,st)
    46815284
    46825285
     
    46875290       
    46885291!
    4689 !--     Level 3 - initialization via building_pars read from file
     5292!--     Level 3 - initialization via building_pars read from file. Note, only
     5293!--     variables that are also defined in the input-standard can be initialized
     5294!--     via file. Other variables will be initialized on level 1 or 2.
    46905295        IF ( building_pars_f%from_file )  THEN
    46915296           DO  m = 1, surf_usm_h%ns
     
    46965301!--           In order to distinguish between ground floor level and
    46975302!--           above-ground-floor level surfaces, set input indices.
    4698               ind_green_frac_r = MERGE( ind_green_frac_r_gfl, ind_green_frac_r_agfl, &
     5303              ind_wall_frac    = MERGE( ind_wall_frac_gfl,                     &
     5304                                        ind_wall_frac_agfl,                    &
    46995305                                        surf_usm_h%ground_level(m) )
    4700               ind_lai_r        = MERGE( ind_lai_r_gfl,        ind_lai_r_agfl,        &
     5306              ind_green_frac_r = MERGE( ind_green_frac_r_gfl,                  &
     5307                                        ind_green_frac_r_agfl,                 &
    47015308                                        surf_usm_h%ground_level(m) )
    4702               ind_z0           = MERGE( ind_z0_gfl,           ind_z0_agfl,           &
     5309              ind_win_frac     = MERGE( ind_win_frac_gfl,                      &
     5310                                        ind_win_frac_agfl,                     &
    47035311                                        surf_usm_h%ground_level(m) )
    4704               ind_z0qh         = MERGE( ind_z0qh_gfl,         ind_z0qh_agfl,         &
     5312              ind_lai_r        = MERGE( ind_lai_r_gfl,                         &
     5313                                        ind_lai_r_agfl,                        &
    47055314                                        surf_usm_h%ground_level(m) )
     5315              ind_z0           = MERGE( ind_z0_gfl,                            &
     5316                                        ind_z0_agfl,                           &
     5317                                        surf_usm_h%ground_level(m) )
     5318              ind_z0qh         = MERGE( ind_z0qh_gfl,                          &
     5319                                        ind_z0qh_agfl,                         &
     5320                                        surf_usm_h%ground_level(m) )
     5321              ind_hc1          = MERGE( ind_hc1_gfl,                           &
     5322                                        ind_hc1_agfl,                          &
     5323                                        surf_usm_h%ground_level(m) )
     5324              ind_hc2          = MERGE( ind_hc2_gfl,                           &
     5325                                        ind_hc2_agfl,                          &
     5326                                        surf_usm_h%ground_level(m) )
     5327              ind_hc3          = MERGE( ind_hc3_gfl,                           &
     5328                                        ind_hc3_agfl,                          &
     5329                                        surf_usm_h%ground_level(m) )
     5330              ind_tc1          = MERGE( ind_tc1_gfl,                           &
     5331                                        ind_tc1_agfl,                          &
     5332                                        surf_usm_h%ground_level(m) )
     5333              ind_tc2          = MERGE( ind_tc2_gfl,                           &
     5334                                        ind_tc2_agfl,                          &
     5335                                        surf_usm_h%ground_level(m) )
     5336              ind_tc3          = MERGE( ind_tc3_gfl,                           &
     5337                                        ind_tc3_agfl,                          &
     5338                                        surf_usm_h%ground_level(m) )
     5339              ind_emis_wall    = MERGE( ind_emis_wall_gfl,                     &
     5340                                        ind_emis_wall_agfl,                    &
     5341                                        surf_usm_h%ground_level(m) )
     5342              ind_emis_green   = MERGE( ind_emis_green_gfl,                    &
     5343                                        ind_emis_green_agfl,                   &
     5344                                        surf_usm_h%ground_level(m) )
     5345              ind_emis_win     = MERGE( ind_emis_win_gfl,                      &
     5346                                        ind_emis_win_agfl,                     &
     5347                                        surf_usm_h%ground_level(m) )
     5348              ind_trans        = MERGE( ind_trans_gfl,                         &
     5349                                        ind_trans_agfl,                        &
     5350                                        surf_usm_h%ground_level(m) )
    47065351
    47075352!
    47085353!--           Initialize relatvie wall- (0), green- (1) and window (2) fractions
    4709               IF ( building_pars_f%pars_xy(ind_wall_frac_r,j,i) /= building_pars_f%fill )    &
    4710                  surf_usm_h%frac(ind_veg_wall,m)  = building_pars_f%pars_xy(ind_wall_frac_r,j,i)   
    4711               IF ( building_pars_f%pars_xy(ind_green_frac_r,j,i) /= building_pars_f%fill )   &
    4712                  surf_usm_h%frac(ind_pav_green,m) = building_pars_f%pars_xy(ind_green_frac_r,j,i)
    4713               IF ( building_pars_f%pars_xy(ind_win_frac_r,j,i) /= building_pars_f%fill )     &
    4714                  surf_usm_h%frac(ind_wat_win,m)   = building_pars_f%pars_xy(ind_win_frac_r,j,i)
    4715 
     5354              IF ( building_pars_f%pars_xy(ind_wall_frac,j,i) /=               &
     5355                   building_pars_f%fill )                                      &
     5356                 surf_usm_h%frac(ind_veg_wall,m)  =                            &
     5357                                    building_pars_f%pars_xy(ind_wall_frac,j,i)   
     5358                 
     5359              IF ( building_pars_f%pars_xy(ind_green_frac_r,j,i) /=            &         
     5360                   building_pars_f%fill )                                      &
     5361                 surf_usm_h%frac(ind_pav_green,m) =                            &
     5362                                    building_pars_f%pars_xy(ind_green_frac_r,j,i)
     5363                 
     5364              IF ( building_pars_f%pars_xy(ind_win_frac,j,i) /=                &
     5365                   building_pars_f%fill )                                      &
     5366                 surf_usm_h%frac(ind_wat_win,m)   =                            &
     5367                                    building_pars_f%pars_xy(ind_win_frac,j,i)
    47165368 
    4717               IF ( building_pars_f%pars_xy(ind_lai_r,j,i) /= building_pars_f%fill )          &
    4718                  surf_usm_h%lai(m)             = building_pars_f%pars_xy(ind_lai_r,j,i)
    4719 
    4720               IF ( building_pars_f%pars_xy(ind_hc1_wall_r,j,i) /= building_pars_f%fill )  THEN
    4721                  surf_usm_h%rho_c_wall(nzb_wall,m)   = building_pars_f%pars_xy(ind_hc1_wall_r,j,i) 
    4722                  surf_usm_h%rho_c_wall(nzb_wall+1,m) = building_pars_f%pars_xy(ind_hc1_wall_r,j,i)
     5369              IF ( building_pars_f%pars_xy(ind_lai_r,j,i) /=                   &
     5370                   building_pars_f%fill )                                      &
     5371                 surf_usm_h%lai(m)  = building_pars_f%pars_xy(ind_lai_r,j,i)
     5372
     5373              IF ( building_pars_f%pars_xy(ind_hc1,j,i) /=                     &
     5374                   building_pars_f%fill )  THEN
     5375                 surf_usm_h%rho_c_wall(nzb_wall,m)   =                         &
     5376                                    building_pars_f%pars_xy(ind_hc1,j,i) 
     5377                 surf_usm_h%rho_c_wall(nzb_wall+1,m) =                         &
     5378                                    building_pars_f%pars_xy(ind_hc1,j,i)
    47235379              ENDIF
    4724               IF ( building_pars_f%pars_xy(ind_hc2_wall_r,j,i) /= building_pars_f%fill )    &
    4725                  surf_usm_h%rho_c_wall(nzb_wall+2,m) = building_pars_f%pars_xy(ind_hc2_wall_r,j,i)
    4726               IF ( building_pars_f%pars_xy(ind_hc3_wall_r,j,i) /= building_pars_f%fill )    &
    4727                  surf_usm_h%rho_c_wall(nzb_wall+3,m) = building_pars_f%pars_xy(ind_hc3_wall_r,j,i)
    4728               IF ( building_pars_f%pars_xy(ind_hc1_wall_r,j,i) /= building_pars_f%fill )  THEN
    4729                  surf_usm_h%rho_c_green(nzb_wall,m)   = rho_c_soil !building_pars_f%pars_xy(ind_hc1_wall_r,j,i) 
    4730                  surf_usm_h%rho_c_green(nzb_wall+1,m) = rho_c_soil !building_pars_f%pars_xy(ind_hc1_wall_r,j,i)
     5380             
     5381             
     5382              IF ( building_pars_f%pars_xy(ind_hc2,j,i) /=                     &
     5383                   building_pars_f%fill )                                      &
     5384                 surf_usm_h%rho_c_wall(nzb_wall+2,m) =                         &
     5385                                    building_pars_f%pars_xy(ind_hc2,j,i)
     5386                 
     5387              IF ( building_pars_f%pars_xy(ind_hc3,j,i) /=                     &
     5388                   building_pars_f%fill )                                      &
     5389                 surf_usm_h%rho_c_wall(nzb_wall+3,m) =                         &
     5390                                    building_pars_f%pars_xy(ind_hc3,j,i)
     5391                 
     5392              IF ( building_pars_f%pars_xy(ind_hc1,j,i) /=                     &
     5393                   building_pars_f%fill )  THEN
     5394                 surf_usm_h%rho_c_green(nzb_wall,m)   =                        &
     5395                                    building_pars_f%pars_xy(ind_hc1,j,i) 
     5396                 surf_usm_h%rho_c_green(nzb_wall+1,m) =                        &
     5397                                    building_pars_f%pars_xy(ind_hc1,j,i)
    47315398              ENDIF
    4732               IF ( building_pars_f%pars_xy(ind_hc2_wall_r,j,i) /= building_pars_f%fill )    &
    4733                  surf_usm_h%rho_c_green(nzb_wall+2,m) = rho_c_soil !building_pars_f%pars_xy(ind_hc2_wall_r,j,i)
    4734               IF ( building_pars_f%pars_xy(ind_hc3_wall_r,j,i) /= building_pars_f%fill )    &
    4735                  surf_usm_h%rho_c_green(nzb_wall+3,m) = rho_c_soil !building_pars_f%pars_xy(ind_hc3_wall_r,j,i)
    4736               IF ( building_pars_f%pars_xy(ind_hc1_win_r,j,i) /= building_pars_f%fill )  THEN
    4737                  surf_usm_h%rho_c_window(nzb_wall,m)   = building_pars_f%pars_xy(ind_hc1_win_r,j,i) 
    4738                  surf_usm_h%rho_c_window(nzb_wall+1,m) = building_pars_f%pars_xy(ind_hc1_win_r,j,i)
     5399              IF ( building_pars_f%pars_xy(ind_hc2,j,i) /=                     &
     5400                   building_pars_f%fill )                                      &
     5401                 surf_usm_h%rho_c_green(nzb_wall+2,m) =                        &
     5402                                    building_pars_f%pars_xy(ind_hc2,j,i)
     5403                 
     5404              IF ( building_pars_f%pars_xy(ind_hc3,j,i) /=                     &
     5405                   building_pars_f%fill )                                      &
     5406                 surf_usm_h%rho_c_green(nzb_wall+3,m) =                        &
     5407                                    building_pars_f%pars_xy(ind_hc3,j,i)
     5408                 
     5409              IF ( building_pars_f%pars_xy(ind_hc1,j,i) /=                     &
     5410                   building_pars_f%fill )  THEN
     5411                 surf_usm_h%rho_c_window(nzb_wall,m)   =                       &
     5412                                    building_pars_f%pars_xy(ind_hc1,j,i) 
     5413                 surf_usm_h%rho_c_window(nzb_wall+1,m) =                       &
     5414                                    building_pars_f%pars_xy(ind_hc1,j,i)
    47395415              ENDIF
    4740               IF ( building_pars_f%pars_xy(ind_hc2_win_r,j,i) /= building_pars_f%fill )    &
    4741                  surf_usm_h%rho_c_window(nzb_wall+2,m) = building_pars_f%pars_xy(ind_hc2_win_r,j,i)
    4742               IF ( building_pars_f%pars_xy(ind_hc3_win_r,j,i) /= building_pars_f%fill )    &
    4743                  surf_usm_h%rho_c_window(nzb_wall+3,m) = building_pars_f%pars_xy(ind_hc3_win_r,j,i)
    4744 
    4745               IF ( building_pars_f%pars_xy(ind_tc1_wall_r,j,i) /= building_pars_f%fill )  THEN
    4746                  surf_usm_h%lambda_h(nzb_wall,m)   = building_pars_f%pars_xy(ind_tc1_wall_r,j,i)         
    4747                  surf_usm_h%lambda_h(nzb_wall+1,m) = building_pars_f%pars_xy(ind_tc1_wall_r,j,i)       
     5416              IF ( building_pars_f%pars_xy(ind_hc2,j,i) /=                     &
     5417                   building_pars_f%fill )                                      &
     5418                 surf_usm_h%rho_c_window(nzb_wall+2,m) =                       &
     5419                                    building_pars_f%pars_xy(ind_hc2,j,i)
     5420                 
     5421              IF ( building_pars_f%pars_xy(ind_hc3,j,i) /=                     &
     5422                   building_pars_f%fill )                                      &
     5423                 surf_usm_h%rho_c_window(nzb_wall+3,m) =                       &
     5424                                    building_pars_f%pars_xy(ind_hc3,j,i)
     5425
     5426              IF ( building_pars_f%pars_xy(ind_tc1,j,i) /=                     &
     5427                   building_pars_f%fill )  THEN
     5428                 surf_usm_h%lambda_h(nzb_wall,m)   =                           &
     5429                                    building_pars_f%pars_xy(ind_tc1,j,i)         
     5430                 surf_usm_h%lambda_h(nzb_wall+1,m) =                           &
     5431                                    building_pars_f%pars_xy(ind_tc1,j,i)       
    47485432              ENDIF
    4749               IF ( building_pars_f%pars_xy(ind_tc2_wall_r,j,i) /= building_pars_f%fill )    &
    4750                  surf_usm_h%lambda_h(nzb_wall+2,m) = building_pars_f%pars_xy(ind_tc2_wall_r,j,i)
    4751               IF ( building_pars_f%pars_xy(ind_tc3_wall_r,j,i) /= building_pars_f%fill )    &
    4752                  surf_usm_h%lambda_h(nzb_wall+3,m) = building_pars_f%pars_xy(ind_tc3_wall_r,j,i)   
    4753               IF ( building_pars_f%pars_xy(ind_tc1_wall_r,j,i) /= building_pars_f%fill )  THEN
    4754                  surf_usm_h%lambda_h_green(nzb_wall,m)   = lambda_h_green_sm !building_pars_f%pars_xy(ind_tc1_wall_r,j,i)         
    4755                  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)       
     5433              IF ( building_pars_f%pars_xy(ind_tc2,j,i) /=                     & 
     5434                   building_pars_f%fill )                                      &
     5435                 surf_usm_h%lambda_h(nzb_wall+2,m) =                           &
     5436                                    building_pars_f%pars_xy(ind_tc2,j,i)
     5437                 
     5438              IF ( building_pars_f%pars_xy(ind_tc3,j,i) /=                     &
     5439                   building_pars_f%fill )                                      &
     5440                 surf_usm_h%lambda_h(nzb_wall+3,m) =                           &
     5441                                    building_pars_f%pars_xy(ind_tc3,j,i)   
     5442                 
     5443              IF ( building_pars_f%pars_xy(ind_tc1,j,i) /=                     &
     5444                   building_pars_f%fill )  THEN
     5445                 surf_usm_h%lambda_h_green(nzb_wall,m)   =                     &
     5446                                     building_pars_f%pars_xy(ind_tc1,j,i)         
     5447                 surf_usm_h%lambda_h_green(nzb_wall+1,m) =                     &
     5448                                     building_pars_f%pars_xy(ind_tc1,j,i)       
    47565449              ENDIF
    4757               IF ( building_pars_f%pars_xy(ind_tc2_wall_r,j,i) /= building_pars_f%fill )    &
    4758                  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)
    4759               IF ( building_pars_f%pars_xy(ind_tc3_wall_r,j,i) /= building_pars_f%fill )    &
    4760                  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)   
    4761               IF ( building_pars_f%pars_xy(ind_tc1_win_r,j,i) /= building_pars_f%fill )  THEN
    4762                  surf_usm_h%lambda_h_window(nzb_wall,m)   = building_pars_f%pars_xy(ind_tc1_win_r,j,i)         
    4763                  surf_usm_h%lambda_h_window(nzb_wall+1,m) = building_pars_f%pars_xy(ind_tc1_win_r,j,i)       
     5450              IF ( building_pars_f%pars_xy(ind_tc2,j,i) /=                     &
     5451                   building_pars_f%fill )                                      &
     5452                 surf_usm_h%lambda_h_green(nzb_wall+2,m) =                     &
     5453                                    building_pars_f%pars_xy(ind_tc2,j,i)
     5454                 
     5455              IF ( building_pars_f%pars_xy(ind_tc3,j,i) /=                     &       
     5456                   building_pars_f%fill )                                      &
     5457                 surf_usm_h%lambda_h_green(nzb_wall+3,m) =                     &
     5458                                    building_pars_f%pars_xy(ind_tc3,j,i)   
     5459                 
     5460              IF ( building_pars_f%pars_xy(ind_tc1_win_r,j,i) /=               &
     5461                   building_pars_f%fill )  THEN
     5462                 surf_usm_h%lambda_h_window(nzb_wall,m)   =                    &
     5463                                     building_pars_f%pars_xy(ind_tc1_win_r,j,i)         
     5464                 surf_usm_h%lambda_h_window(nzb_wall+1,m) =                    &
     5465                                     building_pars_f%pars_xy(ind_tc1_win_r,j,i)       
    47645466              ENDIF
    4765               IF ( building_pars_f%pars_xy(ind_tc2_win_r,j,i) /= building_pars_f%fill )    &
    4766                  surf_usm_h%lambda_h_window(nzb_wall+2,m) = building_pars_f%pars_xy(ind_tc2_win_r,j,i)
    4767               IF ( building_pars_f%pars_xy(ind_tc3_win_r,j,i) /= building_pars_f%fill )    &
    4768                  surf_usm_h%lambda_h_window(nzb_wall+3,m) = building_pars_f%pars_xy(ind_tc3_win_r,j,i)   
    4769 
    4770               IF ( building_pars_f%pars_xy(117,j,i) /= building_pars_f%fill )         &
    4771                  surf_usm_h%target_temp_summer(m)  = building_pars_f%pars_xy(117,j,i)   
    4772               IF ( building_pars_f%pars_xy(118,j,i) /= building_pars_f%fill )         &
    4773                  surf_usm_h%target_temp_winter(m)  = building_pars_f%pars_xy(118,j,i)   
    4774 
    4775               IF ( building_pars_f%pars_xy(ind_emis_wall_r,j,i) /= building_pars_f%fill ) &
    4776                  surf_usm_h%emissivity(ind_veg_wall,m)  = building_pars_f%pars_xy(ind_emis_wall_r,j,i)
    4777               IF ( building_pars_f%pars_xy(ind_emis_green_r,j,i) /= building_pars_f%fill )&
    4778                  surf_usm_h%emissivity(ind_pav_green,m) = building_pars_f%pars_xy(ind_emis_green_r,j,i)
    4779               IF ( building_pars_f%pars_xy(ind_emis_win_r,j,i) /= building_pars_f%fill )  &
    4780                  surf_usm_h%emissivity(ind_wat_win,m)   = building_pars_f%pars_xy(ind_emis_win_r,j,i)
    4781 
    4782               IF ( building_pars_f%pars_xy(ind_trans_r,j,i) /= building_pars_f%fill )    &
    4783                  surf_usm_h%transmissivity(m)      = building_pars_f%pars_xy(ind_trans_r,j,i)
    4784 
    4785               IF ( building_pars_f%pars_xy(ind_z0,j,i) /= building_pars_f%fill )    &
    4786                  surf_usm_h%z0(m)                  = building_pars_f%pars_xy(ind_z0,j,i)
    4787               IF ( building_pars_f%pars_xy(ind_z0qh,j,i) /= building_pars_f%fill )  &
    4788                  surf_usm_h%z0h(m)                 = building_pars_f%pars_xy(ind_z0qh,j,i)
    4789               IF ( building_pars_f%pars_xy(ind_z0qh,j,i) /= building_pars_f%fill )  &
    4790                  surf_usm_h%z0q(m)                 = building_pars_f%pars_xy(ind_z0qh,j,i)
    4791 
    4792               IF ( building_pars_f%pars_xy(ind_alb_wall_r,j,i) /= building_pars_f%fill )    &
    4793                  surf_usm_h%albedo_type(ind_veg_wall,m)  = building_pars_f%pars_xy(ind_alb_wall_r,j,i)
    4794               IF ( building_pars_f%pars_xy(ind_alb_green_r,j,i) /= building_pars_f%fill )   &
    4795                  surf_usm_h%albedo_type(ind_pav_green,m) = building_pars_f%pars_xy(ind_alb_green_r,j,i)
    4796               IF ( building_pars_f%pars_xy(ind_alb_win_r,j,i) /= building_pars_f%fill )     &
    4797                  surf_usm_h%albedo_type(ind_wat_win,m)   = building_pars_f%pars_xy(ind_alb_win_r,j,i)
    4798 
    4799               IF ( building_pars_f%pars_xy(ind_thick_1_wall_r,j,i) /= building_pars_f%fill )    &
    4800                  surf_usm_h%zw(nzb_wall,m) = building_pars_f%pars_xy(ind_thick_1_wall_r,j,i)
    4801               IF ( building_pars_f%pars_xy(ind_thick_2_wall_r,j,i) /= building_pars_f%fill )    &
    4802                  surf_usm_h%zw(nzb_wall+1,m) = building_pars_f%pars_xy(ind_thick_2_wall_r,j,i)
    4803               IF ( building_pars_f%pars_xy(ind_thick_3_wall_r,j,i) /= building_pars_f%fill )    &
    4804                  surf_usm_h%zw(nzb_wall+2,m) = building_pars_f%pars_xy(ind_thick_3_wall_r,j,i)
    4805               IF ( building_pars_f%pars_xy(ind_thick_4_wall_r,j,i) /= building_pars_f%fill )    &
    4806                  surf_usm_h%zw(nzb_wall+3,m) = building_pars_f%pars_xy(ind_thick_4_wall_r,j,i)
    4807               IF ( building_pars_f%pars_xy(ind_thick_1_wall_r,j,i) /= building_pars_f%fill )    &
    4808                  surf_usm_h%zw_green(nzb_wall,m) = building_pars_f%pars_xy(ind_thick_1_wall_r,j,i)
    4809               IF ( building_pars_f%pars_xy(ind_thick_2_wall_r,j,i) /= building_pars_f%fill )    &
    4810                  surf_usm_h%zw_green(nzb_wall+1,m) = building_pars_f%pars_xy(ind_thick_2_wall_r,j,i)
    4811               IF ( building_pars_f%pars_xy(ind_thick_3_wall_r,j,i) /= building_pars_f%fill )    &
    4812                  surf_usm_h%zw_green(nzb_wall+2,m) = building_pars_f%pars_xy(ind_thick_3_wall_r,j,i)
    4813               IF ( building_pars_f%pars_xy(ind_thick_4_wall_r,j,i) /= building_pars_f%fill )    &
    4814                  surf_usm_h%zw_green(nzb_wall+3,m) = building_pars_f%pars_xy(ind_thick_4_wall_r,j,i)
    4815               IF ( building_pars_f%pars_xy(ind_thick_1_win_r,j,i) /= building_pars_f%fill )     &
    4816                  surf_usm_h%zw_window(nzb_wall,m) = building_pars_f%pars_xy(ind_thick_1_win_r,j,i)
    4817               IF ( building_pars_f%pars_xy(ind_thick_2_win_r,j,i) /= building_pars_f%fill )     &
    4818                  surf_usm_h%zw_window(nzb_wall+1,m) = building_pars_f%pars_xy(ind_thick_2_win_r,j,i)
    4819               IF ( building_pars_f%pars_xy(ind_thick_3_win_r,j,i) /= building_pars_f%fill )     &
    4820                  surf_usm_h%zw_window(nzb_wall+2,m) = building_pars_f%pars_xy(ind_thick_3_win_r,j,i)
    4821               IF ( building_pars_f%pars_xy(ind_thick_4_win_r,j,i) /= building_pars_f%fill )     &
    4822                  surf_usm_h%zw_window(nzb_wall+3,m) = building_pars_f%pars_xy(ind_thick_4_win_r,j,i)
    4823 
    4824               IF ( building_pars_f%pars_xy(0,j,i) /= building_pars_f%fill )    &
    4825                  surf_usm_h%c_surface(m)           = building_pars_f%pars_xy(0,j,i)
    4826               IF ( building_pars_f%pars_xy(3,j,i) /= building_pars_f%fill )    &
    4827                  surf_usm_h%lambda_surf(m)         = building_pars_f%pars_xy(3,j,i)
    4828               IF ( building_pars_f%pars_xy(2,j,i) /= building_pars_f%fill )    &
    4829                  surf_usm_h%c_surface_green(m)           = building_pars_f%pars_xy(2,j,i)
    4830               IF ( building_pars_f%pars_xy(5,j,i) /= building_pars_f%fill )    &
    4831                  surf_usm_h%lambda_surf_green(m)         = building_pars_f%pars_xy(5,j,i)
    4832               IF ( building_pars_f%pars_xy(1,j,i) /= building_pars_f%fill )    &
    4833                  surf_usm_h%c_surface_window(m)           = building_pars_f%pars_xy(1,j,i)
    4834               IF ( building_pars_f%pars_xy(4,j,i) /= building_pars_f%fill )    &
    4835                  surf_usm_h%lambda_surf_window(m)         = building_pars_f%pars_xy(4,j,i)
    4836 
    4837               IF ( building_pars_f%pars_xy(ind_green_type_roof,j,i) /= building_pars_f%fill )    &
    4838                  surf_usm_h%green_type_roof(m)            = building_pars_f%pars_xy(ind_green_type_roof,j,i)
     5467              IF ( building_pars_f%pars_xy(ind_tc2,j,i) /=                     &     
     5468                   building_pars_f%fill )                                      &
     5469                 surf_usm_h%lambda_h_window(nzb_wall+2,m) =                    &
     5470                                     building_pars_f%pars_xy(ind_tc2,j,i)
     5471                 
     5472              IF ( building_pars_f%pars_xy(ind_tc3,j,i) /=                     &   
     5473                   building_pars_f%fill )                                      &
     5474                 surf_usm_h%lambda_h_window(nzb_wall+3,m) =                    &
     5475                                    building_pars_f%pars_xy(ind_tc3,j,i)   
     5476
     5477              IF ( building_pars_f%pars_xy(ind_indoor_target_temp_summer,j,i) /=&           
     5478                   building_pars_f%fill )                                      &
     5479                 surf_usm_h%target_temp_summer(m)  =                           &
     5480                      building_pars_f%pars_xy(ind_indoor_target_temp_summer,j,i)   
     5481              IF ( building_pars_f%pars_xy(ind_indoor_target_temp_winter,j,i) /=&           
     5482                   building_pars_f%fill )                                      &
     5483                 surf_usm_h%target_temp_winter(m)  =                           &
     5484                      building_pars_f%pars_xy(ind_indoor_target_temp_winter,j,i)   
     5485
     5486              IF ( building_pars_f%pars_xy(ind_emis_wall,j,i) /=               &   
     5487                   building_pars_f%fill )                                      &
     5488                 surf_usm_h%emissivity(ind_veg_wall,m)  =                      &
     5489                                    building_pars_f%pars_xy(ind_emis_wall,j,i)
     5490                 
     5491              IF ( building_pars_f%pars_xy(ind_emis_green,j,i) /=              &           
     5492                   building_pars_f%fill )                                      &
     5493                 surf_usm_h%emissivity(ind_pav_green,m) =                      &
     5494                                     building_pars_f%pars_xy(ind_emis_green,j,i)
     5495                 
     5496              IF ( building_pars_f%pars_xy(ind_emis_win,j,i) /=                &
     5497                   building_pars_f%fill )                                      &
     5498                 surf_usm_h%emissivity(ind_wat_win,m)   =                      &
     5499                                     building_pars_f%pars_xy(ind_emis_win,j,i)
     5500                 
     5501              IF ( building_pars_f%pars_xy(ind_trans,j,i) /=                   &   
     5502                   building_pars_f%fill )                                      &
     5503                 surf_usm_h%transmissivity(m) =                                &
     5504                                    building_pars_f%pars_xy(ind_trans,j,i)
     5505
     5506              IF ( building_pars_f%pars_xy(ind_z0,j,i) /=                      &         
     5507                   building_pars_f%fill )                                      &
     5508                 surf_usm_h%z0(m) = building_pars_f%pars_xy(ind_z0,j,i)
     5509                 
     5510              IF ( building_pars_f%pars_xy(ind_z0qh,j,i) /=                    &           
     5511                   building_pars_f%fill )                                      &
     5512                 surf_usm_h%z0h(m) = building_pars_f%pars_xy(ind_z0qh,j,i)
     5513              IF ( building_pars_f%pars_xy(ind_z0qh,j,i) /=                    &           
     5514                   building_pars_f%fill )                                      &
     5515                 surf_usm_h%z0q(m) = building_pars_f%pars_xy(ind_z0qh,j,i)
     5516
     5517              IF ( building_pars_f%pars_xy(ind_alb_wall_agfl,j,i) /=           &         
     5518                   building_pars_f%fill )                                      &
     5519                 surf_usm_h%albedo_type(ind_veg_wall,m)  =                     &
     5520                                 building_pars_f%pars_xy(ind_alb_wall_agfl,j,i)
     5521                 
     5522              IF ( building_pars_f%pars_xy(ind_alb_green_agfl,j,i) /=          &           
     5523                   building_pars_f%fill )                                      &
     5524                 surf_usm_h%albedo_type(ind_pav_green,m) =                     &
     5525                                building_pars_f%pars_xy(ind_alb_green_agfl,j,i)
     5526              IF ( building_pars_f%pars_xy(ind_alb_win_agfl,j,i) /=            &         
     5527                   building_pars_f%fill )                                      &
     5528                 surf_usm_h%albedo_type(ind_wat_win,m)   =                     &
     5529                                   building_pars_f%pars_xy(ind_alb_win_agfl,j,i)
     5530
     5531              IF ( building_pars_f%pars_xy(ind_thick_1_agfl,j,i) /=            &         
     5532                   building_pars_f%fill )                                      &
     5533                 surf_usm_h%zw(nzb_wall,m) =                                   &
     5534                                  building_pars_f%pars_xy(ind_thick_1_agfl,j,i)
     5535                 
     5536              IF ( building_pars_f%pars_xy(ind_thick_2_agfl,j,i) /=            &         
     5537                   building_pars_f%fill )                                      &
     5538                 surf_usm_h%zw(nzb_wall+1,m) =                                 &
     5539                                  building_pars_f%pars_xy(ind_thick_2_agfl,j,i)
     5540                 
     5541              IF ( building_pars_f%pars_xy(ind_thick_3_agfl,j,i) /=            &         
     5542                   building_pars_f%fill )                                      &
     5543                 surf_usm_h%zw(nzb_wall+2,m) =                                 &
     5544                                  building_pars_f%pars_xy(ind_thick_3_agfl,j,i)
     5545                 
     5546                 
     5547              IF ( building_pars_f%pars_xy(ind_thick_4_agfl,j,i) /=            &         
     5548                   building_pars_f%fill )                                      &
     5549                 surf_usm_h%zw(nzb_wall+3,m) =                                 &
     5550                                  building_pars_f%pars_xy(ind_thick_4_agfl,j,i)
     5551                 
     5552              IF ( building_pars_f%pars_xy(ind_thick_1_agfl,j,i) /=            &           
     5553                   building_pars_f%fill )                                      &
     5554                 surf_usm_h%zw_green(nzb_wall,m) =                             &
     5555                                  building_pars_f%pars_xy(ind_thick_1_agfl,j,i)
     5556                 
     5557              IF ( building_pars_f%pars_xy(ind_thick_2_agfl,j,i) /=            &         
     5558                   building_pars_f%fill )                                      &
     5559                 surf_usm_h%zw_green(nzb_wall+1,m) =                           &
     5560                                   building_pars_f%pars_xy(ind_thick_2_agfl,j,i)
     5561                 
     5562              IF ( building_pars_f%pars_xy(ind_thick_3_agfl,j,i) /=            &         
     5563                   building_pars_f%fill )                                      &
     5564                 surf_usm_h%zw_green(nzb_wall+2,m) =                           &
     5565                                   building_pars_f%pars_xy(ind_thick_3_agfl,j,i)
     5566                 
     5567              IF ( building_pars_f%pars_xy(ind_thick_4_agfl,j,i) /=            &         
     5568                   building_pars_f%fill )                                      &
     5569                 surf_usm_h%zw_green(nzb_wall+3,m) =                           &
     5570                                   building_pars_f%pars_xy(ind_thick_4_agfl,j,i)
     5571
     5572              IF ( building_pars_f%pars_xy(ind_c_surface,j,i) /=               &       
     5573                   building_pars_f%fill )                                      &
     5574                 surf_usm_h%c_surface(m) =                                     &
     5575                                    building_pars_f%pars_xy(ind_c_surface,j,i)
     5576                 
     5577              IF ( building_pars_f%pars_xy(ind_lambda_surf,j,i) /=             &       
     5578                   building_pars_f%fill )                                      &
     5579                 surf_usm_h%lambda_surf(m) =                                   &
     5580                                    building_pars_f%pars_xy(ind_lambda_surf,j,i)
     5581             
    48395582           ENDDO
    48405583
     
    48475590               
    48485591!
    4849 !--              In order to distinguish between ground floor level and
    4850 !--              above-ground-floor level surfaces, set input indices.
    4851                  ind_alb_green    = MERGE( ind_alb_green_gfl,    ind_alb_green_agfl,    &
    4852                                            surf_usm_v(l)%ground_level(m) )
    4853                  ind_alb_wall     = MERGE( ind_alb_wall_gfl,     ind_alb_wall_agfl,     &
    4854                                            surf_usm_v(l)%ground_level(m) )
    4855                  ind_alb_win      = MERGE( ind_alb_win_gfl,      ind_alb_win_agfl,      &
    4856                                            surf_usm_v(l)%ground_level(m) )
    4857                  ind_wall_frac    = MERGE( ind_wall_frac_gfl,    ind_wall_frac_agfl,    &
    4858                                            surf_usm_v(l)%ground_level(m) )
    4859                  ind_win_frac     = MERGE( ind_win_frac_gfl,     ind_win_frac_agfl,     &
    4860                                            surf_usm_v(l)%ground_level(m) )
    4861                  ind_green_frac_w = MERGE( ind_green_frac_w_gfl, ind_green_frac_w_agfl, &
    4862                                            surf_usm_v(l)%ground_level(m) )
    4863                  ind_green_frac_r = MERGE( ind_green_frac_r_gfl, ind_green_frac_r_agfl, &
    4864                                            surf_usm_v(l)%ground_level(m) )
    4865                  ind_lai_r        = MERGE( ind_lai_r_gfl,        ind_lai_r_agfl,        &
    4866                                            surf_usm_v(l)%ground_level(m) )
    4867                  ind_lai_w        = MERGE( ind_lai_w_gfl,        ind_lai_w_agfl,        &
    4868                                            surf_usm_v(l)%ground_level(m) )
    4869                  ind_hc1          = MERGE( ind_hc1_gfl,          ind_hc1_agfl,          &
    4870                                            surf_usm_v(l)%ground_level(m) )
    4871                  ind_hc1_win      = MERGE( ind_hc1_win_gfl,      ind_hc1_win_agfl,      &
    4872                                            surf_usm_v(l)%ground_level(m) )
    4873                  ind_hc2          = MERGE( ind_hc2_gfl,          ind_hc2_agfl,          &
    4874                                            surf_usm_v(l)%ground_level(m) )
    4875                  ind_hc2_win      = MERGE( ind_hc2_win_gfl,      ind_hc2_win_agfl,      &
    4876                                            surf_usm_v(l)%ground_level(m) )
    4877                  ind_hc3          = MERGE( ind_hc3_gfl,          ind_hc3_agfl,          &
    4878                                            surf_usm_v(l)%ground_level(m) )
    4879                  ind_hc3_win      = MERGE( ind_hc3_win_gfl,      ind_hc3_win_agfl,      &
    4880                                            surf_usm_v(l)%ground_level(m) )
    4881                  ind_tc1          = MERGE( ind_tc1_gfl,          ind_tc1_agfl,          &
    4882                                            surf_usm_v(l)%ground_level(m) )
    4883                  ind_tc1_win      = MERGE( ind_tc1_win_gfl,      ind_tc1_win_agfl,      &
    4884                                            surf_usm_v(l)%ground_level(m) )
    4885                  ind_tc2          = MERGE( ind_tc2_gfl,          ind_tc2_agfl,          &
    4886                                            surf_usm_v(l)%ground_level(m) )
    4887                  ind_tc2_win      = MERGE( ind_tc2_win_gfl,      ind_tc2_win_agfl,      &
    4888                                            surf_usm_v(l)%ground_level(m) )
    4889                  ind_tc3          = MERGE( ind_tc3_gfl,          ind_tc3_agfl,          &
    4890                                            surf_usm_v(l)%ground_level(m) )
    4891                  ind_tc3_win      = MERGE( ind_tc3_win_gfl,      ind_tc3_win_agfl,      &
    4892                                            surf_usm_v(l)%ground_level(m) )
    4893                  ind_thick_1      = MERGE( ind_thick_1_gfl,      ind_thick_1_agfl,      &
    4894                                            surf_usm_v(l)%ground_level(m) )
    4895                  ind_thick_1_win  = MERGE( ind_thick_1_win_gfl,  ind_thick_1_win_agfl,  &
    4896                                            surf_usm_v(l)%ground_level(m) )
    4897                  ind_thick_2      = MERGE( ind_thick_2_gfl,      ind_thick_2_agfl,      &
    4898                                            surf_usm_v(l)%ground_level(m) )
    4899                  ind_thick_2_win  = MERGE( ind_thick_2_win_gfl,  ind_thick_2_win_agfl,  &
    4900                                            surf_usm_v(l)%ground_level(m) )
    4901                  ind_thick_3      = MERGE( ind_thick_3_gfl,      ind_thick_3_agfl,      &
    4902                                            surf_usm_v(l)%ground_level(m) )
    4903                  ind_thick_3_win  = MERGE( ind_thick_3_win_gfl,  ind_thick_3_win_agfl,  &
    4904                                            surf_usm_v(l)%ground_level(m) )
    4905                  ind_thick_4      = MERGE( ind_thick_4_gfl,      ind_thick_4_agfl,      &
    4906                                            surf_usm_v(l)%ground_level(m) )
    4907                  ind_thick_4_win  = MERGE( ind_thick_4_win_gfl,  ind_thick_4_win_agfl,  &
    4908                                            surf_usm_v(l)%ground_level(m) )
    4909                  ind_emis_wall    = MERGE( ind_emis_wall_gfl,    ind_emis_wall_agfl,    &
    4910                                            surf_usm_v(l)%ground_level(m) )
    4911                  ind_emis_green   = MERGE( ind_emis_green_gfl,   ind_emis_green_agfl,   &
    4912                                            surf_usm_v(l)%ground_level(m) )
    4913                  ind_emis_win     = MERGE( ind_emis_win_gfl,     ind_emis_win_agfl,     &
    4914                                            surf_usm_v(l)%ground_level(m) )
    4915                  ind_trans        = MERGE( ind_trans_gfl,       ind_trans_agfl,         &
    4916                                            surf_usm_v(l)%ground_level(m) )
    4917                  ind_z0           = MERGE( ind_z0_gfl,           ind_z0_agfl,           &
    4918                                            surf_usm_v(l)%ground_level(m) )
    4919                  ind_z0qh         = MERGE( ind_z0qh_gfl,         ind_z0qh_agfl,         &
     5592!--                 In order to distinguish between ground floor level and
     5593!--                 above-ground-floor level surfaces, set input indices.
     5594                    ind_wall_frac    = MERGE( ind_wall_frac_gfl,               &
     5595                                              ind_wall_frac_agfl,              &
    49205596                                              surf_usm_v(l)%ground_level(m) )
    4921 
    4922 !
    4923 !--              Initialize relatvie wall- (0), green- (1) and window (2) fractions
    4924                  IF ( building_pars_f%pars_xy(ind_wall_frac,j,i) /=                     &
    4925                       building_pars_f%fill )                                            &
    4926                     surf_usm_v(l)%frac(ind_veg_wall,m)  =                               &
    4927                                       building_pars_f%pars_xy(ind_wall_frac,j,i)   
    4928                  IF ( building_pars_f%pars_xy(ind_green_frac_w,j,i) /=                  &
    4929                       building_pars_f%fill )                                            &
    4930                     surf_usm_v(l)%frac(ind_pav_green,m) =                               &
    4931                                       building_pars_f%pars_xy(ind_green_frac_w,j,i)
    4932                  IF ( building_pars_f%pars_xy(ind_win_frac,j,i) /=                      &
    4933                       building_pars_f%fill )                                            &
    4934                     surf_usm_v(l)%frac(ind_wat_win,m)   =                               &
    4935                                       building_pars_f%pars_xy(ind_win_frac,j,i)
    4936  
    4937                  IF ( building_pars_f%pars_xy(ind_lai_w,j,i) /= building_pars_f%fill )  &
    4938                     surf_usm_v(l)%lai(m) = building_pars_f%pars_xy(ind_lai_w,j,i)
    4939 
    4940                  IF ( building_pars_f%pars_xy(ind_hc1,j,i) /= building_pars_f%fill )    &
    4941                  THEN
    4942                     surf_usm_v(l)%rho_c_wall(nzb_wall,m)   =                            &
    4943                                                     building_pars_f%pars_xy(ind_hc1,j,i) 
    4944                     surf_usm_v(l)%rho_c_wall(nzb_wall+1,m) =                            &
    4945                                                     building_pars_f%pars_xy(ind_hc1,j,i)
    4946                  ENDIF
    4947                  IF ( building_pars_f%pars_xy(ind_hc2,j,i) /= building_pars_f%fill )    &
    4948                     surf_usm_v(l)%rho_c_wall(nzb_wall+2,m) =                            &                           
    4949                                                     building_pars_f%pars_xy(ind_hc2,j,i)
    4950                  IF ( building_pars_f%pars_xy(ind_hc3,j,i) /= building_pars_f%fill )    &
    4951                     surf_usm_v(l)%rho_c_wall(nzb_wall+3,m) =                            &
    4952                                                     building_pars_f%pars_xy(ind_hc3,j,i)
    4953                  IF ( building_pars_f%pars_xy(ind_hc1,j,i) /= building_pars_f%fill )  THEN
    4954                     surf_usm_v(l)%rho_c_green(nzb_wall,m)   =                           &
    4955                                                     rho_c_soil !building_pars_f%pars_xy(ind_hc1,j,i) 
    4956                     surf_usm_v(l)%rho_c_green(nzb_wall+1,m) =                           &
    4957                                                     rho_c_soil !building_pars_f%pars_xy(ind_hc1,j,i)
    4958                  ENDIF
    4959                  IF ( building_pars_f%pars_xy(ind_hc2,j,i) /= building_pars_f%fill )    &
    4960                     surf_usm_v(l)%rho_c_green(nzb_wall+2,m) = rho_c_soil !building_pars_f%pars_xy(ind_hc2,j,i)
    4961                  IF ( building_pars_f%pars_xy(ind_hc3,j,i) /= building_pars_f%fill )    &
    4962                     surf_usm_v(l)%rho_c_green(nzb_wall+3,m) = rho_c_soil !building_pars_f%pars_xy(ind_hc3,j,i)   
    4963                  IF ( building_pars_f%pars_xy(ind_hc1_win,j,i) /= building_pars_f%fill )  THEN
    4964                     surf_usm_v(l)%rho_c_window(nzb_wall,m)   = building_pars_f%pars_xy(ind_hc1_win,j,i) 
    4965                     surf_usm_v(l)%rho_c_window(nzb_wall+1,m) = building_pars_f%pars_xy(ind_hc1_win,j,i)
    4966                  ENDIF
    4967                  IF ( building_pars_f%pars_xy(ind_hc2_win,j,i) /= building_pars_f%fill )    &
    4968                     surf_usm_v(l)%rho_c_window(nzb_wall+2,m) = building_pars_f%pars_xy(ind_hc2_win,j,i)
    4969                  IF ( building_pars_f%pars_xy(ind_hc3_win,j,i) /= building_pars_f%fill )    &
    4970                     surf_usm_v(l)%rho_c_window(nzb_wall+3,m) = building_pars_f%pars_xy(ind_hc3_win,j,i)
    4971 
    4972                  IF ( building_pars_f%pars_xy(ind_tc1,j,i) /= building_pars_f%fill )  THEN
    4973                     surf_usm_v(l)%lambda_h(nzb_wall,m)   = building_pars_f%pars_xy(ind_tc1,j,i)         
    4974                     surf_usm_v(l)%lambda_h(nzb_wall+1,m) = building_pars_f%pars_xy(ind_tc1,j,i)       
    4975                  ENDIF
    4976                  IF ( building_pars_f%pars_xy(ind_tc2,j,i) /= building_pars_f%fill )    &
    4977                     surf_usm_v(l)%lambda_h(nzb_wall+2,m) = building_pars_f%pars_xy(ind_tc2,j,i)
    4978                  IF ( building_pars_f%pars_xy(ind_tc3,j,i) /= building_pars_f%fill )    &
    4979                     surf_usm_v(l)%lambda_h(nzb_wall+3,m) = building_pars_f%pars_xy(ind_tc3,j,i)   
    4980                  IF ( building_pars_f%pars_xy(ind_tc1,j,i) /= building_pars_f%fill )  THEN
    4981                     surf_usm_v(l)%lambda_h_green(nzb_wall,m)   = lambda_h_green_sm !building_pars_f%pars_xy(ind_tc1,j,i)         
    4982                     surf_usm_v(l)%lambda_h_green(nzb_wall+1,m) = lambda_h_green_sm !building_pars_f%pars_xy(ind_tc1,j,i)       
    4983                  ENDIF
    4984                  IF ( building_pars_f%pars_xy(ind_tc2,j,i) /= building_pars_f%fill )    &
    4985                     surf_usm_v(l)%lambda_h_green(nzb_wall+2,m) = lambda_h_green_sm !building_pars_f%pars_xy(ind_tc2,j,i)
    4986                  IF ( building_pars_f%pars_xy(ind_tc3,j,i) /= building_pars_f%fill )    &
    4987                     surf_usm_v(l)%lambda_h_green(nzb_wall+3,m) = lambda_h_green_sm !building_pars_f%pars_xy(ind_tc3,j,i)   
    4988                  IF ( building_pars_f%pars_xy(ind_tc1_win,j,i) /= building_pars_f%fill )  THEN
    4989                     surf_usm_v(l)%lambda_h_window(nzb_wall,m)   = building_pars_f%pars_xy(ind_tc1_win,j,i)         
    4990                     surf_usm_v(l)%lambda_h_window(nzb_wall+1,m) = building_pars_f%pars_xy(ind_tc1_win,j,i)       
    4991                  ENDIF
    4992                  IF ( building_pars_f%pars_xy(ind_tc2_win,j,i) /= building_pars_f%fill )    &
    4993                     surf_usm_v(l)%lambda_h_window(nzb_wall+2,m) = building_pars_f%pars_xy(ind_tc2_win,j,i)
    4994                  IF ( building_pars_f%pars_xy(ind_tc3_win,j,i) /= building_pars_f%fill )    &
    4995                     surf_usm_v(l)%lambda_h_window(nzb_wall+3,m) = building_pars_f%pars_xy(ind_tc3_win,j,i)
    4996 
    4997                  IF ( building_pars_f%pars_xy(117,j,i) /= building_pars_f%fill )         &
    4998                     surf_usm_v(l)%target_temp_summer(m)  = building_pars_f%pars_xy(117,j,i)   
    4999                  IF ( building_pars_f%pars_xy(118,j,i) /= building_pars_f%fill )         &
    5000                     surf_usm_v(l)%target_temp_winter(m)  = building_pars_f%pars_xy(118,j,i)   
    5001 
    5002                  IF ( building_pars_f%pars_xy(ind_emis_wall,j,i) /= building_pars_f%fill ) &
    5003                     surf_usm_v(l)%emissivity(ind_veg_wall,m)  = building_pars_f%pars_xy(ind_emis_wall,j,i)
    5004                  IF ( building_pars_f%pars_xy(ind_emis_green,j,i) /= building_pars_f%fill )&
    5005                     surf_usm_v(l)%emissivity(ind_pav_green,m) = building_pars_f%pars_xy(ind_emis_green,j,i)
    5006                  IF ( building_pars_f%pars_xy(ind_emis_win,j,i) /= building_pars_f%fill )  &
    5007                     surf_usm_v(l)%emissivity(ind_wat_win,m)   = building_pars_f%pars_xy(ind_emis_win,j,i)
    5008 
    5009                  IF ( building_pars_f%pars_xy(ind_trans,j,i) /= building_pars_f%fill )    &
    5010                     surf_usm_v(l)%transmissivity(m)      = building_pars_f%pars_xy(ind_trans,j,i)
    5011 
    5012                  IF ( building_pars_f%pars_xy(ind_z0,j,i) /= building_pars_f%fill )    &
    5013                     surf_usm_v(l)%z0(m)                  = building_pars_f%pars_xy(ind_z0,j,i)
    5014                  IF ( building_pars_f%pars_xy(ind_z0qh,j,i) /= building_pars_f%fill )  &
    5015                     surf_usm_v(l)%z0h(m)                 = building_pars_f%pars_xy(ind_z0qh,j,i)
    5016                  IF ( building_pars_f%pars_xy(ind_z0qh,j,i) /= building_pars_f%fill )  &
    5017                     surf_usm_v(l)%z0q(m)                 = building_pars_f%pars_xy(ind_z0qh,j,i)
    5018 
    5019                  IF ( building_pars_f%pars_xy(ind_alb_wall,j,i) /= building_pars_f%fill )    &
    5020                     surf_usm_v(l)%albedo_type(ind_veg_wall,m)  = building_pars_f%pars_xy(ind_alb_wall,j,i)
    5021                  IF ( building_pars_f%pars_xy(ind_alb_green,j,i) /= building_pars_f%fill )   &
    5022                     surf_usm_v(l)%albedo_type(ind_pav_green,m) = building_pars_f%pars_xy(ind_alb_green,j,i)
    5023                  IF ( building_pars_f%pars_xy(ind_alb_win,j,i) /= building_pars_f%fill )     &
    5024                     surf_usm_v(l)%albedo_type(ind_wat_win,m)   = building_pars_f%pars_xy(ind_alb_win,j,i)
    5025 
    5026                  IF ( building_pars_f%pars_xy(ind_thick_1,j,i) /= building_pars_f%fill )    &
    5027                     surf_usm_v(l)%zw(nzb_wall,m) = building_pars_f%pars_xy(ind_thick_1,j,i)
    5028                  IF ( building_pars_f%pars_xy(ind_thick_2,j,i) /= building_pars_f%fill )    &
    5029                     surf_usm_v(l)%zw(nzb_wall+1,m) = building_pars_f%pars_xy(ind_thick_2,j,i)
    5030                  IF ( building_pars_f%pars_xy(ind_thick_3,j,i) /= building_pars_f%fill )    &
    5031                     surf_usm_v(l)%zw(nzb_wall+2,m) = building_pars_f%pars_xy(ind_thick_3,j,i)
    5032                  IF ( building_pars_f%pars_xy(ind_thick_4,j,i) /= building_pars_f%fill )    &
    5033                     surf_usm_v(l)%zw(nzb_wall+3,m) = building_pars_f%pars_xy(ind_thick_4,j,i)
    5034                  IF ( building_pars_f%pars_xy(ind_thick_1,j,i) /= building_pars_f%fill )    &
    5035                     surf_usm_v(l)%zw_green(nzb_wall,m) = building_pars_f%pars_xy(ind_thick_1,j,i)
    5036                  IF ( building_pars_f%pars_xy(ind_thick_2,j,i) /= building_pars_f%fill )    &
    5037                     surf_usm_v(l)%zw_green(nzb_wall+1,m) = building_pars_f%pars_xy(ind_thick_2,j,i)
    5038                  IF ( building_pars_f%pars_xy(ind_thick_3,j,i) /= building_pars_f%fill )    &
    5039                     surf_usm_v(l)%zw_green(nzb_wall+2,m) = building_pars_f%pars_xy(ind_thick_3,j,i)
    5040                  IF ( building_pars_f%pars_xy(ind_thick_4,j,i) /= building_pars_f%fill )    &
    5041                     surf_usm_v(l)%zw_green(nzb_wall+3,m) = building_pars_f%pars_xy(ind_thick_4,j,i)
    5042                  IF ( building_pars_f%pars_xy(ind_thick_1_win,j,i) /= building_pars_f%fill )    &
    5043                     surf_usm_v(l)%zw_window(nzb_wall,m) = building_pars_f%pars_xy(ind_thick_1_win,j,i)
    5044                  IF ( building_pars_f%pars_xy(ind_thick_2_win,j,i) /= building_pars_f%fill )    &
    5045                     surf_usm_v(l)%zw_window(nzb_wall+1,m) = building_pars_f%pars_xy(ind_thick_2_win,j,i)
    5046                  IF ( building_pars_f%pars_xy(ind_thick_3_win,j,i) /= building_pars_f%fill )    &
    5047                     surf_usm_v(l)%zw_window(nzb_wall+2,m) = building_pars_f%pars_xy(ind_thick_3_win,j,i)
    5048                  IF ( building_pars_f%pars_xy(ind_thick_4_win,j,i) /= building_pars_f%fill )    &
    5049                     surf_usm_v(l)%zw_window(nzb_wall+3,m) = building_pars_f%pars_xy(ind_thick_4_win,j,i)
    5050 
    5051                  IF ( building_pars_f%pars_xy(0,j,i) /= building_pars_f%fill )    &
    5052                     surf_usm_v(l)%c_surface(m)           = building_pars_f%pars_xy(0,j,i)
    5053                  IF ( building_pars_f%pars_xy(3,j,i) /= building_pars_f%fill )    &
    5054                     surf_usm_v(l)%lambda_surf(m)         = building_pars_f%pars_xy(3,j,i)
    5055                  IF ( building_pars_f%pars_xy(2,j,i) /= building_pars_f%fill )    &
    5056                     surf_usm_v(l)%c_surface_green(m)     = building_pars_f%pars_xy(2,j,i)
    5057                  IF ( building_pars_f%pars_xy(5,j,i) /= building_pars_f%fill )    &
    5058                     surf_usm_v(l)%lambda_surf_green(m)   = building_pars_f%pars_xy(5,j,i)
    5059                  IF ( building_pars_f%pars_xy(1,j,i) /= building_pars_f%fill )    &
    5060                     surf_usm_v(l)%c_surface_window(m)    = building_pars_f%pars_xy(1,j,i)
    5061                  IF ( building_pars_f%pars_xy(4,j,i) /= building_pars_f%fill )    &
    5062                     surf_usm_v(l)%lambda_surf_window(m)  = building_pars_f%pars_xy(4,j,i)
    5063 
     5597                    ind_green_frac_w = MERGE( ind_green_frac_w_gfl,            &
     5598                                              ind_green_frac_w_agfl,           &
     5599                                              surf_usm_v(l)%ground_level(m) )
     5600                    ind_win_frac     = MERGE( ind_win_frac_gfl,                &
     5601                                              ind_win_frac_agfl,               &
     5602                                              surf_usm_v(l)%ground_level(m) )
     5603                    ind_lai_w        = MERGE( ind_lai_w_gfl,                   &
     5604                                              ind_lai_w_agfl,                  &
     5605                                              surf_usm_v(l)%ground_level(m) )
     5606                    ind_z0           = MERGE( ind_z0_gfl,                      &
     5607                                              ind_z0_agfl,                     &
     5608                                              surf_usm_v(l)%ground_level(m) )
     5609                    ind_z0qh         = MERGE( ind_z0qh_gfl,                    &
     5610                                              ind_z0qh_agfl,                   &
     5611                                              surf_usm_v(l)%ground_level(m) )
     5612                    ind_hc1          = MERGE( ind_hc1_gfl,                     &
     5613                                              ind_hc1_agfl,                    &
     5614                                              surf_usm_v(l)%ground_level(m) )
     5615                    ind_hc2          = MERGE( ind_hc2_gfl,                     &
     5616                                              ind_hc2_agfl,                    &
     5617                                              surf_usm_v(l)%ground_level(m) )
     5618                    ind_hc3          = MERGE( ind_hc3_gfl,                     &
     5619                                              ind_hc3_agfl,                    &
     5620                                              surf_usm_v(l)%ground_level(m) )
     5621                    ind_tc1          = MERGE( ind_tc1_gfl,                     &
     5622                                              ind_tc1_agfl,                    &
     5623                                              surf_usm_v(l)%ground_level(m) )
     5624                    ind_tc2          = MERGE( ind_tc2_gfl,                     &
     5625                                              ind_tc2_agfl,                    &
     5626                                              surf_usm_v(l)%ground_level(m) )
     5627                    ind_tc3          = MERGE( ind_tc3_gfl,                     &
     5628                                              ind_tc3_agfl,                    &
     5629                                              surf_usm_v(l)%ground_level(m) )
     5630                    ind_emis_wall    = MERGE( ind_emis_wall_gfl,               &
     5631                                              ind_emis_wall_agfl,              &
     5632                                              surf_usm_v(l)%ground_level(m) )
     5633                    ind_emis_green   = MERGE( ind_emis_green_gfl,              &
     5634                                              ind_emis_green_agfl,             &
     5635                                              surf_usm_v(l)%ground_level(m) )
     5636                    ind_emis_win     = MERGE( ind_emis_win_gfl,                &
     5637                                              ind_emis_win_agfl,               &
     5638                                              surf_usm_v(l)%ground_level(m) )
     5639                    ind_trans        = MERGE( ind_trans_gfl,                   &
     5640                                              ind_trans_agfl,                  &
     5641                                              surf_usm_v(l)%ground_level(m) )
     5642                   
     5643!                   
     5644!--                 Initialize relatvie wall- (0), green- (1) and window (2) fractions
     5645                    IF ( building_pars_f%pars_xy(ind_wall_frac,j,i) /=         &
     5646                         building_pars_f%fill )                                &
     5647                       surf_usm_v(l)%frac(ind_veg_wall,m)  =                   &
     5648                                          building_pars_f%pars_xy(ind_wall_frac,j,i)   
     5649                       
     5650                    IF ( building_pars_f%pars_xy(ind_green_frac_w,j,i) /=      &         
     5651                         building_pars_f%fill )                                &
     5652                       surf_usm_v(l)%frac(ind_pav_green,m) =                   &
     5653                                  building_pars_f%pars_xy(ind_green_frac_w,j,i)
     5654                       
     5655                    IF ( building_pars_f%pars_xy(ind_win_frac,j,i) /=          &
     5656                         building_pars_f%fill )                                &
     5657                       surf_usm_v(l)%frac(ind_wat_win,m)   =                   &
     5658                                       building_pars_f%pars_xy(ind_win_frac,j,i)
     5659                   
     5660                    IF ( building_pars_f%pars_xy(ind_lai_w,j,i) /=             &
     5661                         building_pars_f%fill )                                &
     5662                       surf_usm_v(l)%lai(m)  =                                 &
     5663                                       building_pars_f%pars_xy(ind_lai_w,j,i)
     5664                   
     5665                    IF ( building_pars_f%pars_xy(ind_hc1,j,i) /=               &
     5666                         building_pars_f%fill )  THEN
     5667                       surf_usm_v(l)%rho_c_wall(nzb_wall,m)   =                &
     5668                                          building_pars_f%pars_xy(ind_hc1,j,i)
     5669                       surf_usm_v(l)%rho_c_wall(nzb_wall+1,m) =                &
     5670                                          building_pars_f%pars_xy(ind_hc1,j,i)
     5671                    ENDIF
     5672                   
     5673                   
     5674                    IF ( building_pars_f%pars_xy(ind_hc2,j,i) /=               &
     5675                         building_pars_f%fill )                                &
     5676                       surf_usm_v(l)%rho_c_wall(nzb_wall+2,m) =                &
     5677                                          building_pars_f%pars_xy(ind_hc2,j,i)
     5678                       
     5679                    IF ( building_pars_f%pars_xy(ind_hc3,j,i) /=               &         
     5680                         building_pars_f%fill )                                &
     5681                       surf_usm_v(l)%rho_c_wall(nzb_wall+3,m) =                &
     5682                                          building_pars_f%pars_xy(ind_hc3,j,i)
     5683                       
     5684                    IF ( building_pars_f%pars_xy(ind_hc1,j,i) /=               &
     5685                         building_pars_f%fill )  THEN
     5686                       surf_usm_v(l)%rho_c_green(nzb_wall,m)   =               &
     5687                                          building_pars_f%pars_xy(ind_hc1,j,i)
     5688                       surf_usm_v(l)%rho_c_green(nzb_wall+1,m) =               &
     5689                                          building_pars_f%pars_xy(ind_hc1,j,i)
     5690                    ENDIF
     5691                    IF ( building_pars_f%pars_xy(ind_hc2,j,i) /=               &
     5692                         building_pars_f%fill )                                &
     5693                       surf_usm_v(l)%rho_c_green(nzb_wall+2,m) =               &
     5694                                          building_pars_f%pars_xy(ind_hc2,j,i)
     5695                       
     5696                    IF ( building_pars_f%pars_xy(ind_hc3,j,i) /=               &
     5697                         building_pars_f%fill )                                &
     5698                       surf_usm_v(l)%rho_c_green(nzb_wall+3,m) =               &
     5699                                          building_pars_f%pars_xy(ind_hc3,j,i)
     5700                       
     5701                    IF ( building_pars_f%pars_xy(ind_hc1,j,i) /=               &
     5702                         building_pars_f%fill )  THEN
     5703                       surf_usm_v(l)%rho_c_window(nzb_wall,m)   =              &
     5704                                          building_pars_f%pars_xy(ind_hc1,j,i)
     5705                       surf_usm_v(l)%rho_c_window(nzb_wall+1,m) =              &
     5706                                          building_pars_f%pars_xy(ind_hc1,j,i)
     5707                    ENDIF
     5708                    IF ( building_pars_f%pars_xy(ind_hc2,j,i) /=               &
     5709                         building_pars_f%fill )                                &
     5710                       surf_usm_v(l)%rho_c_window(nzb_wall+2,m) =              &
     5711                                          building_pars_f%pars_xy(ind_hc2,j,i)
     5712                       
     5713                    IF ( building_pars_f%pars_xy(ind_hc3,j,i) /=               &
     5714                         building_pars_f%fill )                                &
     5715                       surf_usm_v(l)%rho_c_window(nzb_wall+3,m) =              &
     5716                                          building_pars_f%pars_xy(ind_hc3,j,i)
     5717                   
     5718                    IF ( building_pars_f%pars_xy(ind_tc1,j,i) /=               &
     5719                         building_pars_f%fill )  THEN
     5720                       surf_usm_v(l)%lambda_h(nzb_wall,m)   =                  &
     5721                                          building_pars_f%pars_xy(ind_tc1,j,i)   
     5722                       surf_usm_v(l)%lambda_h(nzb_wall+1,m) =                  &
     5723                                          building_pars_f%pars_xy(ind_tc1,j,i) 
     5724                    ENDIF
     5725                    IF ( building_pars_f%pars_xy(ind_tc2,j,i) /=               & 
     5726                         building_pars_f%fill )                                &
     5727                       surf_usm_v(l)%lambda_h(nzb_wall+2,m) =                  &
     5728                                          building_pars_f%pars_xy(ind_tc2,j,i)
     5729                       
     5730                    IF ( building_pars_f%pars_xy(ind_tc3,j,i) /=               &
     5731                         building_pars_f%fill )                                &
     5732                       surf_usm_v(l)%lambda_h(nzb_wall+3,m) =                  &
     5733                                          building_pars_f%pars_xy(ind_tc3,j,i)
     5734                       
     5735                    IF ( building_pars_f%pars_xy(ind_tc1,j,i) /=               &
     5736                         building_pars_f%fill )  THEN
     5737                       surf_usm_v(l)%lambda_h_green(nzb_wall,m)   =            &
     5738                                           building_pars_f%pars_xy(ind_tc1,j,i)   
     5739                       surf_usm_v(l)%lambda_h_green(nzb_wall+1,m) =            &
     5740                                           building_pars_f%pars_xy(ind_tc1,j,i) 
     5741                    ENDIF
     5742                    IF ( building_pars_f%pars_xy(ind_tc2,j,i) /=               &
     5743                         building_pars_f%fill )                                &
     5744                       surf_usm_v(l)%lambda_h_green(nzb_wall+2,m) =            &
     5745                                          building_pars_f%pars_xy(ind_tc2,j,i)
     5746                       
     5747                    IF ( building_pars_f%pars_xy(ind_tc3,j,i) /=               &       
     5748                         building_pars_f%fill )                                &
     5749                       surf_usm_v(l)%lambda_h_green(nzb_wall+3,m) =            &
     5750                                          building_pars_f%pars_xy(ind_tc3,j,i)
     5751                       
     5752                    IF ( building_pars_f%pars_xy(ind_tc1_win_r,j,i) /=         &
     5753                         building_pars_f%fill )  THEN
     5754                       surf_usm_v(l)%lambda_h_window(nzb_wall,m)   =           &
     5755                                     building_pars_f%pars_xy(ind_tc1_win_r,j,i)         
     5756                       surf_usm_v(l)%lambda_h_window(nzb_wall+1,m) =           &
     5757                                     building_pars_f%pars_xy(ind_tc1_win_r,j,i)       
     5758                    ENDIF
     5759                    IF ( building_pars_f%pars_xy(ind_tc2,j,i) /=               &     
     5760                         building_pars_f%fill )                                &
     5761                       surf_usm_v(l)%lambda_h_window(nzb_wall+2,m) =           &
     5762                                           building_pars_f%pars_xy(ind_tc2,j,i)
     5763                       
     5764                    IF ( building_pars_f%pars_xy(ind_tc3,j,i) /=               &   
     5765                         building_pars_f%fill )                                &
     5766                       surf_usm_v(l)%lambda_h_window(nzb_wall+3,m) =           &
     5767                                          building_pars_f%pars_xy(ind_tc3,j,i)   
     5768                   
     5769                    IF ( building_pars_f%pars_xy(ind_indoor_target_temp_summer,j,i) /=&           
     5770                         building_pars_f%fill )                                &
     5771                       surf_usm_v(l)%target_temp_summer(m)  =                  &
     5772                            building_pars_f%pars_xy(ind_indoor_target_temp_summer,j,i)   
     5773                    IF ( building_pars_f%pars_xy(ind_indoor_target_temp_winter,j,i) /=&           
     5774                         building_pars_f%fill )                                &
     5775                       surf_usm_v(l)%target_temp_winter(m)  =                  &
     5776                            building_pars_f%pars_xy(ind_indoor_target_temp_winter,j,i)   
     5777                   
     5778                    IF ( building_pars_f%pars_xy(ind_emis_wall,j,i) /=         &   
     5779                         building_pars_f%fill )                                &
     5780                       surf_usm_v(l)%emissivity(ind_veg_wall,m)  =             &
     5781                                      building_pars_f%pars_xy(ind_emis_wall,j,i)
     5782                       
     5783                    IF ( building_pars_f%pars_xy(ind_emis_green,j,i) /=        &           
     5784                         building_pars_f%fill )                                &
     5785                       surf_usm_v(l)%emissivity(ind_pav_green,m) =             &
     5786                                      building_pars_f%pars_xy(ind_emis_green,j,i)
     5787                       
     5788                    IF ( building_pars_f%pars_xy(ind_emis_win,j,i) /=          &
     5789                         building_pars_f%fill )                                &
     5790                       surf_usm_v(l)%emissivity(ind_wat_win,m)   =             &
     5791                                      building_pars_f%pars_xy(ind_emis_win,j,i)
     5792                       
     5793                    IF ( building_pars_f%pars_xy(ind_trans,j,i) /=             &   
     5794                         building_pars_f%fill )                                &
     5795                       surf_usm_v(l)%transmissivity(m) =                       &
     5796                                          building_pars_f%pars_xy(ind_trans,j,i)
     5797                   
     5798                    IF ( building_pars_f%pars_xy(ind_z0,j,i) /=                &         
     5799                         building_pars_f%fill )                                &
     5800                       surf_usm_v(l)%z0(m) = building_pars_f%pars_xy(ind_z0,j,i)
     5801                       
     5802                    IF ( building_pars_f%pars_xy(ind_z0qh,j,i) /=              &           
     5803                         building_pars_f%fill )                                &
     5804                       surf_usm_v(l)%z0h(m) =                                  &
     5805                                       building_pars_f%pars_xy(ind_z0qh,j,i)
     5806                    IF ( building_pars_f%pars_xy(ind_z0qh,j,i) /=              &           
     5807                         building_pars_f%fill )                                &
     5808                       surf_usm_v(l)%z0q(m) =                                  &
     5809                                       building_pars_f%pars_xy(ind_z0qh,j,i)
     5810                   
     5811                    IF ( building_pars_f%pars_xy(ind_alb_wall_agfl,j,i) /=     &         
     5812                         building_pars_f%fill )                                &
     5813                       surf_usm_v(l)%albedo_type(ind_veg_wall,m)  =            &
     5814                                 building_pars_f%pars_xy(ind_alb_wall_agfl,j,i)
     5815                       
     5816                    IF ( building_pars_f%pars_xy(ind_alb_green_agfl,j,i) /=    &           
     5817                         building_pars_f%fill )                                &
     5818                       surf_usm_v(l)%albedo_type(ind_pav_green,m) =            &
     5819                                 building_pars_f%pars_xy(ind_alb_green_agfl,j,i)
     5820                    IF ( building_pars_f%pars_xy(ind_alb_win_agfl,j,i) /=      &         
     5821                         building_pars_f%fill )                                &
     5822                       surf_usm_v(l)%albedo_type(ind_wat_win,m)   =            &
     5823                                   building_pars_f%pars_xy(ind_alb_win_agfl,j,i)
     5824                   
     5825                    IF ( building_pars_f%pars_xy(ind_thick_1_agfl,j,i) /=      &         
     5826                         building_pars_f%fill )                                &
     5827                       surf_usm_v(l)%zw(nzb_wall,m) =                          &
     5828                                   building_pars_f%pars_xy(ind_thick_1_agfl,j,i)
     5829                       
     5830                    IF ( building_pars_f%pars_xy(ind_thick_2_agfl,j,i) /=      &         
     5831                         building_pars_f%fill )                                &
     5832                       surf_usm_v(l)%zw(nzb_wall+1,m) =                        &
     5833                                   building_pars_f%pars_xy(ind_thick_2_agfl,j,i)
     5834                       
     5835                    IF ( building_pars_f%pars_xy(ind_thick_3_agfl,j,i) /=      &         
     5836                         building_pars_f%fill )                                &
     5837                       surf_usm_v(l)%zw(nzb_wall+2,m) =                        &
     5838                                   building_pars_f%pars_xy(ind_thick_3_agfl,j,i)
     5839                       
     5840                       
     5841                    IF ( building_pars_f%pars_xy(ind_thick_4_agfl,j,i) /=      &         
     5842                         building_pars_f%fill )                                &
     5843                       surf_usm_v(l)%zw(nzb_wall+3,m) =                        &
     5844                                   building_pars_f%pars_xy(ind_thick_4_agfl,j,i)
     5845                       
     5846                    IF ( building_pars_f%pars_xy(ind_thick_1_agfl,j,i) /=      &           
     5847                         building_pars_f%fill )                                &
     5848                       surf_usm_v(l)%zw_green(nzb_wall,m) =                    &
     5849                                   building_pars_f%pars_xy(ind_thick_1_agfl,j,i)
     5850                       
     5851                    IF ( building_pars_f%pars_xy(ind_thick_2_agfl,j,i) /=      &         
     5852                         building_pars_f%fill )                                &
     5853                       surf_usm_v(l)%zw_green(nzb_wall+1,m) =                  &
     5854                                   building_pars_f%pars_xy(ind_thick_2_agfl,j,i)
     5855                       
     5856                    IF ( building_pars_f%pars_xy(ind_thick_3_agfl,j,i) /=      &         
     5857                         building_pars_f%fill )                                &
     5858                       surf_usm_v(l)%zw_green(nzb_wall+2,m) =                  &
     5859                                   building_pars_f%pars_xy(ind_thick_3_agfl,j,i)
     5860                       
     5861                    IF ( building_pars_f%pars_xy(ind_thick_4_agfl,j,i) /=      &         
     5862                         building_pars_f%fill )                                &
     5863                       surf_usm_v(l)%zw_green(nzb_wall+3,m) =                  &
     5864                                   building_pars_f%pars_xy(ind_thick_4_agfl,j,i)
     5865                   
     5866                    IF ( building_pars_f%pars_xy(ind_c_surface,j,i) /=         &       
     5867                         building_pars_f%fill )                                &
     5868                       surf_usm_v(l)%c_surface(m) =                            &
     5869                                     building_pars_f%pars_xy(ind_c_surface,j,i)
     5870                       
     5871                    IF ( building_pars_f%pars_xy(ind_lambda_surf,j,i) /=       &       
     5872                         building_pars_f%fill )                                &
     5873                       surf_usm_v(l)%lambda_surf(m) =                          &
     5874                                    building_pars_f%pars_xy(ind_lambda_surf,j,i)
     5875                   
    50645876              ENDDO
    50655877           ENDDO
     
    57736585   
    57746586              ENDDO
     6587              lambda_green_temp(nzt_wall+1) = lambda_green_temp(nzt_wall)
    57756588   
    57766589   
Note: See TracChangeset for help on using the changeset viewer.