average_flux.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_flux(PFRAC_TILE,             &
                   psfth_tile, psftq_tile,              &
                   psfts_tile, psfco2_tile,             &
                   psfu_tile, psfv_tile,                &
                   psfth, psftq, psfts, psfco2,         &
                   psfu, psfv                           )  
!     ######################################################################
!
!
!!****  *AVERAGE_FLUX*  
!!
!!    PURPOSE
!!    -------
!      Average the fluxes 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
!!    ------
!!      S. Belair           * Meteo-France *
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    10/03/95 
!!      V.Masson    20/03/96  remove abnormal averages and average TS**4 instead
!!                            of TS
!!      (J.Stein)   27/03/96 use only H and LE in the soil scheme
!!      A. Boone    27/11/02 revised to output ALMA variables, and general applications
!-------------------------------------------------------------------------------
!
!*       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 
!                                              ! a given surface
REAL, DIMENSION(:,:), INTENT(IN) :: psfth_tile ! pot. temp. flux  (mK/s)
REAL, DIMENSION(:,:), INTENT(IN) :: psftq_tile ! water vapor flux (m kg/kg/s)
REAL, DIMENSION(:,:), INTENT(IN) :: psfco2_tile! CO2 flux         (m kg/kg/s)
REAL, DIMENSION(:,:,:),INTENT(IN):: psfts_tile ! scalar flux      (m kg/kg/s)
REAL, DIMENSION(:,:), INTENT(IN) :: psfu_tile  ! zonal momentum flux    (pa)
REAL, DIMENSION(:,:), INTENT(IN) :: psfv_tile  ! meridian momentum flux (pa)
REAL, DIMENSION(:),   INTENT(OUT):: psfth      ! pot. temp. flux  (mK/s)
REAL, DIMENSION(:),   INTENT(OUT):: psftq      ! water vapor flux (m kg/kg/s)
REAL, DIMENSION(:,:), INTENT(OUT):: psfts      ! scalar flux      (m kg/kg/s)
REAL, DIMENSION(:),   INTENT(OUT):: psfco2     ! CO2 flux         (m kg/kg/s)
REAL, DIMENSION(:),   INTENT(OUT):: psfu       ! zonal momentum flux    (pa)
REAL, DIMENSION(:),   INTENT(OUT):: psfv       ! meridian momentum flux (pa)
!
!*      0.2    declarations of local variables
!
INTEGER :: jsv   ! scalar loop counter
INTEGER :: jtile ! tile loop counter
REAL(KIND=JPRB) :: zhook_handle
!-------------------------------------------------------------------------------
!
!       0.     Initialization
!              --------------
!
IF (lhook) CALL dr_hook('AVERAGE_FLUX',0,zhook_handle)
psfth(:)   = 0.
psftq(:)   = 0.
psfco2(:)   = 0.
psfu(:)   = 0.
psfv(:)   = 0.
psfts(:,:) = 0.
!
!       1.     Grid-Box average 1d fluxes
!              --------------------------
!
!
DO jtile = 1, SIZE(psfth_tile,2)
!
!  potential temperature flux:
!
   psfth(:)  = psfth(:) + pfrac_tile(:,jtile) * psfth_tile(:,jtile)
!
!  water vapor flux:
!
   psftq(:)  = psftq(:) + pfrac_tile(:,jtile) * psftq_tile(:,jtile)
!
!  carbon flux:
!
   psfco2(:)  = psfco2(:) + pfrac_tile(:,jtile) * psfco2_tile(:,jtile)
!
!  wind surface friction:
!
   psfu(:)  = psfu(:) + pfrac_tile(:,jtile) * psfu_tile(:,jtile)
   psfv(:)  = psfv(:) + pfrac_tile(:,jtile) * psfv_tile(:,jtile)
!   
END DO
!
!
!
!       2.     Grid-Box average 2d fluxes
!              --------------------------
!
DO jsv = 1, SIZE(psfts_tile,2)
!
   DO jtile = 1, SIZE(psfts_tile,3)
!
!     scalar flux
!
      psfts(:,jsv) = psfts(:,jsv) + pfrac_tile(:,jtile) * psfts_tile(:,jsv,jtile)
!
   END DO
!
END DO
IF (lhook) CALL dr_hook('AVERAGE_FLUX',1,zhook_handle)

!-------------------------------------------------------------------------------
!
END SUBROUTINE average_flux