mode_conv_DOE.F90

Go to the documentation of this file.

!SFX_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
!SFX_LIC This is part of the SURFEX software governed by the CeCILL-C licence
!SFX_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
!SFX_LIC for details. version 1.
!###################
MODULE mode_conv_doe
!###################
!
!!****  *MODE_CONV_DOE* -
!!
!!    PURPOSE
!!    -------
!      
!
!!
!!**  IMPLICIT ARGUMENTS
!!    ------------------
!!       NONE          
!!
!!    REFERENCE
!!    ---------
!!
!!
!!    AUTHOR
!!    ------
!!    G. Pigeon
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    20/08/12 
!
!--------------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!              ------------
!
!--------------------------------------------------------------------------------
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
INTERFACE chtc_vert_doe
        MODULE PROCEDURE chtc_vert_doe_0d
        MODULE PROCEDURE chtc_vert_doe_1d
END INTERFACE
INTERFACE chtc_up_doe
        MODULE PROCEDURE chtc_up_doe_0d
        MODULE PROCEDURE chtc_up_doe_1d
END INTERFACE
INTERFACE chtc_smooth_lee_doe
        MODULE PROCEDURE chtc_smooth_lee_doe_0d
        MODULE PROCEDURE chtc_smooth_lee_doe_1d
END INTERFACE
INTERFACE chtc_smooth_wind_doe
        MODULE PROCEDURE chtc_smooth_wind_doe_0d
        MODULE PROCEDURE chtc_smooth_wind_doe_1d
END INTERFACE
INTERFACE chtc_rough_doe
        MODULE PROCEDURE chtc_rough_doe_0d
        MODULE PROCEDURE chtc_rough_doe_1d
