New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
par_pisces.F90 in branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/NEMO/TOP_SRC/PISCES – NEMO

source: branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/NEMO/TOP_SRC/PISCES/par_pisces.F90 @ 6536

Last change on this file since 6536 was 6536, checked in by cetlod, 8 years ago

Minor bugfix

  • Property svn:keywords set to Id
File size: 14.6 KB
Line 
1MODULE par_pisces
2   !!======================================================================
3   !!                        ***  par_pisces  ***
4   !! TOP :   set the PISCES parameters
5   !!======================================================================
6   !! History :   2.0  !  2007-12  (C. Ethe, G. Madec)  revised architecture
7   !!----------------------------------------------------------------------
8   !! NEMO/TOP 3.3 , NEMO Consortium (2010)
9   !! $Id$
10   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
11   !!----------------------------------------------------------------------
12
13   IMPLICIT NONE
14
15#if defined key_pisces_reduced
16   !!---------------------------------------------------------------------
17   !!   'key_pisces_reduced'   :                                LOBSTER bio-model
18   !!---------------------------------------------------------------------
19   LOGICAL, PUBLIC, PARAMETER ::   lk_pisces     = .TRUE.  !: PISCES flag
20   LOGICAL, PUBLIC, PARAMETER ::   lk_p4z        = .FALSE. !: p4z flag
21   INTEGER, PUBLIC, PARAMETER ::   jp_pisces     =  6      !: number of passive tracers
22   INTEGER, PUBLIC, PARAMETER ::   jp_pisces_2d  =  19     !: additional 2d output
23   INTEGER, PUBLIC, PARAMETER ::   jp_pisces_3d  =   3     !: additional 3d output
24   INTEGER, PUBLIC, PARAMETER ::   jp_pisces_trd =   17    !: number of sms trends for PISCES
25
26   ! assign an index in trc arrays for each LOBSTER prognostic variables
27   INTEGER, PUBLIC, PARAMETER ::   jpdet     =  1        !: detritus                    [mmoleN/m3]
28   INTEGER, PUBLIC, PARAMETER ::   jpzoo     =  2        !: zooplancton concentration   [mmoleN/m3]
29   INTEGER, PUBLIC, PARAMETER ::   jpphy     =  3        !: phytoplancton concentration [mmoleN/m3]
30   INTEGER, PUBLIC, PARAMETER ::   jpno3     =  4        !: nitrate concentration       [mmoleN/m3]
31   INTEGER, PUBLIC, PARAMETER ::   jpnh4     =  5        !: ammonium concentration      [mmoleN/m3]
32   INTEGER, PUBLIC, PARAMETER ::   jpdom     =  6        !: dissolved organic matter    [mmoleN/m3]
33
34   ! productive layer depth
35   INTEGER, PUBLIC, PARAMETER ::   jpkb      = 12        !: first vertical layers where biology is active
36   INTEGER, PUBLIC, PARAMETER ::   jpkbm1    = jpkb - 1  !: first vertical layers where biology is active
37
38#elif defined key_pisces  &&  defined key_kriest
39   !!---------------------------------------------------------------------
40   !!   'key_pisces' & 'key_kriest'                 PISCES bio-model + ???
41   !!---------------------------------------------------------------------
42   LOGICAL, PUBLIC, PARAMETER ::   lk_pisces     = .TRUE.  !: PISCES flag
43   LOGICAL, PUBLIC, PARAMETER ::   lk_p4z        = .TRUE. !: p4z flag
44   LOGICAL, PUBLIC, PARAMETER ::   lk_kriest     = .TRUE.  !: Kriest flag
45   INTEGER, PUBLIC, PARAMETER ::   jp_pisces     =  23     !: number of passive tracers
46   INTEGER, PUBLIC, PARAMETER ::   jp_pisces_2d  =  13     !: additional 2d output
47   INTEGER, PUBLIC, PARAMETER ::   jp_pisces_3d  =  18     !: additional 3d output
48   INTEGER, PUBLIC, PARAMETER ::   jp_pisces_trd =   1     !: number of sms trends for PISCES
49
50   ! assign an index in trc arrays for each LOBSTER prognostic variables
51   !    WARNING: be carefull about the order when reading the restart
52        !   !!gm  this warning should be obsolet with IOM
53   INTEGER, PUBLIC, PARAMETER ::   jpdic =  1    !: dissolved inoganic carbon concentration
54   INTEGER, PUBLIC, PARAMETER ::   jptal =  2    !: total alkalinity
55   INTEGER, PUBLIC, PARAMETER ::   jpoxy =  3    !: oxygen carbon concentration
56   INTEGER, PUBLIC, PARAMETER ::   jpcal =  4    !: calcite  concentration
57   INTEGER, PUBLIC, PARAMETER ::   jppo4 =  5    !: phosphate concentration
58   INTEGER, PUBLIC, PARAMETER ::   jppoc =  6    !: small particulate organic phosphate concentration
59   INTEGER, PUBLIC, PARAMETER ::   jpsil =  7    !: silicate concentration
60   INTEGER, PUBLIC, PARAMETER ::   jpphy =  8    !: phytoplancton concentration
61   INTEGER, PUBLIC, PARAMETER ::   jpzoo =  9    !: zooplancton concentration
62   INTEGER, PUBLIC, PARAMETER ::   jpdoc = 10    !: dissolved organic carbon concentration
63   INTEGER, PUBLIC, PARAMETER ::   jpdia = 11    !: Diatoms Concentration
64   INTEGER, PUBLIC, PARAMETER ::   jpmes = 12    !: Mesozooplankton Concentration
65   INTEGER, PUBLIC, PARAMETER ::   jpdsi = 13    !: Diatoms Silicate Concentration
66   INTEGER, PUBLIC, PARAMETER ::   jpfer = 14    !: Iron Concentration
67   INTEGER, PUBLIC, PARAMETER ::   jpnum = 15    !: Big iron particles Concentration
68   INTEGER, PUBLIC, PARAMETER ::   jpsfe = 16    !: number of particulate organic phosphate concentration
69   INTEGER, PUBLIC, PARAMETER ::   jpdfe = 17    !: Diatoms iron Concentration
70   INTEGER, PUBLIC, PARAMETER ::   jpgsi = 18    !: (big) Silicate Concentration
71   INTEGER, PUBLIC, PARAMETER ::   jpnfe = 19    !: Nano iron Concentration
72   INTEGER, PUBLIC, PARAMETER ::   jpnch = 20    !: Nano Chlorophyll Concentration
73   INTEGER, PUBLIC, PARAMETER ::   jpdch = 21    !: Diatoms Chlorophyll Concentration
74   INTEGER, PUBLIC, PARAMETER ::   jpno3 = 22    !: Nitrates Concentration
75   INTEGER, PUBLIC, PARAMETER ::   jpnh4 = 23    !: Ammonium Concentration
76
77
78#elif defined key_pisces
79   !!---------------------------------------------------------------------
80   !!   'key_pisces'   :                         standard PISCES bio-model
81   !!---------------------------------------------------------------------
82   LOGICAL, PUBLIC, PARAMETER ::   lk_pisces     = .TRUE.  !: PISCES flag
83   LOGICAL, PUBLIC, PARAMETER ::   lk_p4z        = .TRUE.  !: p4z flag
84   LOGICAL, PUBLIC, PARAMETER ::   lk_kriest     = .FALSE. !: Kriest flag
85
86   ! assign an index in trc arrays for each LOBSTER prognostic variables
87   !    WARNING: be carefull about the order when reading the restart
88        !   !!gm  this warning should be obsolet with IOM
89#if ! defined key_n15  && ! defined key_c13
90   LOGICAL, PUBLIC, PARAMETER ::   lk_n15        = .FALSE.  !: N15 flag
91   LOGICAL, PUBLIC, PARAMETER ::   lk_c13        = .FALSE.  !: C13 flag
92   INTEGER, PUBLIC, PARAMETER ::   jp_pisces     = 24      !: number of PISCES passive tracers
93   INTEGER, PUBLIC, PARAMETER ::   jp_pisces_2d  = 20      !: additional 2d output
94   INTEGER, PUBLIC, PARAMETER ::   jp_pisces_3d  = 14      !: additional 3d output
95   INTEGER, PUBLIC, PARAMETER ::   jp_pisces_trd =  1      !: number of sms trends for PISCES
96#endif
97   INTEGER, PUBLIC, PARAMETER ::   jpdic =  1    !: dissolved inoganic carbon concentration
98   INTEGER, PUBLIC, PARAMETER ::   jptal =  2    !: total alkalinity
99   INTEGER, PUBLIC, PARAMETER ::   jpoxy =  3    !: oxygen carbon concentration
100   INTEGER, PUBLIC, PARAMETER ::   jpcal =  4    !: calcite  concentration
101   INTEGER, PUBLIC, PARAMETER ::   jppo4 =  5    !: phosphate concentration
102   INTEGER, PUBLIC, PARAMETER ::   jppoc =  6    !: small particulate organic phosphate concentration
103   INTEGER, PUBLIC, PARAMETER ::   jpsil =  7    !: silicate concentration
104   INTEGER, PUBLIC, PARAMETER ::   jpphy =  8    !: phytoplancton concentration
105   INTEGER, PUBLIC, PARAMETER ::   jpzoo =  9    !: zooplancton concentration
106   INTEGER, PUBLIC, PARAMETER ::   jpdoc = 10    !: dissolved organic carbon concentration
107   INTEGER, PUBLIC, PARAMETER ::   jpdia = 11    !: Diatoms Concentration
108   INTEGER, PUBLIC, PARAMETER ::   jpmes = 12    !: Mesozooplankton Concentration
109   INTEGER, PUBLIC, PARAMETER ::   jpdsi = 13    !: Diatoms Silicate Concentration
110   INTEGER, PUBLIC, PARAMETER ::   jpfer = 14    !: Iron Concentration
111   INTEGER, PUBLIC, PARAMETER ::   jpbfe = 15    !: Big iron particles Concentration
112   INTEGER, PUBLIC, PARAMETER ::   jpgoc = 16    !: big particulate organic phosphate concentration
113   INTEGER, PUBLIC, PARAMETER ::   jpsfe = 17    !: Small iron particles Concentration
114   INTEGER, PUBLIC, PARAMETER ::   jpdfe = 18    !: Diatoms iron Concentration
115   INTEGER, PUBLIC, PARAMETER ::   jpgsi = 19    !: (big) Silicate Concentration
116   INTEGER, PUBLIC, PARAMETER ::   jpnfe = 20    !: Nano iron Concentration
117   INTEGER, PUBLIC, PARAMETER ::   jpnch = 21    !: Nano Chlorophyll Concentration
118   INTEGER, PUBLIC, PARAMETER ::   jpdch = 22    !: Diatoms Chlorophyll Concentration
119   INTEGER, PUBLIC, PARAMETER ::   jpno3 = 23    !: Nitrates Concentration
120   INTEGER, PUBLIC, PARAMETER ::   jpnh4 = 24    !: Ammonium Concentration
121
122#if defined key_n15  && defined key_c13
123   LOGICAL, PUBLIC, PARAMETER ::   lk_n15        = .TRUE.  !: N15 flag
124   LOGICAL, PUBLIC, PARAMETER ::   lk_c13        = .TRUE.  !: C13 flag
125   INTEGER, PUBLIC, PARAMETER ::   jp_pisces     = 42      !: number of PISCES passive tracers
126   INTEGER, PUBLIC, PARAMETER ::   jp_pisces_2d  = 20      !: additional 2d output
127   INTEGER, PUBLIC, PARAMETER ::   jp_pisces_3d  = 14      !: additional 3d output
128   INTEGER, PUBLIC, PARAMETER ::   jp_pisces_trd =  1      !: number of sms trends for PISCES
129
130   INTEGER, PUBLIC, PARAMETER ::   jp15poc = 25  !: 15N small particulate organic concentration
131   INTEGER, PUBLIC, PARAMETER ::   jp15phy = 26  !: 15N phytoplancton concentration
132   INTEGER, PUBLIC, PARAMETER ::   jp15zoo = 27  !: 15N zooplancton concentration
133   INTEGER, PUBLIC, PARAMETER ::   jp15doc = 28  !: 15N dissolved organic carbon concentration
134   INTEGER, PUBLIC, PARAMETER ::   jp15dia = 29  !: 15N Diatoms Concentration
135   INTEGER, PUBLIC, PARAMETER ::   jp15mes = 30  !: 15N Mesozooplankton Concentration
136   INTEGER, PUBLIC, PARAMETER ::   jp15goc = 31  !: 15N big particulate organic concentration
137   INTEGER, PUBLIC, PARAMETER ::   jp15no3 = 32  !: 15N Nitrates Concentration
138   INTEGER, PUBLIC, PARAMETER ::   jp15nh4 = 33  !: 15N Ammonium Concentration
139
140   INTEGER, PUBLIC, PARAMETER ::   jp13dic = 34  !: 13C dissolved inorganic concentration
141   INTEGER, PUBLIC, PARAMETER ::   jp13cal = 35  !: 13C calcite Concentration
142   INTEGER, PUBLIC, PARAMETER ::   jp13poc = 36  !: 13C small particulate organic concentration
143   INTEGER, PUBLIC, PARAMETER ::   jp13phy = 37  !: 13C phytoplancton concentration
144   INTEGER, PUBLIC, PARAMETER ::   jp13zoo = 38  !: 13C zooplancton concentration
145   INTEGER, PUBLIC, PARAMETER ::   jp13doc = 39  !: 13C dissolved organic carbon concentration
146   INTEGER, PUBLIC, PARAMETER ::   jp13dia = 40  !: 13C Diatoms Concentration
147   INTEGER, PUBLIC, PARAMETER ::   jp13mes = 41  !: 13C Mesozooplankton Concentration
148   INTEGER, PUBLIC, PARAMETER ::   jp13goc = 42  !: 13C big particulate organic concentration
149#endif
150#if defined key_n15  && ! defined key_c13
151   LOGICAL, PUBLIC, PARAMETER ::   lk_n15        = .TRUE.  !: N15 flag
152   LOGICAL, PUBLIC, PARAMETER ::   lk_c13        = .FALSE.  !: C13 flag
153   INTEGER, PUBLIC, PARAMETER ::   jp_pisces     = 33      !: number of PISCES passive tracers
154   INTEGER, PUBLIC, PARAMETER ::   jp_pisces_2d  = 20      !: additional 2d output
155   INTEGER, PUBLIC, PARAMETER ::   jp_pisces_3d  = 14      !: additional 3d output
156   INTEGER, PUBLIC, PARAMETER ::   jp_pisces_trd =  1      !: number of sms trends for PISCES
157
158   INTEGER, PUBLIC, PARAMETER ::   jp15poc = 25  !: 15N small particulate organic phosphate concentration
159   INTEGER, PUBLIC, PARAMETER ::   jp15phy = 26  !: 15N phytoplancton concentration
160   INTEGER, PUBLIC, PARAMETER ::   jp15zoo = 27  !: 15N zooplancton concentration
161   INTEGER, PUBLIC, PARAMETER ::   jp15doc = 28  !: 15N dissolved organic carbon concentration
162   INTEGER, PUBLIC, PARAMETER ::   jp15dia = 29  !: 15N Diatoms Concentration
163   INTEGER, PUBLIC, PARAMETER ::   jp15mes = 30  !: 15N Mesozooplankton Concentration
164   INTEGER, PUBLIC, PARAMETER ::   jp15goc = 31  !: 15N big particulate organic phosphate concentration
165   INTEGER, PUBLIC, PARAMETER ::   jp15no3 = 32  !: 15N Nitrates Concentration
166   INTEGER, PUBLIC, PARAMETER ::   jp15nh4 = 33  !: 15N Ammonium Concentration
167#endif
168#if ! defined key_n15  &&  defined key_c13
169   LOGICAL, PUBLIC, PARAMETER ::   lk_n15        = .FALSE.  !: N15 flag
170   LOGICAL, PUBLIC, PARAMETER ::   lk_c13        = .TRUE.  !: C13 flag
171   INTEGER, PUBLIC, PARAMETER ::   jp_pisces     = 33      !: number of PISCES passive tracers
172   INTEGER, PUBLIC, PARAMETER ::   jp_pisces_2d  = 20      !: additional 2d output
173   INTEGER, PUBLIC, PARAMETER ::   jp_pisces_3d  = 14      !: additional 3d output
174   INTEGER, PUBLIC, PARAMETER ::   jp_pisces_trd =  1      !: number of sms trends for PISCES
175
176   INTEGER, PUBLIC, PARAMETER ::   jp13dic = 25  !: 13C dissolved inorganic concentration
177   INTEGER, PUBLIC, PARAMETER ::   jp13cal = 26  !: 13C calcite Concentration
178   INTEGER, PUBLIC, PARAMETER ::   jp13poc = 27  !: 13C small particulate organic concentration
179   INTEGER, PUBLIC, PARAMETER ::   jp13phy = 28  !: 13C phytoplancton concentration
180   INTEGER, PUBLIC, PARAMETER ::   jp13zoo = 29  !: 13C zooplancton concentration
181   INTEGER, PUBLIC, PARAMETER ::   jp13doc = 30  !: 13C dissolved organic carbon concentration
182   INTEGER, PUBLIC, PARAMETER ::   jp13dia = 31  !: 13C Diatoms Concentration
183   INTEGER, PUBLIC, PARAMETER ::   jp13mes = 32  !: 13C Mesozooplankton Concentration
184   INTEGER, PUBLIC, PARAMETER ::   jp13goc = 33  !: 13C big particulate organic concentration
185#endif
186
187#else
188   !!---------------------------------------------------------------------
189   !!   Default                                   No CFC geochemical model
190   !!---------------------------------------------------------------------
191   LOGICAL, PUBLIC, PARAMETER ::   lk_pisces     = .FALSE.  !: PISCES flag
192   LOGICAL, PUBLIC, PARAMETER ::   lk_p4z        = .FALSE.  !: p4z flag
193   INTEGER, PUBLIC, PARAMETER ::   jp_pisces     =  0       !: No CFC tracers
194   INTEGER, PUBLIC, PARAMETER ::   jp_pisces_2d  =  0       !: No CFC additional 2d output arrays
195   INTEGER, PUBLIC, PARAMETER ::   jp_pisces_3d  =  0       !: No CFC additional 3d output arrays
196   INTEGER, PUBLIC, PARAMETER ::   jp_pisces_trd =  0       !: number of sms trends for PISCES
197#endif
198
199   ! Starting/ending PISCES do-loop indices (N.B. no PISCES : jpl_pcs < jpf_pcs the do-loop are never done)
200   INTEGER, PUBLIC, PARAMETER ::   jp_pcs0     = 1                  !: First index of PISCES tracers
201   INTEGER, PUBLIC, PARAMETER ::   jp_pcs1     = jp_pisces          !: Last  index of PISCES tracers
202   INTEGER, PUBLIC, PARAMETER ::   jp_pcs0_2d  = 1               !: First index of 2D diag
203   INTEGER, PUBLIC, PARAMETER ::   jp_pcs1_2d  = jp_pisces_2d    !: Last  index of 2D diag
204   INTEGER, PUBLIC, PARAMETER ::   jp_pcs0_3d  = 1               !: First index of 3D diag
205   INTEGER, PUBLIC, PARAMETER ::   jp_pcs1_3d  = jp_pisces_3d    !: Last  index of 3d diag
206   INTEGER, PUBLIC, PARAMETER ::   jp_pcs0_trd = 1              !: First index of bio diag
207   INTEGER, PUBLIC, PARAMETER ::   jp_pcs1_trd = jp_pisces_trd  !: Last  index of bio diag
208
209
210   !!======================================================================
211END MODULE par_pisces
Note: See TracBrowser for help on using the repository browser.