Changeset 3881 for palm/trunk/SOURCE
- Timestamp:
- Apr 10, 2019 9:31:22 AM (6 years ago)
- Location:
- palm/trunk/SOURCE
- Files:
-
- 6 edited
-
land_surface_model_mod.f90 (modified) (3 diffs)
-
prognostic_equations.f90 (modified) (2 diffs)
-
radiation_model_mod.f90 (modified) (9 diffs)
-
surface_data_output_mod.f90 (modified) (3 diffs)
-
surface_layer_fluxes_mod.f90 (modified) (2 diffs)
-
urban_surface_mod.f90 (modified) (25 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/land_surface_model_mod.f90
r3868 r3881 25 25 ! ----------------- 26 26 ! $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 27 31 ! More strict limitation of roughness length when it is in the order of the 28 32 ! vertical grid spacing … … 4186 4190 pavement_pars_f%fill ) & 4187 4191 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) ) 4189 4193 IF ( pavement_pars_f%pars_xy(ind_p_emis,j,i) /= & 4190 4194 pavement_pars_f%fill ) & 4191 4195 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) 4193 4197 ENDIF 4194 4198 … … 4246 4250 pavement_pars_f%fill ) & 4247 4251 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) ) 4249 4253 4250 4254 IF ( pavement_pars_f%pars_xy(ind_p_emis,j,i) /= & 4251 4255 pavement_pars_f%fill ) & 4252 4256 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) 4254 4258 ENDIF 4255 4259 ENDDO -
palm/trunk/SOURCE/prognostic_equations.f90
r3880 r3881 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Bugfix in OpenMP directive 28 ! 29 ! 3880 2019-04-08 21:43:02Z knoop 27 30 ! Moved wtm_tendencies to module_interface_actions 28 31 ! … … 589 592 ! 590 593 !-- 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) 593 595 594 596 !$ tn = omp_get_thread_num() -
palm/trunk/SOURCE/radiation_model_mod.f90
r3861 r3881 28 28 ! ----------------- 29 29 ! $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 30 34 ! Bugfix: factor of 4.0 required instead of 3.0 in calculation of rad_lw_out_change_0 31 35 ! … … 1397 1401 unit = 'W/m2' 1398 1402 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 1400 1405 IF ( .NOT. radiation ) THEN 1401 1406 message_string = 'output of "' // TRIM( var ) // '" require'& … … 2443 2448 DO ind_type = 0, 2 2444 2449 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 )& 2446 2451 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) 2448 2453 IF ( albedo_pars_f%pars_xy(1,j,i) /= albedo_pars_f%fill )& 2449 2454 surf_lsm_h%aldir(ind_type,m) = & 2450 2455 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) 2451 2459 IF ( albedo_pars_f%pars_xy(2,j,i) /= albedo_pars_f%fill )& 2452 surf_lsm_h%a ldif(ind_type,m) = &2460 surf_lsm_h%asdir(ind_type,m) = & 2453 2461 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 )& 2458 2463 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) 2460 2465 ENDIF 2461 2466 ENDDO … … 2469 2474 j = surf_usm_h%j(m) 2470 2475 ! 2471 !-- Spectralalbedos for wall/green/window surfaces2476 !-- Broadband albedos for wall/green/window surfaces 2472 2477 DO ind_type = 0, 2 2473 2478 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 )& 2475 2480 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) 2489 2482 ENDIF 2490 2483 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 2491 2526 2492 2527 ENDDO … … 2505 2540 DO ind_type = 0, 2 2506 2541 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) /= & 2508 2543 albedo_pars_f%fill ) & 2509 2544 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) 2511 2546 IF ( albedo_pars_f%pars_xy(1,j+joff,i+ioff) /= & 2512 2547 albedo_pars_f%fill ) & 2513 2548 surf_lsm_v(l)%aldir(ind_type,m) = & 2514 2549 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) 2515 2554 IF ( albedo_pars_f%pars_xy(2,j+joff,i+ioff) /= & 2516 2555 albedo_pars_f%fill ) & 2517 surf_lsm_v(l)%a ldif(ind_type,m) = &2556 surf_lsm_v(l)%asdir(ind_type,m) = & 2518 2557 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) /= & 2524 2559 albedo_pars_f%fill ) & 2525 2560 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) 2527 2562 ENDIF 2528 2563 ENDDO … … 2539 2574 j = surf_usm_v(l)%j(m) 2540 2575 ! 2541 !-- Spectralalbedos for wall/green/window surfaces2576 !-- Broadband albedos for wall/green/window surfaces 2542 2577 DO ind_type = 0, 2 2543 2578 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) 2564 2583 ENDIF 2565 2584 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 2567 2633 ENDDO 2568 2634 ENDIF … … 9809 9875 var(1:7) == 'rtm_svf' .OR. var(1:7) == 'rtm_dif' .OR. & 9810 9876 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. & 9811 9878 var == 'rtm_mrt' .OR. var == 'rtm_mrt_sw' .OR. var == 'rtm_mrt_lw' ) THEN 9812 9879 … … 10511 10578 ENDIF 10512 10579 10513 !-- block of RTM output variables10514 !-- variables are intended mainly for debugging and detailed analyse purposes10515 CASE ( 'rtm_skyvf' )10516 !-- sky view factor10517 DO isurf = dirstart(ids), dirend(ids)10518 IF ( surfl(id,isurf) == idsint_u .OR. surfl(id,isurf) == idsint_l ) THEN10519 local_pf(surfl(ix,isurf),surfl(iy,isurf),surfl(iz,isurf)) = skyvf(isurf)10520 ENDIF10521 ENDDO10522 10523 CASE ( 'rtm_skyvft' )10524 !-- sky view factor10525 DO isurf = dirstart(ids), dirend(ids)10526 IF ( surfl(id,isurf) == ids ) THEN10527 local_pf(surfl(ix,isurf),surfl(iy,isurf),surfl(iz,isurf)) = skyvft(isurf)10528 ENDIF10529 ENDDO10530 10531 CASE ( 'rtm_svf', 'rtm_dif' )10532 !-- shape view factors or iradiance factors to selected surface10533 IF ( TRIM(var)=='rtm_svf' ) THEN10534 k = 110535 ELSE10536 k = 210537 ENDIF10538 DO isvf = 1, nsvfl10539 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 ) ) THEN10544 !-- correct source surface10545 local_pf(surfl(ix,isurflt),surfl(iy,isurflt),surfl(iz,isurflt)) = svf(k,isvf)10546 ENDIF10547 ENDDO10548 10549 10580 CASE ( 'rtm_rad_net' ) 10550 10581 !-- array of complete radiation balance … … 10785 10816 ENDIF 10786 10817 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. 10790 10879 10791 10880 END SELECT -
palm/trunk/SOURCE/surface_data_output_mod.f90
r3817 r3881 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Check for zero output timestep (not allowed in parallel NetCDF output mode) 28 ! 29 ! 3817 2019-03-26 13:53:57Z suehring 27 30 ! Correct output coordinates of vertical surface elements 28 31 ! … … 125 128 USE netcdf_interface, & 126 129 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 128 132 129 133 USE pegrid … … 4316 4320 #endif 4317 4321 ! 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 ! 4318 4344 !-- In case of restart runs, check it the number of cores has been changed. 4319 4345 !-- With surface output this is not allowed. -
palm/trunk/SOURCE/surface_layer_fluxes_mod.f90
r3834 r3881 26 26 ! ----------------- 27 27 ! $Id$ 28 ! Assure that Obukhov length does not become zero 29 ! 30 ! 3834 2019-03-28 15:40:15Z forkel 28 31 ! added USE chem_gasphase_mod 29 32 ! … … 1188 1191 CYCLE 1189 1192 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) ) 1190 1197 1191 1198 ENDDO -
palm/trunk/SOURCE/urban_surface_mod.f90
r3832 r3881 28 28 ! ----------------- 29 29 ! $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 30 35 ! instrumented with openmp directives 31 36 ! … … 564 569 REAL(wp) :: roughness_concrete = 0.001_wp !< roughness length of average concrete surface 565 570 ! 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, 574 599 !< 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, 591 603 !< 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, 595 608 !< 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, 597 614 !< ground floor level 598 INTEGER(iwp) :: ind_hc 1_win_r = 102 !< index in input list for heat capacity at firstwindow layer, roof599 INTEGER(iwp) :: ind_hc 2_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, 600 617 !< above ground floor level 601 INTEGER(iwp) :: ind_hc 2_gfl = 26 !< index in input list for heat capacity at second wall layer, ground floor level602 INTEGER(iwp) :: ind_hc 2_wall_r = 90 !< index in input list for heat capacity at second wall layer, roof603 INTEGER(iwp) :: ind_hc 2_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, 604 621 !< above ground floor level 605 INTEGER(iwp) :: ind_hc 2_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, 606 623 !< 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, 609 635 !< 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, 613 640 !< above ground floor level 614 INTEGER(iwp) :: ind_ hc3_win_gfl = 40 !< index in input list for heat capacity at thirdwindow layer,641 INTEGER(iwp) :: ind_tc1_win_gfl = 74 !< index in input list for thermal conductivity at first window layer, 615 642 !< 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, 623 645 !< above ground floor level 624 INTEGER(iwp) :: ind_tc 1_gfl = 28 !< index in input list for thermal conductivity at firstwall layer,646 INTEGER(iwp) :: ind_tc2_gfl = 30 !< index in input list for thermal conductivity at second wall layer, 625 647 !< ground floor level 626 INTEGER(iwp) :: ind_tc 1_wall_r = 92 !< index in input list for thermal conductivity at firstwall layer, roof627 INTEGER(iwp) :: ind_tc 1_win_agfl = 74 !< index in input list for thermal conductivity at firstwindow 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, 628 650 !< above ground floor level 629 INTEGER(iwp) :: ind_tc 1_win_gfl = 41 !< index in input list for thermal conductivity at firstwindow layer,651 INTEGER(iwp) :: ind_tc2_win_gfl = 75 !< index in input list for thermal conductivity at second window layer, 630 652 !< 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, 633 656 !< above ground floor level 634 INTEGER(iwp) :: ind_tc 2_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, 635 658 !< ground floor level 636 INTEGER(iwp) :: ind_tc 2_wall_r = 93 !< index in input list for thermal conductivity at second wall layer, roof637 INTEGER(iwp) :: ind_tc 2_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, 638 661 !< above ground floor level 639 INTEGER(iwp) :: ind_tc 2_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, 640 663 !< 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 691 703 692 704 … … 705 717 'bridges ' & !< type 7 706 718 /) 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) 782 723 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/) ) 981 1664 982 1665 ! … … 2653 3336 CHARACTER(LEN=2) :: ls 2654 3337 CHARACTER(LEN=varnamelength) :: var !< TRIM(variable) 2655 INTEGER(iwp), PARAMETER :: nl1 = 1 6!< number of directional usm variables3338 INTEGER(iwp), PARAMETER :: nl1 = 14 !< number of directional usm variables 2656 3339 CHARACTER(LEN=varnamelength), DIMENSION(nl1) :: varlist1 = & !< list of directional usm variables 2657 3340 (/'usm_wshf ', & … … 2664 3347 'usm_surfwintrans ', & 2665 3348 'usm_surfcat ', & 2666 'usm_surfalb ', &2667 'usm_surfemis ', &2668 3349 'usm_t_surf_wall ', & 2669 3350 'usm_t_surf_window ', & … … 2734 3415 var(1:14) == 'usm_theta_10cm' ) THEN 2735 3416 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 2738 3418 unit = '1' 2739 3419 ELSE … … 2851 3531 LOGICAL, INTENT(OUT) :: found !< 2852 3532 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 procedure3533 REAL(sp), DIMENSION(nzb:nzt+1,nys:nyn,nxl:nxr) :: temp_pf !< temp array for urban surface output procedure 2854 3534 2855 3535 CHARACTER (len=varnamelength) :: var !< trimmed variable name … … 2956 3636 ENDIF 2957 3637 2958 CASE ( 'usm_surfalb' )2959 !2960 !-- surface albedo, weighted average2961 IF ( idsint == iup_u ) THEN2962 DO m = 1, surf_usm_h%ns2963 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 ENDDO2973 ELSE2974 l = idsidx2975 DO m = 1, surf_usm_v(l)%ns2976 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 ENDDO2986 ENDIF2987 2988 CASE ( 'usm_surfemis' )2989 !2990 !-- surface emissivity, weighted average2991 IF ( idsint == iup_u ) THEN2992 DO m = 1, surf_usm_h%ns2993 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 ENDDO3003 ELSE3004 l = idsidx3005 DO m = 1, surf_usm_v(l)%ns3006 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 ENDDO3016 ENDIF3017 3018 3638 CASE ( 'usm_surfwintrans' ) 3019 3639 ! … … 3759 4379 var(1:15) == 'usm_theta_10cm' .OR. & 3760 4380 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. &3762 4381 var(1:16) == 'usm_surfwintrans' .OR. var(1:7) == 'usm_swc' ) THEN 3763 4382 … … 4058 4677 !-- ground-floor-level properties. 4059 4678 surf_usm_h%ground_level = .FALSE. 4060 ! DO m = 1, surf_usm_h%ns4061 ! 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_level4068 ! IF ( building_pars_f%from_file ) THEN4069 ! 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 ! ENDIF4073 ! !4074 ! !-- Determine height of surface element above ground level4075 ! IF ( terrain_height_f%from_file ) THEN4076 ! z_agl = zw(k) - terrain_height_f%var(j,i)4077 ! ELSE4078 ! z_agl = zw(k)4079 ! ENDIF4080 ! !4081 ! !-- Set flag for ground level4082 ! IF ( z_agl <= ground_floor_level_l ) &4083 ! surf_usm_h%ground_level(m) = .TRUE.4084 ! ENDDO4085 4679 4086 4680 DO l = 0, 3 … … 4091 4685 k = surf_usm_v(l)%k(m) 4092 4686 ! 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. 4095 4690 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 4096 4697 IF ( building_pars_f%from_file ) THEN 4097 4698 IF ( building_pars_f%pars_xy(ind_gflh,j,i) /= & 4098 building_pars_f%fill ) &4699 building_pars_f%fill ) & 4099 4700 ground_floor_level_l = building_pars_f%pars_xy(ind_gflh,j,i) 4100 4701 ENDIF … … 4102 4703 !-- Determine height of surface element above ground level. Please 4103 4704 !-- 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 4106 4709 z_agl = zw(k) - terrain_height_f%var(j-surf_usm_v(l)%joff, & 4107 4710 i-surf_usm_v(l)%ioff) … … 4149 4752 ! 4150 4753 !-- 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. 4152 4755 !-- Level 1 - initialization via default attributes 4153 4756 DO m = 1, surf_usm_h%ns … … 4204 4807 surf_usm_h%lambda_h_window(nzb_wall+3,m) = building_pars(ind_tc3_win_r,building_type) 4205 4808 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) 4208 4811 ! 4209 4812 !-- emissivity of wall-, green- and window fraction … … 4238 4841 surf_usm_h%zw_window(nzb_wall+3,m) = building_pars(ind_thick_4_win_r,building_type) 4239 4842 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) 4246 4849 4247 4850 surf_usm_h%green_type_roof(m) = building_pars(ind_green_type_roof,building_type) … … 4367 4970 surf_usm_v(l)%lambda_h_window(nzb_wall+3,m) = building_pars(ind_tc3_win,building_type) 4368 4971 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) 4371 4974 ! 4372 4975 !-- emissivity of wall-, green- and window fraction … … 4400 5003 surf_usm_v(l)%zw_window(nzb_wall+3,m) = building_pars(ind_thick_4_win,building_type) 4401 5004 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) 4408 5011 4409 5012 ENDDO … … 4470 5073 surf_usm_h%lambda_h_window(nzb_wall+3,m) = building_pars(ind_tc3_win_r,st) 4471 5074 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) 4474 5077 ! 4475 5078 !-- emissivity of wall-, green- and window fraction … … 4504 5107 surf_usm_h%zw_window(nzb_wall+3,m) = building_pars(ind_thick_4_win_r,st) 4505 5108 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) 4512 5115 4513 5116 surf_usm_h%green_type_roof(m) = building_pars(ind_green_type_roof,st) … … 4640 5243 surf_usm_v(l)%lambda_h_window(nzb_wall+3,m) = building_pars(ind_tc3_win,st) 4641 5244 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) 4644 5247 ! 4645 5248 !-- emissivity of wall-, green- and window fraction … … 4673 5276 surf_usm_v(l)%zw_window(nzb_wall+3,m) = building_pars(ind_thick_4_win,st) 4674 5277 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) 4681 5284 4682 5285 … … 4687 5290 4688 5291 ! 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. 4690 5295 IF ( building_pars_f%from_file ) THEN 4691 5296 DO m = 1, surf_usm_h%ns … … 4696 5301 !-- In order to distinguish between ground floor level and 4697 5302 !-- 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, & 4699 5305 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, & 4701 5308 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, & 4703 5311 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, & 4705 5314 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) ) 4706 5351 4707 5352 ! 4708 5353 !-- 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) 4716 5368 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) 4723 5379 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) 4731 5398 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) 4739 5415 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) 4748 5432 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) 4756 5449 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) 4764 5466 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 4839 5582 ENDDO 4840 5583 … … 4847 5590 4848 5591 ! 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, & 4920 5596 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 5064 5876 ENDDO 5065 5877 ENDDO … … 5773 6585 5774 6586 ENDDO 6587 lambda_green_temp(nzt_wall+1) = lambda_green_temp(nzt_wall) 5775 6588 5776 6589
Note: See TracChangeset
for help on using the changeset viewer.
