*IF DEF,A03_7A                                                             LLIT7A.2      
C *****************************COPYRIGHT******************************     LLIT7A.3      
C (c) CROWN COPYRIGHT 1997, METEOROLOGICAL OFFICE, All Rights Reserved.    LLIT7A.4      
C                                                                          LLIT7A.5      
C Use, duplication or disclosure of this code is subject to the            LLIT7A.6      
C restrictions as set forth in the contract.                               LLIT7A.7      
C                                                                          LLIT7A.8      
C                Meteorological Office                                     LLIT7A.9      
C                London Road                                               LLIT7A.10     
C                BRACKNELL                                                 LLIT7A.11     
C                Berkshire UK                                              LLIT7A.12     
C                RG12 2SZ                                                  LLIT7A.13     
C                                                                          LLIT7A.14     
C If no contract has been raised with this copy of the code, the use,      LLIT7A.15     
C duplication or disclosure of it is strictly prohibited.  Permission      LLIT7A.16     
C to do so must first be obtained in writing from the Head of Numerical    LLIT7A.17     
C Modelling at the above address.                                          LLIT7A.18     
C ******************************COPYRIGHT******************************    LLIT7A.19     
!***********************************************************************   LLIT7A.20     
! Calculates the leaf turnover rate as a function of temperature and       LLIT7A.21     
! soil water availability                                                  LLIT7A.22     
!***********************************************************************   LLIT7A.23     

      SUBROUTINE LEAF_LIT (LAND_FIELD,VEG_PTS,VEG_INDEX,N,FSMC,TSTAR        1LLIT7A.24     
     &,                    G_LEAF)                                         LLIT7A.25     
                                                                           LLIT7A.26     
      IMPLICIT NONE                                                        LLIT7A.27     
                                                                           LLIT7A.28     
      INTEGER                                                              LLIT7A.29     
     & LAND_FIELD                 ! IN Total number of land points.        LLIT7A.30     
     &,VEG_PTS                    ! IN Number of vegetated points.         LLIT7A.31     
     &,VEG_INDEX(LAND_FIELD)      ! IN Index of vegetated points           LLIT7A.32     
!                                 !    on the land grid.                   LLIT7A.33     
     &,N                          ! IN Plant functional type.              LLIT7A.34     
                                                                           LLIT7A.35     
      REAL                                                                 LLIT7A.36     
     & FSMC(LAND_FIELD)           ! IN Soil moisture availability          LLIT7A.37     
!                                 !    factor.                             LLIT7A.38     
     &,TSTAR(LAND_FIELD)          ! IN Surface temperature (K).            LLIT7A.39     
     &,G_LEAF(LAND_FIELD)         ! OUT Rate of leaf turnover              ABX1F405.932    
!                                 !     (/360days).                        ABX1F405.933    
     &,FM,FT                      ! WORK Soil moisture and leaf            LLIT7A.41     
!                                        temperature amplifiers of         LLIT7A.42     
!                                        leaf turnover.                    LLIT7A.43     
                                                                           LLIT7A.44     
      INTEGER                                                              LLIT7A.45     
     & J,L                        ! Loop counters                          LLIT7A.46     
                                                                           LLIT7A.47     
*CALL NSTYPES                                                              LLIT7A.48     
*CALL TRIF                                                                 LLIT7A.49     
                                                                           LLIT7A.50     
!-----------------------------------------------------------------------   LLIT7A.51     
! Calculate the leaf turnover rate                                         LLIT7A.52     
!-----------------------------------------------------------------------   LLIT7A.53     
      DO J=1,VEG_PTS                                                       LLIT7A.54     
        L = VEG_INDEX(J)                                                   LLIT7A.55     
                                                                           LLIT7A.56     
        FT = 1.0                                                           LLIT7A.57     
        FM = 1.0                                                           LLIT7A.58     
        IF (TSTAR(L) .LT. TLEAF_OF(N)) THEN                                LLIT7A.59     
          FT = 1.0 + DGL_DT(N)*(TLEAF_OF(N)-TSTAR(L))                      LLIT7A.60     
        ELSEIF (FSMC(L) .LT. FSMC_OF(N)) THEN                              LLIT7A.61     
          FM = 1.0 + DGL_DM(N)*(FSMC_OF(N)-FSMC(L))                        LLIT7A.62     
        ENDIF                                                              LLIT7A.63     
                                                                           LLIT7A.64     
        G_LEAF(L) = G_LEAF_0(N)*FT*FM                                      LLIT7A.65     
                                                                           LLIT7A.66     
      ENDDO                                                                LLIT7A.67     
                                                                           LLIT7A.68     
      RETURN                                                               LLIT7A.69     
                                                                           LLIT7A.70     
      END                                                                  LLIT7A.71     
*ENDIF                                                                     LLIT7A.72