6 SUBROUTINE lailoss(PVEG, PSEFOLD, PANMAX, PANDAY, PANFM, PBIOMASS)
60 USE yomhook
,ONLY : lhook, dr_hook
61 USE parkind1
,ONLY : jprb
68 REAL,
DIMENSION(:),
INTENT(IN) :: pveg
69 REAL,
DIMENSION(:),
INTENT(IN) :: psefold
70 REAL,
DIMENSION(:),
INTENT(IN) :: panmax
71 REAL,
DIMENSION(:),
INTENT(IN) :: panday
73 REAL,
DIMENSION(:),
INTENT(INOUT) :: panfm
74 REAL,
DIMENSION(:),
INTENT(INOUT) :: pbiomass
78 REAL,
DIMENSION(SIZE(PSEFOLD)) :: zxsefold, zxm
80 REAL(KIND=JPRB) :: zhook_handle
83 IF (lhook) CALL dr_hook(
'LAILOSS',0,zhook_handle)
85 zbmcoef = xmc/(xmco2*xpcco2)
94 zxsefold(:) = psefold(:)*min(1.0, panfm(:)/panmax(:))/xday
98 zxsefold(:) = max(1.0e-8,zxsefold(:))
102 zxsefold(:) = max(5.,zxsefold(:))
106 zxm(:) = pbiomass(:)*(1.0-exp(-1.0/zxsefold(:)))
110 pbiomass(:) = pbiomass(:) - zxm(:)
115 pbiomass(:) = pbiomass(:) + panday(:)*zbmcoef
123 IF (lhook) CALL dr_hook(
'LAILOSS',1,zhook_handle)
subroutine lailoss(PVEG, PSEFOLD, PANMAX, PANDAY, PANFM, PBIOMASS)