1 |
jahn |
1.1 |
C $Header$ |
2 |
|
|
C $Name$ |
3 |
|
|
|
4 |
|
|
c MONOD_SIZE.h |
5 |
|
|
c description: Size specification for generalized ecosystem model |
6 |
|
|
c Mick Follows, Scott Grant Fall/Winter 2005 |
7 |
|
|
c Stephanie Dutkiewicz Spring/Summer 2006 |
8 |
|
|
c |
9 |
|
|
c |
10 |
|
|
c npmax = no of "functional groups" of phytoplankton |
11 |
|
|
c nzmax = no of "functional groups" of zooplankton |
12 |
|
|
c |
13 |
|
|
c |
14 |
|
|
|
15 |
|
|
INTEGER npmax |
16 |
|
|
INTEGER nzmax |
17 |
|
|
PARAMETER(npmax=78,nzmax=2) |
18 |
|
|
|
19 |
|
|
C iPO4 :: index of PO4 in Ptracer |
20 |
|
|
C iNO3 :: index of NO3 in Ptracer |
21 |
|
|
C iFeT :: index of FeT in Ptracer |
22 |
|
|
C iSi :: index of Si in Ptracer |
23 |
|
|
C iDOP :: index of DOP in Ptracer |
24 |
|
|
C iDON :: index of DON in Ptracer |
25 |
|
|
C iDOFe :: index of DOFe in Ptracer |
26 |
|
|
C iPOP :: index of POP in Ptracer |
27 |
|
|
C iPON :: index of PON in Ptracer |
28 |
|
|
C iPOFe :: index of POFe in Ptracer |
29 |
|
|
C iPOSi :: index of POSi in Ptracer |
30 |
|
|
C iNH4 :: index of NH4 in Ptracer |
31 |
|
|
C iNO2 :: index of NO2 in Ptracer |
32 |
|
|
C iZoo :: index of first zooplankton |
33 |
|
|
C iPhy :: index of first phytoplankton |
34 |
|
|
C iChl :: index of first phytoplankton Chl (if using dynamic chl) |
35 |
|
|
C nDarwin :: total number of ptracers used by DARWIN |
36 |
|
|
C nCompZooMax :: maximum number of components each zooplankton can have (P,N,...) |
37 |
|
|
C strideCompZoo :: increment between components of zooplankton |
38 |
|
|
C :: iZooComp(nz) = iZoo+(nz-1)*strideTypeZoo+(iComp-1)*strideCompZoo |
39 |
|
|
C strideTypeZoo :: increment between types of zooplankton |
40 |
|
|
C :: iZooComp(nz) = iZoo+(nz-1)*strideTypeZoo+(iComp-1)*strideCompZoo |
41 |
|
|
C the only reasonable choices for the above two are: |
42 |
|
|
C strideCompZoo = 1 |
43 |
|
|
C strideTypeZoo = nCompZooMax |
44 |
|
|
C (the traditional choice) or |
45 |
|
|
C strideCompZoo = nzmax |
46 |
|
|
C strideTypeZoo = 1 |
47 |
|
|
C remember to bring the fields in data.ptracers in the right order ! |
48 |
|
|
|
49 |
|
|
INTEGER iPO4 |
50 |
|
|
INTEGER iNO3 |
51 |
|
|
INTEGER iFeT |
52 |
|
|
INTEGER iSi |
53 |
|
|
INTEGER iDOP |
54 |
|
|
INTEGER iDON |
55 |
|
|
INTEGER iDOFe |
56 |
|
|
INTEGER iPOP |
57 |
|
|
INTEGER iPON |
58 |
|
|
INTEGER iPOFe |
59 |
|
|
INTEGER iPOSi |
60 |
|
|
INTEGER iNH4 |
61 |
|
|
INTEGER iNO2 |
62 |
|
|
INTEGER iZoo |
63 |
|
|
INTEGER iPhy |
64 |
|
|
INTEGER nDarwin |
65 |
|
|
INTEGER nCompZooMax |
66 |
|
|
INTEGER strideCompZoo |
67 |
|
|
INTEGER strideTypeZoo |
68 |
|
|
INTEGER iTot |
69 |
|
|
#ifdef ALLOW_CARBON |
70 |
|
|
INTEGER iDIC |
71 |
|
|
INTEGER iDOC |
72 |
|
|
INTEGER iPOC |
73 |
|
|
INTEGER iPIC |
74 |
|
|
INTEGER iALK |
75 |
|
|
INTEGER iO2 |
76 |
|
|
INTEGER iZoC |
77 |
|
|
#endif |
78 |
|
|
#ifdef DYNAMIC_CHL |
79 |
|
|
INTEGER iChl |
80 |
|
|
#endif |
81 |
|
|
PARAMETER (nCompZooMax=4) |
82 |
|
|
PARAMETER (strideCompZoo=1) |
83 |
|
|
PARAMETER (strideTypeZoo=nCompZooMax) |
84 |
|
|
PARAMETER (iPO4 =1) |
85 |
|
|
PARAMETER (iNO3 =2) |
86 |
|
|
PARAMETER (iFeT =3) |
87 |
|
|
PARAMETER (iSi =4) |
88 |
|
|
PARAMETER (iDOP =5) |
89 |
|
|
PARAMETER (iDON =6) |
90 |
|
|
PARAMETER (iDOFe =7) |
91 |
|
|
PARAMETER (iZoo =8) |
92 |
|
|
PARAMETER (iPOP =iZoo+nzmax*nCompZooMax) |
93 |
|
|
PARAMETER (iPON =iPOP +1) |
94 |
|
|
PARAMETER (iPOFe =iPON +1) |
95 |
|
|
PARAMETER (iPOSi =iPOFe +1) |
96 |
|
|
PARAMETER (iNH4 =iPOSi +1) |
97 |
|
|
PARAMETER (iNO2 =iNH4 +1) |
98 |
|
|
PARAMETER (iPhy =iNO2 +1) |
99 |
|
|
#ifdef DYNAMIC_CHL |
100 |
|
|
PARAMETER (iChl =iPhy +npmax) |
101 |
|
|
PARAMETER (iTOT =iPhy +npmax+npmax) |
102 |
|
|
#else |
103 |
|
|
PARAMETER (iTOT =iPhy +npmax) |
104 |
|
|
#endif |
105 |
|
|
#ifdef ALLOW_CARBON |
106 |
|
|
PARAMETER (iDIC =iTOT) |
107 |
|
|
PARAMETER (iDOC =iDIC+1) |
108 |
|
|
PARAMETER (iPOC =iDOC+1) |
109 |
|
|
PARAMETER (iPIC =iPOC+1) |
110 |
|
|
PARAMETER (iALK =iPIC+1) |
111 |
|
|
PARAMETER (iO2 =iALK+1) |
112 |
|
|
PARAMETER (iZoC =iO2+1) |
113 |
|
|
PARAMETER (nDarwin=iZoC+nzmax-1) |
114 |
|
|
#else |
115 |
|
|
PARAMETER (nDarwin=iTOT-1) |
116 |
|
|
#endif |
117 |
|
|
C iZooP(nzmax) :: index of phosphorus content of each zooplankton type |
118 |
|
|
C iZooN(nzmax) :: index of nitrogen content of each zooplankton type |
119 |
|
|
C iZooFe(nzmax) :: index of iron content of each zooplankton type |
120 |
|
|
C iZooSi(nzmax) :: index of silica content of each zooplankton type |
121 |
|
|
C these are computed in darwin_init_fixed from iZoo, strideCompZoo and strideTypeZoo |
122 |
|
|
INTEGER iZooP (nzmax) |
123 |
|
|
INTEGER iZooN (nzmax) |
124 |
|
|
INTEGER iZooFe(nzmax) |
125 |
|
|
INTEGER iZooSi(nzmax) |
126 |
|
|
#ifdef ALLOW_CARBON |
127 |
|
|
INTEGER iZooC(nzmax) |
128 |
|
|
#endif |
129 |
|
|
COMMON/DARWIN_SIZE/ iZooP, iZooN, iZooFe, iZooSi |
130 |
|
|
#ifdef ALLOW_CARBON |
131 |
|
|
& ,iZooC |
132 |
|
|
#endif |
133 |
|
|
|
134 |
|
|
|