average_tsurf.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.
!     #########
      SUBROUTINE average_tsurf(PFRAC_TILE, PTSURF_TILE, PTSURF)  
!     ######################################################
!
!
!!****  *AVERAGE_TSURF*  
!!
!!    PURPOSE
!!    -------
!      Average the surface temperature from the land and water surfaces depending on the
!      fraction of each surface cover type in the mesh area.
!     
!!**  METHOD
!!    ------
!
!!    EXTERNAL
!!    --------
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------ 
!!
!!      
!!    REFERENCE
!!    ---------
!!      
!!    AUTHOR
!!    ------
!!      B. Decharme          * Meteo-France *
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    23/04/2013 
!-------------------------------------------------------------------------------
!
!*       0.     DECLARATIONS
!               ------------
!
!
USE modd_surf_par,   ONLY : xundef
!
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*      0.1    declarations of arguments
!
!
!*      0.1    declarations of arguments
!
REAL, DIMENSION(:,:), INTENT(IN) :: pfrac_tile  ! Fraction in a mesh-area of 
REAL, DIMENSION(:,:), INTENT(IN) :: ptsurf_tile ! surface effective temperature (K)
REAL, DIMENSION(:),   INTENT(OUT):: ptsurf      ! surface effective temperature   (K)
!
!*      0.2    declarations of local variables
!
INTEGER :: ini, inp  ! dimenssion
INTEGER ::  ji, jp    ! loop counter
!
REAL(KIND=JPRB) :: zhook_handle
!-------------------------------------------------------------------------------
!
IF (lhook) CALL dr_hook('AVERAGE_TSURF',0,zhook_handle)
!
ini  = SIZE(pfrac_tile,1)
inp  = SIZE(pfrac_tile,2)
!
ptsurf(:)   = 0.
!
DO jp = 1,inp
  DO ji = 1,ini
     ptsurf(ji) = ptsurf(ji) + pfrac_tile(ji,jp) * ptsurf_tile(ji,jp)
  ENDDO
ENDDO
!
IF (lhook) CALL dr_hook('AVERAGE_TSURF',1,zhook_handle)

!-------------------------------------------------------------------------------
!
END SUBROUTINE average_tsurf