convert_cover_frac.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 convert_cover_frac (DTCO, &
                                        PCOVER, OCOVER,           &
                                         PSEA,PNATURE,PTOWN,PWATER   )  
!     ##############################################################
!
!!**** *CONVERT_COVER* convert surface cover classes into secondary
!!                     physiographic variables
!!
!!    PURPOSE
!!    -------
!!
!!    METHOD
!!    ------
!!   
!!    EXTERNAL
!!    --------
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!    AUTHOR
!!    ------
!!
!!    V. Masson        Meteo-France
!!
!!    MODIFICATION
!!    ------------
!!
!!    Original    01/2004
!     
!----------------------------------------------------------------------------
!
!*    0.     DECLARATION
!            -----------
!
!
USE modd_data_cover_n, ONLY : data_cover_t
!
USE modd_surf_par,   ONLY : xundef
!
USE modi_av_pgd
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*    0.1    Declaration of arguments
!            ------------------------
!
!
TYPE(data_cover_t), INTENT(INOUT) :: DTCO
!
REAL, DIMENSION(:,:), INTENT(IN)    :: PCOVER
LOGICAL, DIMENSION(:), INTENT(IN)   :: OCOVER
REAL, DIMENSION(:),   INTENT(OUT)   :: PSEA
REAL, DIMENSION(:),   INTENT(OUT)   :: PNATURE
REAL, DIMENSION(:),   INTENT(OUT)   :: PTOWN
REAL, DIMENSION(:),   INTENT(OUT)   :: PWATER
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
!
!-------------------------------------------------------------------------------
!
!*    1.      cover main type fractions
!             -------------------------
!
IF (lhook) CALL dr_hook('CONVERT_COVER_FRAC',0,zhook_handle)
 CALL av_pgd(dtco, &
              psea    ,pcover(:,:),dtco%XDATA_SEA    (:),'ALL','ARI',ocover(:))
 CALL av_pgd(dtco, &
              ptown   ,pcover(:,:),dtco%XDATA_TOWN   (:),'ALL','ARI',ocover(:))
 CALL av_pgd(dtco, &
              pnature ,pcover(:,:),dtco%XDATA_NATURE (:),'ALL','ARI',ocover(:))
 CALL av_pgd(dtco, &
              pwater  ,pcover(:,:),dtco%XDATA_WATER  (:),'ALL','ARI',ocover(:))

!
WHERE (psea(:) == xundef) psea(:) = 0.
WHERE (pnature(:) == xundef) pnature(:) = 0.
WHERE (ptown(:) == xundef) ptown(:) = 0.
WHERE (pwater(:) == xundef) pwater(:) = 0.
IF (lhook) CALL dr_hook('CONVERT_COVER_FRAC',1,zhook_handle)
!
!-------------------------------------------------------------------------------
!
END SUBROUTINE convert_cover_frac