END INTERFACE
!
CONTAINS
!----------------------------
!#############################################
FUNCTION chtc_vert_doe_0d(PTS, PTA) RESULT(PCHTC)
!#############################################
!
!!****  *CHTC_VERT_DOE* - 
!!
!!    PURPOSE
!!    -------
!     function to compute convective surface coefficient for
!     a vertical 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, INTENT(IN)                :: PTS     ! Surface temperature (Kelvin)
REAL, INTENT(IN)                :: PTA     ! Air temperature (Kelvin)
REAL                            :: PCHTC   ! Convective heat transfer coefficient [W/(m2.K)]
!
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!-------------------------------------------------------------------------------
IF (lhook) CALL dr_hook('MODE_CONV_DOE:CHTC_VERT_DOE',0,zhook_handle)
!
!*       1.    COMPUTE THE CHTC
!              ----------------
!
pchtc = 1.31 * (abs(pta - pts))**(1./3.) 
!-------------------------------------------------------------------------------
IF (lhook) CALL dr_hook('MODE_CONV_DOE:CHTC_VERT_DOE',1,zhook_handle)
!
END FUNCTION chtc_vert_doe_0d
!#########################
!#############################################
FUNCTION chtc_vert_doe_1d(PTS, PTA) RESULT(PCHTC)
!#############################################
!
!!****  *CHTC_VERT_DOE* - 
!!
!!    PURPOSE
!!    -------
!     function to compute convective surface coefficient for
!     a vertical 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_VERT_DOE',0,zhook_handle)
!
!*       1.    COMPUTE THE CHTC
!              ----------------
!
pchtc = 1.31 * (abs(pta - pts))**(1./3.) 
!-------------------------------------------------------------------------------
IF (lhook) CALL dr_hook('MODE_CONV_DOE:CHTC_VERT_DOE',1,zhook_handle)
!
END FUNCTION chtc_vert_doe_1d
!#########################
!
!#############################################
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
!#######################
!#############################################
FUNCTION chtc_up_doe_0d(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, INTENT(IN)                :: PTS     ! Surface temperature (Kelvin)
REAL, INTENT(IN)                :: PTA     ! Air temperature (Kelvin)
REAL                            :: 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
!              ----------------
!
IF (pta < pts) THEN
   pchtc  = 9.482 * (pts - pta)**(1./3.) / (7.283 - 1)
ELSE 
   pchtc  = 1.810 * (pta - pts)**(1./3.) / (1.382 + 1)  
ENDIF
!
!-------------------------------------------------------------------------------
IF (lhook) CALL dr_hook('MODE_CONV_DOE:CHTC_UP_DOE',1,zhook_handle)
!
END FUNCTION chtc_up_doe_0d
!#######################
!
!#############################################
FUNCTION chtc_down_doe(PTS, PTA) RESULT(PCHTC)
!#############################################
!
!!****  *CHTC_DOWN_DOE* - 
!!
!!    PURPOSE
!!    -------
!     function to compute convective surface coefficient for
!     an downward 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_DOWN_DOE',0,zhook_handle)
!
!*       1.    COMPUTE THE CHTC
!              ----------------
!
WHERE(pta > pts)
   pchtc  = 9.482 * (pta - pts)**(1./3.) / (7.283 - 1)
ELSE WHERE
   pchtc  = 1.810 * (pts - pta)**(1./3.) / (1.382 + 1)  
END WHERE
!
!-------------------------------------------------------------------------------
IF (lhook) CALL dr_hook('MODE_CONV_DOE:CHTC_DOWN_DOE',1,zhook_handle)
!
END FUNCTION chtc_down_doe
!#########################
!
!#############################################
FUNCTION chtc_smooth_lee_doe_0d(PCHTCN, PVMOD) RESULT(PCHTC)
!#############################################
!
!!****  *CHTC_SMOOTH_LEE_DOE* - 
!!
!!    PURPOSE
!!    -------
!     function to compute convective surface coefficient for
!     a leeward smooth surface from the natural convection coef and the
!     wind speed
!
!!**  METHOD
!!    ------
!!
!!    from EnergyPlus Engineering Reference, average the leeward/windward coef 
!!
!!    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,INTENT(IN)                :: PCHTCN  ! Convective heat transfer coefficient for natural conv. [W/(m2.K)]
REAL,INTENT(IN)                :: PVMOD   ! wind speed (m/s)
REAL                           :: PCHTC   ! Convective heat transfer coefficient [W/(m2.K)]
!
!*       0.2   Declarations of local variables
!
REAL(KIND=JPRB) :: ZHOOK_HANDLE
REAL :: ZALEE, ZBLEE ! coef for leeward facade
!----------------------------------------------------------------------
IF (lhook) CALL dr_hook('MODE_CONV_DOE:CHTC_SMOOTH_LEE_DOE',0,zhook_handle)
!
!*       1.    INITIALIZE COEF
!              ----------------
!
zalee = 2.86
zblee = 0.617

!*       2.    COMPUTE THE CHTC
!              ----------------
!
pchtc = sqrt(pchtcn**2+(zalee*pvmod**zblee)**2)
!
!-------------------------------------------------------------------------------
IF (lhook) CALL dr_hook('MODE_CONV_DOE:CHTC_SMOOTH_LEE_DOE',1,zhook_handle)
!
END FUNCTION chtc_smooth_lee_doe_0d
!#########################

!#############################################
FUNCTION chtc_smooth_lee_doe_1d(PCHTCN, PVMOD) RESULT(PCHTC)
!#############################################
!
!!****  *CHTC_SMOOTH_LEE_DOE* - 
!!
!!    PURPOSE
!!    -------
!     function to compute convective surface coefficient for
!     a leeward smooth surface from the natural convection coef and the
!     wind speed
!
!!**  METHOD
!!    ------
!!
!!    from EnergyPlus Engineering Reference, average the leeward/windward coef 
!!
!!    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)                :: PCHTCN  ! Convective heat transfer coefficient for natural conv. [W/(m2.K)]
REAL, DIMENSION(:), INTENT(IN)                :: PVMOD   ! wind speed (m/s)
REAL, DIMENSION(SIZE(PCHTCN))                 :: PCHTC   ! Convective heat transfer coefficient [W/(m2.K)]
!
!*       0.2   Declarations of local variables
!
REAL(KIND=JPRB) :: ZHOOK_HANDLE
REAL :: ZALEE, ZBLEE ! coef for windward/leeward facade
!----------------------------------------------------------------------
IF (lhook) CALL dr_hook('MODE_CONV_DOE:CHTC_SMOOTH_LEE_DOE',0,zhook_handle)
!
!*       1.    INITIALIZE COEF
!              ----------------
!
zalee = 2.86
zblee = 0.617

