*IF DEF,A19_1A,OR,DEF,A19_2A ABX1F405.76
C *****************************COPYRIGHT****************************** ABX1F405.77
C (c) CROWN COPYRIGHT 1997, METEOROLOGICAL OFFICE, All Rights Reserved. INITACC1.4
C INITACC1.5
C Use, duplication or disclosure of this code is subject to the INITACC1.6
C restrictions as set forth in the contract. INITACC1.7
C INITACC1.8
C Meteorological Office INITACC1.9
C London Road INITACC1.10
C BRACKNELL INITACC1.11
C Berkshire UK INITACC1.12
C RG12 2SZ INITACC1.13
C INITACC1.14
C If no contract has been raised with this copy of the code, the use, INITACC1.15
C duplication or disclosure of it is strictly prohibited. Permission INITACC1.16
C to do so must first be obtained in writing from the Head of Numerical INITACC1.17
C Modelling at the above address. INITACC1.18
C ******************************COPYRIGHT****************************** INITACC1.19
! Initialises accumulated carbon fluxes to zero if new calling period INITACC1.20
! INITACC1.21
! Subroutine Interface: INITACC1.22
SUBROUTINE INIT_ACC(LAND_FIELD, 1ABX1F405.78
& NPP_PFT_ACC,G_LEAF_PHEN_PFT_ACC, ABX1F405.79
& RESP_W_PFT_ACC,RESP_S_ACC,ICODE,CMESSAGE) INITACC1.25
INITACC1.26
IMPLICIT NONE INITACC1.27
! INITACC1.28
! Description: INITACC1.29
! Resets accumulation prognostics to zero if a new TRIFFID calling INITACC1.30
! period is starting. This routine is needed when starting an NRUN INITACC1.31
! from an initial dump created in either of the following situations: INITACC1.32
! INITACC1.33
! i) Initial dump created from a non-TRIFFID run INITACC1.34
! INITACC1.35
! ii) Initial dump created in a TRIFFID run mid-way through a TRIFFID INITACC1.36
! calling period. The NRUN may re-start at the same point within INITACC1.37
! this calling period and continue with the accumulation already INITACC1.38
! part-completed in this dump; in this case this routine will not INITACC1.39
! be used. Alternatively, the NRUN may start a new calling INITACC1.40
! period, in which case the accumulation must begin; this routine INITACC1.41
! allows this by re-setting the relevant prognostics to zero. INITACC1.42
! INITACC1.43
! Current Code Owner: Richard Betts INITACC1.44
! INITACC1.45
! History: INITACC1.46
! Version Date Comment INITACC1.47
! ------- ---- ------- INITACC1.48
! 4.4 10/10/97 Original code. Richard Betts INITACC1.49
! INITACC1.50
! Code Description: INITACC1.51
! Language: FORTRAN 77 + common extensions. INITACC1.52
! This code is written to UMDP3 v6 programming standards. INITACC1.53
! INITACC1.54
C Arguments INITACC1.55
INITACC1.56
*CALL NSTYPES 
INITACC1.57
INITACC1.58
INTEGER INITACC1.59
& LAND_FIELD ! IN number of land points ABX1F405.80
INITACC1.61
INITACC1.65
REAL INITACC1.66
& NPP_PFT_ACC(LAND_FIELD,NPFT) !INOUT Accumulated NPP on ABX1F405.81
C ! Plant Functional Types INITACC1.68
&,G_LEAF_PHEN_PFT_ACC(LAND_FIELD,NPFT)!INOUT Accum. phenological ABX1F405.82
C ! leaf turnover rate PFTs INITACC1.72
&,RESP_W_PFT_ACC(LAND_FIELD,NPFT) !INOUT Accumulated wood ABX1F405.83
C ! respiration on PFTs INITACC1.74
&,RESP_S_ACC(LAND_FIELD) !INOUT Accumulated soil resp. ABX1F405.84
INITACC1.76
INITACC1.77
INTEGER INITACC1.78
& L ! Loop counter for land points INITACC1.79
&,N ! Loop counter for plant functional types INITACC1.80
INITACC1.81
INTEGER ICODE ! Work - Internal return code INITACC1.82
CHARACTER*80 CMESSAGE ! Work - Internal error message INITACC1.83
INITACC1.84
INITACC1.85
WRITE (6,*) INITACC1.86
& 'INIT_ACC: setting accumulation prognostics to zero' INITACC1.87
INITACC1.88
DO L=1,LAND_FIELD ABX1F405.85
DO N=1,NPFT ABX1F405.86
NPP_PFT_ACC(L,N) = 0.0 ABX1F405.87
G_LEAF_PHEN_PFT_ACC(L,N) = 0.0 ABX1F405.88
RESP_W_PFT_ACC(L,N) = 0.0 ABX1F405.89
ENDDO ABX1F405.90
RESP_S_ACC(L) = 0.0 ABX1F405.91
ENDDO ABX1F405.92
INITACC1.104
RETURN INITACC1.105
END INITACC1.106
*ENDIF INITACC1.107