albedo_from_nir_vis.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 albedo_from_nir_vis(PSW_BANDS,PALBNIR,PALBVIS,PALBUV,PDIR_ALB,PSCA_ALB)
!     ###########################################################################
!
!!****  *ALBEDO_FROM_NIR_VIS* - routine to initialize albedo for
!!                              any wavelength from near-infra-red,
!!                              visible and UV albedo
!!
!!    PURPOSE
!!    -------
!!
!!**  METHOD
!!    ------
!!
!!    EXTERNAL
!!    --------
!!
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!
!!    AUTHOR
!!    ------
!!      V. Masson   *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    02/2003 
!-------------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!              ------------
!
USE modd_isba_par,   ONLY : xred_edge, xuv_edge
USE modd_surf_par,   ONLY : xundef
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*       0.1   Declarations of arguments
!              -------------------------
!
REAL, DIMENSION(:),   INTENT(IN) :: psw_bands ! mean wavelength of each shortwave band (m)
REAL, DIMENSION(:),   INTENT(IN) :: palbnir   ! near infra-red albedo
REAL, DIMENSION(:),   INTENT(IN) :: palbvis   ! visible albedo
REAL, DIMENSION(:),   INTENT(IN) :: palbuv    ! UV albedo
REAL, DIMENSION(:,:), INTENT(OUT):: pdir_alb  ! direct  albedo for each wavelength
REAL, DIMENSION(:,:), INTENT(OUT):: psca_alb  ! diffuse albedo for each wavelength
!
!*       0.2   Declarations of local variables
!              -------------------------------
!
INTEGER :: iswb ! number of SW spectral bands
INTEGER :: jswb ! loop counter on number of SW spectral bands
REAL(KIND=JPRB) :: zhook_handle
!-------------------------------------------------------------------------------
!
IF (lhook) CALL dr_hook('ALBEDO_FROM_NIR_VIS',0,zhook_handle)
iswb = SIZE(psw_bands)
!
pdir_alb(:,:) = xundef
psca_alb(:,:) = xundef
!
IF (iswb==1) THEN
  WHERE(palbnir(:)/= xundef) pdir_alb(:,1) = 0.5*(palbnir(:)+palbvis(:))
ELSE 
  DO jswb=1,iswb
    IF (psw_bands(jswb)>xred_edge ) THEN      ! XRED_EDGE=0.7 micro-m
      pdir_alb(:,jswb) = palbnir(:)
    ELSE IF (psw_bands(jswb)<xuv_edge ) THEN  ! XUV_EDGE=0.25 micro-m
      pdir_alb(:,jswb) = palbuv(:)
    ELSE
      pdir_alb(:,jswb) = palbvis(:)
    END IF
  END DO
END IF
!
psca_alb(:,:) = pdir_alb(:,:)
IF (lhook) CALL dr_hook('ALBEDO_FROM_NIR_VIS',1,zhook_handle)
!-------------------------------------------------------------------------------
!
END SUBROUTINE albedo_from_nir_vis