Changeset 3502
- Timestamp:
- Nov 7, 2018 2:45:23 PM (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/urban_surface_mod.f90
r3469 r3502 28 28 ! ----------------- 29 29 ! $Id$ 30 ! Disable initialization of building roofs with ground-floor-level properties, 31 ! since this causes strong oscillations of surface temperature during the 32 ! spinup. 33 ! 34 ! 3469 2018-10-30 20:05:07Z kanani 30 35 ! Add missing PUBLIC variables for new indoor model 31 36 ! … … 1016 1021 REAL(wp) :: window_inner_temperature = 295.0_wp !< temperature of the inner window surface (~22 degrees C) (K) 1017 1022 1018 1023 REAL(wp) :: m_total = 0.0_wp !< weighted total water content of the soil (m3/m3) 1019 1024 INTEGER(iwp) :: soil_type 1020 1025 … … 4942 4947 !-- use terrain height array from file, if available. This flag is later used 4943 4948 !-- to control initialization of surface attributes. 4949 !-- Todo: for the moment disable initialization of building roofs with 4950 !-- ground-floor-level properties. 4944 4951 surf_usm_h%ground_level = .FALSE. 4945 DO m = 1, surf_usm_h%ns4946 i = surf_usm_h%i(m)4947 j = surf_usm_h%j(m)4948 k = surf_usm_h%k(m)4949 ! 4950 ! -- Get local ground level. If no ground level is given in input file,4951 ! -- use default value.4952 ground_floor_level_l = ground_floor_level4953 IF ( building_pars_f%from_file ) THEN4954 IF ( building_pars_f%pars_xy(ind_gflh,j,i) /= &4955 building_pars_f%fill ) &4956 ground_floor_level_l = building_pars_f%pars_xy(ind_gflh,j,i)4957 ENDIF4958 ! 4959 ! -- Determine height of surface element above ground level4960 IF ( terrain_height_f%from_file ) THEN4961 z_agl = zw(k) - terrain_height_f%var(j,i)4962 ELSE4963 z_agl = zw(k)4964 ENDIF4965 ! 4966 ! -- Set flag for ground level4967 IF ( z_agl <= ground_floor_level_l ) &4968 surf_usm_h%ground_level(m) = .TRUE.4969 ENDDO4952 ! DO m = 1, surf_usm_h%ns 4953 ! i = surf_usm_h%i(m) 4954 ! j = surf_usm_h%j(m) 4955 ! k = surf_usm_h%k(m) 4956 ! ! 4957 ! !-- Get local ground level. If no ground level is given in input file, 4958 ! !-- use default value. 4959 ! ground_floor_level_l = ground_floor_level 4960 ! IF ( building_pars_f%from_file ) THEN 4961 ! IF ( building_pars_f%pars_xy(ind_gflh,j,i) /= & 4962 ! building_pars_f%fill ) & 4963 ! ground_floor_level_l = building_pars_f%pars_xy(ind_gflh,j,i) 4964 ! ENDIF 4965 ! ! 4966 ! !-- Determine height of surface element above ground level 4967 ! IF ( terrain_height_f%from_file ) THEN 4968 ! z_agl = zw(k) - terrain_height_f%var(j,i) 4969 ! ELSE 4970 ! z_agl = zw(k) 4971 ! ENDIF 4972 ! ! 4973 ! !-- Set flag for ground level 4974 ! IF ( z_agl <= ground_floor_level_l ) & 4975 ! surf_usm_h%ground_level(m) = .TRUE. 4976 ! ENDDO 4970 4977 4971 4978 DO l = 0, 3
Note: See TracChangeset
for help on using the changeset viewer.