|
SURFEX v7.3
General documentation of Surfex
|
|
SURFEX v7.3
General documentation of Surfex
|
00001 !################### 00002 MODULE MODE_CONV_DOE 00003 !################### 00004 ! 00005 !!**** *MODE_CONV_DOE* - 00006 !! 00007 !! PURPOSE 00008 !! ------- 00009 ! 00010 ! 00011 !! 00012 !!** IMPLICIT ARGUMENTS 00013 !! ------------------ 00014 !! NONE 00015 !! 00016 !! REFERENCE 00017 !! --------- 00018 !! 00019 !! 00020 !! AUTHOR 00021 !! ------ 00022 !! G. Pigeon 00023 !! 00024 !! MODIFICATIONS 00025 !! ------------- 00026 !! Original 20/08/12 00027 ! 00028 !-------------------------------------------------------------------------------- 00029 ! 00030 !* 0. DECLARATIONS 00031 ! ------------ 00032 ! 00033 !-------------------------------------------------------------------------------- 00034 USE YOMHOOK ,ONLY : LHOOK, DR_HOOK 00035 USE PARKIND1 ,ONLY : JPRB 00036 INTERFACE CHTC_VERT_DOE 00037 MODULE PROCEDURE CHTC_VERT_DOE_0D 00038 MODULE PROCEDURE CHTC_VERT_DOE_1D 00039 END INTERFACE 00040 INTERFACE CHTC_UP_DOE 00041 MODULE PROCEDURE CHTC_UP_DOE_0D 00042 MODULE PROCEDURE CHTC_UP_DOE_1D 00043 END INTERFACE 00044 INTERFACE CHTC_SMOOTH_LEE_DOE 00045 MODULE PROCEDURE CHTC_SMOOTH_LEE_DOE_0D 00046 MODULE PROCEDURE CHTC_SMOOTH_LEE_DOE_1D 00047 END INTERFACE 00048 INTERFACE CHTC_SMOOTH_WIND_DOE 00049 MODULE PROCEDURE CHTC_SMOOTH_WIND_DOE_0D 00050 MODULE PROCEDURE CHTC_SMOOTH_WIND_DOE_1D 00051 END INTERFACE 00052 INTERFACE CHTC_ROUGH_DOE 00053 MODULE PROCEDURE CHTC_ROUGH_DOE_0D 00054 MODULE PROCEDURE CHTC_ROUGH_DOE_1D 00055 END INTERFACE 00056 ! 00057 CONTAINS 00058 !---------------------------- 00059 !############################################# 00060 FUNCTION CHTC_VERT_DOE_0D(PTS, PTA) RESULT(PCHTC) 00061 !############################################# 00062 ! 00063 !!**** *CHTC_VERT_DOE* - 00064 !! 00065 !! PURPOSE 00066 !! ------- 00067 ! function to compute convective surface coefficient for 00068 ! a vertical surface from surface temperature and air 00069 ! temperature 00070 ! 00071 !!** METHOD 00072 !! ------ 00073 !! 00074 !! EXTERNAL 00075 !! -------- 00076 !! NONE 00077 !! 00078 !! IMPLICIT ARGUMENTS 00079 !! ------------------ 00080 !! 00081 !! REFERENCE 00082 !! --------- 00083 !! EnergyPlus, Engineering Reference, DOE-2 model for convection on outside 00084 !! surfaces, p68 00085 !! 00086 !! AUTHOR 00087 !! ------ 00088 !! G. Pigeon * Meteo France * 00089 !! 00090 !! MODIFICATIONS 00091 !! ------------- 00092 !! Original 20/00/2012 00093 ! 00094 !------------------------------------------------------------------------------- 00095 !* 0. DECLARATIONS 00096 ! ------------ 00097 ! 00098 IMPLICIT NONE 00099 ! 00100 !* 0.1 Declarations of arguments and results 00101 ! 00102 ! 00103 REAL, INTENT(IN) :: PTS ! Surface temperature (Kelvin) 00104 REAL, INTENT(IN) :: PTA ! Air temperature (Kelvin) 00105 REAL :: PCHTC ! Convective heat transfer coefficient [W/(m²K)] 00106 ! 00107 REAL(KIND=JPRB) :: ZHOOK_HANDLE 00108 !------------------------------------------------------------------------------- 00109 IF (LHOOK) CALL DR_HOOK('MODE_CONV_DOE:CHTC_VERT_DOE',0,ZHOOK_HANDLE) 00110 ! 00111 !* 1. COMPUTE THE CHTC 00112 ! ---------------- 00113 ! 00114 PCHTC = 1.31 * (ABS(PTA - PTS))**(1./3.) 00115 !------------------------------------------------------------------------------- 00116 IF (LHOOK) CALL DR_HOOK('MODE_CONV_DOE:CHTC_VERT_DOE',1,ZHOOK_HANDLE) 00117 ! 00118 END FUNCTION CHTC_VERT_DOE_0D 00119 !######################### 00120 !############################################# 00121 FUNCTION CHTC_VERT_DOE_1D(PTS, PTA) RESULT(PCHTC) 00122 !############################################# 00123 ! 00124 !!**** *CHTC_VERT_DOE* - 00125 !! 00126 !! PURPOSE 00127 !! ------- 00128 ! function to compute convective surface coefficient for 00129 ! a vertical surface from surface temperature and air 00130 ! temperature 00131 ! 00132 !!** METHOD 00133 !! ------ 00134 !! 00135 !! EXTERNAL 00136 !! -------- 00137 !! NONE 00138 !! 00139 !! IMPLICIT ARGUMENTS 00140 !! ------------------ 00141 !! 00142 !! REFERENCE 00143 !! --------- 00144 !! EnergyPlus, Engineering Reference, DOE-2 model for convection on outside 00145 !! surfaces, p68 00146 !! 00147 !! AUTHOR 00148 !! ------ 00149 !! G. Pigeon * Meteo France * 00150 !! 00151 !! MODIFICATIONS 00152 !! ------------- 00153 !! Original 20/00/2012 00154 ! 00155 !------------------------------------------------------------------------------- 00156 !* 0. DECLARATIONS 00157 ! ------------ 00158 ! 00159 IMPLICIT NONE 00160 ! 00161 !* 0.1 Declarations of arguments and results 00162 ! 00163 ! 00164 REAL, DIMENSION(:), INTENT(IN) :: PTS ! Surface temperature (Kelvin) 00165 REAL, DIMENSION(:), INTENT(IN) :: PTA ! Air temperature (Kelvin) 00166 REAL, DIMENSION(SIZE(PTS)) :: PCHTC ! Convective heat transfer coefficient [W/(m²K)] 00167 ! 00168 REAL(KIND=JPRB) :: ZHOOK_HANDLE 00169 !------------------------------------------------------------------------------- 00170 IF (LHOOK) CALL DR_HOOK('MODE_CONV_DOE:CHTC_VERT_DOE',0,ZHOOK_HANDLE) 00171 ! 00172 !* 1. COMPUTE THE CHTC 00173 ! ---------------- 00174 ! 00175 PCHTC = 1.31 * (ABS(PTA - PTS))**(1./3.) 00176 !------------------------------------------------------------------------------- 00177 IF (LHOOK) CALL DR_HOOK('MODE_CONV_DOE:CHTC_VERT_DOE',1,ZHOOK_HANDLE) 00178 ! 00179 END FUNCTION CHTC_VERT_DOE_1D 00180 !######################### 00181 ! 00182 !############################################# 00183 FUNCTION CHTC_UP_DOE_1D(PTS, PTA) RESULT(PCHTC) 00184 !############################################# 00185 ! 00186 !!**** *CHTC_UP_DOE* - 00187 !! 00188 !! PURPOSE 00189 !! ------- 00190 ! function to compute convective surface coefficient for 00191 ! an upward surface from surface temperature and air 00192 ! temperature 00193 ! 00194 !!** METHOD 00195 !! ------ 00196 !! 00197 !! EXTERNAL 00198 !! -------- 00199 !! NONE 00200 !! 00201 !! IMPLICIT ARGUMENTS 00202 !! ------------------ 00203 !! 00204 !! REFERENCE 00205 !! --------- 00206 !! EnergyPlus, Engineering Reference, DOE-2 model for convection on outside 00207 !! surfaces, p68 00208 !! 00209 !! AUTHOR 00210 !! ------ 00211 !! G. Pigeon * Meteo France * 00212 !! 00213 !! MODIFICATIONS 00214 !! ------------- 00215 !! Original 20/00/2012 00216 ! 00217 !------------------------------------------------------------------------------- 00218 !* 0. DECLARATIONS 00219 ! ------------ 00220 ! 00221 IMPLICIT NONE 00222 ! 00223 !* 0.1 Declarations of arguments and results 00224 ! 00225 REAL, DIMENSION(:), INTENT(IN) :: PTS ! Surface temperature (Kelvin) 00226 REAL, DIMENSION(:), INTENT(IN) :: PTA ! Air temperature (Kelvin) 00227 REAL, DIMENSION(SIZE(PTS)) :: PCHTC ! Convective heat transfer coefficient [W/(m²K)] 00228 ! 00229 REAL(KIND=JPRB) :: ZHOOK_HANDLE 00230 !------------------------------------------------------------------------------- 00231 IF (LHOOK) CALL DR_HOOK('MODE_CONV_DOE:CHTC_UP_DOE',0,ZHOOK_HANDLE) 00232 ! 00233 !* 1. COMPUTE THE CHTC 00234 ! ---------------- 00235 ! 00236 WHERE(PTA < PTS) 00237 PCHTC = 9.482 * (PTS - PTA)**(1./3.) / (7.283 - 1) 00238 ELSE WHERE 00239 PCHTC = 1.810 * (PTA - PTS)**(1./3.) / (1.382 + 1) 00240 END WHERE 00241 ! 00242 !------------------------------------------------------------------------------- 00243 IF (LHOOK) CALL DR_HOOK('MODE_CONV_DOE:CHTC_UP_DOE',1,ZHOOK_HANDLE) 00244 ! 00245 END FUNCTION CHTC_UP_DOE_1D 00246 !####################### 00247 !############################################# 00248 FUNCTION CHTC_UP_DOE_0D(PTS, PTA) RESULT(PCHTC) 00249 !############################################# 00250 ! 00251 !!**** *CHTC_UP_DOE* - 00252 !! 00253 !! PURPOSE 00254 !! ------- 00255 ! function to compute convective surface coefficient for 00256 ! an upward surface from surface temperature and air 00257 ! temperature 00258 ! 00259 !!** METHOD 00260 !! ------ 00261 !! 00262 !! EXTERNAL 00263 !! -------- 00264 !! NONE 00265 !! 00266 !! IMPLICIT ARGUMENTS 00267 !! ------------------ 00268 !! 00269 !! REFERENCE 00270 !! --------- 00271 !! EnergyPlus, Engineering Reference, DOE-2 model for convection on outside 00272 !! surfaces, p68 00273 !! 00274 !! AUTHOR 00275 !! ------ 00276 !! G. Pigeon * Meteo France * 00277 !! 00278 !! MODIFICATIONS 00279 !! ------------- 00280 !! Original 20/00/2012 00281 ! 00282 !------------------------------------------------------------------------------- 00283 !* 0. DECLARATIONS 00284 ! ------------ 00285 ! 00286 IMPLICIT NONE 00287 ! 00288 !* 0.1 Declarations of arguments and results 00289 ! 00290 REAL, INTENT(IN) :: PTS ! Surface temperature (Kelvin) 00291 REAL, INTENT(IN) :: PTA ! Air temperature (Kelvin) 00292 REAL :: PCHTC ! Convective heat transfer coefficient [W/(m²K)] 00293 ! 00294 REAL(KIND=JPRB) :: ZHOOK_HANDLE 00295 !------------------------------------------------------------------------------- 00296 IF (LHOOK) CALL DR_HOOK('MODE_CONV_DOE:CHTC_UP_DOE',0,ZHOOK_HANDLE) 00297 ! 00298 !* 1. COMPUTE THE CHTC 00299 ! ---------------- 00300 ! 00301 IF (PTA < PTS) THEN 00302 PCHTC = 9.482 * (PTS - PTA)**(1./3.) / (7.283 - 1) 00303 ELSE 00304 PCHTC = 1.810 * (PTA - PTS)**(1./3.) / (1.382 + 1) 00305 ENDIF 00306 ! 00307 !------------------------------------------------------------------------------- 00308 IF (LHOOK) CALL DR_HOOK('MODE_CONV_DOE:CHTC_UP_DOE',1,ZHOOK_HANDLE) 00309 ! 00310 END FUNCTION CHTC_UP_DOE_0D 00311 !####################### 00312 ! 00313 !############################################# 00314 FUNCTION CHTC_DOWN_DOE(PTS, PTA) RESULT(PCHTC) 00315 !############################################# 00316 ! 00317 !!**** *CHTC_DOWN_DOE* - 00318 !! 00319 !! PURPOSE 00320 !! ------- 00321 ! function to compute convective surface coefficient for 00322 ! an downward surface from surface temperature and air 00323 ! temperature 00324 ! 00325 !!** METHOD 00326 !! ------ 00327 !! 00328 !! EXTERNAL 00329 !! -------- 00330 !! NONE 00331 !! 00332 !! IMPLICIT ARGUMENTS 00333 !! ------------------ 00334 !! 00335 !! REFERENCE 00336 !! --------- 00337 !! EnergyPlus, Engineering Reference, DOE-2 model for convection on outside 00338 !! surfaces, p68 00339 !! 00340 !! AUTHOR 00341 !! ------ 00342 !! G. Pigeon * Meteo France * 00343 !! 00344 !! MODIFICATIONS 00345 !! ------------- 00346 !! Original 20/00/2012 00347 ! 00348 !------------------------------------------------------------------------------- 00349 !* 0. DECLARATIONS 00350 ! ------------ 00351 ! 00352 IMPLICIT NONE 00353 ! 00354 !* 0.1 Declarations of arguments and results 00355 ! 00356 REAL, DIMENSION(:), INTENT(IN) :: PTS ! Surface temperature (Kelvin) 00357 REAL, DIMENSION(:), INTENT(IN) :: PTA ! Air temperature (Kelvin) 00358 REAL, DIMENSION(SIZE(PTS)) :: PCHTC ! Convective heat transfer coefficient [W/(m²K)] 00359 ! 00360 REAL(KIND=JPRB) :: ZHOOK_HANDLE 00361 !------------------------------------------------------------------------------- 00362 IF (LHOOK) CALL DR_HOOK('MODE_CONV_DOE:CHTC_DOWN_DOE',0,ZHOOK_HANDLE) 00363 ! 00364 !* 1. COMPUTE THE CHTC 00365 ! ---------------- 00366 ! 00367 WHERE(PTA > PTS) 00368 PCHTC = 9.482 * (PTA - PTS)**(1./3.) / (7.283 - 1) 00369 ELSE WHERE 00370 PCHTC = 1.810 * (PTS - PTA)**(1./3.) / (1.382 + 1) 00371 END WHERE 00372 ! 00373 !------------------------------------------------------------------------------- 00374 IF (LHOOK) CALL DR_HOOK('MODE_CONV_DOE:CHTC_DOWN_DOE',1,ZHOOK_HANDLE) 00375 ! 00376 END FUNCTION CHTC_DOWN_DOE 00377 !######################### 00378 ! 00379 !############################################# 00380 FUNCTION CHTC_SMOOTH_LEE_DOE_0D(PCHTCN, PVMOD) RESULT(PCHTC) 00381 !############################################# 00382 ! 00383 !!**** *CHTC_SMOOTH_LEE_DOE* - 00384 !! 00385 !! PURPOSE 00386 !! ------- 00387 ! function to compute convective surface coefficient for 00388 ! a leeward smooth surface from the natural convection coef and the 00389 ! wind speed 00390 ! 00391 !!** METHOD 00392 !! ------ 00393 !! 00394 !! from EnergyPlus Engineering Reference, average the leeward/windward coef 00395 !! 00396 !! EXTERNAL 00397 !! -------- 00398 !! NONE 00399 !! 00400 !! IMPLICIT ARGUMENTS 00401 !! ------------------ 00402 !! 00403 !! REFERENCE 00404 !! --------- 00405 !! EnergyPlus, Engineering Reference, DOE-2 model for convection on outside 00406 !! surfaces, p68 00407 !! 00408 !! AUTHOR 00409 !! ------ 00410 !! G. Pigeon * Meteo France * 00411 !! 00412 !! MODIFICATIONS 00413 !! ------------- 00414 !! Original 20/00/2012 00415 ! 00416 !------------------------------------------------------------------------------- 00417 !* 0. DECLARATIONS 00418 ! ------------ 00419 ! 00420 IMPLICIT NONE 00421 ! 00422 !* 0.1 Declarations of arguments and results 00423 ! 00424 REAL,INTENT(IN) :: PCHTCN ! Convective heat transfer coefficient for natural conv. [W/(m²K)] 00425 REAL,INTENT(IN) :: PVMOD ! wind speed (m/s) 00426 REAL :: PCHTC ! Convective heat transfer coefficient [W/(m²K)] 00427 ! 00428 !* 0.2 Declarations of local variables 00429 ! 00430 REAL(KIND=JPRB) :: ZHOOK_HANDLE 00431 REAL :: ZALEE, ZBLEE ! coef for leeward facade 00432 !---------------------------------------------------------------------- 00433 IF (LHOOK) CALL DR_HOOK('MODE_CONV_DOE:CHTC_SMOOTH_LEE_DOE',0,ZHOOK_HANDLE) 00434 ! 00435 !* 1. INITIALIZE COEF 00436 ! ---------------- 00437 ! 00438 ZALEE = 2.86 00439 ZBLEE = 0.617 00440 00441 !* 2. COMPUTE THE CHTC 00442 ! ---------------- 00443 ! 00444 PCHTC = SQRT(PCHTCN**2+(ZALEE*PVMOD**ZBLEE)**2) 00445 ! 00446 !------------------------------------------------------------------------------- 00447 IF (LHOOK) CALL DR_HOOK('MODE_CONV_DOE:CHTC_SMOOTH_LEE_DOE',1,ZHOOK_HANDLE) 00448 ! 00449 END FUNCTION CHTC_SMOOTH_LEE_DOE_0D 00450 !######################### 00451 00452 !############################################# 00453 FUNCTION CHTC_SMOOTH_LEE_DOE_1D(PCHTCN, PVMOD) RESULT(PCHTC) 00454 !############################################# 00455 ! 00456 !!**** *CHTC_SMOOTH_LEE_DOE* - 00457 !! 00458 !! PURPOSE 00459 !! ------- 00460 ! function to compute convective surface coefficient for 00461 ! a leeward smooth surface from the natural convection coef and the 00462 ! wind speed 00463 ! 00464 !!** METHOD 00465 !! ------ 00466 !! 00467 !! from EnergyPlus Engineering Reference, average the leeward/windward coef 00468 !! 00469 !! EXTERNAL 00470 !! -------- 00471 !! NONE 00472 !! 00473 !! IMPLICIT ARGUMENTS 00474 !! ------------------ 00475 !! 00476 !! REFERENCE 00477 !! --------- 00478 !! EnergyPlus, Engineering Reference, DOE-2 model for convection on outside 00479 !! surfaces, p68 00480 !! 00481 !! AUTHOR 00482 !! ------ 00483 !! G. Pigeon * Meteo France * 00484 !! 00485 !! MODIFICATIONS 00486 !! ------------- 00487 !! Original 20/00/2012 00488 ! 00489 !------------------------------------------------------------------------------- 00490 !* 0. DECLARATIONS 00491 ! ------------ 00492 ! 00493 IMPLICIT NONE 00494 ! 00495 !* 0.1 Declarations of arguments and results 00496 ! 00497 REAL, DIMENSION(:), INTENT(IN) :: PCHTCN ! Convective heat transfer coefficient for natural conv. [W/(m²K)] 00498 REAL, DIMENSION(:), INTENT(IN) :: PVMOD ! wind speed (m/s) 00499 REAL, DIMENSION(SIZE(PCHTCN)) :: PCHTC ! Convective heat transfer coefficient [W/(m²K)] 00500 ! 00501 !* 0.2 Declarations of local variables 00502 ! 00503 REAL(KIND=JPRB) :: ZHOOK_HANDLE 00504 REAL :: ZALEE, ZBLEE ! coef for windward/leeward facade 00505 !---------------------------------------------------------------------- 00506 IF (LHOOK) CALL DR_HOOK('MODE_CONV_DOE:CHTC_SMOOTH_LEE_DOE',0,ZHOOK_HANDLE) 00507 ! 00508 !* 1. INITIALIZE COEF 00509 ! ---------------- 00510 ! 00511 ZALEE = 2.86 00512 ZBLEE = 0.617 00513 00514 !* 2. COMPUTE THE CHTC 00515 ! ---------------- 00516 ! 00517 PCHTC = SQRT(PCHTCN**2+(ZALEE*PVMOD**ZBLEE)**2) 00518 ! 00519 !------------------------------------------------------------------------------- 00520 IF (LHOOK) CALL DR_HOOK('MODE_CONV_DOE:CHTC_SMOOTH_LEE_DOE',1,ZHOOK_HANDLE) 00521 ! 00522 END FUNCTION CHTC_SMOOTH_LEE_DOE_1D 00523 !######################### 00524 !############################################# 00525 FUNCTION CHTC_SMOOTH_WIND_DOE_0D(PCHTCN, PVMOD) RESULT(PCHTC) 00526 !############################################# 00527 ! 00528 !!**** *CHTC_SMOOTH_WIND_DOE* - 00529 !! 00530 !! PURPOSE 00531 !! ------- 00532 ! function to compute convective surface coefficient for 00533 ! a smooth surface from the natural convection coef and the 00534 ! wind speed 00535 ! 00536 !!** METHOD 00537 !! ------ 00538 !! 00539 !! from EnergyPlus Engineering Reference, average the leeward/windward coef 00540 !! 00541 !! EXTERNAL 00542 !! -------- 00543 !! NONE 00544 !! 00545 !! IMPLICIT ARGUMENTS 00546 !! ------------------ 00547 !! 00548 !! REFERENCE 00549 !! --------- 00550 !! EnergyPlus, Engineering Reference, DOE-2 model for convection on outside 00551 !! surfaces, p68 00552 !! 00553 !! AUTHOR 00554 !! ------ 00555 !! G. Pigeon * Meteo France * 00556 !! 00557 !! MODIFICATIONS 00558 !! ------------- 00559 !! Original 20/00/2012 00560 ! 00561 !------------------------------------------------------------------------------- 00562 !* 0. DECLARATIONS 00563 ! ------------ 00564 ! 00565 IMPLICIT NONE 00566 ! 00567 !* 0.1 Declarations of arguments and results 00568 ! 00569 REAL,INTENT(IN) :: PCHTCN ! Convective heat transfer coefficient for natural conv. [W/(m²K)] 00570 REAL,INTENT(IN) :: PVMOD ! wind speed (m/s) 00571 REAL :: PCHTC ! Convective heat transfer coefficient [W/(m²K)] 00572 ! 00573 !* 0.2 Declarations of local variables 00574 ! 00575 REAL(KIND=JPRB) :: ZHOOK_HANDLE 00576 REAL :: ZAWIND, ZBWIND ! coef for windward/leeward facade 00577 !---------------------------------------------------------------------- 00578 IF (LHOOK) CALL DR_HOOK('MODE_CONV_DOE:CHTC_SMOOTH_WIND_DOE',0,ZHOOK_HANDLE) 00579 ! 00580 !* 1. INITIALIZE COEF 00581 ! ---------------- 00582 ! 00583 ZAWIND = 2.38 00584 ZBWIND = 0.89 00585 00586 !* 2. COMPUTE THE CHTC 00587 ! ---------------- 00588 ! 00589 PCHTC = SQRT(PCHTCN**2+(ZAWIND*PVMOD**ZBWIND)**2) 00590 ! 00591 !------------------------------------------------------------------------------- 00592 IF (LHOOK) CALL DR_HOOK('MODE_CONV_DOE:CHTC_SMOOTH_WIND_DOE',1,ZHOOK_HANDLE) 00593 ! 00594 END FUNCTION CHTC_SMOOTH_WIND_DOE_0D 00595 !######################### 00596 00597 !############################################# 00598 FUNCTION CHTC_SMOOTH_WIND_DOE_1D(PCHTCN, PVMOD) RESULT(PCHTC) 00599 !############################################# 00600 ! 00601 !!**** *CHTC_SMOOTH_WIND_DOE* - 00602 !! 00603 !! PURPOSE 00604 !! ------- 00605 ! function to compute convective surface coefficient for 00606 ! a windward smooth surface from the natural convection coef and the 00607 ! wind speed 00608 ! 00609 !!** METHOD 00610 !! ------ 00611 !! 00612 !! from EnergyPlus Engineering Reference, average the leeward/windward coef 00613 !! 00614 !! EXTERNAL 00615 !! -------- 00616 !! NONE 00617 !! 00618 !! IMPLICIT ARGUMENTS 00619 !! ------------------ 00620 !! 00621 !! REFERENCE 00622 !! --------- 00623 !! EnergyPlus, Engineering Reference, DOE-2 model for convection on outside 00624 !! surfaces, p68 00625 !! 00626 !! AUTHOR 00627 !! ------ 00628 !! G. Pigeon * Meteo France * 00629 !! 00630 !! MODIFICATIONS 00631 !! ------------- 00632 !! Original 20/00/2012 00633 ! 00634 !------------------------------------------------------------------------------- 00635 !* 0. DECLARATIONS 00636 ! ------------ 00637 ! 00638 IMPLICIT NONE 00639 ! 00640 !* 0.1 Declarations of arguments and results 00641 ! 00642 REAL, DIMENSION(:), INTENT(IN) :: PCHTCN ! Convective heat transfer coefficient for natural conv. [W/(m²K)] 00643 REAL, DIMENSION(:), INTENT(IN) :: PVMOD ! wind speed (m/s) 00644 REAL, DIMENSION(SIZE(PCHTCN)) :: PCHTC ! Convective heat transfer coefficient [W/(m²K)] 00645 ! 00646 !* 0.2 Declarations of local variables 00647 ! 00648 REAL(KIND=JPRB) :: ZHOOK_HANDLE 00649 REAL :: ZAWIND, ZBWIND ! coef for windward facade 00650 !---------------------------------------------------------------------- 00651 IF (LHOOK) CALL DR_HOOK('MODE_CONV_DOE:CHTC_SMOOTH_WIND_DOE',0,ZHOOK_HANDLE) 00652 ! 00653 !* 1. INITIALIZE COEF 00654 ! ---------------- 00655 ! 00656 ZAWIND = 2.38 00657 ZBWIND = 0.89 00658 00659 !* 2. COMPUTE THE CHTC 00660 ! ---------------- 00661 ! 00662 PCHTC = SQRT(PCHTCN**2+(ZAWIND*PVMOD**ZBWIND)**2) 00663 ! 00664 !------------------------------------------------------------------------------- 00665 IF (LHOOK) CALL DR_HOOK('MODE_CONV_DOE:CHTC_SMOOTH_WIND_DOE',1,ZHOOK_HANDLE) 00666 ! 00667 END FUNCTION CHTC_SMOOTH_WIND_DOE_1D 00668 !######################### 00669 !############################################# 00670 FUNCTION CHTC_ROUGH_DOE_0D(PCHTCN, PCHTCS, PROUGH) RESULT(PCHTC) 00671 !############################################# 00672 ! 00673 !!**** *CHTC_ROUGH_DOE* - 00674 !! 00675 !! PURPOSE 00676 !! ------- 00677 ! function to compute convective surface coefficient for 00678 ! a rough surface from the natural convection coef, the 00679 ! smooth surface convective coef and the roughness coef 00680 ! 00681 !!** METHOD 00682 !! ------ 00683 !! 00684 !! 00685 !! EXTERNAL 00686 !! -------- 00687 !! NONE 00688 !! 00689 !! IMPLICIT ARGUMENTS 00690 !! ------------------ 00691 !! 00692 !! REFERENCE 00693 !! --------- 00694 !! EnergyPlus, Engineering Reference, DOE-2 model for convection on outside 00695 !! surfaces, p68 00696 !! 00697 !! AUTHOR 00698 !! ------ 00699 !! G. Pigeon * Meteo France * 00700 !! 00701 !! MODIFICATIONS 00702 !! ------------- 00703 !! Original 20/00/2012 00704 ! 00705 !------------------------------------------------------------------------------- 00706 !* 0. DECLARATIONS 00707 ! ------------ 00708 ! 00709 IMPLICIT NONE 00710 ! 00711 !* 0.1 Declarations of arguments and results 00712 ! 00713 REAL, INTENT(IN) :: PCHTCN ! Convective heat transfer coefficient for natural conv. [W/(m²K)] 00714 REAL, INTENT(IN) :: PCHTCS ! Convective heat transfer coefficient for a smooth surface [W/(m²K)] 00715 REAL, INTENT(IN) :: PROUGH ! Convective heat transfer coefficient for a smooth surface [W/(m²K)] 00716 REAL :: PCHTC ! Convective heat transfer coefficient [W/(m²K)] 00717 ! 00718 !* 0.2 Declarations of local variables 00719 ! 00720 REAL(KIND=JPRB) :: ZHOOK_HANDLE 00721 !---------------------------------------------------------------------- 00722 IF (LHOOK) CALL DR_HOOK('MODE_CONV_DOE:CHTC_ROUGH_DOE',0,ZHOOK_HANDLE) 00723 ! 00724 !* 1. COMPUTE THE CHTC 00725 ! ---------------- 00726 ! 00727 PCHTC = PCHTCN + PROUGH * (PCHTCS - PCHTCN) 00728 ! 00729 !------------------------------------------------------------------------------- 00730 IF (LHOOK) CALL DR_HOOK('MODE_CONV_DOE:CHTC_ROUGH_DOE',1,ZHOOK_HANDLE) 00731 ! 00732 END FUNCTION CHTC_ROUGH_DOE_0D 00733 !######################### 00734 !############################################# 00735 FUNCTION CHTC_ROUGH_DOE_1D(PCHTCN, PCHTCS, PROUGH) RESULT(PCHTC) 00736 !############################################# 00737 ! 00738 !!**** *CHTC_ROUGH_DOE* - 00739 !! 00740 !! PURPOSE 00741 !! ------- 00742 ! function to compute convective surface coefficient for 00743 ! a rough surface from the natural convection coef, the 00744 ! smooth surface convective coef and the roughness coef 00745 ! 00746 !!** METHOD 00747 !! ------ 00748 !! 00749 !! 00750 !! EXTERNAL 00751 !! -------- 00752 !! NONE 00753 !! 00754 !! IMPLICIT ARGUMENTS 00755 !! ------------------ 00756 !! 00757 !! REFERENCE 00758 !! --------- 00759 !! EnergyPlus, Engineering Reference, DOE-2 model for convection on outside 00760 !! surfaces, p68 00761 !! 00762 !! AUTHOR 00763 !! ------ 00764 !! G. Pigeon * Meteo France * 00765 !! 00766 !! MODIFICATIONS 00767 !! ------------- 00768 !! Original 20/00/2012 00769 ! 00770 !------------------------------------------------------------------------------- 00771 !* 0. DECLARATIONS 00772 ! ------------ 00773 ! 00774 IMPLICIT NONE 00775 ! 00776 !* 0.1 Declarations of arguments and results 00777 ! 00778 REAL, DIMENSION(:), INTENT(IN) :: PCHTCN ! Convective heat transfer coefficient for natural conv. [W/(m²K)] 00779 REAL, DIMENSION(:), INTENT(IN) :: PCHTCS ! Convective heat transfer coefficient for a smooth surface [W/(m²K)] 00780 REAL, DIMENSION(:), INTENT(IN) :: PROUGH ! Convective heat transfer coefficient for a smooth surface [W/(m²K)] 00781 REAL, DIMENSION(SIZE(PCHTCN)) :: PCHTC ! Convective heat transfer coefficient [W/(m²K)] 00782 ! 00783 !* 0.2 Declarations of local variables 00784 ! 00785 REAL(KIND=JPRB) :: ZHOOK_HANDLE 00786 !---------------------------------------------------------------------- 00787 IF (LHOOK) CALL DR_HOOK('MODE_CONV_DOE:CHTC_ROUGH_DOE',0,ZHOOK_HANDLE) 00788 ! 00789 !* 1. COMPUTE THE CHTC 00790 ! ---------------- 00791 ! 00792 PCHTC = PCHTCN + PROUGH * (PCHTCS - PCHTCN) 00793 ! 00794 !------------------------------------------------------------------------------- 00795 IF (LHOOK) CALL DR_HOOK('MODE_CONV_DOE:CHTC_ROUGH_DOE',1,ZHOOK_HANDLE) 00796 ! 00797 END FUNCTION CHTC_ROUGH_DOE_1D 00798 00799 END MODULE MODE_CONV_DOE
1.8.0