!*       2.    COMPUTE THE CHTC
!              ----------------
!
pchtc = sqrt(pchtcn**2+(zalee*pvmod**zblee)**2)
!
!-------------------------------------------------------------------------------
IF (lhook) CALL dr_hook('MODE_CONV_DOE:CHTC_SMOOTH_LEE_DOE',1,zhook_handle)
!
END FUNCTION chtc_smooth_lee_doe_1d
!#########################
!#############################################
FUNCTION chtc_smooth_wind_doe_0d(PCHTCN, PVMOD) RESULT(PCHTC)
!#############################################
!
!!****  *CHTC_SMOOTH_WIND_DOE* - 
!!
!!    PURPOSE
!!    -------
!     function to compute convective surface coefficient for
!     a smooth surface from the natural convection coef and the
!     wind speed
!
!!**  METHOD
!!    ------
!!
!!    from EnergyPlus Engineering Reference, average the leeward/windward coef 
!!
!!    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,INTENT(IN)                :: PCHTCN  ! Convective heat transfer coefficient for natural conv. [W/(m2.K)]
REAL,INTENT(IN)                :: PVMOD   ! wind speed (m/s)
REAL                           :: PCHTC   ! Convective heat transfer coefficient [W/(m2.K)]
!
!*       0.2   Declarations of local variables
!
REAL(KIND=JPRB) :: ZHOOK_HANDLE
REAL :: ZAWIND, ZBWIND ! coef for windward/leeward facade
!----------------------------------------------------------------------
IF (lhook) CALL dr_hook('MODE_CONV_DOE:CHTC_SMOOTH_WIND_DOE',0,zhook_handle)
!
!*       1.    INITIALIZE COEF
!              ----------------
!
zawind = 2.38
zbwind = 0.89

!*       2.    COMPUTE THE CHTC
!              ----------------
!
pchtc = sqrt(pchtcn**2+(zawind*pvmod**zbwind)**2)
!
!-------------------------------------------------------------------------------
IF (lhook) CALL dr_hook('MODE_CONV_DOE:CHTC_SMOOTH_WIND_DOE',1,zhook_handle)
!
END FUNCTION chtc_smooth_wind_doe_0d
!#########################

!#############################################
FUNCTION chtc_smooth_wind_doe_1d(PCHTCN, PVMOD) RESULT(PCHTC)
!#############################################
!
!!****  *CHTC_SMOOTH_WIND_DOE* - 
!!
!!    PURPOSE
!!    -------
!     function to compute convective surface coefficient for
!     a windward smooth surface from the natural convection coef and the
!     wind speed
!
!!**  METHOD
!!    ------
!!
!!    from EnergyPlus Engineering Reference, average the leeward/windward coef 
!!
!!    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)                :: PCHTCN  ! Convective heat transfer coefficient for natural conv. [W/(m2.K)]
REAL, DIMENSION(:), INTENT(IN)                :: PVMOD   ! wind speed (m/s)
REAL, DIMENSION(SIZE(PCHTCN))                 :: PCHTC   ! Convective heat transfer coefficient [W/(m2.K)]
!
!*       0.2   Declarations of local variables
!
REAL(KIND=JPRB) :: ZHOOK_HANDLE
REAL :: ZAWIND, ZBWIND ! coef for windward facade
!----------------------------------------------------------------------
IF (lhook) CALL dr_hook('MODE_CONV_DOE:CHTC_SMOOTH_WIND_DOE',0,zhook_handle)
!
!*       1.    INITIALIZE COEF
!              ----------------
!
zawind = 2.38
zbwind = 0.89

