1 | MODULE sms_medusa |
---|
2 | !!---------------------------------------------------------------------- |
---|
3 | !! *** sms_medusa.F90 *** |
---|
4 | !! TOP : MEDUSA Source Minus Sink variables |
---|
5 | !!---------------------------------------------------------------------- |
---|
6 | !! History : - ! 1999-09 (M. Levy) original code |
---|
7 | !! - ! 2000-12 (O. Aumont, E. Kestenare) add sediment |
---|
8 | !! 1.0 ! 2005-10 (C. Ethe) F90 |
---|
9 | !! 1.0 ! 2005-03 (A-S Kremeur) add fphylab, fzoolab, fdetlab, fdbod |
---|
10 | !! - ! 2005-06 (A-S Kremeur) add sedpocb, sedpocn, sedpoca |
---|
11 | !! 2.0 ! 2007-04 (C. Deltel, G. Madec) Free form and modules |
---|
12 | !! - ! 2008-08 (K. Popova) adaptation for MEDUSA |
---|
13 | !! - ! 2008-11 (A. Yool) continuing adaptation for MEDUSA |
---|
14 | !! - ! 2010-03 (A. Yool) updated for branch inclusion |
---|
15 | !! - ! 2011-04 (A. Yool) updated for ROAM project |
---|
16 | !!---------------------------------------------------------------------- |
---|
17 | |
---|
18 | #if defined key_medusa |
---|
19 | !!---------------------------------------------------------------------- |
---|
20 | !! 'key_medusa' MEDUSA model |
---|
21 | !!---------------------------------------------------------------------- |
---|
22 | USE par_oce |
---|
23 | USE par_trc |
---|
24 | |
---|
25 | IMPLICIT NONE |
---|
26 | PUBLIC |
---|
27 | |
---|
28 | !!---------------------------------------------------------------------- |
---|
29 | !! NEMO/TOP 1.0 , LOCEAN-IPSL (2005) |
---|
30 | !! $Id$ |
---|
31 | !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) |
---|
32 | !!---------------------------------------------------------------------- |
---|
33 | |
---|
34 | INTEGER :: numnatp_ref = -1 !! Logical units for namelist medusa |
---|
35 | INTEGER :: numnatp_cfg = -1 !! Logical units for namelist medusa |
---|
36 | INTEGER :: numonp = -1 !! Logical unit for namelist medusa output |
---|
37 | |
---|
38 | !!---------------------------------------------------------------------- |
---|
39 | !! Biological parameters |
---|
40 | !!---------------------------------------------------------------------- |
---|
41 | !! |
---|
42 | !! Primary production and chl related quantities |
---|
43 | REAL(wp) :: xxi !: conversion factor from gC to mmolN |
---|
44 | REAL(wp) :: xaln !: Chl-a specific initial slope of P-I curve for non-diatoms |
---|
45 | REAL(wp) :: xald !: Chl-a specific initial slope of P-I curve for diatoms |
---|
46 | INTEGER :: jphy !: phytoplankton T-dependent growth switch |
---|
47 | REAL(wp) :: xvpn !: maximum growth rate for non-diatoms |
---|
48 | REAL(wp) :: xvpd !: maximum growth rate for diatoms |
---|
49 | REAL(wp) :: xthetam !: maximum Chl to C ratio for non-diatoms |
---|
50 | REAL(wp) :: xthetamd !: maximum Chl to C ratio for diatoms |
---|
51 | REAL(wp) :: jq10 !: specific Q10 value (jphy==2) |
---|
52 | !! |
---|
53 | !! Diatom silicon parameters |
---|
54 | REAL(wp) :: xsin0 !: minimum diatom Si:N ratio |
---|
55 | REAL(wp) :: xnsi0 !: minimum diatom N:Si ratio |
---|
56 | REAL(wp) :: xuif !: hypothetical growth ratio at infinite Si:N ratio |
---|
57 | !! |
---|
58 | !! Nutrient limitation |
---|
59 | INTEGER :: jliebig !: Liebig nutrient uptake switch |
---|
60 | REAL(wp) :: xnln !: half-sat constant for DIN uptake by non-diatoms |
---|
61 | REAL(wp) :: xnld !: half-sat constant for DIN uptake by diatoms |
---|
62 | REAL(wp) :: xsld !: half-sat constant for Si uptake by diatoms |
---|
63 | REAL(wp) :: xfln !: half-sat constant for Fe uptake by non-diatoms |
---|
64 | REAL(wp) :: xfld !: half-sat constant for Fe uptake by diatoms |
---|
65 | !! |
---|
66 | !! Grazing |
---|
67 | REAL(wp) :: xgmi !: microzoo maximum growth rate |
---|
68 | REAL(wp) :: xgme !: mesozoo maximum growth rate |
---|
69 | REAL(wp) :: xkmi !: microzoo grazing half-sat parameter |
---|
70 | REAL(wp) :: xkme !: mesozoo grazing half-sat parameter |
---|
71 | REAL(wp) :: xphi !: micro/mesozoo grazing inefficiency |
---|
72 | REAL(wp) :: xbetan !: micro/mesozoo N assimilation efficiency |
---|
73 | REAL(wp) :: xbetac !: micro/mesozoo C assimilation efficiency |
---|
74 | REAL(wp) :: xkc !: micro/mesozoo net C growth efficiency |
---|
75 | REAL(wp) :: xpmipn !: grazing preference of microzoo for non-diatoms |
---|
76 | REAL(wp) :: xpmid !: grazing preference of microzoo for diatoms |
---|
77 | REAL(wp) :: xpmepn !: grazing preference of mesozoo for non-diatoms |
---|
78 | REAL(wp) :: xpmepd !: grazing preference of mesozoo for diatoms |
---|
79 | REAL(wp) :: xpmezmi !: grazing preference of mesozoo for microzoo |
---|
80 | REAL(wp) :: xpmed !: grazing preference of mesozoo for detritus |
---|
81 | !! |
---|
82 | !! Metabolic losses |
---|
83 | REAL(wp) :: xmetapn !: non-diatom metabolic loss rate |
---|
84 | REAL(wp) :: xmetapd !: diatom metabolic loss rate |
---|
85 | REAL(wp) :: xmetazmi !: microzoo metabolic loss rate |
---|
86 | REAL(wp) :: xmetazme !: mesozoo metabolic loss rate |
---|
87 | !! |
---|
88 | !! Mortality losses |
---|
89 | INTEGER :: jmpn !: non-diatom mortality functional form |
---|
90 | REAL(wp) :: xmpn !: non-diatom mortality rate |
---|
91 | REAL(wp) :: xkphn !: non-diatom mortality half-sat constant |
---|
92 | INTEGER :: jmpd !: diatom mortality functional form |
---|
93 | REAL(wp) :: xmpd !: diatom mortality rate |
---|
94 | REAL(wp) :: xkphd !: diatom mortality half-sat constant |
---|
95 | INTEGER :: jmzmi !: microzoo mortality functional form |
---|
96 | REAL(wp) :: xmzmi !: microzoo mortality rate |
---|
97 | REAL(wp) :: xkzmi !: microzoo mortality half-sat constant |
---|
98 | INTEGER :: jmzme !: mesozoo mortality functional form |
---|
99 | REAL(wp) :: xmzme !: mesozoo mortality rate |
---|
100 | REAL(wp) :: xkzme !: mesozoo mortality half-sat constant |
---|
101 | !! |
---|
102 | !! Remineralisation |
---|
103 | INTEGER :: jmd !: detritus T-dependent remineralisation switch |
---|
104 | INTEGER :: jsfd !: accelerate seafloor detritus remin. switch |
---|
105 | REAL(wp) :: xmd !: detrital nitrogen remineralisation rate |
---|
106 | REAL(wp) :: xmdc !: detrital carbon remineralisation rate |
---|
107 | !! |
---|
108 | !! Stochiometric ratios |
---|
109 | REAL(wp) :: xthetapn !: non-diatom C:N ratio |
---|
110 | REAL(wp) :: xthetapd !: diatom C:N ratio |
---|
111 | REAL(wp) :: xthetazmi !: microzoo C:N ratio |
---|
112 | REAL(wp) :: xthetazme !: mesozoo C:N ratio |
---|
113 | REAL(wp) :: xthetad !: detritus C:N ratio |
---|
114 | REAL(wp) :: xrfn !: phytoplankton Fe:N ratio |
---|
115 | REAL(wp) :: xrsn !: diatom Si:N ratio (NOT USED HERE; RETAINED FOR LOBSTER) |
---|
116 | !! |
---|
117 | !! Iron parameters |
---|
118 | INTEGER :: jiron !: iron scavenging submodel switch |
---|
119 | REAL(wp) :: xfe_mass !: iron atomic mass |
---|
120 | REAL(wp) :: xfe_sol !: aeolian iron solubility |
---|
121 | REAL(wp) :: xfe_sed !: sediment iron input |
---|
122 | REAL(wp) :: xLgT !: total ligand concentration (umol/m3) |
---|
123 | REAL(wp) :: xk_FeL !: dissociation constant for (Fe + L) |
---|
124 | REAL(wp) :: xk_sc_Fe !: scavenging rate of "free" iron |
---|
125 | !! |
---|
126 | !! Gravitational sinking |
---|
127 | REAL(wp) :: vsed !: detritus gravitational sinking rate |
---|
128 | REAL(wp) :: xhr !: coefficient for Martin et al. (1987) remineralisation |
---|
129 | !! |
---|
130 | !! Fast-sinking detritus parameters |
---|
131 | INTEGER :: jexport !: fast detritus remineralisation switch |
---|
132 | INTEGER :: jfdfate !: fate of fast detritus at seafloor switch |
---|
133 | INTEGER :: jrratio !: rain ratio switch |
---|
134 | INTEGER :: jocalccd !: CCD switch |
---|
135 | REAL(wp) :: xridg_r0 !: Ridgwell rain ratio coefficient |
---|
136 | REAL(wp) :: xfdfrac1 !: fast-sinking fraction of diatom nat. mort. losses |
---|
137 | REAL(wp) :: xfdfrac2 !: fast-sinking fraction of mesozooplankton mort. losses |
---|
138 | REAL(wp) :: xfdfrac3 !: fast-sinking fraction of diatom silicon grazing losses |
---|
139 | REAL(wp) :: xcaco3a !: polar (high latitude) CaCO3 fraction |
---|
140 | REAL(wp) :: xcaco3b !: equatorial (low latitude) CaCO3 fraction |
---|
141 | REAL(wp) :: xmassc !: organic C mass:mole ratio, C106 H175 O40 N16 P1 |
---|
142 | REAL(wp) :: xmassca !: calcium carbonate mass:mole ratio, CaCO3 |
---|
143 | REAL(wp) :: xmasssi !: biogenic silicon mass:mole ratio, (H2SiO3)n |
---|
144 | REAL(wp) :: xprotca !: calcium carbonate protection ratio |
---|
145 | REAL(wp) :: xprotsi !: biogenic silicon protection ratio |
---|
146 | REAL(wp) :: xfastc !: organic C remineralisation length scale |
---|
147 | REAL(wp) :: xfastca !: calcium carbonate dissolution length scale |
---|
148 | REAL(wp) :: xfastsi !: biogenic silicon dissolution length scale |
---|
149 | !! |
---|
150 | !! Benthos parameters |
---|
151 | INTEGER :: jorgben !: does organic detritus go to the benthos? |
---|
152 | INTEGER :: jinorgben !: does inorganic detritus go to the benthos? |
---|
153 | !! |
---|
154 | REAL(wp) :: xsedn !: organic nitrogen sediment remineralisation rate |
---|
155 | REAL(wp) :: xsedfe !: organic iron sediment remineralisation rate |
---|
156 | REAL(wp) :: xsedsi !: inorganic silicon sediment dissolution rate |
---|
157 | REAL(wp) :: xsedc !: organic carbon sediment remineralisation rate |
---|
158 | REAL(wp) :: xsedca !: inorganic carbon sediment dissolution rate |
---|
159 | REAL(wp) :: xburial !: burial rate of seafloor detritus |
---|
160 | !! |
---|
161 | !! River parameters |
---|
162 | INTEGER :: jriver_n !: riverine nitrogen? 0 = no, 1 = conc, 2 = flux |
---|
163 | INTEGER :: jriver_si !: riverine silicon? 0 = no, 1 = conc, 2 = flux |
---|
164 | INTEGER :: jriver_c !: riverine carbon? 0 = no, 1 = conc, 2 = flux |
---|
165 | INTEGER :: jriver_alk!: riverine alkalinity? 0 = no, 1 = conc, 2 = flux |
---|
166 | INTEGER :: jriver_dep!: depth river input added to? 1 = surface, >1 possible |
---|
167 | !! |
---|
168 | !! Miscellaneous |
---|
169 | REAL(wp) :: xsdiss !: diatom frustule dissolution rate |
---|
170 | !! |
---|
171 | !! Additional parameters |
---|
172 | INTEGER :: jpkb !: vertical layer for diagnostic of the vertical flux |
---|
173 | !! |
---|
174 | !! UKESM diagnostics |
---|
175 | INTEGER :: jdms !: include DMS diagnostics ? Jpalm (27-08-2014) |
---|
176 | INTEGER :: jdms_input !: use instant (0) or diel-average (1) inputs (AXY, 08/07/2015) |
---|
177 | INTEGER :: jdms_model !: choice of DMS model passed to atmosphere |
---|
178 | !! 1 = ANDR, 2 = SIMO, 3 = ARAN, 4 = HALL |
---|
179 | !! JPALM --19-12-2017 -- UM people need to tune the Anderson DMS scheme |
---|
180 | REAL(wp) :: dmsmin !: Anderson DMS scheme - DMS minimum value |
---|
181 | REAL(wp) :: dmscut !: Anderson DMS scheme - DMS cutoff value |
---|
182 | REAL(wp) :: dmsslp !: Anderson DMS scheme - DMS slope |
---|
183 | !! FOR UKESM |
---|
184 | REAL(wp) :: scl_chl !: scaling factor for tuned Chl passed to the UM |
---|
185 | INTEGER :: chl_out !: select Chl field exported and scaled for UM: |
---|
186 | !: 1- Surface Chl ; 2- MLD Chl |
---|
187 | !! |
---|
188 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: remdmp !: depth dependent damping coefficient of passive tracers |
---|
189 | !! |
---|
190 | !! AXY (27/07/10): add in indices for depth horizons (for sinking flux |
---|
191 | !! and seafloor iron inputs) |
---|
192 | INTEGER :: i0100, i0150, i0200, i0500, i1000, i1100 |
---|
193 | #if defined key_roam |
---|
194 | !! |
---|
195 | !! ROAM carbon, alkalinity and oxygen cycle parameters |
---|
196 | REAL(wp) :: xthetaphy !: oxygen evolution/consumption by phytoplankton |
---|
197 | REAL(wp) :: xthetazoo !: oxygen consumption by zooplankton |
---|
198 | REAL(wp) :: xthetanit !: oxygen consumption by nitrogen remineralisation |
---|
199 | REAL(wp) :: xthetarem !: oxygen consumption by carbon remineralisation |
---|
200 | REAL(wp) :: xo2min !: oxygen minimum concentration |
---|
201 | !! |
---|
202 | !! 3D fields of carbonate system parameters |
---|
203 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: f3_pH !: 3D pH |
---|
204 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: f3_h2co3 !: 3D carbonic acid |
---|
205 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: f3_hco3 !: 3D bicarbonate |
---|
206 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: f3_co3 !: 3D carbonate |
---|
207 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: f3_omcal !: 3D omega calcite |
---|
208 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: f3_omarg !: 3D omega aragonite |
---|
209 | !! |
---|
210 | !! 2D fields of calcium carbonate compensation depth |
---|
211 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: f2_ccd_cal !: 2D calcite CCD depth |
---|
212 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: f2_ccd_arg !: 2D aragonite CCD depth |
---|
213 | !! |
---|
214 | !! 2D fields of organic and inorganic material sedimented on the seafloor |
---|
215 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: zb_sed_n !: 2D organic nitrogen (before) |
---|
216 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: zn_sed_n !: 2D organic nitrogen (now) |
---|
217 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: za_sed_n !: 2D organic nitrogen (after) |
---|
218 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: zb_sed_fe !: 2D organic iron (before) |
---|
219 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: zn_sed_fe !: 2D organic iron (now) |
---|
220 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: za_sed_fe !: 2D organic iron (after) |
---|
221 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: zb_sed_si !: 2D inorganic silicon (before) |
---|
222 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: zn_sed_si !: 2D inorganic silicon (now) |
---|
223 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: za_sed_si !: 2D inorganic silicon (after) |
---|
224 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: zb_sed_c !: 2D organic carbon (before) |
---|
225 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: zn_sed_c !: 2D organic carbon (now) |
---|
226 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: za_sed_c !: 2D organic carbon (after) |
---|
227 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: zb_sed_ca !: 2D inorganic carbon (before) |
---|
228 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: zn_sed_ca !: 2D inorganic carbon (now) |
---|
229 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: za_sed_ca !: 2D inorganic carbon (after) |
---|
230 | !! |
---|
231 | !! 2D fields of temporally averaged properties for DMS calculations (AXY, 07/07/15) |
---|
232 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: zb_dms_chn !: 2D avg CHN (before) |
---|
233 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: zn_dms_chn !: 2D avg CHN (now) |
---|
234 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: za_dms_chn !: 2D avg CHN (after) |
---|
235 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: zb_dms_chd !: 2D avg CHD (before) |
---|
236 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: zn_dms_chd !: 2D avg CHD (now) |
---|
237 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: za_dms_chd !: 2D avg CHD (after) |
---|
238 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: zb_dms_mld !: 2D avg MLD (before) |
---|
239 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: zn_dms_mld !: 2D avg MLD (now) |
---|
240 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: za_dms_mld !: 2D avg MLD (after) |
---|
241 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: zb_dms_qsr !: 2D avg QSR (before) |
---|
242 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: zn_dms_qsr !: 2D avg QSR (now) |
---|
243 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: za_dms_qsr !: 2D avg QSR (after) |
---|
244 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: zb_dms_din !: 2D avg DIN (before) |
---|
245 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: zn_dms_din !: 2D avg DIN (now) |
---|
246 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: za_dms_din !: 2D avg DIN (after) |
---|
247 | !! |
---|
248 | !! 2D fields needing to be knows at first tstp for coupling with atm - UKEMS(Jpalm,14-06-2016) |
---|
249 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: zb_co2_flx !: 2D avg fx co2 (before) |
---|
250 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: zn_co2_flx !: 2D avg fx co2 (now) |
---|
251 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: za_co2_flx !: 2D avg fx co2 (after) |
---|
252 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: zb_dms_srf !: 2D avg sfr dms (before) |
---|
253 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: zn_dms_srf !: 2D avg sfr dms (now) |
---|
254 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: za_dms_srf !: 2D avg srf dms (after) |
---|
255 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: zn_chl_srf !: 2D avg srf chl (now) |
---|
256 | |
---|
257 | #endif |
---|
258 | |
---|
259 | !!---------------------------------------------------------------------- |
---|
260 | !! CCD parameter |
---|
261 | !!---------------------------------------------------------------------- |
---|
262 | !! |
---|
263 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: ocal_ccd !: CCD depth |
---|
264 | |
---|
265 | !!---------------------------------------------------------------------- |
---|
266 | !! Dust parameters |
---|
267 | !!---------------------------------------------------------------------- |
---|
268 | !! |
---|
269 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: dust !: dust parameter 1 |
---|
270 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: zirondep !! Fe deposition |
---|
271 | |
---|
272 | !!---------------------------------------------------------------------- |
---|
273 | !! River parameters |
---|
274 | !!---------------------------------------------------------------------- |
---|
275 | !! |
---|
276 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: riv_n !: riverine N |
---|
277 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: riv_si !: riverine Si |
---|
278 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: riv_c !: riverine C |
---|
279 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: riv_alk !: riverine alkalinity |
---|
280 | !! AXY (19/07/12): add this to permit river fluxes to be added below top box |
---|
281 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: friver_dep !: where river fluxes added |
---|
282 | |
---|
283 | #if defined key_roam |
---|
284 | !!---------------------------------------------------------------------- |
---|
285 | !! JPALM -- change hist_pco2 init |
---|
286 | !!---------------------------------------------------------------------- |
---|
287 | !! |
---|
288 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:) :: hist_pco2 !: pCO2 |
---|
289 | INTEGER :: co2_rec |
---|
290 | REAL(wp) :: co2_yinit, co2_yend |
---|
291 | #endif |
---|
292 | |
---|
293 | !!---------------------------------------------------------------------- |
---|
294 | !! JPALM -- PI CO2 key |
---|
295 | !!---------------------------------------------------------------------- |
---|
296 | !! |
---|
297 | #if defined key_axy_pi_co2 |
---|
298 | LOGICAL , PUBLIC :: lk_pi_co2 = .TRUE. !: PI xCO2 used |
---|
299 | #else |
---|
300 | LOGICAL , PUBLIC :: lk_pi_co2 = .FALSE. !: PI xCO2 unused |
---|
301 | #endif |
---|
302 | |
---|
303 | !!---------------------------------------------------------------------- |
---|
304 | !! Optical parameters |
---|
305 | !!---------------------------------------------------------------------- |
---|
306 | !! |
---|
307 | REAL(wp) :: xkr0 !: water coefficient absorption in red (NAMELIST) |
---|
308 | REAL(wp) :: xkg0 !: water coefficient absorption in green (NAMELIST) |
---|
309 | REAL(wp) :: xkrp !: pigment coefficient absorption in red (NAMELIST) |
---|
310 | REAL(wp) :: xkgp !: pigment coefficient absorption in green (NAMELIST) |
---|
311 | REAL(wp) :: xlr !: exposant for pigment absorption in red (NAMELIST) |
---|
312 | REAL(wp) :: xlg !: exposant for pigment absorption in green (NAMELIST) |
---|
313 | REAL(wp) :: rpig !: chla/chla+phea ratio (NAMELIST) |
---|
314 | |
---|
315 | INTEGER , ALLOCATABLE, SAVE, DIMENSION(:,:) :: neln !: number of levels in the euphotic layer |
---|
316 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: xze !: euphotic layer depth |
---|
317 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: xpar !: par (photosynthetic available radiation) |
---|
318 | |
---|
319 | !!---------------------------------------------------------------------- |
---|
320 | !! Sediment parameters |
---|
321 | !! |
---|
322 | !! AXY (16/01/12): these parameters were originally part of the pre- |
---|
323 | !! cursor model on which MEDUSA's code was grounded; |
---|
324 | !! they do not relate to the sediment/benthos submodel |
---|
325 | !! added as part of the ROAM project; they have only |
---|
326 | !! been retained because they are distributed through |
---|
327 | !! MEDUSA and require a proper clean-up to purge |
---|
328 | !!---------------------------------------------------------------------- |
---|
329 | !! |
---|
330 | REAL(wp) :: sedlam !: time coefficient of POC remineralization in sediments |
---|
331 | REAL(wp) :: sedlostpoc !: ??? |
---|
332 | REAL(wp) :: areacot !: ??? |
---|
333 | |
---|
334 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: dminl !: fraction of sinking POC released in sediments |
---|
335 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: dmin3 !: fraction of sinking POC released at each level |
---|
336 | |
---|
337 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: sedpocb !: mass of POC in sediments |
---|
338 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: sedpocn !: mass of POC in sediments |
---|
339 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: sedpoca !: mass of POC in sediments |
---|
340 | |
---|
341 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: fbodf !: rapid sinking particles |
---|
342 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: fbods !: rapid sinking particles |
---|
343 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: fbodn !: rapid sinking particles |
---|
344 | |
---|
345 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: ffln !: |
---|
346 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: fflf !: |
---|
347 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: ffls !: |
---|
348 | |
---|
349 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: cmask !: ??? |
---|
350 | |
---|
351 | !!---------------------------------------------------------------------- |
---|
352 | !! NEMO/TOP 3.3 , NEMO Consortium (2010) |
---|
353 | !! $Id$ |
---|
354 | !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) |
---|
355 | !!---------------------------------------------------------------------- |
---|
356 | CONTAINS |
---|
357 | |
---|
358 | INTEGER FUNCTION sms_medusa_alloc() |
---|
359 | !!---------------------------------------------------------------------- |
---|
360 | !! *** ROUTINE sms_medusa_alloc *** |
---|
361 | !!---------------------------------------------------------------------- |
---|
362 | USE lib_mpp , ONLY: ctl_warn |
---|
363 | INTEGER :: ierr(8) ! Local variables |
---|
364 | !!---------------------------------------------------------------------- |
---|
365 | ierr(:) = 0 |
---|
366 | ! |
---|
367 | #if defined key_medusa |
---|
368 | !* depth-dependent damping coefficient |
---|
369 | ALLOCATE( remdmp(jpk,jp_medusa), STAT=ierr(1) ) |
---|
370 | # if defined key_roam |
---|
371 | !* 2D and 3D fields of carbonate system parameters |
---|
372 | ALLOCATE( f2_ccd_cal(jpi,jpj) , f2_ccd_arg(jpi,jpj) , & |
---|
373 | & f3_pH(jpi,jpj,jpk) , f3_h2co3(jpi,jpj,jpk), & |
---|
374 | & f3_hco3(jpi,jpj,jpk) , f3_co3(jpi,jpj,jpk) , & |
---|
375 | & f3_omcal(jpi,jpj,jpk), f3_omarg(jpi,jpj,jpk), STAT=ierr(2) ) |
---|
376 | !* 2D fields of organic and inorganic material sedimented on the seafloor |
---|
377 | ALLOCATE( zb_sed_n(jpi,jpj) , zn_sed_n(jpi,jpj) , & |
---|
378 | & za_sed_n(jpi,jpj) , & |
---|
379 | & zb_sed_fe(jpi,jpj) , zn_sed_fe(jpi,jpj) , & |
---|
380 | & za_sed_fe(jpi,jpj) , & |
---|
381 | & zb_sed_si(jpi,jpj) , zn_sed_si(jpi,jpj) , & |
---|
382 | & za_sed_si(jpi,jpj) , & |
---|
383 | & zb_sed_c(jpi,jpj) , zn_sed_c(jpi,jpj) , & |
---|
384 | & za_sed_c(jpi,jpj) , & |
---|
385 | & zb_sed_ca(jpi,jpj) , zn_sed_ca(jpi,jpj) , & |
---|
386 | & za_sed_ca(jpi,jpj) , STAT=ierr(3) ) |
---|
387 | !* 2D fields of temporally averaged properties for DMS calculations (AXY, 07/07/15) |
---|
388 | ALLOCATE( zb_dms_chn(jpi,jpj) , zn_dms_chn(jpi,jpj) , & |
---|
389 | & za_dms_chn(jpi,jpj) , & |
---|
390 | & zb_dms_chd(jpi,jpj) , zn_dms_chd(jpi,jpj) , & |
---|
391 | & za_dms_chd(jpi,jpj) , & |
---|
392 | & zb_dms_mld(jpi,jpj) , zn_dms_mld(jpi,jpj) , & |
---|
393 | & za_dms_mld(jpi,jpj) , & |
---|
394 | & zb_dms_qsr(jpi,jpj) , zn_dms_qsr(jpi,jpj) , & |
---|
395 | & za_dms_qsr(jpi,jpj) , & |
---|
396 | & zb_dms_din(jpi,jpj) , zn_dms_din(jpi,jpj) , & |
---|
397 | & za_dms_din(jpi,jpj) , STAT=ierr(4) ) |
---|
398 | !* 2D fields needing to be knows at first tstp for coupling with atm - |
---|
399 | !UKEMSi (Jpalm,14-06-2016) |
---|
400 | ALLOCATE( zb_co2_flx(jpi,jpj) , zn_co2_flx(jpi,jpj) , & |
---|
401 | & za_co2_flx(jpi,jpj) , & |
---|
402 | & zb_dms_srf(jpi,jpj) , zn_dms_srf(jpi,jpj) , & |
---|
403 | & za_dms_srf(jpi,jpj) , zn_chl_srf(jpi,jpj) , STAT=ierr(5) ) |
---|
404 | # endif |
---|
405 | !* 2D fields of miscellaneous parameters |
---|
406 | ALLOCATE( ocal_ccd(jpi,jpj) , dust(jpi,jpj) , & |
---|
407 | & zirondep(jpi,jpj) , & |
---|
408 | & riv_n(jpi,jpj) , & |
---|
409 | & riv_si(jpi,jpj) , riv_c(jpi,jpj) , & |
---|
410 | & riv_alk(jpi,jpj) , friver_dep(jpk,jpk) , STAT=ierr(6) ) |
---|
411 | !* 2D and 3D fields of light parameters |
---|
412 | ALLOCATE( neln(jpi,jpj) , xze(jpi,jpj) , & |
---|
413 | & xpar(jpi,jpj,jpk) , STAT=ierr(7) ) |
---|
414 | !* 2D and 3D fields of sediment-associated parameters |
---|
415 | ALLOCATE( dminl(jpi,jpj) , dmin3(jpi,jpj,jpk) , & |
---|
416 | & sedpocb(jpi,jpj) , sedpocn(jpi,jpj) , & |
---|
417 | & sedpoca(jpi,jpj) , fbodn(jpi,jpj) , & |
---|
418 | & fbodf(jpi,jpj) , fbods(jpi,jpj) , & |
---|
419 | & ffln(jpi,jpj,jpk) , fflf(jpi,jpj,jpk) , & |
---|
420 | & ffls(jpi,jpj,jpk) , cmask(jpi,jpj) , STAT=ierr(8) ) |
---|
421 | #endif |
---|
422 | ! |
---|
423 | sms_medusa_alloc = MAXVAL( ierr ) |
---|
424 | ! |
---|
425 | IF( sms_medusa_alloc /= 0 ) CALL ctl_warn('sms_medusa_alloc: failed to allocate arrays') |
---|
426 | ! |
---|
427 | END FUNCTION sms_medusa_alloc |
---|
428 | |
---|
429 | #else |
---|
430 | !!---------------------------------------------------------------------- |
---|
431 | !! Empty module : NO MEDUSA model |
---|
432 | !!---------------------------------------------------------------------- |
---|
433 | #endif |
---|
434 | |
---|
435 | !!====================================================================== |
---|
436 | END MODULE sms_medusa |
---|