convert_cover_ch_isba.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_ch_isba (DTCO, &
                                           pcover,ocover,psoilrc_so2,psoilrc_o3)
!     ##############################################################
!
!!**** *CONVERT_COVER* convert surface cover classes into secondary 
!!                     physiographic variables for ISBA
!!
!!    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_data_cover,     ONLY : xdata_soilrc_so2, xdata_soilrc_o3 
!
USE modd_data_cover_par, ONLY : nvegtype, jpcover
!
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)   :: psoilrc_so2
REAL, DIMENSION(:,:),   INTENT(OUT)   :: psoilrc_o3
REAL(KIND=JPRB) :: zhook_handle
!
!
!*    0.2    Declaration of local variables
!            ------------------------------
!
!-------------------------------------------------------------------------------
!
IF (lhook) CALL dr_hook('CONVERT_COVER_CH_ISBA',0,zhook_handle)
!
IF (ASSOCIATED(dtco%XDATA_WEIGHT)) DEALLOCATE(dtco%XDATA_WEIGHT)
!
 CALL av_pgd(dtco, &
              psoilrc_so2 ,pcover ,xdata_soilrc_so2(:,:) ,'NAT','ARI',ocover,kdecade=1)
 CALL av_pgd(dtco, &
              psoilrc_o3  ,pcover ,xdata_soilrc_o3(:,:) ,'NAT','ARI',ocover,kdecade=1)
!
IF (ASSOCIATED(dtco%XDATA_WEIGHT)) DEALLOCATE(dtco%XDATA_WEIGHT)
!
IF (lhook) CALL dr_hook('CONVERT_COVER_CH_ISBA',1,zhook_handle)
!
!-------------------------------------------------------------------------------
!
END SUBROUTINE convert_cover_ch_isba