!*       2.    COMPUTE THE CHTC
!              ----------------
!
pchtc = sqrt(pchtcn**2+(zawind*pvmod**zbwind)**2) 
!
!-------------------------------------------------------------------------------
IF (lhook) CALL dr_hook('MODE_CONV_DOE:CHTC_SMOOTH_WIND_DOE',1,zhook_handle)
!
END FUNCTION chtc_smooth_wind_doe_1d
!#########################
!#############################################
FUNCTION chtc_rough_doe_0d(PCHTCN, PCHTCS, PROUGH) RESULT(PCHTC)
!#############################################
!
!!****  *CHTC_ROUGH_DOE* - 
!!
!!    PURPOSE
!!    -------
!     function to compute convective surface coefficient for
!     a rough surface from the natural convection coef, the
!     smooth surface convective coef and the roughness coef
!
!!**  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, INTENT(IN)                :: PCHTCN  ! Convective heat transfer coefficient for natural conv. [W/(m2.K)]
REAL, INTENT(IN)                :: PCHTCS  ! Convective heat transfer coefficient for a smooth surface [W/(m2.K)]
REAL, INTENT(IN)                :: PROUGH  ! Convective heat transfer coefficient for a smooth surface [W/(m2.K)]
REAL                            :: PCHTC   ! Convective heat transfer coefficient [W/(m2.K)]
!
!*       0.2   Declarations of local variables
!
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!----------------------------------------------------------------------
IF (lhook) CALL dr_hook('MODE_CONV_DOE:CHTC_ROUGH_DOE',0,zhook_handle)
!
!*       1.    COMPUTE THE CHTC
!              ----------------
!
pchtc = pchtcn + prough * (pchtcs - pchtcn)
!
!-------------------------------------------------------------------------------
IF (lhook) CALL dr_hook('MODE_CONV_DOE:CHTC_ROUGH_DOE',1,zhook_handle)
!
END FUNCTION chtc_rough_doe_0d
!#########################
!#############################################
FUNCTION chtc_rough_doe_1d(PCHTCN, PCHTCS, PROUGH) RESULT(PCHTC)
!#############################################
!
!!****  *CHTC_ROUGH_DOE* - 
!!
!!    PURPOSE
!!    -------
!     function to compute convective surface coefficient for
!     a rough surface from the natural convection coef, the
!     smooth surface convective coef and the roughness coef
!
!!**  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)                :: PCHTCN  ! Convective heat transfer coefficient for natural conv. [W/(m2.K)]
REAL, DIMENSION(:), INTENT(IN)                :: PCHTCS  ! Convective heat transfer coefficient for a smooth surface [W/(m2.K)]
REAL, DIMENSION(:), INTENT(IN)                :: PROUGH  ! Convective heat transfer coefficient for a smooth surface [W/(m2.K)]
REAL, DIMENSION(SIZE(PCHTCN))                 :: PCHTC   ! Convective heat transfer coefficient [W/(m2.K)]
!
!*       0.2   Declarations of local variables
!
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!----------------------------------------------------------------------
IF (lhook) CALL dr_hook('MODE_CONV_DOE:CHTC_ROUGH_DOE',0,zhook_handle)
!
!*       1.    COMPUTE THE CHTC
!              ----------------
!
pchtc = pchtcn + prough * (pchtcs - pchtcn)
!
!-------------------------------------------------------------------------------
IF (lhook) CALL dr_hook('MODE_CONV_DOE:CHTC_ROUGH_DOE',1,zhook_handle)
!
END FUNCTION chtc_rough_doe_1d

END MODULE mode_conv_doe