SURFEX v7.3: /home/dasprezs/EXPORT_v7_3/src/SURFEX/albedo

Go to the documentation of this file.
00001 !     #########
00002       FUNCTION ALBEDO_TA96(PZENITH) RESULT(PDIR_ALB)
00003 !     ##################################################################
00004 !
00005 !!****  *ALBEDO_TA96*  
00006 !!
00007 !!    PURPOSE
00008 !!    -------
00009 !       computes the direct albedo over open water
00010 !
00011 !     
00012 !!**  METHOD
00013 !!    ------
00014 !
00015 !!    EXTERNAL
00016 !!    --------
00017 !!
00018 !!    IMPLICIT ARGUMENTS
00019 !!    ------------------ 
00020 !!
00021 !!      
00022 !!    REFERENCE
00023 !!    ---------
00024 !!
00025 !!      
00026 !!    AUTHOR
00027 !!    ------
00028 !!      V. Masson           * Meteo-France *
00029 !!
00030 !!    MODIFICATIONS
00031 !!    -------------
00032 !!      Original    19/09/03
00033 !       
00034 !-------------------------------------------------------------------------------
00035 !
00036 !*       0.     DECLARATIONS
00037 !               ------------
00038 !
00039 USE MODD_WATER_PAR,  ONLY : XALBCOEF_TA96
00040 !
00041 !
00042 USE YOMHOOK   ,ONLY : LHOOK,   DR_HOOK
00043 USE PARKIND1  ,ONLY : JPRB
00044 !
00045 IMPLICIT NONE
00046 !
00047 !*      0.1    declarations of arguments
00048 !              -------------------------
00049 !
00050 REAL, DIMENSION(:), INTENT(IN)  :: PZENITH     ! zenithal angle (radian)
00051 !
00052 REAL, DIMENSION(SIZE(PZENITH))  :: PDIR_ALB    ! direct albedo on water
00053 REAL(KIND=JPRB) :: ZHOOK_HANDLE
00054 !
00055 !-------------------------------------------------------------------------------
00056 !
00057 IF (LHOOK) CALL DR_HOOK('ALBEDO_TA96',0,ZHOOK_HANDLE)
00058 PDIR_ALB(:) = XALBCOEF_TA96/(1.1*(MAX(COS(PZENITH(:)),0.))**1.4+0.15)
00059 IF (LHOOK) CALL DR_HOOK('ALBEDO_TA96',1,ZHOOK_HANDLE)
00060 !
00061 !-------------------------------------------------------------------------------
00062 !
00063 END FUNCTION ALBEDO_TA96