Changeset 3593
- Timestamp:
- Dec 3, 2018 1:51:13 PM (5 years ago)
- Location:
- palm/trunk/SOURCE
- Files:
-
- 5 edited
-
biometeorology_mod.f90 (modified) (32 diffs)
-
indoor_model_mod.f90 (modified) (6 diffs)
-
multi_agent_system_mod.f90 (modified) (4 diffs)
-
surface_mod.f90 (modified) (2 diffs)
-
wind_turbine_model_mod.f90 (modified) (4 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/biometeorology_mod.f90
r3582 r3593 27 27 ! ----------------- 28 28 ! $Id$ 29 ! Bugfix: additional tmrt_grid allocation in case bio_mrt not selected as ouput, 30 ! replace degree symbol by degree_C 31 ! 32 ! 3582 2018-11-29 19:16:36Z suehring 29 33 ! Consistently use bio_fill_value everywhere, 30 34 ! move allocation and initialization of output variables to bio_check_data_output … … 128 132 ! 129 133 !-- Declare all global variables within the module (alphabetical order) 130 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: tmrt_grid !< tmrt results ( °C)131 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: perct !< PT results ( °C)132 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: utci !< UTCI results ( °C)133 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: pet !< PET results ( °C)134 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: tmrt_grid !< tmrt results (degree_C) 135 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: perct !< PT results (degree_C) 136 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: utci !< UTCI results (degree_C) 137 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: pet !< PET results (degree_C) 134 138 ! 135 139 !-- Grids for averaged thermal indices 136 140 REAL(wp), DIMENSION(:), ALLOCATABLE :: mrt_av_grid !< time average mean 137 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: perct_av !< PT results (aver. input) ( °C)138 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: utci_av !< UTCI results (aver. input) ( °C)139 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: pet_av !< PET results (aver. input) ( °C)141 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: perct_av !< PT results (aver. input) (degree_C) 142 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: utci_av !< UTCI results (aver. input) (degree_C) 143 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: pet_av !< PET results (aver. input) (degree_C) 140 144 141 145 … … 636 640 ! 637 641 !-- Allocate a temporary array with the desired output dimensions. 638 CASE ( 'bio_mrt' )642 CASE ( 'bio_mrt' ) 639 643 unit = 'degree_C' 640 644 IF ( .NOT. ALLOCATED( tmrt_grid ) ) THEN … … 1185 1189 1186 1190 ! 1187 !-- Calculate biometeorology MRT from local radiation 1188 ! fluxes calculated by RTM and assign into 2D grid 1189 tmrt_grid = bio_fill_value 1191 !-- Calculate biometeorology MRT from local radiation fluxes calculated by RTM and assign 1192 !-- into 2D grid. Depending on selected output quantities, tmrt_grid might not have been 1193 !-- allocated in bio_check_data_output yet. 1194 IF ( .NOT. ALLOCATED( tmrt_grid ) ) THEN 1195 ALLOCATE( tmrt_grid (nys:nyn,nxl:nxr) ) 1196 ENDIF 1197 tmrt_grid = REAL( bio_fill_value, KIND = wp ) 1198 1190 1199 DO l = 1, nmrtbl 1191 1200 i = mrtbl(ix,l) … … 1224 1233 ! 1225 1234 !-- Output parameters 1226 REAL(wp), INTENT ( OUT ) :: tmrt !< Mean radiant temperature ( °C)1227 REAL(wp), INTENT ( OUT ) :: ta !< Air temperature ( °C)1235 REAL(wp), INTENT ( OUT ) :: tmrt !< Mean radiant temperature (degree_C) 1236 REAL(wp), INTENT ( OUT ) :: ta !< Air temperature (degree_C) 1228 1237 REAL(wp), INTENT ( OUT ) :: vp !< Vapour pressure (hPa) 1229 1238 REAL(wp), INTENT ( OUT ) :: ws !< Wind speed (local level) (m/s) … … 1322 1331 1323 1332 REAL(wp) :: clo !< Clothing index (no dimension) 1324 REAL(wp) :: ta !< Air temperature ( °C)1333 REAL(wp) :: ta !< Air temperature (degree_C) 1325 1334 REAL(wp) :: vp !< Vapour pressure (hPa) 1326 1335 REAL(wp) :: ws !< Wind speed (local level) (m/s) 1327 1336 REAL(wp) :: pair !< Air pressure (hPa) 1328 REAL(wp) :: perct_ij !< Perceived temperature ( °C)1329 REAL(wp) :: utci_ij !< Universal thermal climate index ( °C)1330 REAL(wp) :: pet_ij !< Physiologically equivalent temperature ( °C)1331 REAL(wp) :: tmrt_ij !< Mean radiant temperature ( °C)1337 REAL(wp) :: perct_ij !< Perceived temperature (degree_C) 1338 REAL(wp) :: utci_ij !< Universal thermal climate index (degree_C) 1339 REAL(wp) :: pet_ij !< Physiologically equivalent temperature (degree_C) 1340 REAL(wp) :: tmrt_ij !< Mean radiant temperature (degree_C) 1332 1341 1333 1342 ! … … 1423 1432 ! 1424 1433 !-- Input parameters 1425 REAL(wp), INTENT ( IN ) :: ta !< Air temperature ( °C)1434 REAL(wp), INTENT ( IN ) :: ta !< Air temperature (degree_C) 1426 1435 REAL(wp), INTENT ( IN ) :: vp !< Vapour pressure (hPa) 1427 1436 REAL(wp), INTENT ( IN ) :: ws !< Wind speed (local level) (m/s) 1428 1437 REAL(wp), INTENT ( IN ) :: pair !< Air pressure (hPa) 1429 REAL(wp), INTENT ( IN ) :: tmrt !< Mean radiant temperature ( °C)1438 REAL(wp), INTENT ( IN ) :: tmrt !< Mean radiant temperature (degree_C) 1430 1439 REAL(wp), INTENT ( IN ) :: dt !< Time past since last calculation (s) 1431 1440 REAL(wp), INTENT ( IN ) :: age !< Age of agent (y) … … 1438 1447 !-- Both, input and output parameters 1439 1448 Real(wp), INTENT ( INOUT ) :: energy_storage !< Energy storage (W/m²) 1440 Real(wp), INTENT ( INOUT ) :: t_clo !< Clothing temperature ( °C)1449 Real(wp), INTENT ( INOUT ) :: t_clo !< Clothing temperature (degree_C) 1441 1450 Real(wp), INTENT ( INOUT ) :: clo !< Current clothing in sulation 1442 1451 Real(wp), INTENT ( INOUT ) :: actlev !< Individuals activity level … … 1444 1453 ! 1445 1454 !-- Output parameters 1446 REAL(wp), INTENT ( OUT ) :: ipt !< Instationary perceived temp. ( °C)1455 REAL(wp), INTENT ( OUT ) :: ipt !< Instationary perceived temp. (degree_C) 1447 1456 ! 1448 1457 !-- If clo equals the initial value, this is the initial call … … 1484 1493 ! 1485 1494 !-- Type of input of the argument list 1486 REAL(WP), INTENT ( IN ) :: ta_in !< Local air temperature ( °C)1495 REAL(WP), INTENT ( IN ) :: ta_in !< Local air temperature (degree_C) 1487 1496 REAL(WP), INTENT ( IN ) :: vp !< Loacl vapour pressure (hPa) 1488 1497 REAL(WP), INTENT ( IN ) :: ws_hag !< Incident wind speed (m/s) 1489 REAL(WP), INTENT ( IN ) :: tmrt !< Local mean radiant temperature ( °C)1498 REAL(WP), INTENT ( IN ) :: tmrt !< Local mean radiant temperature (degree_C) 1490 1499 REAL(WP), INTENT ( IN ) :: hag !< Height of wind speed input (m) 1491 1500 ! 1492 1501 !-- Type of output of the argument list 1493 REAL(wp), INTENT ( OUT ) :: utci_ij !< Universal Thermal Climate Index ( °C)1494 1495 REAL(WP) :: ta !< air temperature modified by offset ( °C)1502 REAL(wp), INTENT ( OUT ) :: utci_ij !< Universal Thermal Climate Index (degree_C) 1503 1504 REAL(WP) :: ta !< air temperature modified by offset (degree_C) 1496 1505 REAL(WP) :: pa !< air pressure in kPa (kPa) 1497 REAL(WP) :: d_tmrt !< delta-tmrt ( °C)1506 REAL(WP) :: d_tmrt !< delta-tmrt (degree_C) 1498 1507 REAL(WP) :: va !< wind speed at 10 m above ground level (m/s) 1499 REAL(WP) :: offset !< utci deviation by ta cond. exceeded ( °C)1508 REAL(WP) :: offset !< utci deviation by ta cond. exceeded (degree_C) 1500 1509 REAL(WP) :: part_ta !< Air temperature related part of the regression 1501 1510 REAL(WP) :: ta2 !< 2 times ta … … 1823 1832 ! Description: 1824 1833 ! ------------ 1825 !> calculate_perct_static: Estimation of perceived temperature (PT, deg C)1834 !> calculate_perct_static: Estimation of perceived temperature (PT, degree_C) 1826 1835 !> Value of perct is the Perceived Temperature, degree centigrade 1827 1836 !------------------------------------------------------------------------------! … … 1860 1869 REAL(wp) :: pmv_s !< Fangers predicted mean vote for summer clothing 1861 1870 REAL(wp) :: pmva !< adjusted predicted mean vote 1862 REAL(wp) :: ptc !< perceived temp. for cold conditions ( °C)1871 REAL(wp) :: ptc !< perceived temp. for cold conditions (degree_C) 1863 1872 REAL(wp) :: d_std !< factor to threshold for sultriness 1864 1873 REAL(wp) :: pmvs !< pred. mean vote considering sultrieness 1865 REAL(wp) :: top !< Gagge's operative temperatures ( °C)1874 REAL(wp) :: top !< Gagge's operative temperatures (degree_C) 1866 1875 1867 1876 INTEGER(iwp) :: ncount !< running index … … 2244 2253 REAL(wp) :: heat_convection !< energy loss by autocnvection (W) 2245 2254 REAL(wp) :: activity !< persons activity (must stay == actlev, W) 2246 REAL(wp) :: t_skin_aver !< average skin temperature ( °C)2255 REAL(wp) :: t_skin_aver !< average skin temperature (degree_C) 2247 2256 REAL(wp) :: bc !< preliminary result storage 2248 2257 REAL(wp) :: cc !< preliminary result storage … … 2250 2259 REAL(wp) :: ec !< preliminary result storage 2251 2260 REAL(wp) :: gc !< preliminary result storage 2252 REAL(wp) :: t_clothing !< clothing temperature ( °C)2261 REAL(wp) :: t_clothing !< clothing temperature (degree_C) 2253 2262 REAL(wp) :: hr !< radiational heat resistence 2254 2263 REAL(wp) :: clo !< clothing insulation index (clo) … … 2892 2901 REAL(wp), INTENT(in) :: height !< Persons height (m) 2893 2902 REAL(wp), INTENT(in) :: work !< Current workload (W) 2894 REAL(wp), INTENT(in) :: ta !< Air Temperature ( °C)2903 REAL(wp), INTENT(in) :: ta !< Air Temperature (degree_C) 2895 2904 REAL(wp), INTENT(in) :: vp !< Vapor pressure (hPa) 2896 2905 REAL(wp), INTENT(in) :: ws !< Wind speed in approx. 1.1m (m/s) 2897 REAL(wp), INTENT(in) :: tmrt !< Mean radiant temperature ( °C)2906 REAL(wp), INTENT(in) :: tmrt !< Mean radiant temperature (degree_C) 2898 2907 REAL(wp), INTENT(in) :: pair !< Air pressure (hPa) 2899 2908 REAL(wp), INTENT(in) :: dt !< Timestep (s) … … 3084 3093 !> SUBROUTINE ipt_cycle 3085 3094 !> Calculates one timestep for the instationary version of perceived 3086 !> temperature (iPT, °C) for3095 !> temperature (iPT, degree_C) for 3087 3096 !> - standard measured/predicted meteorological values and TMRT 3088 3097 !> as input; … … 3099 3108 ! 3100 3109 !-- Type of input of the argument list 3101 REAL(wp), INTENT ( IN ) :: ta !< Air temperature ( °C)3110 REAL(wp), INTENT ( IN ) :: ta !< Air temperature (degree_C) 3102 3111 REAL(wp), INTENT ( IN ) :: vp !< Vapor pressure (hPa) 3103 REAL(wp), INTENT ( IN ) :: tmrt !< Mean radiant temperature ( °C)3112 REAL(wp), INTENT ( IN ) :: tmrt !< Mean radiant temperature (degree_C) 3104 3113 REAL(wp), INTENT ( IN ) :: ws !< Wind speed (m/s) 3105 3114 REAL(wp), INTENT ( IN ) :: pair !< Air pressure (hPa) … … 3111 3120 !-- In and output parameters 3112 3121 REAL(wp), INTENT (INOUT) :: storage !< Heat storage (W) 3113 REAL(wp), INTENT (INOUT) :: t_clothing !< Clothig temperature ( °C)3122 REAL(wp), INTENT (INOUT) :: t_clothing !< Clothig temperature (degree_C) 3114 3123 ! 3115 3124 !-- Type of output of the argument list 3116 REAL(wp), INTENT ( OUT ) :: ipt !< Instationary perceived temperature ( °C)3125 REAL(wp), INTENT ( OUT ) :: ipt !< Instationary perceived temperature (degree_C) 3117 3126 ! 3118 3127 !-- Type of internal variables … … 3191 3200 ! 3192 3201 !-- Input argument types 3193 REAL(wp), INTENT ( IN ) :: ta !< Air temperature ( °C)3194 REAL(wp), INTENT ( IN ) :: tmrt !< Mean radiant temperature ( °C)3202 REAL(wp), INTENT ( IN ) :: ta !< Air temperature (degree_C) 3203 REAL(wp), INTENT ( IN ) :: tmrt !< Mean radiant temperature (degree_C) 3195 3204 REAL(wp), INTENT ( IN ) :: pa !< Vapour pressure (hPa) 3196 3205 REAL(wp), INTENT ( IN ) :: pair !< Air pressure (hPa) … … 3205 3214 3206 3215 REAL(wp), INTENT (INOUT) :: s !< storage var. of energy balance (W/m2) 3207 REAL(wp), INTENT (INOUT) :: t_cloth !< clothing temperature ( °C)3216 REAL(wp), INTENT (INOUT) :: t_cloth !< clothing temperature (degree_C) 3208 3217 ! 3209 3218 !-- Internal variables … … 3212 3221 REAL(wp) :: f_cl !< Increase in surface due to clothing (factor) 3213 3222 REAL(wp) :: heat_convection !< energy loss by autocnvection (W) 3214 REAL(wp) :: t_skin_aver !< average skin temperature ( °C)3223 REAL(wp) :: t_skin_aver !< average skin temperature (degree_C) 3215 3224 REAL(wp) :: bc !< preliminary result storage 3216 3225 REAL(wp) :: cc !< preliminary result storage … … 3218 3227 REAL(wp) :: ec !< preliminary result storage 3219 3228 REAL(wp) :: gc !< preliminary result storage 3220 REAL(wp) :: t_clothing !< clothing temperature ( °C)3229 REAL(wp) :: t_clothing !< clothing temperature (degree_C) 3221 3230 ! REAL(wp) :: hr !< radiational heat resistence 3222 3231 REAL(wp) :: clo !< clothing insulation index (clo) … … 3338 3347 ! 3339 3348 !-- Input arguments: 3340 REAL(wp), INTENT( IN ) :: ta !< Air temperature ( °C)3341 REAL(wp), INTENT( IN ) :: tmrt !< Mean radiant temperature ( °C)3349 REAL(wp), INTENT( IN ) :: ta !< Air temperature (degree_C) 3350 REAL(wp), INTENT( IN ) :: tmrt !< Mean radiant temperature (degree_C) 3342 3351 REAL(wp), INTENT( IN ) :: v !< Wind speed (m/s) 3343 3352 REAL(wp), INTENT( IN ) :: vpa !< Vapor pressure (hPa) … … 3345 3354 ! 3346 3355 !-- Output arguments: 3347 REAL(wp), INTENT ( OUT ) :: pet_ij !< PET ( °C)3356 REAL(wp), INTENT ( OUT ) :: pet_ij !< PET (degree_C) 3348 3357 ! 3349 3358 !-- Internal variables: … … 3360 3369 REAL(wp) :: rtv 3361 3370 REAL(wp) :: vpts !< Sat. vapor pressure over skin (hPa) 3362 REAL(wp) :: tsk !< Skin temperature ( °C)3363 REAL(wp) :: tcl !< Clothing temperature ( °C)3371 REAL(wp) :: tsk !< Skin temperature (degree_C) 3372 REAL(wp) :: tcl !< Clothing temperature (degree_C) 3364 3373 REAL(wp) :: wetsk !< Fraction of wet skin (factor) 3365 3374 ! … … 3412 3421 !-- Input arguments: 3413 3422 REAL(wp), INTENT( IN ) :: pair !< air pressure (hPa) 3414 REAL(wp), INTENT( IN ) :: ta !< air temperature ( °C)3423 REAL(wp), INTENT( IN ) :: ta !< air temperature (degree_C) 3415 3424 REAL(wp), INTENT( IN ) :: vpa !< vapor pressure (hPa) 3416 3425 REAL(wp), INTENT( IN ) :: age !< Persons age (a) … … 3474 3483 REAL(wp), INTENT( IN ) :: int_heat !< internal heat production (W) 3475 3484 REAL(wp), INTENT( IN ) :: pair !< Air pressure (hPa) 3476 REAL(wp), INTENT( IN ) :: ta !< Air temperature ( °C)3477 REAL(wp), INTENT( IN ) :: tmrt !< Mean radiant temperature ( °C)3485 REAL(wp), INTENT( IN ) :: ta !< Air temperature (degree_C) 3486 REAL(wp), INTENT( IN ) :: tmrt !< Mean radiant temperature (degree_C) 3478 3487 REAL(wp), INTENT( IN ) :: v !< Wind speed (m/s) 3479 3488 REAL(wp), INTENT( IN ) :: vpa !< Vapor pressure (hPa) … … 3492 3501 REAL(wp), INTENT( OUT ) :: rdcl !< Diffusion resistence of clothing (factor) 3493 3502 REAL(wp), INTENT( OUT ) :: rdsk !< Diffusion resistence of skin (factor) 3494 REAL(wp), INTENT( OUT ) :: tcl !< Clothing temperature ( °C)3495 REAL(wp), INTENT( OUT ) :: tsk !< Skin temperature ( °C)3503 REAL(wp), INTENT( OUT ) :: tcl !< Clothing temperature (degree_C) 3504 REAL(wp), INTENT( OUT ) :: tsk !< Skin temperature (degree_C) 3496 3505 REAL(wp), INTENT( OUT ) :: vpts !< Sat. vapor pressure over skin (hPa) 3497 3506 REAL(wp), INTENT( OUT ) :: wetsk !< Fraction of wet skin (dimensionless) … … 3508 3517 ! 3509 3518 !-- Internal variables 3510 REAL(wp) :: c(0:10) !< Core temperature array ( °C)3519 REAL(wp) :: c(0:10) !< Core temperature array (degree_C) 3511 3520 REAL(wp) :: cbare !< Convection through bare skin 3512 3521 REAL(wp) :: cclo !< Convection through clothing … … 3753 3762 REAL(wp), INTENT( IN ) :: rdsk !< diffusion resistence of skin (factor) 3754 3763 REAL(wp), INTENT( IN ) :: rtv !< respiratory volume 3755 REAL(wp), INTENT( IN ) :: ta !< air temperature ( °C)3756 REAL(wp), INTENT( IN ) :: tcl !< clothing temperature ( °C)3757 REAL(wp), INTENT( IN ) :: tsk !< skin temperature ( °C)3764 REAL(wp), INTENT( IN ) :: ta !< air temperature (degree_C) 3765 REAL(wp), INTENT( IN ) :: tcl !< clothing temperature (degree_C) 3766 REAL(wp), INTENT( IN ) :: tsk !< skin temperature (degree_C) 3758 3767 REAL(wp), INTENT( IN ) :: vpts !< sat. vapor pressure over skin (hPa) 3759 3768 REAL(wp), INTENT( IN ) :: wetsk !< fraction of wet skin (dimensionless) … … 3761 3770 !-- Output arguments: 3762 3771 REAL(wp), INTENT( OUT ) :: aeff !< effective surface area (m²) 3763 REAL(wp), INTENT( OUT ) :: pet_ij !< PET ( °C)3772 REAL(wp), INTENT( OUT ) :: pet_ij !< PET (degree_C) 3764 3773 ! 3765 3774 !-- Cconstants: … … 3782 3791 REAL ( wp ) :: rclo !< Radiational loss of clothing (W/m²) 3783 3792 REAL ( wp ) :: rsum !< Radiational loss or gain (W/m²) 3784 REAL ( wp ) :: tex !< Temperat. of exhaled air ( °C)3793 REAL ( wp ) :: tex !< Temperat. of exhaled air (degree_C) 3785 3794 REAL ( wp ) :: vpex !< Vapor pressure of exhaled air (hPa) 3786 3795 REAL ( wp ) :: xx !< Delta PET per iteration (K) -
palm/trunk/SOURCE/indoor_model_mod.f90
r3524 r3593 26 26 ! ----------------- 27 27 ! $Id$ 28 ! Replace degree symbol by degree_C 29 ! 30 ! 3524 2018-11-14 13:36:44Z raasch 28 31 ! working precision added to make code Fortran 2008 conform 29 32 ! … … 182 185 REAL(wp) :: schedule_d !< activation for internal loads (low or high + low) 183 186 REAL(wp) :: skip_time_do_indoor = 0.0_wp !< [s] Indoor model is not called before this time 184 REAL(wp) :: theta_air !<! [ °C] air temperature of the RC-node185 REAL(wp) :: theta_air_0 !<! [ oC] air temperature of the RC-node in equilibrium186 REAL(wp) :: theta_air_10 !<! [ oC] air temperature of the RC-node from a heating capacity of 10 W/m²187 REAL(wp) :: theta_air_ac !< [ oC] actual room temperature after heating/cooling188 REAL(wp) :: theta_air_set !< [ oC] Setpoint_temperature for the room189 REAL(wp) :: theta_int_c_set !< [ oC] Max. Setpoint temperature (summer)190 REAL(wp) :: theta_int_h_set !< [ oC] Max. Setpoint temperature (winter)191 REAL(wp) :: theta_m !<! [ oC} inner temperature of the RC-node192 REAL(wp) :: theta_m_t !<! [ oC] (Fictive) component temperature timestep193 REAL(wp) :: theta_m_t_prev !< [ oC] (Fictive) component temperature previous timestep (do not change)194 REAL(wp) :: theta_op !< [ oC] operative temperature195 REAL(wp) :: theta_s !<! [ oC] surface temperature of the RC-node187 REAL(wp) :: theta_air !<! [degree_C] air temperature of the RC-node 188 REAL(wp) :: theta_air_0 !<! [degree_C] air temperature of the RC-node in equilibrium 189 REAL(wp) :: theta_air_10 !<! [degree_C] air temperature of the RC-node from a heating capacity of 10 W/m² 190 REAL(wp) :: theta_air_ac !< [degree_C] actual room temperature after heating/cooling 191 REAL(wp) :: theta_air_set !< [degree_C] Setpoint_temperature for the room 192 REAL(wp) :: theta_int_c_set !< [degree_C] Max. Setpoint temperature (summer) 193 REAL(wp) :: theta_int_h_set !< [degree_C] Max. Setpoint temperature (winter) 194 REAL(wp) :: theta_m !<! [degree_C} inner temperature of the RC-node 195 REAL(wp) :: theta_m_t !<! [degree_C] (Fictive) component temperature timestep 196 REAL(wp) :: theta_m_t_prev !< [degree_C] (Fictive) component temperature previous timestep (do not change) 197 REAL(wp) :: theta_op !< [degree_C] operative temperature 198 REAL(wp) :: theta_s !<! [degree_C] surface temperature of the RC-node 196 199 REAL(wp) :: time_indoor = 0.0_wp !< [s] time since last call of indoor model 197 200 REAL(wp) :: time_utc_hour !< Time in hours per day (UTC) … … 279 282 + near_facade_temperature ) & 280 283 ) / h_tr_2 & 281 ) !< [ oC] Eq. (C.5)284 ) !< [degree_C] Eq. (C.5) 282 285 283 286 !< Calculation of component temperature at factual timestep … … 286 289 ) & 287 290 / ( ( c_m / 3600 ) + 0.5 * ( h_tr_3 + h_tr_em ) ) & 288 ) !< [ oC] Eq. (C.4)291 ) !< [degree_C] Eq. (C.4) 289 292 290 293 !< Calculation of mean inner temperature for the RC-node in actual timestep 291 theta_m = ( theta_m_t + theta_m_t_prev ) * 0.5 !< [ oC] Eq. (C.9)294 theta_m = ( theta_m_t + theta_m_t_prev ) * 0.5 !< [degree_C] Eq. (C.9) 292 295 293 296 !< Calculation of mean surface temperature of the RC-node in actual timestep … … 296 299 ) & 297 300 / ( h_tr_ms + h_tr_w + h_tr_1 ) & 298 ) !< [ oC] Eq. (C.10)301 ) !< [degree_C] Eq. (C.10) 299 302 300 303 !< Calculation of the air temperature of the RC-node 301 304 theta_air = ( h_tr_is * theta_s + h_ve * near_facade_temperature & 302 + phi_ia + phi_hc_nd_dummy ) / ( h_tr_is + h_ve ) !< [ oC] Eq. (C.11)305 + phi_ia + phi_hc_nd_dummy ) / ( h_tr_is + h_ve ) !< [degree_C] Eq. (C.11) 303 306 304 307 END SUBROUTINE im_calc_temperatures … … 1100 1103 !-- Calculate the operating temperature with weighted mean temperature of air and mean solar temperature 1101 1104 !-- Will be used for thermal comfort calculations 1102 theta_op = 0.3 * theta_air_ac + 0.7 * theta_s !< [ °C] operative Temperature Eq. (C.12)1105 theta_op = 0.3 * theta_air_ac + 0.7 * theta_s !< [degree_C] operative Temperature Eq. (C.12) 1103 1106 ! 1104 1107 !-- Heat flux into the wall. Value needed in urban_surface_mod to -
palm/trunk/SOURCE/multi_agent_system_mod.f90
r3587 r3593 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Replace degree symbol by degree_C/degrees 28 ! 29 ! 3587 2018-11-30 13:52:19Z sward 27 30 ! Added output of agent substep time 28 31 ! … … 237 240 REAL(wp) :: speed_x !< speed of agent in x 238 241 REAL(wp) :: speed_y !< speed of agent in y 239 REAL(wp) :: ipt !< instationary thermal index iPT ( °C)242 REAL(wp) :: ipt !< instationary thermal index iPT (degree_C) 240 243 REAL(wp) :: windspeed !< absolute value of windspeed at agent position 241 244 REAL(wp) :: x !< x-position … … 363 366 INTEGER(iwp) :: a !< agent iterator 364 367 !-- local meteorological conditions 365 REAL(wp) :: tmrt !< mean radiant temperature ( °C)366 REAL(wp) :: ta !< air temperature ( °C)368 REAL(wp) :: tmrt !< mean radiant temperature (degree_C) 369 REAL(wp) :: ta !< air temperature (degree_C) 367 370 REAL(wp) :: vp !< vapour pressure (hPa) 368 371 REAL(wp) :: v !< wind speed (local level) (m/s) … … 4574 4577 4575 4578 ! !-- forces that are located outside of a sight radius of 4576 ! !-- 200 °(-> COS(100./180.*pi) = COS(.555*pi)) of4579 ! !-- 200 degrees (-> COS(100./180.*pi) = COS(.555*pi)) of 4577 4580 ! !-- current agent are considered to have an effect of 50% 4578 4581 ! IF ( force_d_x * agents(nl)%speed_e_x + & -
palm/trunk/SOURCE/surface_mod.f90
r3572 r3593 26 26 ! ----------------- 27 27 ! $Id$ 28 ! Replace degree symbol by 'degrees' 29 ! 30 ! 3572 2018-11-28 11:40:28Z suehring 28 31 ! Define short- and longwave radiation flux arrays (e.g. diffuse, direct, 29 32 ! reflected, resedual) for all surfaces. This is required to surface outputs … … 376 379 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: frac !< relative surface fraction (LSM: vegetation, water, pavement; USM: wall, green, window) 377 380 378 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: aldif !< albedo for longwave diffusive radiation, solar angle of 60 °379 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: aldir !< albedo for longwave direct radiation, solar angle of 60 °380 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: asdif !< albedo for shortwave diffusive radiation, solar angle of 60 °381 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: asdir !< albedo for shortwave direct radiation, solar angle of 60 °382 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: rrtm_aldif !< albedo for longwave diffusive radiation, solar angle of 60 °383 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: rrtm_aldir !< albedo for longwave direct radiation, solar angle of 60 °384 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: rrtm_asdif !< albedo for shortwave diffusive radiation, solar angle of 60 °385 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: rrtm_asdir !< albedo for shortwave direct radiation, solar angle of 60 °381 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: aldif !< albedo for longwave diffusive radiation, solar angle of 60 degrees 382 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: aldir !< albedo for longwave direct radiation, solar angle of 60 degrees 383 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: asdif !< albedo for shortwave diffusive radiation, solar angle of 60 degrees 384 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: asdir !< albedo for shortwave direct radiation, solar angle of 60 degrees 385 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: rrtm_aldif !< albedo for longwave diffusive radiation, solar angle of 60 degrees 386 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: rrtm_aldir !< albedo for longwave direct radiation, solar angle of 60 degrees 387 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: rrtm_asdif !< albedo for shortwave diffusive radiation, solar angle of 60 degrees 388 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: rrtm_asdir !< albedo for shortwave direct radiation, solar angle of 60 degrees 386 389 387 390 REAL(wp), DIMENSION(:), ALLOCATABLE :: q_surface !< skin-surface mixing ratio -
palm/trunk/SOURCE/wind_turbine_model_mod.f90
r3274 r3593 26 26 ! ----------------- 27 27 ! $Id$ 28 ! Replace degree symbol by 'degrees' 29 ! 30 ! 3274 2018-09-24 15:42:55Z knoop 28 31 ! Modularization of all bulk cloud physics code components 29 32 ! … … 1536 1539 turb_cl_sel2 = 0.0_wp 1537 1540 1538 turb_cd_tab(1,iir) = 0.0_wp ! For -180 °(iialpha=1) the values1541 turb_cd_tab(1,iir) = 0.0_wp ! For -180 degrees (iialpha=1) the values 1539 1542 turb_cl_tab(1,iir) = 0.0_wp ! for each radius has to be set 1540 1543 ! explicitly … … 1545 1548 turb_cl_sel2 = turb_cl_table(:,t2) 1546 1549 ! 1547 !-- For -180 °(iialpha=1) the values for each radius has to be set1550 !-- For -180 degrees (iialpha=1) the values for each radius has to be set 1548 1551 !-- explicitly 1549 1552 turb_cd_tab(1,iir) = ( weight_a * turb_cd_table(1,t1) + weight_b & … … 2544 2547 !-- The yaw controller computes a 30s running mean of the wind direction. 2545 2548 !-- If the difference between turbine alignment and wind direction exceeds 2546 !-- 5 °, the turbine is yawed. The mechanism stops as soon as the 2s-running2547 !-- mean of the missalignment is smaller than 0.5 °.2549 !-- 5 degrees, the turbine is yawed. The mechanism stops as soon as the 2s-running 2550 !-- mean of the missalignment is smaller than 0.5 degrees. 2548 2551 !-- Attention: If the timestep during the simulation changes significantly 2549 2552 !-- the lengths of the running means change and it does not correspond to
Note: See TracChangeset
for help on using the changeset viewer.
