source: branches/NERC/dev_r5518_GO6_under_ice_relax/NEMOGCM/NEMO/TOP_SRC/MEDUSA/bio_medusa_init.F90 @ 10047

Last change on this file since 10047 was 10047, checked in by jpalmier, 3 years ago

merge with GO6_package_branch 9385-10020 ; plus debug OMIP_DIC

File size: 35.8 KB
Line 
1MODULE bio_medusa_init_mod
2   !!======================================================================
3   !!                         ***  MODULE bio_medusa_init  ***
4   !! Initialisation for TRC_BIO_MEDUSA
5   !!======================================================================
6   !! History :
7   !!   -   ! 2017-04 (M. Stringer)        Code taken from trcbio_medusa.F90
8   !!   -   ! 2017-08 (A. Yool)            Add slow-sinking detrius variables
9   !!   -   ! 2018-08 (A. Yool)            add OMIP preindustrial DIC
10   !!----------------------------------------------------------------------
11#if defined key_medusa
12   !!----------------------------------------------------------------------
13   !!                                                   MEDUSA bio-model
14   !!----------------------------------------------------------------------
15
16   IMPLICIT NONE
17   PRIVATE
18     
19   PUBLIC   bio_medusa_init    ! Called in trcbio_medusa.F90
20
21   !!----------------------------------------------------------------------
22   !! NEMO/TOP 2.0 , LOCEAN-IPSL (2007)
23   !! $Id$
24   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
25   !!----------------------------------------------------------------------
26
27CONTAINS
28
29   SUBROUTINE bio_medusa_init( kt )
30      !!---------------------------------------------------------------------
31      !!                     ***  ROUTINE bio_medusa_init  ***
32      !! This called from TRC_BIO_MEDUSA and
33      !!  - zeroes arrays used in TRC_BIO_MEDUSA
34      !!  - sets up the arrays needed for the diagnostics
35      !!----------------------------------------------------------------------
36      USE bio_medusa_mod
37      USE par_oce,           ONLY: jpi, jpj, jpk
38      USE sms_medusa,        ONLY: jdms
39      USE trc,               ONLY: ln_diatrc, med_diag, nittrc000 
40      USE in_out_manager,    ONLY: lwp, numout
41
42      USE iom,               ONLY: lk_iomput
43      USE trcnam_medusa,     ONLY: trc_nam_iom_medusa
44
45      !! time (integer timestep)
46      INTEGER, INTENT( in ) ::    kt
47
48      IF( ln_diatrc ) THEN
49         IF (lwp) write (numout,*) 'Diagnostics are now ALL through XIOS (key_xios)'
50         IF (lwp) write (numout,*) 'No more key_diatrc anymore.'
51      ENDIF
52
53      !!----------------------------------------------------------------------
54      !! Zero fast-sinking detritus 2D fields
55      !!----------------------------------------------------------------------
56      !!
57      ffastn(:,:)  = 0.0        !! organic nitrogen
58      ffastsi(:,:) = 0.0        !! biogenic silicon
59      ffastfe(:,:) = 0.0        !! organic iron
60      ffastc(:,:)  = 0.0        !! organic carbon
61      ffastca(:,:) = 0.0        !! biogenic calcium carbonate
62      !!
63      fsedn(:,:)   = 0.0        !! Seafloor flux of N
64      fsedsi(:,:)  = 0.0        !! Seafloor flux of Si
65      fsedfe(:,:)  = 0.0        !! Seafloor flux of Fe
66      fsedc(:,:)   = 0.0        !! Seafloor flux of C
67      fsedca(:,:)  = 0.0        !! Seafloor flux of CaCO3
68      !!
69      fregenfast(:,:)   = 0.0   !! integrated  N regeneration (fast detritus)
70      fregenfastsi(:,:) = 0.0   !! integrated Si regeneration (fast detritus)
71# if defined key_roam
72      fregenfastc(:,:)  = 0.0   !! integrated  C regeneration (fast detritus)
73# endif
74      !!
75      fccd(:,:)    = 0.0        !! last depth level before CCD
76
77      !!----------------------------------------------------------------------
78      !! blank nutrient/flux inventories
79      !!----------------------------------------------------------------------
80      !!
81      fflx_n(:,:)  = 0.0        !! nitrogen flux total
82      fflx_si(:,:) = 0.0        !! silicon  flux total
83      fflx_fe(:,:) = 0.0        !! iron     flux total
84      fifd_n(:,:)  = 0.0        !! nitrogen fast detritus production
85      fifd_si(:,:) = 0.0        !! silicon  fast detritus production
86      fifd_fe(:,:) = 0.0        !! iron     fast detritus production
87      fofd_n(:,:)  = 0.0        !! nitrogen fast detritus remineralisation
88      fofd_si(:,:) = 0.0        !! silicon  fast detritus remineralisation
89      fofd_fe(:,:) = 0.0        !! iron     fast detritus remineralisation
90# if defined key_roam
91      fflx_c(:,:)  = 0.0        !! carbon     flux total
92      fflx_a(:,:)  = 0.0        !! alkalinity flux total
93      fflx_o2(:,:) = 0.0        !! oxygen     flux total
94      ftot_c(:,:)  = 0.0        !! carbon     inventory
95      ftot_a(:,:)  = 0.0        !! alkalinity inventory
96      ftot_o2(:,:) = 0.0        !! oxygen     inventory
97      fifd_c(:,:)  = 0.0        !! carbon     fast detritus production
98! I don't think fifd_a or fifd_o2 are used - marc 11/5/17
99!      fifd_a(:,:)  = 0.0        !! alkalinity fast detritus production
100!      fifd_o2(:,:) = 0.0        !! oxygen     fast detritus production
101      fofd_c(:,:)  = 0.0        !! carbon     fast detritus remineralisation
102! I don't think fofd_a or fofd_o2 are used - marc 11/5/17
103!      fofd_a(:,:)  = 0.0        !! alkalinity fast detritus remineralisation
104!      fofd_o2(:,:) = 0.0        !! oxygen     fast detritus remineralisation
105      !!
106      fnit_prod(:,:) = 0.0      !! (organic)   nitrogen production
107      fnit_cons(:,:) = 0.0      !! (organic)   nitrogen consumption
108      fsil_prod(:,:) = 0.0      !! (inorganic) silicon production
109      fsil_cons(:,:) = 0.0      !! (inorganic) silicon consumption
110      fcar_prod(:,:) = 0.0      !! (organic)   carbon production
111      fcar_cons(:,:) = 0.0      !! (organic)   carbon consumption
112      !!
113      foxy_prod(:,:) = 0.0      !! oxygen production
114      foxy_cons(:,:) = 0.0      !! oxygen consumption
115      foxy_anox(:,:) = 0.0      !! unrealised oxygen consumption
116      !!
117# endif
118      ftot_n(:,:)   = 0.0       !! N inventory
119      ftot_si(:,:)  = 0.0       !! Si inventory
120      ftot_fe(:,:)  = 0.0       !! Fe inventory
121      ftot_pn(:,:)  = 0.0       !! integrated non-diatom phytoplankton
122      ftot_pd(:,:)  = 0.0       !! integrated diatom     phytoplankton
123      ftot_zmi(:,:) = 0.0       !! integrated microzooplankton
124      ftot_zme(:,:) = 0.0       !! integrated mesozooplankton
125      ftot_det(:,:) = 0.0       !! integrated slow detritus, nitrogen
126      ftot_dtc(:,:) = 0.0       !! integrated slow detritus, carbon
127      !!
128      fzmi_i(:,:)  = 0.0        !! material grazed by microzooplankton
129      fzmi_o(:,:)  = 0.0        !! ... sum of fate of this material
130      fzme_i(:,:)  = 0.0        !! material grazed by  mesozooplankton
131      fzme_o(:,:)  = 0.0        !! ... sum of fate of this material
132      !!
133      f_sbenin_n(:,:)  = 0.0    !! slow detritus N  -> benthic pool
134      f_sbenin_fe(:,:) = 0.0    !! slow detritus Fe -> benthic pool
135      f_sbenin_c(:,:)  = 0.0    !! slow detritus C  -> benthic pool
136      f_fbenin_n(:,:)  = 0.0    !! fast detritus N  -> benthic pool
137      f_fbenin_fe(:,:) = 0.0    !! fast detritus Fe -> benthic pool
138      f_fbenin_si(:,:) = 0.0    !! fast detritus Si -> benthic pool
139      f_fbenin_c(:,:)  = 0.0    !! fast detritus C  -> benthic pool
140      f_fbenin_ca(:,:) = 0.0    !! fast detritus Ca -> benthic pool
141      !!
142      f_benout_n(:,:)  = 0.0    !! benthic N  pool  -> dissolved
143      f_benout_fe(:,:) = 0.0    !! benthic Fe pool  -> dissolved
144      f_benout_si(:,:) = 0.0    !! benthic Si pool  -> dissolved
145      f_benout_c(:,:)  = 0.0    !! benthic C  pool  -> dissolved
146      f_benout_ca(:,:) = 0.0    !! benthic Ca pool  -> dissolved
147      !!
148      f_benout_lyso_ca(:,:) = 0.0 !! benthic Ca pool  -> dissolved (when it shouldn't!)
149      !!
150      f_runoff(:,:)  = 0.0      !! riverine runoff
151      f_riv_n(:,:)   = 0.0      !! riverine N   input
152      f_riv_si(:,:)  = 0.0      !! riverine Si  input
153      f_riv_c(:,:)   = 0.0      !! riverine C   input
154      f_riv_alk(:,:) = 0.0      !! riverine alk input
155      !!
156      !! Jpalm -- 06-03-2017 -- Forgotten var to init
157      f_omarg(:,:) = 0.0        !!
158      f_omcal(:,:) = 0.0 
159      xFree(:,:) = 0.0          !! state variables for iron-ligand system
160      fcomm_resp(:,:) = 0.0 
161      fprn_ml(:,:) = 0.0        !! mixed layer PP diagnostics
162      fprd_ml(:,:) = 0.0        !! mixed layer PP diagnostics
163      !! AXY (16/08/17)
164      fchl_ml(:,:) = 0.0   !! mixed layer chlorophyll diagnostics
165      !!
166      fslownflux(:,:) = 0.0 
167      fslowcflux(:,:) = 0.0 
168      !!
169      !! JPALM -- 21-09-2017 -- needed to debug air-sea carb
170      f_xco2a(:,:)  = 0.0
171      f_pco2w(:,:)  = 0.0
172      f_ph(:,:)     = 0.0
173      f_kw660(:,:)  = 0.0
174      ztmp(:,:)  = 0.0
175      zsal(:,:)  = 0.0
176      zalk(:,:)  = 0.0
177      zdic(:,:)  = 0.0
178      zsil(:,:)  = 0.0
179# if defined key_mocsy
180      ! zpho is only defined if key_mocsy
181      ! is active, so we must protect this
182      ! initialisation accordingly.
183      zpho(:,:)  = 0.0
184# endif
185      f_co2flux(:,:)  = 0.0 
186      f_pco2atm(:,:)  = 0.0
187      f_h2co3(:,:)    = 0.0
188      f_hco3(:,:)     = 0.0
189      f_co3(:,:)      = 0.0
190      f_omarg(:,:)    = 0.0
191      f_omcal(:,:)    = 0.0
192# if defined key_omip_dic
193      zomd(:,:)          = 0.0
194      f_pi_xco2a(:,:)    = 0.0
195      f_pi_pco2w(:,:)    = 0.0
196      f_pi_ph(:,:)       = 0.0
197      f_pi_co2flux(:,:)  = 0.0 
198      f_pi_pco2atm(:,:)  = 0.0
199      f_pi_h2co3(:,:)    = 0.0
200      f_pi_hco3(:,:)     = 0.0
201      f_pi_co3(:,:)      = 0.0
202      f_pi_omarg(:,:)    = 0.0
203      f_pi_omcal(:,:)    = 0.0
204# endif     
205      !!
206      !! AXY (08/08/17): zero slow detritus fluxes
207      fslowsink(:,:)  = 0.0
208# if defined key_roam
209      fslowsinkc(:,:) = 0.0
210# endif     
211      !!
212      !! allocate and initiate 2D diag
213      !! -----------------------------
214      !! Juju :: add kt condition !!
215      IF ( lk_iomput ) THEN 
216
217         !! initialise iom_use test
218         IF ( kt == nittrc000 )   CALL trc_nam_iom_medusa 
219
220         ALLOCATE( zw2d(1:jpi, 1:jpj) )
221         zw2d(:,:)      = 0.0   !!
222         IF ( med_diag%PRN%dgsave ) THEN
223            ALLOCATE( fprn2d(1:jpi, 1:jpj) )
224            fprn2d(:,:)      = 0.0 !!
225         ENDIF
226         IF ( med_diag%MPN%dgsave ) THEN
227            ALLOCATE( fdpn2d(1:jpi, 1:jpj) )
228            fdpn2d(:,:)      = 0.0 !!
229         ENDIF
230         IF ( med_diag%PRD%dgsave ) THEN
231            ALLOCATE( fprd2d(1:jpi, 1:jpj) )
232            fprd2d(:,:)      = 0.0 !!
233         ENDIF
234         IF( med_diag%MPD%dgsave ) THEN
235            ALLOCATE( fdpd2d(1:jpi, 1:jpj) )
236            fdpd2d(:,:)      = 0.0 !!
237         ENDIF
238         IF( med_diag%OPAL%dgsave ) THEN
239            ALLOCATE( fprds2d(1:jpi, 1:jpj) )
240            fprds2d(:,:)      = 0.0 !!
241         ENDIF
242         IF( med_diag%OPALDISS%dgsave ) THEN
243            ALLOCATE( fsdiss2d(1:jpi, 1:jpj) )
244            fsdiss2d(:,:)      = 0.0 !!
245         ENDIF
246         IF( med_diag%GMIPn%dgsave ) THEN
247            ALLOCATE( fgmipn2d(1:jpi, 1:jpj) )
248            fgmipn2d(:,:)      = 0.0 !!
249         ENDIF
250         IF( med_diag%GMID%dgsave ) THEN
251            ALLOCATE( fgmid2d(1:jpi, 1:jpj) )
252            fgmid2d(:,:)      = 0.0 !!
253         ENDIF
254         IF( med_diag%MZMI%dgsave ) THEN
255            ALLOCATE( fdzmi2d(1:jpi, 1:jpj) )
256            fdzmi2d(:,:)      = 0.0 !!
257         ENDIF
258         IF( med_diag%GMEPN%dgsave ) THEN
259            ALLOCATE( fgmepn2d(1:jpi, 1:jpj) )
260            fgmepn2d(:,:)      = 0.0 !!
261         ENDIF
262         IF( med_diag%GMEPD%dgsave ) THEN
263            ALLOCATE( fgmepd2d(1:jpi, 1:jpj) )
264            fgmepd2d(:,:)      = 0.0 !!
265         ENDIF
266         IF( med_diag%GMEZMI%dgsave ) THEN
267            ALLOCATE( fgmezmi2d(1:jpi, 1:jpj) )
268            fgmezmi2d(:,:)      = 0.0 !!
269         ENDIF
270         IF( med_diag%GMED%dgsave ) THEN
271            ALLOCATE( fgmed2d(1:jpi, 1:jpj) )
272            fgmed2d(:,:)      = 0.0 !!
273         ENDIF
274         IF( med_diag%MZME%dgsave ) THEN
275            ALLOCATE( fdzme2d(1:jpi, 1:jpj) )
276            fdzme2d(:,:)      = 0.0 !!
277         ENDIF
278         IF( med_diag%DETN%dgsave ) THEN
279            ALLOCATE( fslown2d(1:jpi, 1:jpj) )
280            fslown2d(:,:)      = 0.0 !!
281         ENDIF
282         IF( med_diag%MDET%dgsave ) THEN
283            ALLOCATE( fdd2d(1:jpi, 1:jpj) )
284            fdd2d(:,:)      = 0.0 !!
285         ENDIF     
286         IF( med_diag%AEOLIAN%dgsave ) THEN
287            ALLOCATE( ffetop2d(1:jpi, 1:jpj) )
288            ffetop2d(:,:)      = 0.0 !!
289         ENDIF
290         IF( med_diag%BENTHIC%dgsave ) THEN
291            ALLOCATE( ffebot2d(1:jpi, 1:jpj) )
292            ffebot2d(:,:)      = 0.0 !!
293         ENDIF
294         IF( med_diag%SCAVENGE%dgsave ) THEN
295            ALLOCATE( ffescav2d(1:jpi, 1:jpj) )
296            ffescav2d(:,:)      = 0.0 !!
297         ENDIF
298         IF( med_diag%PN_JLIM%dgsave ) THEN
299            ALLOCATE( fjln2d(1:jpi, 1:jpj) )
300            fjln2d(:,:)      = 0.0 !!
301         ENDIF
302         IF( med_diag%PN_NLIM%dgsave ) THEN
303            ALLOCATE( fnln2d(1:jpi, 1:jpj) )
304            fnln2d(:,:)      = 0.0 !!
305         ENDIF
306         IF( med_diag%PN_FELIM%dgsave ) THEN
307            ALLOCATE( ffln2d(1:jpi, 1:jpj) )
308            ffln2d(:,:)      = 0.0 !!
309         ENDIF
310         IF( med_diag%PD_JLIM%dgsave ) THEN
311            ALLOCATE( fjld2d(1:jpi, 1:jpj) )
312            fjld2d(:,:)      = 0.0 !!
313         ENDIF
314         IF( med_diag%PD_NLIM%dgsave ) THEN
315            ALLOCATE( fnld2d(1:jpi, 1:jpj) )
316            fnld2d(:,:)      = 0.0 !!
317         ENDIF
318         IF( med_diag%PD_FELIM%dgsave ) THEN
319            ALLOCATE( ffld2d(1:jpi, 1:jpj) )
320            ffld2d(:,:)      = 0.0 !!
321         ENDIF
322         IF( med_diag%PD_SILIM%dgsave ) THEN
323            ALLOCATE( fsld2d2(1:jpi, 1:jpj) )
324            fsld2d2(:,:)      = 0.0 !!
325         ENDIF
326         IF( med_diag%PDSILIM2%dgsave ) THEN
327            ALLOCATE( fsld2d(1:jpi, 1:jpj) )
328            fsld2d(:,:)      = 0.0 !!
329         ENDIF
330!!
331!! skip SDT_XXXX diagnostics here
332!!
333         IF( med_diag%TOTREG_N%dgsave ) THEN
334            ALLOCATE( fregen2d(1:jpi, 1:jpj) )
335            fregen2d(:,:)      = 0.0 !!
336         ENDIF
337         IF( med_diag%TOTRG_SI%dgsave ) THEN
338            ALLOCATE( fregensi2d(1:jpi, 1:jpj) )
339            fregensi2d(:,:)      = 0.0 !!
340         ENDIF
341!!
342!! skip REG_XXXX diagnostics here
343!!
344         IF( med_diag%FASTN%dgsave ) THEN
345            ALLOCATE( ftempn2d(1:jpi, 1:jpj) )
346            ftempn2d(:,:)      = 0.0 !!
347         ENDIF
348         IF( med_diag%FASTSI%dgsave ) THEN
349            ALLOCATE( ftempsi2d(1:jpi, 1:jpj) )
350            ftempsi2d(:,:)      = 0.0 !!
351         ENDIF
352         IF( med_diag%FASTFE%dgsave ) THEN
353            ALLOCATE( ftempfe2d(1:jpi, 1:jpj) )
354            ftempfe2d(:,:)      = 0.0 !!
355         ENDIF
356         IF( med_diag%FASTC%dgsave ) THEN
357            ALLOCATE( ftempc2d(1:jpi, 1:jpj) )
358            ftempc2d(:,:)      = 0.0 !!
359         ENDIF
360         IF( med_diag%FASTCA%dgsave ) THEN
361            ALLOCATE( ftempca2d(1:jpi, 1:jpj) )
362            ftempca2d(:,:)      = 0.0 !!
363         ENDIF     
364!!
365!! skip FDT_XXXX, RG_XXXXF, FDS_XXXX, RGS_XXXXF diagnostics here
366!!
367         IF( med_diag%REMINN%dgsave ) THEN
368            ALLOCATE( freminn2d(1:jpi, 1:jpj) )
369            freminn2d(:,:)      = 0.0 !!
370         ENDIF
371         IF( med_diag%REMINSI%dgsave ) THEN
372            ALLOCATE( freminsi2d(1:jpi, 1:jpj) )
373            freminsi2d(:,:)      = 0.0 !!
374         ENDIF
375         IF( med_diag%REMINFE%dgsave ) THEN
376            ALLOCATE( freminfe2d(1:jpi, 1:jpj) )
377            freminfe2d(:,:)      = 0.0 !!
378         ENDIF
379         IF( med_diag%REMINC%dgsave ) THEN
380            ALLOCATE( freminc2d(1:jpi, 1:jpj) )
381            freminc2d(:,:)      = 0.0 !!
382         ENDIF
383         IF( med_diag%REMINCA%dgsave ) THEN
384            ALLOCATE( freminca2d(1:jpi, 1:jpj) )
385            freminca2d(:,:)      = 0.0 !!
386         ENDIF
387
388# if defined key_roam
389!!
390!! skip SEAFLRXX, MED_XXXX, INTFLX_XX, INT_XX, ML_XXX, OCAL_XXX, FE_XXXX, MED_XZE, WIND diagnostics here
391!!
392         IF( med_diag%RR_0100%dgsave ) THEN
393            ALLOCATE( ffastca2d(1:jpi, 1:jpj) )
394            ffastca2d(:,:)      = 0.0 !!
395         ENDIF
396
397         IF( med_diag%ATM_PCO2%dgsave ) THEN
398            ALLOCATE( f_pco2a2d(1:jpi, 1:jpj) )
399            f_pco2a2d(:,:)      = 0.0 !!
400         ENDIF
401#  if defined key_omip_dic         
402         IF( med_diag%PI_ATM_PCO2%dgsave ) THEN
403            ALLOCATE( f_pi_pco2a_2d(1:jpi, 1:jpj) )
404            f_pi_pco2a_2d(:,:)      = 0.0 !!
405         ENDIF
406#  endif
407!!
408!! skip OCN_PH diagnostic here
409!!
410#  if defined key_omip_dic         
411         IF( med_diag%PI_OCN_PH%dgsave ) THEN
412            ALLOCATE( f_pi_ph_2d(1:jpi, 1:jpj) )
413            f_pi_ph_2d(:,:)   = 0.0 !!
414         ENDIF
415#  endif
416         IF( med_diag%OCN_PCO2%dgsave ) THEN
417            ALLOCATE( f_pco2w2d(1:jpi, 1:jpj) )
418            f_pco2w2d(:,:)      = 0.0 !!
419         ENDIF
420#  if defined key_omip_dic         
421         IF( med_diag%PI_OCN_PCO2%dgsave ) THEN
422            ALLOCATE( f_pi_pco2w_2d(1:jpi, 1:jpj) )
423            f_pi_pco2w_2d(:,:)   = 0.0 !!
424         ENDIF
425#  endif
426!!
427!! skip OCNH2CO3, OCN_HCO3, OCN_CO3 diagnostics here
428!!
429#  if defined key_omip_dic         
430         IF( med_diag%PI_OCNH2CO3%dgsave ) THEN
431            ALLOCATE( f_pi_h2co3_2d(1:jpi, 1:jpj) )
432            f_pi_h2co3_2d(:,:)   = 0.0 !!
433         ENDIF
434         IF( med_diag%PI_OCN_HCO3%dgsave ) THEN
435            ALLOCATE( f_pi_hco3_2d(1:jpi, 1:jpj) )
436            f_pi_hco3_2d(:,:)   = 0.0 !!
437         ENDIF
438         IF( med_diag%PI_OCN_CO3%dgsave ) THEN
439            ALLOCATE( f_pi_co3_2d(1:jpi, 1:jpj) )
440            f_pi_co3_2d(:,:)   = 0.0 !!
441         ENDIF
442#  endif
443
444         IF( med_diag%CO2FLUX%dgsave ) THEN
445            ALLOCATE( f_co2flux2d(1:jpi, 1:jpj) )
446            f_co2flux2d(:,:)      = 0.0 !!
447         ENDIF
448#  if defined key_omip_dic         
449         IF( med_diag%PI_CO2FLUX%dgsave ) THEN
450            ALLOCATE( f_pi_co2flux2d(1:jpi, 1:jpj) )
451            f_pi_co2flux2d(:,:)   = 0.0 !!
452         ENDIF
453#  endif
454!!
455!! skip OM_XXX diagnostics here
456!!
457#  if defined key_omip_dic         
458         IF( med_diag%PI_OM_CAL%dgsave ) THEN
459            ALLOCATE( f_pi_omcal_2d(1:jpi, 1:jpj) )
460            f_pi_omcal_2d(:,:)   = 0.0 !!
461         ENDIF
462         IF( med_diag%PI_OM_ARG%dgsave ) THEN
463            ALLOCATE( f_pi_omarg_2d(1:jpi, 1:jpj) )
464            f_pi_omarg_2d(:,:)   = 0.0 !!
465         ENDIF
466#  endif
467         IF( med_diag%TCO2%dgsave ) THEN
468            ALLOCATE( f_TDIC2d(1:jpi, 1:jpj) )
469            f_TDIC2d(:,:)      = 0.0 !!
470         ENDIF
471#  if defined key_omip_dic         
472         IF( med_diag%PI_TCO2%dgsave ) THEN
473            ALLOCATE( f_pi_TDIC2d(1:jpi, 1:jpj) )
474            f_pi_TDIC2d(:,:)      = 0.0 !!
475         ENDIF
476#  endif
477         IF( med_diag%TALK%dgsave ) THEN
478            ALLOCATE( f_TALK2d(1:jpi, 1:jpj) )
479            f_TALK2d(:,:)      = 0.0 !!
480         ENDIF
481         IF( med_diag%KW660%dgsave ) THEN
482            ALLOCATE( f_kw6602d(1:jpi, 1:jpj) )
483            f_kw6602d(:,:)      = 0.0 !!
484         ENDIF
485         IF( med_diag%ATM_PP0%dgsave ) THEN
486            ALLOCATE( f_pp02d(1:jpi, 1:jpj) )
487            f_pp02d(:,:)      = 0.0 !!
488         ENDIF
489         IF( med_diag%O2FLUX%dgsave ) THEN
490            ALLOCATE( f_o2flux2d(1:jpi, 1:jpj) )
491            f_o2flux2d(:,:)      = 0.0 !!
492         ENDIF
493         IF( med_diag%O2SAT%dgsave ) THEN
494            ALLOCATE( f_o2sat2d(1:jpi, 1:jpj) )
495            f_o2sat2d(:,:)      = 0.0 !!
496         ENDIF 
497!!
498!! skip XXX_CCD diagnostics here
499!!
500         IF( med_diag%SFR_OCAL%dgsave ) THEN
501            ALLOCATE( sfr_ocal2d(1:jpi, 1:jpj) )
502            sfr_ocal2d(:,:)      = 0.0 !!
503         ENDIF
504         IF( med_diag%SFR_OARG%dgsave ) THEN
505            ALLOCATE( sfr_oarg2d(1:jpi, 1:jpj) )
506            sfr_oarg2d(:,:)      = 0.0 !!
507         ENDIF
508!!
509!! skip XX_PROD, XX_CONS, O2_ANOX, RR_XXXX diagnostics here
510!!
511         IF( med_diag%IBEN_N%dgsave ) THEN
512            ALLOCATE( iben_n2d(1:jpi, 1:jpj) )
513            iben_n2d(:,:)      = 0.0 !!
514         ENDIF
515         IF( med_diag%IBEN_FE%dgsave ) THEN
516            ALLOCATE( iben_fe2d(1:jpi, 1:jpj) )
517            iben_fe2d(:,:)      = 0.0 !!
518         ENDIF
519         IF( med_diag%IBEN_C%dgsave ) THEN
520            ALLOCATE( iben_c2d(1:jpi, 1:jpj) )
521            iben_c2d(:,:)      = 0.0 !!
522         ENDIF
523         IF( med_diag%IBEN_SI%dgsave ) THEN
524            ALLOCATE( iben_si2d(1:jpi, 1:jpj) )
525            iben_si2d(:,:)      = 0.0 !!
526         ENDIF
527         IF( med_diag%IBEN_CA%dgsave ) THEN
528            ALLOCATE( iben_ca2d(1:jpi, 1:jpj) )
529            iben_ca2d(:,:)      = 0.0 !!
530         ENDIF
531         IF( med_diag%OBEN_N%dgsave ) THEN
532            ALLOCATE( oben_n2d(1:jpi, 1:jpj) )
533            oben_n2d(:,:)      = 0.0 !!
534         ENDIF
535         IF( med_diag%OBEN_FE%dgsave ) THEN
536            ALLOCATE( oben_fe2d(1:jpi, 1:jpj) )
537            oben_fe2d(:,:)      = 0.0 !!
538         ENDIF
539         IF( med_diag%OBEN_C%dgsave ) THEN
540            ALLOCATE( oben_c2d(1:jpi, 1:jpj) )
541            oben_c2d(:,:)      = 0.0 !!
542         ENDIF
543         IF( med_diag%OBEN_SI%dgsave ) THEN
544            ALLOCATE( oben_si2d(1:jpi, 1:jpj) )
545            oben_si2d(:,:)      = 0.0 !!
546         ENDIF
547         IF( med_diag%OBEN_CA%dgsave ) THEN
548            ALLOCATE( oben_ca2d(1:jpi, 1:jpj) )
549            oben_ca2d(:,:)      = 0.0 !!
550         ENDIF
551!!
552!! skip BEN_XX diagnostics here
553!!
554
555         IF( med_diag%RIV_N%dgsave ) THEN
556            ALLOCATE( rivn2d(1:jpi, 1:jpj) )
557            rivn2d(:,:)      = 0.0 !!
558         ENDIF
559         IF( med_diag%RIV_SI%dgsave ) THEN
560            ALLOCATE( rivsi2d(1:jpi, 1:jpj) )
561            rivsi2d(:,:)      = 0.0 !!
562         ENDIF
563         IF( med_diag%RIV_C%dgsave ) THEN
564            ALLOCATE( rivc2d(1:jpi, 1:jpj) )
565            rivc2d(:,:)      = 0.0 !!
566         ENDIF
567         IF( med_diag%RIV_ALK%dgsave ) THEN
568            ALLOCATE( rivalk2d(1:jpi, 1:jpj) )
569            rivalk2d(:,:)      = 0.0 !!
570         ENDIF
571         IF( med_diag%DETC%dgsave ) THEN
572            ALLOCATE( fslowc2d(1:jpi, 1:jpj) )
573            fslowc2d(:,:)      = 0.0 !!
574         ENDIF 
575!!
576!! skip SDC_XXXX, INVTXXX diagnostics here
577!!
578         IF( med_diag%LYSO_CA%dgsave ) THEN
579            ALLOCATE( lyso_ca2d(1:jpi, 1:jpj) )
580            lyso_ca2d(:,:)      = 0.0 !!
581         ENDIF
582!!
583!! skip COM_RESP diagnostic here
584!!
585         IF( med_diag%PN_LLOSS%dgsave ) THEN
586            ALLOCATE( fdpn22d(1:jpi, 1:jpj) )
587            fdpn22d(:,:)      = 0.0 !!
588         ENDIF
589         IF( med_diag%PD_LLOSS%dgsave ) THEN
590            ALLOCATE( fdpd22d(1:jpi, 1:jpj) )
591            fdpd22d(:,:)      = 0.0 !!
592         ENDIF
593         IF( med_diag%ZI_LLOSS%dgsave ) THEN
594            ALLOCATE( fdzmi22d(1:jpi, 1:jpj) )
595            fdzmi22d(:,:)      = 0.0 !!
596         ENDIF
597         IF( med_diag%ZE_LLOSS%dgsave ) THEN
598            ALLOCATE( fdzme22d(1:jpi, 1:jpj) )
599            fdzme22d(:,:)      = 0.0 !!
600         ENDIF
601         IF( med_diag%ZI_MES_N%dgsave ) THEN   
602            ALLOCATE( zimesn2d(1:jpi, 1:jpj) )
603            zimesn2d(:,:)      = 0.0 !!
604         ENDIF
605         IF( med_diag%ZI_MES_D%dgsave ) THEN
606            ALLOCATE( zimesd2d(1:jpi, 1:jpj) )
607            zimesd2d(:,:)      = 0.0 !!
608         ENDIF
609         IF( med_diag%ZI_MES_C%dgsave ) THEN
610            ALLOCATE( zimesc2d(1:jpi, 1:jpj) )
611            zimesc2d(:,:)      = 0.0 !!
612         ENDIF
613         IF( med_diag%ZI_MESDC%dgsave ) THEN
614            ALLOCATE( zimesdc2d(1:jpi, 1:jpj) )
615            zimesdc2d(:,:)      = 0.0 !!
616         ENDIF
617         IF( med_diag%ZI_EXCR%dgsave ) THEN
618            ALLOCATE( ziexcr2d(1:jpi, 1:jpj) )
619            ziexcr2d(:,:)      = 0.0 !!
620         ENDIF
621         IF( med_diag%ZI_RESP%dgsave ) THEN
622            ALLOCATE( ziresp2d(1:jpi, 1:jpj) )
623            ziresp2d(:,:)      = 0.0 !!
624         ENDIF
625         IF( med_diag%ZI_GROW%dgsave ) THEN
626            ALLOCATE( zigrow2d(1:jpi, 1:jpj) )
627            zigrow2d(:,:)      = 0.0 !!
628         ENDIF
629         IF( med_diag%ZE_MES_N%dgsave ) THEN
630            ALLOCATE( zemesn2d(1:jpi, 1:jpj) )
631            zemesn2d(:,:)      = 0.0 !!
632         ENDIF
633         IF( med_diag%ZE_MES_D%dgsave ) THEN
634            ALLOCATE( zemesd2d(1:jpi, 1:jpj) )
635            zemesd2d(:,:)      = 0.0 !!
636         ENDIF
637         IF( med_diag%ZE_MES_C%dgsave ) THEN
638            ALLOCATE( zemesc2d(1:jpi, 1:jpj) )
639            zemesc2d(:,:)      = 0.0 !!
640         ENDIF
641         IF( med_diag%ZE_MESDC%dgsave ) THEN
642            ALLOCATE( zemesdc2d(1:jpi, 1:jpj) )
643            zemesdc2d(:,:)      = 0.0 !!
644         ENDIF
645         IF( med_diag%ZE_EXCR%dgsave ) THEN
646            ALLOCATE( zeexcr2d(1:jpi, 1:jpj) )
647            zeexcr2d(:,:)      = 0.0 !!
648         ENDIF                 
649         IF( med_diag%ZE_RESP%dgsave ) THEN
650            ALLOCATE( zeresp2d(1:jpi, 1:jpj) )
651            zeresp2d(:,:)      = 0.0 !!
652         ENDIF
653         IF( med_diag%ZE_GROW%dgsave ) THEN
654            ALLOCATE( zegrow2d(1:jpi, 1:jpj) )
655            zegrow2d(:,:)      = 0.0 !!
656         ENDIF
657         IF( med_diag%MDETC%dgsave ) THEN
658            ALLOCATE( mdetc2d(1:jpi, 1:jpj) )
659            mdetc2d(:,:)      = 0.0 !!
660         ENDIF
661         IF( med_diag%GMIDC%dgsave ) THEN
662            ALLOCATE( gmidc2d(1:jpi, 1:jpj) )
663            gmidc2d(:,:)      = 0.0 !!
664         ENDIF
665         IF( med_diag%GMEDC%dgsave ) THEN
666            ALLOCATE( gmedc2d(1:jpi, 1:jpj) )
667            gmedc2d(:,:)      = 0.0 !!
668         ENDIF
669!!
670!! skip INT_XXX diagnostics here
671!!
672         IF (jdms .eq. 1) THEN
673            IF( med_diag%DMS_SURF%dgsave ) THEN
674               ALLOCATE( dms_surf2d(1:jpi, 1:jpj) )
675               dms_surf2d(:,:)      = 0.0 !!
676            ENDIF
677            IF( med_diag%DMS_ANDR%dgsave ) THEN
678               ALLOCATE( dms_andr2d(1:jpi, 1:jpj) )
679               dms_andr2d(:,:)      = 0.0 !!
680            ENDIF
681            IF( med_diag%DMS_SIMO%dgsave ) THEN
682               ALLOCATE( dms_simo2d(1:jpi, 1:jpj) )
683               dms_simo2d(:,:)      = 0.0 !!
684            ENDIF
685            IF( med_diag%DMS_ARAN%dgsave ) THEN
686               ALLOCATE( dms_aran2d(1:jpi, 1:jpj) )
687               dms_aran2d(:,:)      = 0.0 !!
688            ENDIF
689            IF( med_diag%DMS_HALL%dgsave ) THEN
690               ALLOCATE( dms_hall2d(1:jpi, 1:jpj) )
691               dms_hall2d(:,:)      = 0.0 !!
692            ENDIF
693            IF( med_diag%DMS_ANDM%dgsave ) THEN
694               ALLOCATE( dms_andm2d(1:jpi, 1:jpj) )
695               dms_andm2d(:,:)      = 0.0 !!
696            ENDIF
697         ENDIF   
698         !!
699         !! AXY (24/11/16): extra MOCSY diagnostics, 2D
700         IF( med_diag%ATM_XCO2%dgsave ) THEN
701            ALLOCATE( f_xco2a_2d(1:jpi, 1:jpj) )
702            f_xco2a_2d(:,:)      = 0.0 !!
703         ENDIF
704         IF( med_diag%OCN_FCO2%dgsave ) THEN
705            ALLOCATE( f_fco2w_2d(1:jpi, 1:jpj) )
706            f_fco2w_2d(:,:)      = 0.0 !!
707         ENDIF
708         IF( med_diag%ATM_FCO2%dgsave ) THEN
709            ALLOCATE( f_fco2a_2d(1:jpi, 1:jpj) )
710            f_fco2a_2d(:,:)      = 0.0 !!
711         ENDIF
712         IF( med_diag%OCN_RHOSW%dgsave ) THEN
713            ALLOCATE( f_ocnrhosw_2d(1:jpi, 1:jpj) )
714            f_ocnrhosw_2d(:,:)      = 0.0 !!
715         ENDIF
716         IF( med_diag%OCN_SCHCO2%dgsave ) THEN
717            ALLOCATE( f_ocnschco2_2d(1:jpi, 1:jpj) )
718            f_ocnschco2_2d(:,:)      = 0.0 !!
719         ENDIF
720         IF( med_diag%OCN_KWCO2%dgsave ) THEN
721            ALLOCATE( f_ocnkwco2_2d(1:jpi, 1:jpj) )
722            f_ocnkwco2_2d(:,:)      = 0.0 !!
723         ENDIF
724         IF( med_diag%OCN_K0%dgsave ) THEN
725            ALLOCATE( f_ocnk0_2d(1:jpi, 1:jpj) )
726            f_ocnk0_2d(:,:)      = 0.0 !!
727         ENDIF
728         IF( med_diag%CO2STARAIR%dgsave ) THEN
729            ALLOCATE( f_co2starair_2d(1:jpi, 1:jpj) )
730            f_co2starair_2d(:,:)      = 0.0 !!
731         ENDIF
732         IF( med_diag%OCN_DPCO2%dgsave ) THEN
733            ALLOCATE( f_ocndpco2_2d(1:jpi, 1:jpj) )
734            f_ocndpco2_2d(:,:)      = 0.0 !!
735         ENDIF
736#  if defined key_omip_dic         
737         IF( med_diag%PI_ATM_XCO2%dgsave ) THEN
738            ALLOCATE( f_pi_xco2a_2d(1:jpi, 1:jpj) )
739            f_pi_xco2a_2d(:,:)      = 0.0 !!
740         ENDIF
741         IF( med_diag%PI_OCN_FCO2%dgsave ) THEN
742            ALLOCATE( f_pi_fco2w_2d(1:jpi, 1:jpj) )
743            f_pi_fco2w_2d(:,:)      = 0.0 !!
744         ENDIF
745         IF( med_diag%PI_ATM_FCO2%dgsave ) THEN
746            ALLOCATE( f_pi_fco2a_2d(1:jpi, 1:jpj) )
747            f_pi_fco2a_2d(:,:)      = 0.0 !!
748         ENDIF
749         IF( med_diag%PI_CO2STARAIR%dgsave ) THEN
750            ALLOCATE( f_pi_co2starair_2d(1:jpi, 1:jpj) )
751            f_pi_co2starair_2d(:,:)      = 0.0 !!
752         ENDIF
753         IF( med_diag%OCN_PI_DPCO2%dgsave ) THEN
754            ALLOCATE( f_pi_ocndpco2_2d(1:jpi, 1:jpj) )
755            f_pi_ocndpco2_2d(:,:)      = 0.0 !!
756         ENDIF
757         IF( med_diag%PI_FGCO2%dgsave ) THEN
758            ALLOCATE( pi_fgco2(1:jpi, 1:jpj) )
759            pi_fgco2(:,:)      = 0.0 !!
760         ENDIF
761#  endif
762# endif 
763         IF( med_diag%TPP3%dgsave ) THEN
764            ALLOCATE( tpp3d(1:jpi, 1:jpj, 1:jpk) )
765            tpp3d(:,:,:)      = 0.0 !!
766         ENDIF
767         IF( med_diag%DETFLUX3%dgsave ) THEN
768            ALLOCATE( detflux3d(1:jpi, 1:jpj, 1:jpk) )
769            detflux3d(:,:,:)      = 0.0 !!
770         ENDIF
771         IF( med_diag%REMIN3N%dgsave ) THEN
772             ALLOCATE( remin3dn(1:jpi, 1:jpj, 1:jpk) )
773             remin3dn(:,:,:)      = 0.0 !!
774         ENDIF
775         !!
776         !! AXY (10/11/16): CMIP6 diagnostics, 2D
777         !! JPALM -- 17-11-16 -- put fgco2 alloc out of diag request
778         !!                   needed for coupling/passed through restart
779         !! IF( med_diag%FGCO2%dgsave ) THEN
780            ALLOCATE( fgco2(1:jpi, 1:jpj) )
781            fgco2(:,:)      = 0.0 !!
782         !! ENDIF
783         IF( med_diag%INTDISSIC%dgsave ) THEN
784            ALLOCATE( intdissic(1:jpi, 1:jpj) )
785            intdissic(:,:)  = 0.0 !!
786         ENDIF         
787         IF( med_diag%INTDISSIN%dgsave ) THEN
788            ALLOCATE( intdissin(1:jpi, 1:jpj) )
789            intdissin(:,:)  = 0.0 !!
790         ENDIF         
791         IF( med_diag%INTDISSISI%dgsave ) THEN
792            ALLOCATE( intdissisi(1:jpi, 1:jpj) )
793            intdissisi(:,:)  = 0.0 !!
794         ENDIF         
795         IF( med_diag%INTTALK%dgsave ) THEN
796            ALLOCATE( inttalk(1:jpi, 1:jpj) )
797            inttalk(:,:)  = 0.0 !!
798         ENDIF         
799         IF( med_diag%O2min%dgsave ) THEN
800            ALLOCATE( o2min(1:jpi, 1:jpj) )
801            o2min(:,:)  = 1.e3 !! set to high value as we're looking for min(o2)
802         ENDIF         
803         IF( med_diag%ZO2min%dgsave ) THEN
804            ALLOCATE( zo2min(1:jpi, 1:jpj) )
805            zo2min(:,:)  = 0.0 !!
806         ENDIF         
807         IF( med_diag%FBDDTALK%dgsave  ) THEN
808            ALLOCATE( fbddtalk(1:jpi, 1:jpj) )
809            fbddtalk(:,:)  = 0.0 !!
810         ENDIF
811         IF( med_diag%FBDDTDIC%dgsave  ) THEN
812            ALLOCATE( fbddtdic(1:jpi, 1:jpj) )
813            fbddtdic(:,:)  = 0.0 !!
814         ENDIF
815         IF( med_diag%FBDDTDIFE%dgsave ) THEN
816            ALLOCATE( fbddtdife(1:jpi, 1:jpj) )
817            fbddtdife(:,:) = 0.0 !!
818         ENDIF
819         IF( med_diag%FBDDTDIN%dgsave  ) THEN
820            ALLOCATE( fbddtdin(1:jpi, 1:jpj) )
821            fbddtdin(:,:)  = 0.0 !!
822         ENDIF
823         IF( med_diag%FBDDTDISI%dgsave ) THEN
824            ALLOCATE( fbddtdisi(1:jpi, 1:jpj) )
825            fbddtdisi(:,:) = 0.0 !!
826         ENDIF
827         !!
828         !! AXY (10/11/16): CMIP6 diagnostics, 3D
829         IF( med_diag%TPPD3%dgsave     ) THEN
830            ALLOCATE( tppd3(1:jpi, 1:jpj, 1:jpk) )
831            tppd3(:,:,:)     = 0.0 !!
832         ENDIF
833         IF( med_diag%BDDTALK3%dgsave  ) THEN
834            ALLOCATE( bddtalk3(1:jpi, 1:jpj, 1:jpk) )
835            bddtalk3(:,:,:)  = 0.0 !!
836         ENDIF
837         IF( med_diag%BDDTDIC3%dgsave  ) THEN
838            ALLOCATE( bddtdic3(1:jpi, 1:jpj, 1:jpk) )
839            bddtdic3(:,:,:)  = 0.0 !!
840         ENDIF
841         IF( med_diag%BDDTDIFE3%dgsave ) THEN
842            ALLOCATE( bddtdife3(1:jpi, 1:jpj, 1:jpk) )
843            bddtdife3(:,:,:) = 0.0 !!
844         ENDIF
845         IF( med_diag%BDDTDIN3%dgsave  ) THEN
846            ALLOCATE( bddtdin3(1:jpi, 1:jpj, 1:jpk) )
847            bddtdin3(:,:,:)  = 0.0 !!
848         ENDIF
849         IF( med_diag%BDDTDISI3%dgsave ) THEN
850            ALLOCATE( bddtdisi3(1:jpi, 1:jpj, 1:jpk) )
851            bddtdisi3(:,:,:) = 0.0 !!
852         ENDIF
853         IF( med_diag%FD_NIT3%dgsave   ) THEN
854            ALLOCATE( fd_nit3(1:jpi, 1:jpj, 1:jpk) )
855            fd_nit3(:,:,:)   = 0.0 !!
856         ENDIF
857         IF( med_diag%FD_SIL3%dgsave   ) THEN
858            ALLOCATE( fd_sil3(1:jpi, 1:jpj, 1:jpk) )
859            fd_sil3(:,:,:)   = 0.0 !!
860         ENDIF
861         IF( med_diag%FD_CAR3%dgsave   ) THEN
862            ALLOCATE( fd_car3(1:jpi, 1:jpj, 1:jpk) )
863            fd_car3(:,:,:)   = 0.0 !!
864         ENDIF
865         IF( med_diag%FD_CAL3%dgsave   ) THEN
866            ALLOCATE( fd_cal3(1:jpi, 1:jpj, 1:jpk) )
867            fd_cal3(:,:,:)   = 0.0 !!
868         ENDIF
869         IF( med_diag%DCALC3%dgsave    ) THEN
870            ALLOCATE( dcalc3(1:jpi, 1:jpj, 1:jpk) )
871            dcalc3(:,:,: )   = 0.0 !!
872         ENDIF
873         IF( med_diag%EXPC3%dgsave     ) THEN
874            ALLOCATE( expc3(1:jpi, 1:jpj, 1:jpk) )
875            expc3(:,:,: )    = 0.0 !!
876         ENDIF
877         IF( med_diag%EXPN3%dgsave     ) THEN
878            ALLOCATE( expn3(1:jpi, 1:jpj, 1:jpk) )
879            expn3(:,:,: )    = 0.0 !!
880         ENDIF
881         IF( med_diag%FEDISS3%dgsave   ) THEN
882            ALLOCATE( fediss3(1:jpi, 1:jpj, 1:jpk) )
883            fediss3(:,:,: )  = 0.0 !!
884         ENDIF
885         IF( med_diag%FESCAV3%dgsave   ) THEN
886            ALLOCATE( fescav3(1:jpi, 1:jpj, 1:jpk) )
887            fescav3(:,:,: )  = 0.0 !!
888         ENDIF
889         IF( med_diag%MIGRAZP3%dgsave   ) THEN
890            ALLOCATE( migrazp3(1:jpi, 1:jpj, 1:jpk) )
891            migrazp3(:,:,: )  = 0.0 !!
892         ENDIF
893         IF( med_diag%MIGRAZD3%dgsave   ) THEN
894            ALLOCATE( migrazd3(1:jpi, 1:jpj, 1:jpk) )
895            migrazd3(:,:,: )  = 0.0 !!
896         ENDIF
897         IF( med_diag%MEGRAZP3%dgsave   ) THEN
898            ALLOCATE( megrazp3(1:jpi, 1:jpj, 1:jpk) )
899            megrazp3(:,:,: )  = 0.0 !!
900         ENDIF
901         IF( med_diag%MEGRAZD3%dgsave   ) THEN
902            ALLOCATE( megrazd3(1:jpi, 1:jpj, 1:jpk) )
903            megrazd3(:,:,: )  = 0.0 !!
904         ENDIF
905         IF( med_diag%MEGRAZZ3%dgsave   ) THEN
906            ALLOCATE( megrazz3(1:jpi, 1:jpj, 1:jpk) )
907            megrazz3(:,:,: )  = 0.0 !!
908         ENDIF
909         IF( med_diag%O2SAT3%dgsave     ) THEN
910            ALLOCATE( o2sat3(1:jpi, 1:jpj, 1:jpk) )
911            o2sat3(:,:,: )    = 0.0 !!
912         ENDIF
913         IF( med_diag%PBSI3%dgsave      ) THEN
914            ALLOCATE( pbsi3(1:jpi, 1:jpj, 1:jpk) )
915            pbsi3(:,:,: )     = 0.0 !!
916         ENDIF
917         IF( med_diag%PCAL3%dgsave      ) THEN
918            ALLOCATE( pcal3(1:jpi, 1:jpj, 1:jpk) )
919            pcal3(:,:,: )     = 0.0 !!
920         ENDIF
921         IF( med_diag%REMOC3%dgsave     ) THEN
922            ALLOCATE( remoc3(1:jpi, 1:jpj, 1:jpk) )
923            remoc3(:,:,: )    = 0.0 !!
924         ENDIF
925         IF( med_diag%PNLIMJ3%dgsave    ) THEN
926            ALLOCATE( pnlimj3(1:jpi, 1:jpj, 1:jpk) )
927            pnlimj3(:,:,: )   = 0.0 !!
928         ENDIF
929         IF( med_diag%PNLIMN3%dgsave    ) THEN
930            ALLOCATE( pnlimn3(1:jpi, 1:jpj, 1:jpk) )
931            pnlimn3(:,:,: )   = 0.0 !!
932         ENDIF
933         IF( med_diag%PNLIMFE3%dgsave   ) THEN
934            ALLOCATE( pnlimfe3(1:jpi, 1:jpj, 1:jpk) )
935            pnlimfe3(:,:,: )  = 0.0 !!
936         ENDIF
937         IF( med_diag%PDLIMJ3%dgsave    ) THEN
938            ALLOCATE( pdlimj3(1:jpi, 1:jpj, 1:jpk) )
939            pdlimj3(:,:,: )   = 0.0 !!
940         ENDIF
941         IF( med_diag%PDLIMN3%dgsave    ) THEN
942            ALLOCATE( pdlimn3(1:jpi, 1:jpj, 1:jpk) )
943            pdlimn3(:,:,: )   = 0.0 !!
944         ENDIF
945         IF( med_diag%PDLIMFE3%dgsave   ) THEN
946            ALLOCATE( pdlimfe3(1:jpi, 1:jpj, 1:jpk) )
947            pdlimfe3(:,:,: )  = 0.0 !!
948         ENDIF
949         IF( med_diag%PDLIMSI3%dgsave   ) THEN
950            ALLOCATE( pdlimsi3(1:jpi, 1:jpj, 1:jpk) )
951            pdlimsi3(:,:,: )  = 0.0 !!
952         ENDIF
953      ENDIF
954      !! lk_iomput
955
956   END SUBROUTINE bio_medusa_init
957
958#else
959   !!======================================================================
960   !!  Dummy module :                                   No MEDUSA bio-model
961   !!======================================================================
962CONTAINS
963   SUBROUTINE bio_medusa_init( )                   ! Empty routine
964      IMPLICIT NONE
965      WRITE(*,*) 'bio_medusa_init: You should not have seen this print! error?'
966   END SUBROUTINE bio_medusa_init
967#endif 
968
969   !!======================================================================
970END MODULE bio_medusa_init_mod
Note: See TracBrowser for help on using the repository browser.