Changeset 3593 for palm/trunk/SOURCE/biometeorology_mod.f90
- Timestamp:
- Dec 3, 2018 1:51:13 PM (5 years ago)
- File:
-
- 1 edited
-
palm/trunk/SOURCE/biometeorology_mod.f90 (modified) (32 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)
Note: See TracChangeset
for help on using the changeset viewer.
