FUNCTION CHTC_UP_DOE_1D(PTS, PTA) RESULT(PCHTC)
!#############################################
!
!!**** *CHTC_UP_DOE* -
!!
!! PURPOSE
!! -------
! function to compute convective surface coefficient for
! an upward surface from surface temperature and air
! temperature
!
!!** METHOD
!! ------
!!
!! EXTERNAL
!! --------
!! NONE
!!
!! IMPLICIT ARGUMENTS
!! ------------------
!!
!! REFERENCE
!! ---------
!! EnergyPlus, Engineering Reference, DOE-2 model for convection on outside
!! surfaces, p68
!!
!! AUTHOR
!! ------
!! G. Pigeon * Meteo France *
!!
!! MODIFICATIONS
!! -------------
!! Original 20/00/2012
!
!-------------------------------------------------------------------------------
!* 0. DECLARATIONS
! ------------
!
IMPLICIT NONE
!
!* 0.1 Declarations of arguments and results
!
REAL, DIMENSION(:), INTENT(IN) :: PTS ! Surface temperature (Kelvin)
REAL, DIMENSION(:), INTENT(IN) :: PTA ! Air temperature (Kelvin)
REAL, DIMENSION(SIZE(PTS)) :: PCHTC ! Convective heat transfer coefficient [W/(m2.K)]
!
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!-------------------------------------------------------------------------------
IF (LHOOK) CALL DR_HOOK('MODE_CONV_DOE:CHTC_UP_DOE',0,ZHOOK_HANDLE)
!
!* 1. COMPUTE THE CHTC
! ----------------
!
WHERE(PTA < PTS)
PCHTC = 9.482 * (PTS - PTA)**(1./3.) / (7.283 - 1)
ELSE WHERE
PCHTC = 1.810 * (PTA - PTS)**(1./3.) / (1.382 + 1)
END WHERE
!
!-------------------------------------------------------------------------------
IF (LHOOK) CALL DR_HOOK('MODE_CONV_DOE:CHTC_UP_DOE',1,ZHOOK_HANDLE)
!
END FUNCTION CHTC_UP_DOE_1D