source: branches/UKMO/dev_r5518_GO6_package_asm_surf_bgc/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcnam_medusa.F90 @ 8457

Last change on this file since 8457 was 8457, checked in by dford, 3 years ago

Option to read in a value for atmospheric xCO2.

File size: 90.4 KB
Line 
1MODULE trcnam_medusa
2   !!======================================================================
3   !!                      ***  MODULE trcnam_medusa  ***
4   !! TOP :   initialisation of some run parameters for MEDUSA bio-model
5   !!======================================================================
6   !! History :   2.0  !  2007-12  (C. Ethe, G. Madec) Original code
7   !!              -   !  2008-08  (K. Popova) adaptation for MEDUSA
8   !!              -   !  2008-11  (A. Yool) continuing adaptation for MEDUSA
9   !!              -   !  2010-03  (A. Yool) updated for branch inclusion
10   !!              -   !  2011-04  (A. Yool) updated for ROAM project
11   !!              -   !  2013-05  (A. Yool) renamed (from trclsm) for v3.5
12   !!              -   !  2015-11  (J. Palmieri) added iom_use for diags
13   !!              -   !  2016-11  (A. Yool) updated diags for CMIP6
14   !!----------------------------------------------------------------------
15#if defined key_medusa
16   !!----------------------------------------------------------------------
17   !!   'key_medusa'   :                                       MEDUSA model
18   !!----------------------------------------------------------------------
19   !! trc_nam_medusa      : MEDUSA model initialisation
20   !!----------------------------------------------------------------------
21   USE oce_trc         ! Ocean variables
22   USE par_trc         ! TOP parameters
23   USE trc             ! TOP variables
24   USE sms_medusa      ! sms trends
25   USE iom             ! I/O manager
26   !!USE trc_nam_dia     ! JPALM 13-11-2015 -- if iom_use for diag
27
28   !! AXY (04/02/14): necessary to find NaNs on HECTOR
29   USE, INTRINSIC :: ieee_arithmetic 
30
31   IMPLICIT NONE
32   PRIVATE
33
34   PUBLIC   trc_nam_medusa       ! called by trcnam.F90 module
35   PUBLIC   trc_nam_iom_medusa   ! called by trcnam.F90 module
36
37   !!* Substitution
38#  include "domzgr_substitute.h90"
39   !!----------------------------------------------------------------------
40   !! NEMO/TOP 2.0 , LOCEAN-IPSL (2007)
41   !! $Id$
42   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
43   !!----------------------------------------------------------------------
44
45CONTAINS
46
47   SUBROUTINE trc_nam_medusa
48      !!----------------------------------------------------------------------
49      !!                     ***  trc_nam_medusa  *** 
50      !!
51      !! ** Purpose :   read MEDUSA namelist
52      !!
53      !! ** input   :   file 'namelist.trc.sms' containing the following
54      !!             namelist: natbio, natopt, and natdbi ("key_trc_diabio")
55      !!
56      !! ekp: namelist nabio contains ALL parameters of the ecosystem
57      !!      point sourses and sinks PLUS sediment exchange
58      !!      dia_bio - used by Lobster to output all point terms
59      !!                (sourses and sinks of bio)
60      !!      dia_add - additional diagnostics for biology such as
61      !!                primary production (2d depth integrated field or 3d)
62      !!----------------------------------------------------------------------
63      !!
64      INTEGER            :: ji,jj,jk
65      REAL(wp)           :: fthk, fdep, fdep1
66      REAL(wp)           :: q1, q2, q3
67      !
68      NAMELIST/natbio/ xxi,xaln,xald,jphy,xvpn,xvpd,          &
69      &    xsin0,xnsi0,xuif,jliebig, jq10,                    &
70      &    xthetam,xthetamd,xnln,xnld,xsld,xfln,xfld,         &
71      &  xgmi,xgme,xkmi,xkme,xphi,xbetan,xbetac,xkc,          &
72      &    xpmipn,xpmid,xpmepn,xpmepd,xpmezmi,xpmed,          &
73      &  xmetapn,xmetapd,xmetazmi,xmetazme,                   &
74      &  jmpn,xmpn,xkphn,jmpd,xmpd,xkphd,jmzmi,xmzmi,xkzmi,   &
75      &    jmzme,xmzme,xkzme,jmd,jsfd,xmd,xmdc,               &
76      &  xthetapn,xthetapd,xthetazmi,xthetazme,xthetad,       &
77      &    xrfn,xrsn,vsed,xhr,                                &
78      &  jiron,xfe_mass,xfe_sol,xfe_sed,xLgT,xk_FeL,xk_sc_Fe, &
79      &  jexport,jfdfate,jrratio,jocalccd,xridg_r0,           &
80      &    xfdfrac1,xfdfrac2,xfdfrac3,                        &
81      &    xcaco3a,xcaco3b,xmassc,xmassca,xmasssi,xprotca,    &
82      &    xprotsi,xfastc,xfastca,xfastsi,                    &
83      &  jorgben,jinorgben,xsedn,xsedfe,xsedsi,xsedc,xsedca,  &
84      &    xburial,                                           &
85      &  jriver_n,jriver_si,jriver_c,jriver_alk,jriver_dep,   &
86      &  xsdiss,                                              &
87      &  sedlam,sedlostpoc,jpkb,jdms,jdms_input,jdms_model
88#if defined key_roam
89      NAMELIST/natroam/ xthetaphy,xthetazoo,xthetanit,        &
90# if defined key_foam_medusa
91      &    xobs_xco2a,                                        &
92# endif
93      &    xthetarem,xo2min 
94#endif
95      NAMELIST/natopt/xkg0,xkr0,xkgp,xkrp,xlg,xlr,rpig
96      INTEGER :: jl, jn
97      INTEGER :: ios                 ! Local integer output status for namelist read
98      TYPE(DIAG), DIMENSION(jp_medusa_2d)  :: meddia2d
99      TYPE(DIAG), DIMENSION(jp_medusa_3d)  :: meddia3d
100      TYPE(DIAG), DIMENSION(jp_medusa_trd) :: meddiabio
101      CHARACTER(LEN=32)   ::   clname
102      !!
103      NAMELIST/nammeddia/ meddia3d, meddia2d     ! additional diagnostics
104
105      !!----------------------------------------------------------------------
106
107      IF(lwp) WRITE(numout,*)
108      clname = 'namelist_medusa'
109      IF(lwp) WRITE(numout,*) ' trc_nam_medusa: read MEDUSA namelist'
110      IF(lwp) WRITE(numout,*) ' ~~~~~~~~~~~~~~'
111# if defined key_debug_medusa
112      CALL flush(numout)
113# endif
114
115
116      CALL ctl_opn( numnatp_ref, TRIM( clname )//'_ref', 'OLD', 'FORMATTED', 'SEQUENTIAL', -1, numout, .FALSE. )
117      CALL ctl_opn( numnatp_cfg, TRIM( clname )//'_cfg', 'OLD', 'FORMATTED', 'SEQUENTIAL', -1, numout, .FALSE. )
118      IF(lwm) CALL ctl_opn( numonp     , 'output.namelist.pis' , 'UNKNOWN', 'FORMATTED', 'SEQUENTIAL', -1, numout, .FALSE. )
119
120# if defined key_debug_medusa
121      CALL flush(numout)
122      IF (lwp) write (numout,*) '------------------------------'
123      IF (lwp) write (numout,*) 'Jpalm - debug'
124      IF (lwp) write (numout,*) 'open namelist_medusa -- OK'
125      IF (lwp) write (numout,*) 'Now, read namilists inside :'
126      IF (lwp) write (numout,*) ' '
127# endif
128      !
129# if defined key_debug_medusa
130      CALL flush(numout)
131# endif
132      !
133# if defined key_debug_medusa
134      IF (lwp) write (numout,*) '------------------------------'
135      IF (lwp) write (numout,*) 'Jpalm - debug'
136      IF (lwp) write (numout,*) 'Just before reading namelist_medusa :: nammeddia'
137      IF (lwp) write (numout,*) ' '
138      CALL flush(numout)
139# endif
140
141     IF( ( .NOT.lk_iomput .AND. ln_diatrc ) .OR. ( ln_diatrc .AND. lk_medusa ) ) THEN
142         !
143         ! Namelist nammeddia
144         ! -------------------
145         REWIND( numnatp_ref )              ! Namelist nammeddia in reference namelist : MEDUSA diagnostics
146         READ  ( numnatp_ref, nammeddia, IOSTAT = ios, ERR = 901)
147901      IF( ios /= 0 ) CALL ctl_nam ( ios , 'nammeddia in reference namelist', lwp )
148
149         REWIND( numnatp_cfg )              ! Namelist nammeddia in configuration namelist : MEDUSA diagnostics
150         READ  ( numnatp_cfg, nammeddia, IOSTAT = ios, ERR = 902 )
151902      IF( ios /= 0 ) CALL ctl_nam ( ios , 'nammeddia in configuration namelist', lwp )
152         IF(lwm) WRITE ( numonp, nammeddia )
153
154# if defined key_debug_medusa
155         IF (lwp) write (numout,*) '------------------------------'
156         IF (lwp) write (numout,*) 'Jpalm - debug'
157         IF (lwp) write (numout,*) 'reading namelist_medusa :: nammeddia OK'
158         IF (lwp) write (numout,*) 'Check number of variable in nammeddia:'
159         IF (lwp) write (numout,*) 'jp_medusa_2d: ',jp_medusa_2d ,'jp_medusa_3d: ',jp_medusa_3d
160         IF (lwp) write (numout,*) ' '
161         CALL flush(numout)
162# endif
163         DO jl = 1, jp_medusa_2d
164            jn = jp_msa0_2d + jl - 1
165# if defined key_debug_medusa
166            IF (lwp) write (numout,*) 'Check what is readden in nammeddia: -- 2D'
167            IF (lwp) write (numout,*) jl,'meddia2d-sname: ',meddia2d(jl)%sname 
168            IF (lwp) write (numout,*) jl,'meddia2d-lname: ',meddia2d(jl)%lname 
169            IF (lwp) write (numout,*) jl,'meddia2d-units: ',meddia2d(jl)%units 
170            CALL flush(numout)
171# endif
172            ctrc2d(jn) = meddia2d(jl)%sname
173            ctrc2l(jn) = meddia2d(jl)%lname
174            ctrc2u(jn) = meddia2d(jl)%units
175         END DO
176
177         DO jl = 1, jp_medusa_3d
178            jn = jp_msa0_3d + jl - 1
179# if defined key_debug_medusa
180            IF (lwp) write (numout,*) 'Check what is readden in nammeddia: -- 3D'
181            IF (lwp) write (numout,*) jl,'meddia3d-sname: ',meddia3d(jl)%sname 
182            IF (lwp) write (numout,*) jl,'meddia3d-lname: ',meddia3d(jl)%lname
183            IF (lwp) write (numout,*) jl,'meddia3d-units: ',meddia3d(jl)%units
184            CALL flush(numout)
185# endif
186            ctrc3d(jn) = meddia3d(jl)%sname
187            ctrc3l(jn) = meddia3d(jl)%lname
188            ctrc3u(jn) = meddia3d(jl)%units
189         END DO
190
191         IF(lwp) THEN                   ! control print
192# if defined key_debug_medusa
193            IF (lwp) write (numout,*) '------------------------------'
194            IF (lwp) write (numout,*) 'Jpalm - debug'
195            IF (lwp) write (numout,*) 'Var name assignation OK'
196            IF (lwp) write (numout,*) 'next check var names'
197            IF (lwp) write (numout,*) ' '
198            CALL flush(numout)
199# endif
200            WRITE(numout,*)
201            WRITE(numout,*) ' Namelist : natadd'
202            DO jl = 1, jp_medusa_3d
203               jn = jp_msa0_3d + jl - 1
204               WRITE(numout,*) '  3d diag nb : ', jn, '    short name : ', ctrc3d(jn), &
205                 &             '  long name  : ', ctrc3l(jn), '   unit : ', ctrc3u(jn)
206            END DO
207            WRITE(numout,*) ' '
208
209            DO jl = 1, jp_medusa_2d
210               jn = jp_msa0_2d + jl - 1
211               WRITE(numout,*) '  2d diag nb : ', jn, '    short name : ', ctrc2d(jn), &
212                 &             '  long name  : ', ctrc2l(jn), '   unit : ', ctrc2u(jn)
213            END DO
214            WRITE(numout,*) ' '
215         ENDIF
216         !
217      ENDIF   
218         !
219# if defined key_debug_medusa
220            CALL flush(numout)
221# endif
222
223      ! 1.4 namelist natbio : biological parameters
224      ! -------------------------------------------
225     
226      xxi         = 0.
227      xaln        = 0.
228      xald        = 0.
229      jphy        = 0
230      xvpn        = 0.
231      xvpd        = 0.
232      xthetam     = 0.
233      xthetamd    = 0.
234!!
235      xsin0       = 0.
236      xnsi0       = 0.
237      xuif        = 0.
238!!
239      jliebig     = 0
240      jq10        = 0.
241      xnln        = 0.
242      xnld        = 0.
243      xsld        = 0.
244      xfln        = 0.
245      xfld        = 0.
246!!
247      xgmi        = 0.
248      xgme        = 0.
249      xkmi        = 0.
250      xkme        = 0.
251      xphi    = 0.
252      xbetan      = 0.
253      xbetac      = 0.
254      xkc         = 0.
255      xpmipn      = 0.
256      xpmid       = 0.
257      xpmepn      = 0.
258      xpmepd      = 0.
259      xpmezmi     = 0.
260      xpmed       = 0.
261!!
262      xmetapn     = 0.
263      xmetapd     = 0.
264      xmetazmi    = 0.
265      xmetazme    = 0.
266!!
267      jmpn        = 0
268      xmpn        = 0.
269      xkphn       = 0.
270      jmpd        = 0
271      xmpd        = 0.
272      xkphd       = 0.
273      jmzmi       = 0
274      xmzmi       = 0.
275      xkzmi       = 0.
276      jmzme       = 0
277      xmzme       = 0.
278      xkzme       = 0.
279!!
280      jmd         = 0
281      jsfd        = 0
282      xmd         = 0.
283      xmdc        = 0.
284!!
285      xthetapn    = 0.
286      xthetapd    = 0.
287      xthetazmi   = 0.
288      xthetazme   = 0.
289      xthetad     = 0.
290      xrfn        = 0.
291      xrsn        = 0.  !: (NOT USED HERE; RETAINED FOR LOBSTER)
292!!
293      jiron       = 0
294      xfe_mass    = 0.
295      xfe_sol     = 0.
296      xfe_sed     = 0.
297      xLgT        = 0.
298      xk_FeL     = 0.
299      xk_sc_Fe    = 0.
300!!
301      jexport     = 0
302      jfdfate     = 0
303      jrratio     = 0
304      jocalccd    = 0
305      xridg_r0    = 0.
306      xfdfrac1   = 0.
307      xfdfrac2   = 0.
308      xfdfrac3   = 0.
309      xcaco3a    = 0.
310      xcaco3b    = 0.
311      xmassc     = 0.
312      xmassca    = 0.
313      xmasssi    = 0.
314      xprotca    = 0.
315      xprotsi    = 0.
316      xfastc     = 0.
317      xfastca    = 0.
318      xfastsi    = 0.
319!!
320      jorgben     = 0
321      jinorgben   = 0
322      xsedn       = 0.
323      xsedfe      = 0.
324      xsedsi      = 0.
325      xsedc       = 0.
326      xsedca      = 0.
327      xburial     = 0.
328!!
329      jriver_n    = 0
330      jriver_si   = 0
331      jriver_c    = 0
332      jriver_alk  = 0
333      jriver_dep  = 1
334!!
335      xsdiss     = 0.
336!!
337      vsed        = 0.
338      xhr         = 0.
339!!
340      sedlam     = 0.
341      sedlostpoc  = 0.
342      jpkb    = 0.
343      jdms        = 0
344      jdms_input  = 0
345      jdms_model  = 0
346           
347      !REWIND(numnatm)
348      !READ(numnatm,natbio)
349         ! Namelist natbio
350         ! -------------------
351         REWIND( numnatp_ref )              ! Namelist natbio in reference namelist : MEDUSA diagnostics
352         READ  ( numnatp_ref, natbio, IOSTAT = ios, ERR = 903)
353903      IF( ios /= 0 ) CALL ctl_nam ( ios , 'natbio in reference namelist', lwp )
354
355         REWIND( numnatp_cfg )              ! Namelist natbio in configuration namelist : MEDUSA diagnostics
356         READ  ( numnatp_cfg, natbio, IOSTAT = ios, ERR = 904 )
357904      IF( ios /= 0 ) CALL ctl_nam ( ios , 'natbio in configuration namelist', lwp )
358         IF(lwm) WRITE ( numonp, natbio )
359
360!! Primary production and chl related quantities
361!!       xxi         :  conversion factor from gC to mmolN
362!!       xaln        :  Chl-a specific initial slope of P-I curve for non-diatoms
363!!       xald        :  Chl-a specific initial slope of P-I curve for diatoms
364!!       jphy        :  phytoplankton T-dependent growth switch
365!!       xvpn        :  maximum growth rate for non-diatoms
366!!       xvpd        :  maximum growth rate for diatoms
367!!       xthetam     :  maximum Chl to C ratio for non-diatoms     
368!!       xthetamd    :  maximum Chl to C ratio for diatoms     
369!!
370!! Diatom silicon parameters
371!!       xsin0       :  minimum diatom Si:N ratio
372!!       xnsi0       :  minimum diatom N:Si ratio
373!!       xuif        :  hypothetical growth ratio at infinite Si:N ratio
374!!
375!! Nutrient limitation
376!!       jliebig     :  Liebig nutrient uptake switch
377!!       xnln        :  half-sat constant for DIN uptake by non-diatoms
378!!       xnld        :  half-sat constant for DIN uptake by diatoms
379!!       xsl         :  half-sat constant for Si uptake by diatoms
380!!       xfld        :  half-sat constant for Fe uptake by diatoms 
381!!       xfln        :  half-sat constant for Fe uptake by non-datoms
382!!
383!! Grazing
384!!       xgmi        :  microzoo maximum growth rate
385!!       xgme        :  mesozoo maximum growth rate
386!!       xkmi        :  microzoo grazing half-sat parameter
387!!       xkme        :  mesozoo grazing half-sat parameter
388!!       xphi        :  micro/mesozoo grazing inefficiency
389!!       xbetan      :  micro/mesozoo N assimilation efficiency
390!!       xbetac      :  micro/mesozoo C assimilation efficiency
391!!       xkc         :  micro/mesozoo net C growth efficiency
392!!       xpmipn      :  grazing preference of microzoo for non-diatoms
393!!       xpmid       :  grazing preference of microzoo for diatoms
394!!       xpmepn      :  grazing preference of mesozoo for non-diatoms
395!!       xpmepd      :  grazing preference of mesozoo for diatoms
396!!       xpmezmi     :  grazing preference of mesozoo for microzoo
397!!       xpmed       :  grazing preference of mesozoo for detritus
398!!
399!! Metabolic losses
400!!       xmetapn     :  non-diatom metabolic loss rate
401!!       xmetapd     :  diatom     metabolic loss rate
402!!       xmetazmi    :  microzoo   metabolic loss rate
403!!       xmetazme    :  mesozoo    metabolic loss rate
404!!
405!! Mortality/Remineralisation
406!!       jmpn        :  non-diatom mortality functional form
407!!       xmpn        :  non-diatom mortality rate
408!!       xkphn       :  non-diatom mortality half-sat constant
409!!       jmpd        :  diatom     mortality functional form
410!!       xmpd        :  diatom mortality rate
411!!       xkphd       :  diatom mortality half-sat constant
412!!       jmzmi       :  microzoo   mortality functional form
413!!       xmzmi       :  microzoo mortality rate
414!!       xkzmi       :  microzoo mortality half-sat constant
415!!       jmzme       :  mesozoo    mortality functional form
416!!       xmzme       :  mesozoo mortality rate
417!!       xkzme       :  mesozoo mortality half-sat constant
418!!
419!! Remineralisation
420!!       jmd         :  detritus T-dependent remineralisation switch
421!!       jsfd        :  accelerate seafloor detritus remin. switch
422!!       xmd         :  detrital nitrogen remineralisation rate
423!!       xmdc        :  detrital carbon remineralisation rate
424!!
425!! Stochiometric ratios
426!!       xthetapn    :  non-diatom C:N ratio
427!!       xthetapd    :  diatom C:N ratio
428!!       xthetazmi   :  microzoo C:N ratio
429!!       xthetazme   :  mesozoo C:N ratio
430!!       xthetad     :  detritus C:N ratio
431!!       xrfn        :  phytoplankton Fe:N ratio
432!!  xrsn        :  diatom Si:N ratio (*NOT* used)
433!!
434!! Iron parameters
435!!       jiron       :  iron scavenging submodel switch
436!!       xfe_mass    :  iron atomic mass
437!!  xfe_sol     :  aeolian iron solubility
438!!  xfe_sed     :  sediment iron input
439!!  xLgT      :  total ligand concentration (umol/m3)
440!!  xk_FeL       :  dissociation constant for (Fe + L)
441!!  xk_sc_Fe    :  scavenging rate of "free" iron
442!! 
443!! Fast-sinking detritus parameters
444!!       jexport     :  fast detritus remineralisation switch
445!!       jfdfate     :  fate of fast detritus at seafloor switch
446!!       jrratio     :  rain ratio switch
447!!       jocalccd    :  CCD switch
448!!       xridg_r0    :  Ridgwell rain ratio coefficient
449!!       xfdfrac1    :  fast-sinking fraction of diatom nat. mort. losses
450!!       xfdfrac2    :  fast-sinking fraction of meszooplankton mort. losses
451!!       xfdfrac3    :  fast-sinking fraction of diatom silicon grazing losses
452!!       xcaco3a     :  polar (high latitude) CaCO3 fraction
453!!       xcaco3b     :  equatorial (low latitude) CaCO3 fraction
454!!       xmassc      :  organic C mass:mole ratio, C106 H175 O40 N16 P1
455!!       xmassca     :  calcium carbonate mass:mole ratio, CaCO3
456!!       xmasssi     :  biogenic silicon mass:mole ratio, (H2SiO3)n
457!!       xprotca     :  calcium carbonate protection ratio
458!!       xprotsi     :  biogenic silicon protection ratio
459!!       xfastc      :  organic C remineralisation length scale
460!!       xfastca     :  calcium carbonate dissolution length scale
461!!       xfastsi     :  biogenic silicon dissolution length scale
462!!
463!! Benthic
464!!       jorgben     :  does   organic detritus go to the benthos?
465!!       jinorgben   :  does inorganic detritus go to the benthos?
466!!       xsedn       :  organic   nitrogen sediment remineralisation rate
467!!       xsedfe      :  organic   iron     sediment remineralisation rate
468!!       xsedsi      :  inorganic silicon  sediment dissolution      rate
469!!       xsedc       :  organic   carbon   sediment remineralisation rate
470!!       xsedca      :  inorganic carbon   sediment dissolution      rate
471!!       xburial     :  burial rate of seafloor detritus
472!!
473!! Riverine inputs
474!!       jriver_n    :  riverine N          input?
475!!       jriver_si   :  riverine Si         input?
476!!       jriver_c    :  riverine C          input?
477!!       jriver_alk  :  riverine alkalinity input?
478!!       jriver_dep  :  depth of riverine   input?
479!!
480!! Miscellaneous
481!!       xsdiss      :  diatom frustule dissolution rate
482!!
483!! Gravitational sinking     
484!!       vsed        :  detritus gravitational sinking rate
485!!       xhr         :  coeff for Martin's remineralisation profile
486!!
487!! Additional parameters
488!!       sedlam      :  time coeff of POC in sediments
489!!      sedlostpoc   :  sediment geol loss for POC
490!!       jpkb        :  vertical layer for diagnostic of the vertical flux
491!!                      NOTE that in LOBSTER it is a first vertical layers where
492!!                      biology is active 
493!!
494!! UKESM1 - new diagnostics  !! Jpalm
495!!       jdms        :  include dms diagnostics
496!!  jdms_input  :  use instant (0) or diel-avg (1) inputs
497!!       jdms_model  :  choice of DMS model passed to atmosphere
498!!                      1 = ANDR, 2 = SIMO, 3 = ARAN, 4 = HALL, 5 = ANDM
499!!
500      IF(lwp) THEN
501!!
502!! AXY (08/11/13): compilation key notification
503         WRITE(numout,*) '=== Compilation keys'
504#if defined key_roam
505         WRITE(numout,*)     &
506         &   ' key_roam                                                               = ACTIVE'
507#else
508         WRITE(numout,*)     &
509         &   ' key_roam                                                               = INACTIVE'
510#endif       
511#if defined key_axy_carbchem
512         WRITE(numout,*)     &
513         &   ' key_axy_carbchem                                                       = ACTIVE'
514#else
515         WRITE(numout,*)     &
516         &   ' key_axy_carbchem                                                       = INACTIVE'
517#endif       
518#if defined key_mocsy
519         WRITE(numout,*)     &
520         &   ' key_mocsy                                                              = ACTIVE'
521#else
522         WRITE(numout,*)     &
523         &   ' key_mocsy                                                              = INACTIVE'
524#endif       
525#if defined key_avgqsr_medusa
526         WRITE(numout,*)     &
527         &   ' key_avgqsr_medusa                                                      = ACTIVE'
528#else
529         WRITE(numout,*)     &
530         &   ' key_avgqsr_medusa                                                      = INACTIVE'
531#endif       
532#if defined key_bs_axy_zforce
533         WRITE(numout,*)     &
534         &   ' key_bs_axy_zforce                                                      = ACTIVE'
535#else
536         WRITE(numout,*)     &
537         &   ' key_bs_axy_zforce                                                      = INACTIVE'
538#endif       
539#if defined key_bs_axy_yrlen
540         WRITE(numout,*)     &
541         &   ' key_bs_axy_yrlen                                                       = ACTIVE'
542#else
543         WRITE(numout,*)     &
544         &   ' key_bs_axy_yrlen                                                       = INACTIVE'
545#endif       
546#if defined key_deep_fe_fix
547         WRITE(numout,*)     &
548         &   ' key_deep_fe_fix                                                        = ACTIVE'
549#else
550         WRITE(numout,*)     &
551         &   ' key_deep_fe_fix                                                        = INACTIVE'
552#endif       
553#if defined key_axy_nancheck
554         WRITE(numout,*)     &
555         &   ' key_axy_nancheck                                                       = ACTIVE'
556#else
557         WRITE(numout,*)     &
558         &   ' key_axy_nancheck                                                       = INACTIVE'
559#endif       
560# if defined key_axy_pi_co2
561         WRITE(numout,*)     &
562         &   ' key_axy_pi_co2                                                         = ACTIVE'
563#else
564         WRITE(numout,*)     &
565         &   ' key_axy_pi_co2                                                         = INACTIVE'
566# endif
567# if defined key_debug_medusa
568         WRITE(numout,*)     &
569         &   ' key_debug_medusa                                                       = ACTIVE'
570#else
571         WRITE(numout,*)     &
572         &   ' key_debug_medusa                                                       = INACTIVE'
573# endif
574#if defined key_foam_medusa
575         WRITE(numout,*)     &
576         &   ' key_foam_medusa                                                        = ACTIVE'
577#else
578         WRITE(numout,*)     &
579         &   ' key_foam_medusa                                                        = INACTIVE'
580#endif
581         WRITE(numout,*) ' '
582
583         WRITE(numout,*) 'natbio'
584         WRITE(numout,*) ' '
585!!
586!! Primary production and chl related quantities
587         WRITE(numout,*) '=== Primary production'
588         WRITE(numout,*)     &
589         &   ' conversion factor from gC to mmolN,                        xxi         =', xxi
590         WRITE(numout,*)     &
591         &   ' Chl-a specific initial slope of P-I curve for non-diatoms, xaln        = ', xaln
592         WRITE(numout,*)     &
593         &   ' Chl-a specific initial slope of P-I curve for diatoms,     xald        = ', xald
594         if (jphy.eq.1) then
595            WRITE(numout,*) &
596            &   ' phytoplankton growth is *temperature-dependent*            jphy        = ', jphy
597         elseif (jphy.eq.2) then
598            WRITE(numout,*) &
599            &   ' phytoplankton growth is *temperature-dependent(Q10)*       jphy        = ', jphy
600         elseif (jphy.eq.0) then
601            WRITE(numout,*) &
602            &   ' phytoplankton growth is *temperature-independent*          jphy        = ', jphy
603         endif
604         WRITE(numout,*)     &
605         &   ' maximum growth rate for non-diatoms,                       xvpn        = ', xvpn
606         WRITE(numout,*)     &
607         &   ' maximum growth rate for diatoms,                           xvpn        = ', xvpd
608         WRITE(numout,*)     &
609         &   ' maximum Chl to C ratio for non-diatoms,                    xthetam     = ', xthetam
610         WRITE(numout,*)     &
611         &   ' maximum Chl to C ratio for diatoms,                        xthetamd    = ', xthetamd
612         WRITE(numout,*)     &
613         &   ' specific Q10 value (jphy==2),                                  jq10    = ', jq10
614!!
615!! Diatom silicon parameters
616         WRITE(numout,*) '=== Diatom silicon parameters'
617         WRITE(numout,*)     &
618         &   ' minimum diatom Si:N ratio,                                 xsin0       = ', xsin0
619         WRITE(numout,*)     &
620         &   ' minimum diatom N:Si ratio,                                 xnsi0       = ', xnsi0
621         WRITE(numout,*)     &
622         &   ' hypothetical growth ratio at infinite Si:N ratio,          xuif        = ', xuif
623!!
624!! Nutrient limitation
625         WRITE(numout,*) '=== Nutrient limitation'
626         if (jliebig.eq.1) then
627            WRITE(numout,*) &
628            &   ' nutrient uptake is a Liebig Law (= most limiting) function jliebig     = ', jliebig
629         elseif (jliebig.eq.0) then
630            WRITE(numout,*) &
631            &   ' nutrient uptake is a multiplicative function               jliebig     = ', jliebig
632         endif
633         WRITE(numout,*)     &
634         &   ' half-sat constant for DIN uptake by non-diatoms,           xnln        = ', xnln
635         WRITE(numout,*)     &
636         &   ' half-sat constant for DIN uptake by diatoms,               xnld        = ', xnld
637         WRITE(numout,*)     &
638         &   ' half-sat constant for Si uptake by diatoms,                xsld        = ', xsld
639         WRITE(numout,*)     &
640         &   ' half-sat constant for Fe uptake by non-diatoms,            xfln        = ', xfln
641         WRITE(numout,*)     &
642         &   ' half-sat constant for Fe uptake by diatoms,                xfld        = ', xfld
643!!
644!! Grazing
645         WRITE(numout,*) '=== Zooplankton grazing'
646         WRITE(numout,*)     &
647         &   ' microzoo maximum growth rate,                              xgmi        = ', xgmi
648         WRITE(numout,*)     &
649         &   ' mesozoo maximum growth rate,                               xgme        = ', xgme
650         WRITE(numout,*)     &
651         &   ' microzoo grazing half-sat parameter,                       xkmi        = ', xkmi
652         WRITE(numout,*)     &
653         &   ' mesozoo grazing half-sat parameter,                        xkme        = ', xkme
654         WRITE(numout,*)     &
655         &   ' micro/mesozoo grazing inefficiency,                        xphi        = ', xphi
656         WRITE(numout,*)     &
657         &   ' micro/mesozoo N assimilation efficiency,                   xbetan      = ', xbetan
658         WRITE(numout,*)     &
659         &   ' micro/mesozoo C assimilation efficiency,                   xbetac      = ', xbetan
660         WRITE(numout,*)     &
661         &   ' micro/mesozoo net C growth efficiency,                     xkc         = ', xkc
662         WRITE(numout,*)     &
663         &   ' grazing preference of microzoo for non-diatoms,            xpmipn      = ', xpmipn
664         WRITE(numout,*)     &
665         &   ' grazing preference of microzoo for detritus,               xpmid       = ', xpmid
666         WRITE(numout,*)     &
667         &   ' grazing preference of mesozoo for non-diatoms,             xpmepn      = ', xpmepn
668         WRITE(numout,*)     &
669         &   ' grazing preference of mesozoo for diatoms,                 xpmepd      = ', xpmepd
670         WRITE(numout,*)     &
671         &   ' grazing preference of mesozoo for microzoo,                xpmezmi     = ', xpmezmi
672         WRITE(numout,*)     &
673         &   ' grazing preference of mesozoo for detritus,                xpmed       = ', xpmed
674!!
675!! Metabolic losses
676         WRITE(numout,*) '=== Metabolic losses'
677         WRITE(numout,*)     &
678         &   ' non-diatom metabolic loss rate,                            xmetapn     = ', xmetapn
679         WRITE(numout,*)     &
680         &   ' diatom     metabolic loss rate,                            xmetapd     = ', xmetapd
681         WRITE(numout,*)     &
682         &   ' microzoo   metabolic loss rate,                            xmetazmi    = ', xmetazmi
683         WRITE(numout,*)     &
684         &   ' mesozoo    metabolic loss rate,                            xmetazme    = ', xmetazme
685!!
686!! Mortality losses
687         WRITE(numout,*) '=== Mortality losses'
688         if (jmpn.eq.1) then
689            WRITE(numout,*)     &
690            &   ' non-diatom mortality functional form,            LINEAR    jmpn        = ', jmpn
691         elseif (jmpn.eq.2) then
692            WRITE(numout,*)     &
693            &   ' non-diatom mortality functional form,         QUADRATIC    jmpn        = ', jmpn
694         elseif (jmpn.eq.3) then
695            WRITE(numout,*)     &
696            &   ' non-diatom mortality functional form,        HYPERBOLIC    jmpn        = ', jmpn
697         elseif (jmpn.eq.4) then
698            WRITE(numout,*)     &
699            &   ' non-diatom mortality functional form,           SIGMOID    jmpn        = ', jmpn
700         endif
701         WRITE(numout,*)     &
702         &   ' non-diatom mortality rate,                                 xmpn        = ', xmpn
703         WRITE(numout,*)     &
704         &   ' non-diatom mortality half-sat constant                     xkphn       = ', xkphn
705         if (jmpd.eq.1) then
706            WRITE(numout,*)     &
707            &   ' diatom mortality functional form,                LINEAR    jmpd        = ', jmpd
708         elseif (jmpd.eq.2) then
709            WRITE(numout,*)     &
710            &   ' diatom mortality functional form,             QUADRATIC    jmpd        = ', jmpd
711         elseif (jmpd.eq.3) then
712            WRITE(numout,*)     &
713            &   ' diatom mortality functional form,            HYPERBOLIC    jmpd        = ', jmpd
714         elseif (jmpd.eq.4) then
715            WRITE(numout,*)     &
716            &   ' diatom mortality functional form,               SIGMOID    jmpd        = ', jmpd
717         endif
718         WRITE(numout,*)     &
719         &   ' diatom mortality rate,                                     xmpd        = ', xmpd
720         WRITE(numout,*)     &
721         &   ' diatom mortality half-sat constant                         xkphd       = ', xkphd
722         if (jmzmi.eq.1) then
723            WRITE(numout,*)     &
724            &   ' microzoo mortality functional form,              LINEAR    jmzmi       = ', jmzmi
725         elseif (jmzmi.eq.2) then
726            WRITE(numout,*)     &
727            &   ' microzoo mortality functional form,           QUADRATIC    jmzmi       = ', jmzmi
728         elseif (jmzmi.eq.3) then
729            WRITE(numout,*)     &
730            &   ' microzoo mortality functional form,          HYPERBOLIC    jmzmi       = ', jmzmi
731         elseif (jmzmi.eq.4) then
732            WRITE(numout,*)     &
733            &   ' microzoo mortality functional form,             SIGMOID    jmzmi       = ', jmzmi
734         endif
735         WRITE(numout,*)     &
736         &   ' microzoo mortality rate,                                   xmzmi       = ', xmzmi
737         WRITE(numout,*)     &
738         &   ' mesozoo mortality half-sat constant,                       xkzmi       = ', xkzmi
739         if (jmzme.eq.1) then
740            WRITE(numout,*)     &
741            &   ' mesozoo mortality functional form,               LINEAR    jmzme       = ', jmzme
742         elseif (jmzme.eq.2) then
743            WRITE(numout,*)     &
744            &   ' mesozoo mortality functional form,            QUADRATIC    jmzme       = ', jmzme
745         elseif (jmzme.eq.3) then
746            WRITE(numout,*)     &
747            &   ' mesozoo mortality functional form,           HYPERBOLIC    jmzme       = ', jmzme
748         elseif (jmzme.eq.4) then
749            WRITE(numout,*)     &
750            &   ' mesozoo mortality functional form,              SIGMOID    jmzme       = ', jmzme
751         endif
752         WRITE(numout,*)     &
753         &   ' mesozoo mortality rate,                                    xmzme       = ', xmzme
754         WRITE(numout,*)     &
755         &   ' mesozoo mortality half-sat constant,                       xkzme       = ', xkzme
756!!
757!! Remineralisation
758         WRITE(numout,*) '=== Remineralisation'
759         if (jmd.eq.1) then
760            WRITE(numout,*) &
761            &   ' detritus remineralisation is *temperature-dependent*       jmd         = ', jmd
762         elseif (jmd.eq.2) then
763            WRITE(numout,*) &
764            &   ' detritus remineralisation is *temperature-dependent(Q10)*  jmd         = ', jmd
765         elseif (jmd.eq.0) then
766            WRITE(numout,*) &
767            &   ' detritus remineralisation is *temperature-independent*     jmd         = ', jmd
768         endif
769         if (jsfd.eq.1) then
770            WRITE(numout,*) &
771            &   ' detritus seafloor remineralisation is *accelerated*        jsfd        = ', jsfd
772         else
773            WRITE(numout,*) &
774            &   ' detritus seafloor remineralisation occurs at same rate     jsfd        = ', jsfd
775         endif
776         WRITE(numout,*)     &
777         &   ' detrital nitrogen remineralisation rate,                   xmd         = ', xmd
778         WRITE(numout,*)     &
779         &   ' detrital carbon remineralisation rate,                     xmdc        = ', xmdc
780!!
781!! Stochiometric ratios
782         WRITE(numout,*) '=== Stoichiometric ratios'
783         WRITE(numout,*)     &
784         &   ' non-diatom C:N ratio,                                      xthetapn    = ', xthetapn
785         WRITE(numout,*)     &
786         &   ' diatom C:N ratio,                                          xthetapd    = ', xthetapd
787         WRITE(numout,*)     &
788         &   ' microzoo C:N ratio,                                        xthetazmi   = ', xthetazmi
789         WRITE(numout,*)     &
790         &   ' mesozoo C:N ratio,                                         xthetazme   = ', xthetazme
791         WRITE(numout,*)     &
792         &   ' detritus C:N ratio,                                        xthetad     = ', xthetad
793         WRITE(numout,*)     &
794         &   ' phytoplankton Fe:N ratio,                                  xrfn        = ', xrfn
795         WRITE(numout,*)     &
796         &   ' diatom Si:N ratio,                                         xrsn        = ', xrsn
797!!   
798!! Iron parameters
799         WRITE(numout,*) '=== Iron parameters'
800         if (jiron.eq.1) then
801            WRITE(numout,*)     &
802            &   ' Dutkiewicz et al. (2005) iron scavenging                   jiron       = ', jiron
803         elseif (jiron.eq.2) then
804            WRITE(numout,*)     &
805            &   ' Moore et al. (2004) iron scavenging                        jiron       = ', jiron
806         elseif (jiron.eq.3) then
807            WRITE(numout,*)     &
808            &   ' Moore et al. (2008) iron scavenging                        jiron       = ', jiron
809         elseif (jiron.eq.4) then
810            WRITE(numout,*)     &
811            &   ' Galbraith et al. (2010) iron scavenging                    jiron       = ', jiron
812         else
813            WRITE(numout,*)     &
814            &   ' There is **no** iron scavenging                            jiron       = ', jiron
815         endif
816         WRITE(numout,*)     &
817         &   ' iron atomic mass,                                          xfe_mass    = ', xfe_mass
818         WRITE(numout,*)     &
819         &   ' aeolian iron solubility,                                   xfe_sol     = ', xfe_sol
820         WRITE(numout,*)     &
821         &   ' sediment iron input,                                       xfe_sed     = ', xfe_sed
822         WRITE(numout,*)     &
823         &   ' total ligand concentration (umol/m3),                      xLgT        = ', xLgT
824         WRITE(numout,*)     &
825         &   ' dissociation constant for (Fe + L),                        xk_FeL      = ', xk_FeL
826         WRITE(numout,*)     &
827         &   ' scavenging rate for free iron,                             xk_sc_Fe    = ', xk_sc_Fe
828!!
829!! Fast-sinking detritus parameters
830         WRITE(numout,*) '=== Fast-sinking detritus'
831         if (jexport.eq.1) then
832            WRITE(numout,*) &
833            &   ' fast-detritus remin. uses Dunne et al. (2007; ballast)     jexport     = ', jexport
834         elseif (jexport.eq.2) then
835            WRITE(numout,*) &
836            &   ' fast-detritus remin. uses Martin et al. (1987)             jexport     = ', jexport
837         elseif (jexport.eq.2) then
838            WRITE(numout,*) &
839            &   ' fast-detritus remin. uses Henson et al. (2011)             jexport     = ', jexport
840         endif
841         if (jfdfate.eq.1) then
842            WRITE(numout,*) &
843            &   ' fast-detritus reaching seafloor becomes slow-detritus      jfdfate     = ', jfdfate
844         elseif (jfdfate.eq.0) then
845            WRITE(numout,*) &
846            &   ' fast-detritus reaching seafloor instantly remineralised    jfdfate     = ', jfdfate
847         endif
848#if defined key_roam
849         if (jrratio.eq.0) then
850            WRITE(numout,*) &
851            &   ' Dunne et al. (2005) rain ratio submodel                    jrratio     = ', jrratio
852         elseif (jrratio.eq.1) then
853            WRITE(numout,*) &
854            &   ' Ridgwell et al. (2007) rain ratio submodel (surface omega) jrratio     = ', jrratio
855         elseif (jrratio.eq.2) then
856            WRITE(numout,*) &
857            &   ' Ridgwell et al. (2007) rain ratio submodel (3D omega)      jrratio     = ', jrratio
858         endif
859#else         
860         jrratio = 0
861         WRITE(numout,*) &
862         &   ' Dunne et al. (2005) rain ratio submodel                    jrratio     = ', jrratio
863#endif         
864#if defined key_roam
865         if (jocalccd.eq.0) then
866            WRITE(numout,*) &
867            &   ' Default, fixed CCD used                                    jocalccd    = ', jocalccd
868         elseif (jocalccd.eq.1) then
869            WRITE(numout,*) &
870            &   ' Calculated, dynamic CCD used                               jocalccd    = ', jocalccd
871         endif
872#else         
873         jocalccd = 0
874         WRITE(numout,*) &
875         &   ' Default, fixed CCD used                                    jocalccd    = ', jocalccd
876#endif
877         WRITE(numout,*)     &
878         &   ' Ridgwell rain ratio coefficient,                           xridg_r0    = ', xridg_r0
879         WRITE(numout,*)     &
880         &   ' fast-sinking fraction of diatom nat. mort. losses,         xfdfrac1    = ', xfdfrac1
881         WRITE(numout,*)     &
882         &   ' fast-sinking fraction of mesozooplankton mort. losses,     xfdfrac2    = ', xfdfrac2
883         WRITE(numout,*)     &
884         &   ' fast-sinking fraction of diatom silicon grazing losses,    xfdfrac3    = ', xfdfrac3
885         WRITE(numout,*)     &
886         &   ' polar (high latitude) CaCO3 fraction,                      xcaco3a     = ', xcaco3a
887         WRITE(numout,*)     &
888         &   ' equatorial (low latitude) CaCO3 fraction,                  xcaco3b     = ', xcaco3b
889         WRITE(numout,*)     &
890         &   ' organic C mass:mole ratio, C106 H175 O40 N16 P1,           xmassc      = ', xmassc
891         WRITE(numout,*)     &
892         &   ' calcium carbonate mass:mole ratio, CaCO3,                  xmassca     = ', xmassca
893         WRITE(numout,*)     &
894         &   ' biogenic silicon mass:mole ratio, (H2SiO3)n,               xmasssi     = ', xmasssi
895         WRITE(numout,*)     &
896         &   ' calcium carbonate protection ratio,                        xprotca     = ', xprotca
897         WRITE(numout,*)     &
898         &   ' biogenic silicon protection ratio,                         xprotsi     = ', xprotsi
899         WRITE(numout,*)     &
900         &   ' organic C remineralisation length scale,                   xfastc      = ', xfastc
901         WRITE(numout,*)     &
902         &   ' calcium carbonate dissolution length scale,                xfastca     = ', xfastca
903         WRITE(numout,*)     &
904         &   ' biogenic silicon dissolution length scale,                 xfastsi     = ', xfastsi
905!!
906!! Benthos parameters
907         WRITE(numout,*) '=== Benthos parameters'
908         WRITE(numout,*)     &
909         &   ' does   organic detritus go to the benthos?,                jorgben     = ', jorgben
910         WRITE(numout,*)     &
911         &   ' does inorganic detritus go to the benthos?,                jinorgben   = ', jinorgben
912!!
913!! Some checks on parameters related to benthos parameters
914         if (jorgben.eq.1 .and. jsfd.eq.1) then
915            !! slow detritus going to benthos at an accelerated rate
916            WRITE(numout,*) '  === WARNING! ==='
917            WRITE(numout,*) '  jsfd *and* jorgben are active - please check that you wish this'
918            WRITE(numout,*) '  === WARNING! ==='
919         endif
920         if (jorgben.eq.1 .and. jfdfate.eq.1) then
921            !! fast detritus going to benthos but via slow detritus
922            WRITE(numout,*) '  === WARNING! ==='
923            WRITE(numout,*) '  jfdfate *and* jorgben are active - please check that you wish this'
924            WRITE(numout,*) '  === WARNING! ==='
925         endif
926         if (jorgben.eq.0 .and. jinorgben.eq.1) then
927            !! inorganic fast detritus going to benthos but organic fast detritus is not
928            WRITE(numout,*) '  === WARNING! ==='
929            WRITE(numout,*) '  jinorgben is active but jorgben is not - please check that you wish this'
930            WRITE(numout,*) '  === WARNING! ==='
931         endif
932         WRITE(numout,*)     &
933         &   ' organic   nitrogen sediment remineralisation rate,         xsedn       = ', xsedn
934         WRITE(numout,*)     &
935         &   ' organic   iron     sediment remineralisation rate,         xsedfe      = ', xsedfe
936         WRITE(numout,*)     &
937         &   ' inorganic silicon  sediment remineralisation rate,         xsedsi      = ', xsedsi
938         WRITE(numout,*)     &
939         &   ' organic   carbon   sediment remineralisation rate,         xsedc       = ', xsedc
940         WRITE(numout,*)     &
941         &   ' inorganic carbon   sediment remineralisation rate,         xsedca      = ', xsedca
942         WRITE(numout,*)     &
943         &   ' burial rate of seafloor detritus,                          xburial     = ', xburial
944!!
945!! Riverine inputs
946         WRITE(numout,*) '=== Riverine inputs'
947         if (jriver_n.eq.0) then
948            WRITE(numout,*)     &
949            &   ' *no* riverine N input,                                     jriver_n    = ', jriver_n
950         elseif (jriver_n.eq.1) then
951            WRITE(numout,*)     &
952            &   ' riverine N concentrations specified,                       jriver_n    = ', jriver_n
953         elseif (jriver_n.eq.2) then
954            WRITE(numout,*)     &
955            &   ' riverine N inputs specified,                               jriver_n    = ', jriver_n
956         endif
957         if (jriver_si.eq.0) then
958            WRITE(numout,*)     &
959            &   ' *no* riverine Si input,                                    jriver_si   = ', jriver_si
960         elseif (jriver_si.eq.1) then
961            WRITE(numout,*)     &
962            &   ' riverine Si concentrations specified,                      jriver_si   = ', jriver_si
963         elseif (jriver_si.eq.2) then
964            WRITE(numout,*)     &
965            &   ' riverine Si inputs specified,                              jriver_si   = ', jriver_si
966         endif
967         if (jriver_c.eq.0) then
968            WRITE(numout,*)     &
969            &   ' *no* riverine C input,                                     jriver_c    = ', jriver_c
970         elseif (jriver_c.eq.1) then
971            WRITE(numout,*)     &
972            &   ' riverine C concentrations specified,                       jriver_c    = ', jriver_c
973         elseif (jriver_c.eq.2) then
974            WRITE(numout,*)     &
975            &   ' riverine C inputs specified,                               jriver_c    = ', jriver_c
976         endif
977         if (jriver_alk.eq.0) then
978            WRITE(numout,*)     &
979            &   ' *no* riverine alkalinity input,                            jriver_alk  = ', jriver_alk
980         elseif (jriver_alk.eq.1) then
981            WRITE(numout,*)     &
982            &   ' riverine alkalinity concentrations specified,              jriver_alk  = ', jriver_alk
983         elseif (jriver_alk.eq.2) then
984            WRITE(numout,*)     &
985            &   ' riverine alkalinity inputs specified,                      jriver_alk  = ', jriver_alk
986         endif
987         !! AXY (19/07/12): prevent (gross) stupidity on part of user
988         if (jriver_dep.lt.1.or.jriver_dep.ge.jpk) then
989            jriver_dep = 1
990         endif
991         WRITE(numout,*)     &
992         &   ' riverine input applied to down to depth k = ...            jriver_dep  = ', jriver_dep
993!!
994!! Miscellaneous
995         WRITE(numout,*) '=== Miscellaneous'
996         WRITE(numout,*)     &
997         &   ' diatom frustule dissolution rate,                          xsdiss      = ', xsdiss
998!!
999!! Gravitational sinking     
1000         WRITE(numout,*) '=== Gravitational sinking'
1001         WRITE(numout,*)     &
1002         &   ' detritus gravitational sinking rate,                       vsed        = ', vsed
1003         WRITE(numout,*)     & 
1004         &   ' coefficient for Martin et al. (1987) remineralisation,     xhr         = ', xhr
1005!!
1006!! Non-Medusa parameters
1007         WRITE(numout,*) '=== Non-Medusa parameters'
1008         WRITE(numout,*)     & 
1009         &   ' time coeff of POC in sediments,                            sedlam      = ', sedlam
1010         WRITE(numout,*)     &
1011         &   ' Sediment geol loss for POC,                                sedlostpoc  = ', sedlostpoc
1012         WRITE(numout,*)     &
1013         &   ' Vert layer for diagnostic of vertical flux,                jpkp        = ', jpkb
1014!!
1015!! UKESM1 - new diagnostics  !! Jpalm; AXY (08/07/15)
1016         WRITE(numout,*) '=== UKESM1-related parameters'
1017         WRITE(numout,*)     &
1018         &   ' include DMS diagnostic?,                                   jdms        = ', jdms
1019         if (jdms_input .eq. 0) then
1020            WRITE(numout,*)     &
1021            &   ' use instant (0) or diel-avg (1) inputs,                    jdms_input  = instantaneous'
1022         else
1023            WRITE(numout,*)     &
1024            &   ' use instant (0) or diel-avg (1) inputs,                    jdms_input  = diel-average'
1025         endif
1026    if (jdms_model .eq. 1) then
1027            WRITE(numout,*)     &
1028            &   ' choice of DMS model passed to atmosphere,                  jdms_model  = Anderson et al. (2001)'
1029    elseif (jdms_model .eq. 2) then
1030            WRITE(numout,*)     &
1031            &   ' choice of DMS model passed to atmosphere,                  jdms_model  = Simo & Dachs (2002)'
1032    elseif (jdms_model .eq. 3) then
1033            WRITE(numout,*)     &
1034            &   ' choice of DMS model passed to atmosphere,                  jdms_model  = Aranami & Tsunogai (2004)'
1035    elseif (jdms_model .eq. 4) then
1036            WRITE(numout,*)     &
1037            &   ' choice of DMS model passed to atmosphere,                  jdms_model  = Halloran et al. (2010)'
1038    elseif (jdms_model .eq. 5) then
1039            WRITE(numout,*)     &
1040            &   ' choice of DMS model passed to atmosphere,                  jdms_model  = Anderson et al. (2001; modified)'
1041         endif
1042!!
1043      ENDIF
1044!!
1045!! Key depth positions (with thanks to Andrew Coward for bug-fixing this bit)
1046      DO jk = 1,jpk
1047         !! level thickness
1048         fthk  = e3t_1d(jk)
1049         !! level depth (top of level)
1050         fdep  = gdepw_1d(jk)
1051         !! level depth (bottom of level)
1052         fdep1 = fdep + fthk
1053         !!
1054         if (fdep.lt.100.0.AND.fdep1.gt.100.0) then        !  100 m
1055            i0100 = jk
1056         elseif (fdep.lt.150.0.AND.fdep1.gt.150.0) then    !  150 m (for BASIN)
1057            i0150 = jk
1058         elseif (fdep.lt.200.0.AND.fdep1.gt.200.0) then    !  200 m
1059            i0200 = jk
1060         elseif (fdep.lt.500.0.AND.fdep1.gt.500.0) then    !  500 m
1061            i0500 = jk
1062         elseif (fdep.lt.1000.0.AND.fdep1.gt.1000.0) then  ! 1000 m
1063            i1000 = jk
1064         elseif (fdep1.lt.1100.0) then                     ! 1100 m (for Moore et al. sedimentary iron)
1065            i1100 = jk
1066         endif
1067      enddo
1068      !!
1069      IF(lwp) THEN
1070          WRITE(numout,*) '=== Position of key depths'
1071          WRITE(numout,*)     & 
1072          &   ' jk position of  100 m horizon                              i0100       = ', i0100
1073          WRITE(numout,*)     &
1074          &   ' jk position of  150 m horizon                              i0150       = ', i0150
1075          WRITE(numout,*)     & 
1076          &   ' jk position of  200 m horizon                              i0200       = ', i0200
1077          WRITE(numout,*)     & 
1078          &   ' jk position of  500 m horizon                              i0500       = ', i0500
1079          WRITE(numout,*)     & 
1080          &   ' jk position of 1000 m horizon                              i1000       = ', i1000
1081          WRITE(numout,*)     & 
1082          &   ' jk position of 1100 m horizon [*]                          i1100       = ', i1100
1083          WRITE(numout,*) 'Got here ' , SIZE(friver_dep)
1084          CALL flush(numout)
1085      ENDIF
1086
1087#if defined key_roam
1088
1089      ! 1.4b namelist natroam : ROAM parameters
1090      ! ---------------------------------------
1091     
1092      xthetaphy = 0.
1093      xthetazoo = 0.
1094      xthetanit = 0.
1095      xthetarem = 0.
1096      xo2min    = 0.
1097# if defined key_foam_medusa
1098      xobs_xco2a = 0.
1099# endif
1100
1101      !READ(numnatm,natroam)
1102         ! Namelist natroam
1103         ! -------------------
1104         REWIND( numnatp_ref )              ! Namelist natroam in reference namelist : MEDUSA diagnostics
1105         READ  ( numnatp_ref, natroam, IOSTAT = ios, ERR = 905)
1106905      IF( ios /= 0 ) CALL ctl_nam ( ios , 'natroam in reference namelist', lwp )
1107
1108         REWIND( numnatp_cfg )              ! Namelist natroam in configuration namelist : MEDUSA diagnostics
1109         READ  ( numnatp_cfg, natroam, IOSTAT = ios, ERR = 906 )
1110906      IF( ios /= 0 ) CALL ctl_nam ( ios , 'natroam in configuration namelist', lwp )
1111         IF(lwm) WRITE ( numonp, natroam )
1112
1113!! ROAM carbon, alkalinity and oxygen cycle parameters
1114!!       xthetaphy :  oxygen evolution/consumption by phytoplankton
1115!!       xthetazoo :  oxygen consumption by zooplankton
1116!!       xthetanit :  oxygen consumption by nitrogen remineralisation
1117!!       xthetarem :  oxygen consumption by carbon remineralisation
1118!!       xo2min    :  oxygen minimum concentration
1119# if defined key_foam_medusa
1120!!       xobs_xco2a : observed atmospheric xCO2
1121# endif
1122
1123      IF(lwp) THEN
1124          WRITE(numout,*) 'natroam'
1125          WRITE(numout,*) ' '
1126!!
1127!! ROAM carbon, alkalinity and oxygen cycle parameters
1128          WRITE(numout,*) '=== ROAM carbon, alkalinity and oxygen cycle parameters'
1129          WRITE(numout,*)     &
1130          &   ' oxygen evolution/consumption by phytoplankton              xthetaphy   = ', xthetaphy
1131          WRITE(numout,*)     &
1132          &   ' oxygen consumption by zooplankton                          xthetazoo   = ', xthetazoo
1133          WRITE(numout,*)     &
1134          &   ' oxygen consumption by nitrogen remineralisation            xthetanit   = ', xthetanit
1135          WRITE(numout,*)     &
1136          &   ' oxygen consumption by carbon remineralisation              xthetarem   = ', xthetarem
1137          WRITE(numout,*)     &
1138          &   ' oxygen minimum concentration                               xo2min      = ', xo2min
1139# if defined key_foam_medusa
1140          WRITE(numout,*)     &
1141          &   ' observed atmospheric xCO2                                  xobs_xco2a  = ', xobs_xco2a
1142# endif
1143       ENDIF
1144
1145#endif
1146
1147      CALL flush(numout)
1148
1149      ! 1.5 namelist natopt : parameters for optic
1150      ! ------------------------------------------
1151
1152      xkg0  = 0.
1153      xkr0  = 0.
1154      xkgp  = 0.
1155      xkrp  = 0.
1156      xlg   = 0.
1157      xlr   = 0.
1158      rpig  = 0.
1159
1160      !READ(numnatm,natopt)
1161         ! Namelist natopt
1162         ! -------------------
1163         REWIND( numnatp_ref )              ! Namelist natopt in reference namelist : MEDUSA diagnostics
1164         READ  ( numnatp_ref, natopt, IOSTAT = ios, ERR = 907)
1165907      IF( ios /= 0 ) CALL ctl_nam ( ios , 'natopt in reference namelist', lwp )
1166
1167         REWIND( numnatp_cfg )              ! Namelist natopt in configuration namelist : MEDUSA diagnostics
1168         READ  ( numnatp_cfg, natopt, IOSTAT = ios, ERR = 908 )
1169908      IF( ios /= 0 ) CALL ctl_nam ( ios , 'natopt in configuration namelist', lwp )
1170         IF(lwm) WRITE ( numonp, natopt )
1171
1172      IF(lwp) THEN
1173         WRITE(numout,*) 'natopt'
1174         WRITE(numout,*) ' '
1175         WRITE(numout,*) ' green   water absorption coeff  xkg0  = ',xkg0
1176         WRITE(numout,*) ' red water absorption coeff      xkr0  = ',xkr0
1177         WRITE(numout,*) ' pigment red absorption coeff    xkrp  = ',xkrp
1178         WRITE(numout,*) ' pigment green absorption coeff  xkgp  = ',xkgp
1179         WRITE(numout,*) ' green chl exposant              xlg   = ',xlg
1180         WRITE(numout,*) ' red   chl exposant              xlr   = ',xlr
1181         WRITE(numout,*) ' chla/chla+phea ratio            rpig  = ',rpig
1182         WRITE(numout,*) ' '
1183
1184      ENDIF
1185
1186      IF(lwp) THEN
1187         WRITE(numout,*) 'NaN check'
1188         WRITE(numout,*) ' '
1189         q1 = -1.
1190         q2 = 0.
1191         q3 = log(q1)
1192         write (numout,*) 'q3 = ', q3
1193         if ( ieee_is_nan( q3 ) ) then
1194            write (numout,*) 'NaN detected'
1195         else
1196            write (numout,*) 'NaN not detected'
1197         endif
1198         WRITE(numout,*) ' '
1199       ENDIF
1200
1201   END SUBROUTINE trc_nam_medusa
1202   
1203   SUBROUTINE trc_nam_iom_medusa
1204      !!---------------------------------------------------------------------
1205      !!                     ***  ROUTINE trc_nam_iom_medusa  ***
1206      !!
1207      !! ** Purpose : read all diag requested in iodef file through iom_use
1208      !!              So it is done only once
1209      !!            ** All diagnostic MEDUSA could asked are registered in
1210      !!            the med_diag type with a boolean value
1211      !!            So if required, one diagnostic will be true.
1212      !!
1213      !!---------------------------------------------------------------------
1214      !!
1215      !!
1216      !!----------------------------------------------------------------------           
1217      !! Variable conventions
1218      !!----------------------------------------------------------------------
1219      !!
1220      IF (iom_use("INVTN")) THEN
1221          med_diag%INVTN%dgsave = .TRUE.
1222      ELSE
1223          med_diag%INVTN%dgsave = .FALSE.
1224      ENDIF
1225      IF  (iom_use("INVTSI")) THEN
1226          med_diag%INVTSI%dgsave = .TRUE.
1227      ELSE
1228          med_diag%INVTSI%dgsave = .FALSE.
1229      ENDIF
1230      IF  (iom_use("INVTFE")) THEN
1231          med_diag%INVTFE%dgsave = .TRUE.
1232      ELSE
1233          med_diag%INVTFE%dgsave = .FALSE.
1234      ENDIF
1235      IF  (iom_use("PRN")) THEN
1236          med_diag%PRN%dgsave = .TRUE.
1237      ELSE
1238          med_diag%PRN%dgsave = .FALSE.
1239      ENDIF
1240      IF  (iom_use("MPN")) THEN
1241          med_diag%MPN%dgsave = .TRUE.
1242      ELSE
1243          med_diag%MPN%dgsave = .FALSE.
1244      ENDIF
1245      IF  (iom_use("PRD")) THEN
1246          med_diag%PRD%dgsave = .TRUE.
1247      ELSE
1248          med_diag%PRD%dgsave = .FALSE.
1249      ENDIF
1250      IF  (iom_use("MPD")) THEN
1251          med_diag%MPD%dgsave = .TRUE.
1252      ELSE
1253          med_diag%MPD%dgsave = .FALSE.
1254      ENDIF
1255      IF  (iom_use("DSED")) THEN
1256          med_diag%DSED%dgsave = .TRUE.
1257      ELSE
1258          med_diag%DSED%dgsave = .FALSE.
1259      ENDIF
1260      IF  (iom_use("OPAL")) THEN
1261          med_diag%OPAL%dgsave = .TRUE.
1262      ELSE
1263          med_diag%OPAL%dgsave = .FALSE.
1264      ENDIF
1265      IF  (iom_use("OPALDISS")) THEN
1266          med_diag%OPALDISS%dgsave = .TRUE.
1267      ELSE
1268          med_diag%OPALDISS%dgsave = .FALSE.
1269      ENDIF
1270      IF  (iom_use("GMIPn")) THEN
1271          med_diag%GMIPn%dgsave = .TRUE.
1272      ELSE
1273          med_diag%GMIPn%dgsave = .FALSE.
1274      ENDIF
1275      IF  (iom_use("GMID")) THEN
1276          med_diag%GMID%dgsave = .TRUE.
1277      ELSE
1278          med_diag%GMID%dgsave = .FALSE.
1279      ENDIF
1280      IF  (iom_use("MZMI")) THEN
1281          med_diag%MZMI%dgsave = .TRUE.
1282      ELSE
1283          med_diag%MZMI%dgsave = .FALSE.
1284      ENDIF
1285      IF  (iom_use("GMEPN")) THEN
1286          med_diag%GMEPN%dgsave = .TRUE.
1287      ELSE
1288          med_diag%GMEPN%dgsave = .FALSE.
1289      ENDIF
1290      IF  (iom_use("GMEPD")) THEN
1291          med_diag%GMEPD%dgsave = .TRUE.
1292      ELSE
1293          med_diag%GMEPD%dgsave = .FALSE.
1294      ENDIF
1295      IF  (iom_use("GMEZMI")) THEN
1296          med_diag%GMEZMI%dgsave = .TRUE.
1297      ELSE
1298          med_diag%GMEZMI%dgsave = .FALSE.
1299      ENDIF
1300      IF  (iom_use("GMED")) THEN
1301          med_diag%GMED%dgsave = .TRUE.
1302      ELSE
1303          med_diag%GMED%dgsave = .FALSE.
1304      ENDIF
1305      IF  (iom_use("MZME")) THEN
1306          med_diag%MZME%dgsave = .TRUE.
1307      ELSE
1308          med_diag%MZME%dgsave = .FALSE.
1309      ENDIF
1310      IF  (iom_use("DEXP")) THEN
1311          med_diag%DEXP%dgsave = .TRUE.
1312      ELSE
1313          med_diag%DEXP%dgsave = .FALSE.
1314      ENDIF
1315      IF  (iom_use("DETN")) THEN
1316          med_diag%DETN%dgsave = .TRUE.
1317      ELSE
1318          med_diag%DETN%dgsave = .FALSE.
1319      ENDIF
1320      IF  (iom_use("MDET")) THEN
1321          med_diag%MDET%dgsave = .TRUE.
1322      ELSE
1323          med_diag%MDET%dgsave = .FALSE.
1324      ENDIF
1325      IF  (iom_use("AEOLIAN")) THEN
1326          med_diag%AEOLIAN%dgsave = .TRUE.
1327      ELSE
1328          med_diag%AEOLIAN%dgsave = .FALSE.
1329      ENDIF
1330      IF  (iom_use("BENTHIC")) THEN
1331          med_diag%BENTHIC%dgsave = .TRUE.
1332      ELSE
1333          med_diag%BENTHIC%dgsave = .FALSE.
1334      ENDIF
1335      IF  (iom_use("SCAVENGE")) THEN
1336          med_diag%SCAVENGE%dgsave = .TRUE.
1337      ELSE
1338          med_diag%SCAVENGE%dgsave = .FALSE.
1339      ENDIF
1340      IF  (iom_use("PN_JLIM")) THEN
1341          med_diag%PN_JLIM%dgsave = .TRUE.
1342      ELSE
1343          med_diag%PN_JLIM%dgsave = .FALSE.
1344      ENDIF
1345      IF  (iom_use("PN_NLIM")) THEN
1346          med_diag%PN_NLIM%dgsave = .TRUE.
1347      ELSE
1348          med_diag%PN_NLIM%dgsave = .FALSE.
1349      ENDIF
1350      IF  (iom_use("PN_FELIM")) THEN
1351          med_diag%PN_FELIM%dgsave = .TRUE.
1352      ELSE
1353          med_diag%PN_FELIM%dgsave = .FALSE.
1354      ENDIF
1355      IF  (iom_use("PD_JLIM")) THEN
1356          med_diag%PD_JLIM%dgsave = .TRUE.
1357      ELSE
1358          med_diag%PD_JLIM%dgsave = .FALSE.
1359      ENDIF
1360      IF  (iom_use("PD_NLIM")) THEN
1361          med_diag%PD_NLIM%dgsave = .TRUE.
1362      ELSE
1363          med_diag%PD_NLIM%dgsave = .FALSE.
1364      ENDIF
1365      IF  (iom_use("PD_FELIM")) THEN
1366          med_diag%PD_FELIM%dgsave = .TRUE.
1367      ELSE
1368          med_diag%PD_FELIM%dgsave = .FALSE.
1369      ENDIF
1370      IF  (iom_use("PD_SILIM")) THEN
1371          med_diag%PD_SILIM%dgsave = .TRUE.
1372      ELSE
1373          med_diag%PD_SILIM%dgsave = .FALSE.
1374      ENDIF
1375      IF  (iom_use("PDSILIM2")) THEN
1376          med_diag%PDSILIM2%dgsave = .TRUE.
1377      ELSE
1378          med_diag%PDSILIM2%dgsave = .FALSE.
1379      ENDIF
1380      IF  (iom_use("SDT__100")) THEN
1381          med_diag%SDT__100%dgsave = .TRUE.
1382      ELSE
1383          med_diag%SDT__100%dgsave = .FALSE.
1384      ENDIF
1385      IF  (iom_use("SDT__200")) THEN
1386          med_diag%SDT__200%dgsave = .TRUE.
1387      ELSE
1388          med_diag%SDT__200%dgsave = .FALSE.
1389      ENDIF
1390      IF  (iom_use("SDT__500")) THEN
1391          med_diag%SDT__500%dgsave = .TRUE.
1392      ELSE
1393          med_diag%SDT__500%dgsave = .FALSE.
1394      ENDIF
1395      IF  (iom_use("SDT_1000")) THEN
1396          med_diag%SDT_1000%dgsave = .TRUE.
1397      ELSE
1398          med_diag%SDT_1000%dgsave = .FALSE.
1399      ENDIF
1400      IF  (iom_use("TOTREG_N")) THEN
1401          med_diag%TOTREG_N%dgsave = .TRUE.
1402      ELSE
1403          med_diag%TOTREG_N%dgsave = .FALSE.
1404      ENDIF
1405      IF  (iom_use("TOTRG_SI")) THEN
1406          med_diag%TOTRG_SI%dgsave = .TRUE.
1407      ELSE
1408          med_diag%TOTRG_SI%dgsave = .FALSE.
1409      ENDIF
1410      IF  (iom_use("REG__100")) THEN
1411          med_diag%REG__100%dgsave = .TRUE.
1412      ELSE
1413          med_diag%REG__100%dgsave = .FALSE.
1414      ENDIF
1415      IF  (iom_use("REG__200")) THEN
1416          med_diag%REG__200%dgsave = .TRUE.
1417      ELSE
1418          med_diag%REG__200%dgsave = .FALSE.
1419      ENDIF
1420      IF  (iom_use("REG__500")) THEN
1421          med_diag%REG__500%dgsave = .TRUE.
1422      ELSE
1423          med_diag%REG__500%dgsave = .FALSE.
1424      ENDIF
1425      IF  (iom_use("REG_1000")) THEN
1426          med_diag%REG_1000%dgsave = .TRUE.
1427      ELSE
1428          med_diag%REG_1000%dgsave = .FALSE.
1429      ENDIF
1430      IF  (iom_use("FASTN")) THEN
1431          med_diag%FASTN%dgsave = .TRUE.
1432      ELSE
1433          med_diag%FASTN%dgsave = .FALSE.
1434      ENDIF
1435      IF  (iom_use("FASTSI")) THEN
1436          med_diag%FASTSI%dgsave = .TRUE.
1437      ELSE
1438          med_diag%FASTSI%dgsave = .FALSE.
1439      ENDIF
1440      IF  (iom_use("FASTFE")) THEN
1441          med_diag%FASTFE%dgsave = .TRUE.
1442      ELSE
1443          med_diag%FASTFE%dgsave = .FALSE.
1444      ENDIF
1445      IF  (iom_use("FASTC")) THEN
1446          med_diag%FASTC%dgsave = .TRUE.
1447      ELSE
1448          med_diag%FASTC%dgsave = .FALSE.
1449      ENDIF
1450      IF  (iom_use("FASTCA")) THEN
1451          med_diag%FASTCA%dgsave = .TRUE.
1452      ELSE
1453          med_diag%FASTCA%dgsave = .FALSE.
1454      ENDIF
1455      IF  (iom_use("FDT__100")) THEN
1456          med_diag%FDT__100%dgsave = .TRUE.
1457      ELSE
1458          med_diag%FDT__100%dgsave = .FALSE.
1459      ENDIF
1460      IF  (iom_use("FDT__200")) THEN
1461          med_diag%FDT__200%dgsave = .TRUE.
1462      ELSE
1463          med_diag%FDT__200%dgsave = .FALSE.
1464      ENDIF
1465      IF  (iom_use("FDT__500")) THEN
1466          med_diag%FDT__500%dgsave = .TRUE.
1467      ELSE
1468          med_diag%FDT__500%dgsave = .FALSE.
1469      ENDIF
1470      IF  (iom_use("FDT_1000")) THEN
1471          med_diag%FDT_1000%dgsave = .TRUE.
1472      ELSE
1473          med_diag%FDT_1000%dgsave = .FALSE.
1474      ENDIF
1475      IF  (iom_use("RG__100F")) THEN
1476          med_diag%RG__100F%dgsave = .TRUE.
1477      ELSE
1478          med_diag%RG__100F%dgsave = .FALSE.
1479      ENDIF
1480      IF  (iom_use("RG__200F")) THEN
1481          med_diag%RG__200F%dgsave = .TRUE.
1482      ELSE
1483          med_diag%RG__200F%dgsave = .FALSE.
1484      ENDIF
1485      IF  (iom_use("RG__500F")) THEN
1486          med_diag%RG__500F%dgsave = .TRUE.
1487      ELSE
1488          med_diag%RG__500F%dgsave = .FALSE.
1489      ENDIF
1490      IF  (iom_use("RG_1000F")) THEN
1491          med_diag%RG_1000F%dgsave = .TRUE.
1492      ELSE
1493          med_diag%RG_1000F%dgsave = .FALSE.
1494      ENDIF
1495      IF  (iom_use("FDS__100")) THEN
1496          med_diag%FDS__100%dgsave = .TRUE.
1497      ELSE
1498          med_diag%FDS__100%dgsave = .FALSE.
1499      ENDIF
1500      IF  (iom_use("FDS__200")) THEN
1501          med_diag%FDS__200%dgsave = .TRUE.
1502      ELSE
1503          med_diag%FDS__200%dgsave = .FALSE.
1504      ENDIF
1505      IF  (iom_use("FDS__500")) THEN
1506          med_diag%FDS__500%dgsave = .TRUE.
1507      ELSE
1508          med_diag%FDS__500%dgsave = .FALSE.
1509      ENDIF
1510      IF  (iom_use("FDS_1000")) THEN
1511          med_diag%FDS_1000%dgsave = .TRUE.
1512      ELSE
1513          med_diag%FDS_1000%dgsave = .FALSE.
1514      ENDIF
1515      IF  (iom_use("RGS_100F")) THEN
1516          med_diag%RGS_100F%dgsave = .TRUE.
1517      ELSE
1518          med_diag%RGS_100F%dgsave = .FALSE.
1519      ENDIF
1520      IF  (iom_use("RGS_200F")) THEN
1521          med_diag%RGS_200F%dgsave = .TRUE.
1522      ELSE
1523          med_diag%RGS_200F%dgsave = .FALSE.
1524      ENDIF
1525      IF  (iom_use("RGS_500F")) THEN
1526          med_diag%RGS_500F%dgsave = .TRUE.
1527      ELSE
1528          med_diag%RGS_500F%dgsave = .FALSE.
1529      ENDIF
1530      IF  (iom_use("RGS1000F")) THEN
1531          med_diag%RGS1000F%dgsave = .TRUE.
1532      ELSE
1533          med_diag%RGS1000F%dgsave = .FALSE.
1534      ENDIF
1535      IF  (iom_use("REMINN")) THEN
1536          med_diag%REMINN%dgsave = .TRUE.
1537      ELSE
1538          med_diag%REMINN%dgsave = .FALSE.
1539      ENDIF
1540      IF  (iom_use("REMINSI")) THEN
1541          med_diag%REMINSI%dgsave = .TRUE.
1542      ELSE
1543          med_diag%REMINSI%dgsave = .FALSE.
1544      ENDIF
1545      IF  (iom_use("REMINFE")) THEN
1546          med_diag%REMINFE%dgsave = .TRUE.
1547      ELSE
1548          med_diag%REMINFE%dgsave = .FALSE.
1549      ENDIF
1550      IF  (iom_use("REMINC")) THEN
1551          med_diag%REMINC%dgsave = .TRUE.
1552      ELSE
1553          med_diag%REMINC%dgsave = .FALSE.
1554      ENDIF
1555      IF  (iom_use("REMINCA")) THEN
1556          med_diag%REMINCA%dgsave = .TRUE.
1557      ELSE
1558          med_diag%REMINCA%dgsave = .FALSE.
1559      ENDIF
1560      IF  (iom_use("SEAFLRN")) THEN
1561          med_diag%SEAFLRN%dgsave = .TRUE.
1562      ELSE
1563          med_diag%SEAFLRN%dgsave = .FALSE.
1564      ENDIF
1565      IF  (iom_use("SEAFLRSI")) THEN
1566          med_diag%SEAFLRSI%dgsave = .TRUE.
1567      ELSE
1568          med_diag%SEAFLRSI%dgsave = .FALSE.
1569      ENDIF
1570      IF  (iom_use("SEAFLRFE")) THEN
1571          med_diag%SEAFLRFE%dgsave = .TRUE.
1572      ELSE
1573          med_diag%SEAFLRFE%dgsave = .FALSE.
1574      ENDIF
1575      IF  (iom_use("SEAFLRC")) THEN
1576          med_diag%SEAFLRC%dgsave = .TRUE.
1577      ELSE
1578          med_diag%SEAFLRC%dgsave = .FALSE.
1579      ENDIF
1580      IF  (iom_use("SEAFLRCA")) THEN
1581          med_diag%SEAFLRCA%dgsave = .TRUE.
1582      ELSE
1583          med_diag%SEAFLRCA%dgsave = .FALSE.
1584      ENDIF
1585      IF  (iom_use("MED_QSR")) THEN
1586          med_diag%MED_QSR%dgsave = .TRUE.
1587      ELSE
1588          med_diag%MED_QSR%dgsave = .FALSE.
1589      ENDIF
1590      IF  (iom_use("MED_XPAR")) THEN
1591          med_diag%MED_XPAR%dgsave = .TRUE.
1592      ELSE
1593          med_diag%MED_XPAR%dgsave = .FALSE.
1594      ENDIF
1595      IF  (iom_use("INTFLX_N")) THEN
1596          med_diag%INTFLX_N%dgsave = .TRUE.
1597      ELSE
1598          med_diag%INTFLX_N%dgsave = .FALSE.
1599      ENDIF
1600      IF  (iom_use("INTFLX_SI")) THEN
1601          med_diag%INTFLX_SI%dgsave = .TRUE.
1602      ELSE
1603          med_diag%INTFLX_SI%dgsave = .FALSE.
1604      ENDIF
1605      IF  (iom_use("INTFLX_FE")) THEN
1606          med_diag%INTFLX_FE%dgsave = .TRUE.
1607      ELSE
1608          med_diag%INTFLX_FE%dgsave = .FALSE.
1609      ENDIF
1610      IF  (iom_use("INT_PN")) THEN
1611          med_diag%INT_PN%dgsave = .TRUE.
1612      ELSE
1613          med_diag%INT_PN%dgsave = .FALSE.
1614      ENDIF
1615      IF  (iom_use("INT_PD")) THEN
1616          med_diag%INT_PD%dgsave = .TRUE.
1617      ELSE
1618          med_diag%INT_PD%dgsave = .FALSE.
1619      ENDIF
1620      IF  (iom_use("ML_PRN")) THEN
1621          med_diag%ML_PRN%dgsave = .TRUE.
1622      ELSE
1623          med_diag%ML_PRN%dgsave = .FALSE.
1624      ENDIF
1625      IF  (iom_use("ML_PRD")) THEN
1626          med_diag%ML_PRD%dgsave = .TRUE.
1627      ELSE
1628          med_diag%ML_PRD%dgsave = .FALSE.
1629      ENDIF
1630      IF  (iom_use("OCAL_CCD")) THEN
1631          med_diag%OCAL_CCD%dgsave = .TRUE.
1632      ELSE
1633          med_diag%OCAL_CCD%dgsave = .FALSE.
1634      ENDIF
1635      IF  (iom_use("OCAL_LVL")) THEN
1636          med_diag%OCAL_LVL%dgsave = .TRUE.
1637      ELSE
1638          med_diag%OCAL_LVL%dgsave = .FALSE.
1639      ENDIF
1640      IF  (iom_use("FE_0000")) THEN
1641          med_diag%FE_0000%dgsave = .TRUE.
1642      ELSE
1643          med_diag%FE_0000%dgsave = .FALSE.
1644      ENDIF
1645      IF  (iom_use("FE_0100")) THEN
1646          med_diag%FE_0100%dgsave = .TRUE.
1647      ELSE
1648          med_diag%FE_0100%dgsave = .FALSE.
1649      ENDIF
1650      IF  (iom_use("FE_0200")) THEN
1651          med_diag%FE_0200%dgsave = .TRUE.
1652      ELSE
1653          med_diag%FE_0200%dgsave = .FALSE.
1654      ENDIF
1655      IF  (iom_use("FE_0500")) THEN
1656          med_diag%FE_0500%dgsave = .TRUE.
1657      ELSE
1658          med_diag%FE_0500%dgsave = .FALSE.
1659      ENDIF
1660      IF  (iom_use("FE_1000")) THEN
1661          med_diag%FE_1000%dgsave = .TRUE.
1662      ELSE
1663          med_diag%FE_1000%dgsave = .FALSE.
1664      ENDIF
1665      IF  (iom_use("MED_XZE")) THEN
1666          med_diag%MED_XZE%dgsave = .TRUE.
1667      ELSE
1668          med_diag%MED_XZE%dgsave = .FALSE.
1669      ENDIF
1670      IF  (iom_use("WIND")) THEN
1671          med_diag%WIND%dgsave = .TRUE.
1672      ELSE
1673          med_diag%WIND%dgsave = .FALSE.
1674      ENDIF
1675      IF  (iom_use("ATM_PCO2")) THEN
1676          med_diag%ATM_PCO2%dgsave = .TRUE.
1677      ELSE
1678          med_diag%ATM_PCO2%dgsave = .FALSE.
1679      ENDIF
1680      IF  (iom_use("OCN_PH")) THEN
1681          med_diag%OCN_PH%dgsave = .TRUE.
1682      ELSE
1683          med_diag%OCN_PH%dgsave = .FALSE.
1684      ENDIF
1685      IF  (iom_use("OCN_PCO2")) THEN
1686          med_diag%OCN_PCO2%dgsave = .TRUE.
1687      ELSE
1688          med_diag%OCN_PCO2%dgsave = .FALSE.
1689      ENDIF
1690      IF  (iom_use("OCNH2CO3")) THEN
1691          med_diag%OCNH2CO3%dgsave = .TRUE.
1692      ELSE
1693          med_diag%OCNH2CO3%dgsave = .FALSE.
1694      ENDIF
1695      IF  (iom_use("OCN_HCO3")) THEN
1696          med_diag%OCN_HCO3%dgsave = .TRUE.
1697      ELSE
1698          med_diag%OCN_HCO3%dgsave = .FALSE.
1699      ENDIF
1700      IF  (iom_use("OCN_CO3")) THEN
1701          med_diag%OCN_CO3%dgsave = .TRUE.
1702      ELSE
1703          med_diag%OCN_CO3%dgsave = .FALSE.
1704      ENDIF
1705      IF  (iom_use("CO2FLUX")) THEN
1706          med_diag%CO2FLUX%dgsave = .TRUE.
1707      ELSE
1708          med_diag%CO2FLUX%dgsave = .FALSE.
1709      ENDIF
1710      IF  (iom_use("OM_CAL")) THEN
1711          med_diag%OM_CAL%dgsave = .TRUE.
1712      ELSE
1713          med_diag%OM_CAL%dgsave = .FALSE.
1714      ENDIF
1715      IF  (iom_use("OM_ARG")) THEN
1716          med_diag%OM_ARG%dgsave = .TRUE.
1717      ELSE
1718          med_diag%OM_ARG%dgsave = .FALSE.
1719      ENDIF
1720      IF  (iom_use("TCO2")) THEN
1721          med_diag%TCO2%dgsave = .TRUE.
1722      ELSE
1723          med_diag%TCO2%dgsave = .FALSE.
1724      ENDIF
1725      IF  (iom_use("TALK")) THEN
1726          med_diag%TALK%dgsave = .TRUE.
1727      ELSE
1728          med_diag%TALK%dgsave = .FALSE.
1729      ENDIF
1730      IF  (iom_use("KW660")) THEN
1731          med_diag%KW660%dgsave = .TRUE.
1732      ELSE
1733          med_diag%KW660%dgsave = .FALSE.
1734      ENDIF
1735      IF  (iom_use("ATM_PP0")) THEN
1736          med_diag%ATM_PP0%dgsave = .TRUE.
1737      ELSE
1738          med_diag%ATM_PP0%dgsave = .FALSE.
1739      ENDIF
1740      IF  (iom_use("O2FLUX")) THEN
1741          med_diag%O2FLUX%dgsave = .TRUE.
1742      ELSE
1743          med_diag%O2FLUX%dgsave = .FALSE.
1744      ENDIF
1745      IF  (iom_use("O2SAT")) THEN
1746          med_diag%O2SAT%dgsave = .TRUE.
1747      ELSE
1748          med_diag%O2SAT%dgsave = .FALSE.
1749      ENDIF
1750      IF  (iom_use("CAL_CCD")) THEN
1751          med_diag%CAL_CCD%dgsave = .TRUE.
1752      ELSE
1753          med_diag%CAL_CCD%dgsave = .FALSE.
1754      ENDIF
1755      IF  (iom_use("ARG_CCD")) THEN
1756          med_diag%ARG_CCD%dgsave = .TRUE.
1757      ELSE
1758          med_diag%ARG_CCD%dgsave = .FALSE.
1759      ENDIF
1760      IF  (iom_use("SFR_OCAL")) THEN
1761          med_diag%SFR_OCAL%dgsave = .TRUE.
1762      ELSE
1763          med_diag%SFR_OCAL%dgsave = .FALSE.
1764      ENDIF
1765      IF  (iom_use("SFR_OARG")) THEN
1766          med_diag%SFR_OARG%dgsave = .TRUE.
1767      ELSE
1768          med_diag%SFR_OARG%dgsave = .FALSE.
1769      ENDIF
1770      IF  (iom_use("N_PROD")) THEN
1771          med_diag%N_PROD%dgsave = .TRUE.
1772      ELSE
1773          med_diag%N_PROD%dgsave = .FALSE.
1774      ENDIF
1775      IF  (iom_use("N_CONS")) THEN
1776          med_diag%N_CONS%dgsave = .TRUE.
1777      ELSE
1778          med_diag%N_CONS%dgsave = .FALSE.
1779      ENDIF
1780      IF  (iom_use("C_PROD")) THEN
1781          med_diag%C_PROD%dgsave = .TRUE.
1782      ELSE
1783          med_diag%C_PROD%dgsave = .FALSE.
1784      ENDIF
1785      IF  (iom_use("C_CONS")) THEN
1786          med_diag%C_CONS%dgsave = .TRUE.
1787      ELSE
1788          med_diag%C_CONS%dgsave = .FALSE.
1789      ENDIF
1790      IF  (iom_use("O2_PROD")) THEN
1791          med_diag%O2_PROD%dgsave = .TRUE.
1792      ELSE
1793          med_diag%O2_PROD%dgsave = .FALSE.
1794      ENDIF
1795      IF  (iom_use("O2_CONS")) THEN
1796          med_diag%O2_CONS%dgsave = .TRUE.
1797      ELSE
1798          med_diag%O2_CONS%dgsave = .FALSE.
1799      ENDIF
1800      IF  (iom_use("O2_ANOX")) THEN
1801          med_diag%O2_ANOX%dgsave = .TRUE.
1802      ELSE
1803          med_diag%O2_ANOX%dgsave = .FALSE.
1804      ENDIF
1805      IF  (iom_use("RR_0100")) THEN
1806          med_diag%RR_0100%dgsave = .TRUE.
1807      ELSE
1808          med_diag%RR_0100%dgsave = .FALSE.
1809      ENDIF
1810      IF  (iom_use("RR_0500")) THEN
1811          med_diag%RR_0500%dgsave = .TRUE.
1812      ELSE
1813          med_diag%RR_0500%dgsave = .FALSE.
1814      ENDIF
1815      IF  (iom_use("RR_1000")) THEN
1816          med_diag%RR_1000%dgsave = .TRUE.
1817      ELSE
1818          med_diag%RR_1000%dgsave = .FALSE.
1819      ENDIF
1820      IF  (iom_use("IBEN_N")) THEN
1821          med_diag%IBEN_N%dgsave = .TRUE.
1822      ELSE
1823          med_diag%IBEN_N%dgsave = .FALSE.
1824      ENDIF
1825      IF  (iom_use("IBEN_FE")) THEN
1826          med_diag%IBEN_FE%dgsave = .TRUE.
1827      ELSE
1828          med_diag%IBEN_FE%dgsave = .FALSE.
1829      ENDIF
1830      IF  (iom_use("IBEN_C")) THEN
1831          med_diag%IBEN_C%dgsave = .TRUE.
1832      ELSE
1833          med_diag%IBEN_C%dgsave = .FALSE.
1834      ENDIF
1835      IF  (iom_use("IBEN_SI")) THEN
1836          med_diag%IBEN_SI%dgsave = .TRUE.
1837      ELSE
1838          med_diag%IBEN_SI%dgsave = .FALSE.
1839      ENDIF
1840      IF  (iom_use("IBEN_CA")) THEN
1841          med_diag%IBEN_CA%dgsave = .TRUE.
1842      ELSE
1843          med_diag%IBEN_CA%dgsave = .FALSE.
1844      ENDIF
1845      IF  (iom_use("OBEN_N")) THEN
1846          med_diag%OBEN_N%dgsave = .TRUE.
1847      ELSE
1848          med_diag%OBEN_N%dgsave = .FALSE.
1849      ENDIF
1850      IF  (iom_use("OBEN_FE")) THEN
1851          med_diag%OBEN_FE%dgsave = .TRUE.
1852      ELSE
1853          med_diag%OBEN_FE%dgsave = .FALSE.
1854      ENDIF
1855      IF  (iom_use("OBEN_C")) THEN
1856          med_diag%OBEN_C%dgsave = .TRUE.
1857      ELSE
1858          med_diag%OBEN_C%dgsave = .FALSE.
1859      ENDIF
1860      IF  (iom_use("OBEN_SI")) THEN
1861          med_diag%OBEN_SI%dgsave = .TRUE.
1862      ELSE
1863          med_diag%OBEN_SI%dgsave = .FALSE.
1864      ENDIF
1865      IF  (iom_use("OBEN_CA")) THEN
1866          med_diag%OBEN_CA%dgsave = .TRUE.
1867      ELSE
1868          med_diag%OBEN_CA%dgsave = .FALSE.
1869      ENDIF
1870      IF  (iom_use("BEN_N")) THEN
1871          med_diag%BEN_N%dgsave = .TRUE.
1872      ELSE
1873          med_diag%BEN_N%dgsave = .FALSE.
1874      ENDIF
1875      IF  (iom_use("BEN_FE")) THEN
1876          med_diag%BEN_FE%dgsave = .TRUE.
1877      ELSE
1878          med_diag%BEN_FE%dgsave = .FALSE.
1879      ENDIF
1880      IF  (iom_use("BEN_C")) THEN
1881          med_diag%BEN_C%dgsave = .TRUE.
1882      ELSE
1883          med_diag%BEN_C%dgsave = .FALSE.
1884      ENDIF
1885      IF  (iom_use("BEN_SI")) THEN
1886          med_diag%BEN_SI%dgsave = .TRUE.
1887      ELSE
1888          med_diag%BEN_SI%dgsave = .FALSE.
1889      ENDIF
1890      IF  (iom_use("BEN_CA")) THEN
1891          med_diag%BEN_CA%dgsave = .TRUE.
1892      ELSE
1893          med_diag%BEN_CA%dgsave = .FALSE.
1894      ENDIF
1895      IF  (iom_use("RUNOFF")) THEN
1896          med_diag%RUNOFF%dgsave = .TRUE.
1897      ELSE
1898          med_diag%RUNOFF%dgsave = .FALSE.
1899      ENDIF
1900      IF  (iom_use("RIV_N")) THEN
1901          med_diag%RIV_N%dgsave = .TRUE.
1902      ELSE
1903          med_diag%RIV_N%dgsave = .FALSE.
1904      ENDIF
1905      IF  (iom_use("RIV_SI")) THEN
1906          med_diag%RIV_SI%dgsave = .TRUE.
1907      ELSE
1908          med_diag%RIV_SI%dgsave = .FALSE.
1909      ENDIF
1910      IF  (iom_use("RIV_C")) THEN
1911          med_diag%RIV_C%dgsave = .TRUE.
1912      ELSE
1913          med_diag%RIV_C%dgsave = .FALSE.
1914      ENDIF
1915      IF  (iom_use("RIV_ALK")) THEN
1916          med_diag%RIV_ALK%dgsave = .TRUE.
1917      ELSE
1918          med_diag%RIV_ALK%dgsave = .FALSE.
1919      ENDIF
1920      IF  (iom_use("DETC")) THEN
1921          med_diag%DETC%dgsave = .TRUE.
1922      ELSE
1923          med_diag%DETC%dgsave = .FALSE.
1924      ENDIF
1925      IF  (iom_use("SDC__100")) THEN
1926          med_diag%SDC__100%dgsave = .TRUE.
1927      ELSE
1928          med_diag%SDC__100%dgsave = .FALSE.
1929      ENDIF
1930      IF  (iom_use("SDC__200")) THEN
1931          med_diag%SDC__200%dgsave = .TRUE.
1932      ELSE
1933          med_diag%SDC__200%dgsave = .FALSE.
1934      ENDIF
1935      IF  (iom_use("SDC__500")) THEN
1936          med_diag%SDC__500%dgsave = .TRUE.
1937      ELSE
1938          med_diag%SDC__500%dgsave = .FALSE.
1939      ENDIF
1940      IF  (iom_use("SDC_1000")) THEN
1941          med_diag%SDC_1000%dgsave = .TRUE.
1942      ELSE
1943          med_diag%SDC_1000%dgsave = .FALSE.
1944      ENDIF
1945      IF  (iom_use("INVTC")) THEN
1946          med_diag%INVTC%dgsave = .TRUE.
1947      ELSE
1948          med_diag%INVTC%dgsave = .FALSE.
1949      ENDIF
1950      IF  (iom_use("INVTALK")) THEN
1951          med_diag%INVTALK%dgsave = .TRUE.
1952      ELSE
1953          med_diag%INVTALK%dgsave = .FALSE.
1954      ENDIF
1955      IF  (iom_use("INVTO2")) THEN
1956          med_diag%INVTO2%dgsave = .TRUE.
1957      ELSE
1958          med_diag%INVTO2%dgsave = .FALSE.
1959      ENDIF
1960      IF  (iom_use("LYSO_CA")) THEN
1961          med_diag%LYSO_CA%dgsave = .TRUE.
1962      ELSE
1963          med_diag%LYSO_CA%dgsave = .FALSE.
1964      ENDIF
1965      IF  (iom_use("COM_RESP")) THEN
1966          med_diag%COM_RESP%dgsave = .TRUE.
1967      ELSE
1968          med_diag%COM_RESP%dgsave = .FALSE.
1969      ENDIF
1970      IF  (iom_use("PN_LLOSS")) THEN
1971          med_diag%PN_LLOSS%dgsave = .TRUE.
1972      ELSE
1973          med_diag%PN_LLOSS%dgsave = .FALSE.
1974      ENDIF
1975      IF  (iom_use("PD_LLOSS")) THEN
1976          med_diag%PD_LLOSS%dgsave = .TRUE.
1977      ELSE
1978          med_diag%PD_LLOSS%dgsave = .FALSE.
1979      ENDIF
1980      IF  (iom_use("ZI_LLOSS")) THEN
1981          med_diag%ZI_LLOSS%dgsave = .TRUE.
1982      ELSE
1983          med_diag%ZI_LLOSS%dgsave = .FALSE.
1984      ENDIF
1985      IF  (iom_use("ZE_LLOSS")) THEN
1986          med_diag%ZE_LLOSS%dgsave = .TRUE.
1987      ELSE
1988          med_diag%ZE_LLOSS%dgsave = .FALSE.
1989      ENDIF
1990      IF  (iom_use("ZI_MES_N")) THEN
1991          med_diag%ZI_MES_N%dgsave = .TRUE.
1992      ELSE
1993          med_diag%ZI_MES_N%dgsave = .FALSE.
1994      ENDIF
1995      IF  (iom_use("ZI_MES_D")) THEN
1996          med_diag%ZI_MES_D%dgsave = .TRUE.
1997      ELSE
1998          med_diag%ZI_MES_D%dgsave = .FALSE.
1999      ENDIF
2000      IF  (iom_use("ZI_MES_C")) THEN
2001          med_diag%ZI_MES_C%dgsave = .TRUE.
2002      ELSE
2003          med_diag%ZI_MES_C%dgsave = .FALSE.
2004      ENDIF
2005      IF  (iom_use("ZI_MESDC")) THEN
2006          med_diag%ZI_MESDC%dgsave = .TRUE.
2007      ELSE
2008          med_diag%ZI_MESDC%dgsave = .FALSE.
2009      ENDIF
2010      IF  (iom_use("ZI_EXCR")) THEN
2011          med_diag%ZI_EXCR%dgsave = .TRUE.
2012      ELSE
2013          med_diag%ZI_EXCR%dgsave = .FALSE.
2014      ENDIF
2015      IF  (iom_use("ZI_RESP")) THEN
2016          med_diag%ZI_RESP%dgsave = .TRUE.
2017      ELSE
2018          med_diag%ZI_RESP%dgsave = .FALSE.
2019      ENDIF
2020      IF  (iom_use("ZI_GROW")) THEN
2021          med_diag%ZI_GROW%dgsave = .TRUE.
2022      ELSE
2023          med_diag%ZI_GROW%dgsave = .FALSE.
2024      ENDIF
2025      IF  (iom_use("ZE_MES_N")) THEN
2026          med_diag%ZE_MES_N%dgsave = .TRUE.
2027      ELSE
2028          med_diag%ZE_MES_N%dgsave = .FALSE.
2029      ENDIF
2030      IF  (iom_use("ZE_MES_D")) THEN
2031          med_diag%ZE_MES_D%dgsave = .TRUE.
2032      ELSE
2033          med_diag%ZE_MES_D%dgsave = .FALSE.
2034      ENDIF
2035      IF  (iom_use("ZE_MES_C")) THEN
2036          med_diag%ZE_MES_C%dgsave = .TRUE.
2037      ELSE
2038          med_diag%ZE_MES_C%dgsave = .FALSE.
2039      ENDIF
2040      IF  (iom_use("ZE_MESDC")) THEN
2041          med_diag%ZE_MESDC%dgsave = .TRUE.
2042      ELSE
2043          med_diag%ZE_MESDC%dgsave = .FALSE.
2044      ENDIF
2045      IF  (iom_use("ZE_EXCR")) THEN
2046          med_diag%ZE_EXCR%dgsave = .TRUE.
2047      ELSE
2048          med_diag%ZE_EXCR%dgsave = .FALSE.
2049      ENDIF
2050      IF  (iom_use("ZE_RESP")) THEN
2051          med_diag%ZE_RESP%dgsave = .TRUE.
2052      ELSE
2053          med_diag%ZE_RESP%dgsave = .FALSE.
2054      ENDIF
2055      IF  (iom_use("ZE_GROW")) THEN
2056          med_diag%ZE_GROW%dgsave = .TRUE.
2057      ELSE
2058          med_diag%ZE_GROW%dgsave = .FALSE.
2059      ENDIF
2060      IF  (iom_use("MDETC")) THEN
2061          med_diag%MDETC%dgsave = .TRUE.
2062      ELSE
2063          med_diag%MDETC%dgsave = .FALSE.
2064      ENDIF
2065      IF  (iom_use("GMIDC")) THEN
2066          med_diag%GMIDC%dgsave = .TRUE.
2067      ELSE
2068          med_diag%GMIDC%dgsave = .FALSE.
2069      ENDIF
2070      IF  (iom_use("GMEDC")) THEN
2071          med_diag%GMEDC%dgsave = .TRUE.
2072      ELSE
2073          med_diag%GMEDC%dgsave = .FALSE.
2074      ENDIF
2075      IF  (iom_use("INT_ZMI")) THEN
2076          med_diag%INT_ZMI%dgsave = .TRUE.
2077      ELSE
2078          med_diag%INT_ZMI%dgsave = .FALSE.
2079      ENDIF
2080      IF  (iom_use("INT_ZME")) THEN
2081          med_diag%INT_ZME%dgsave = .TRUE.
2082      ELSE
2083          med_diag%INT_ZME%dgsave = .FALSE.
2084      ENDIF
2085      IF  (iom_use("INT_DET")) THEN
2086          med_diag%INT_DET%dgsave = .TRUE.
2087      ELSE
2088          med_diag%INT_DET%dgsave = .FALSE.
2089      ENDIF
2090      IF  (iom_use("INT_DTC")) THEN
2091          med_diag%INT_DTC%dgsave = .TRUE.
2092      ELSE
2093          med_diag%INT_DTC%dgsave = .FALSE.
2094      ENDIF
2095      IF  (iom_use("DMS_SURF")) THEN
2096          med_diag%DMS_SURF%dgsave = .TRUE.
2097      ELSE
2098          med_diag%DMS_SURF%dgsave = .FALSE.
2099      ENDIF
2100      IF  (iom_use("DMS_ANDR")) THEN
2101          med_diag%DMS_ANDR%dgsave = .TRUE.
2102      ELSE
2103          med_diag%DMS_ANDR%dgsave = .FALSE.
2104      ENDIF
2105      IF  (iom_use("DMS_SIMO")) THEN
2106          med_diag%DMS_SIMO%dgsave = .TRUE.
2107      ELSE
2108          med_diag%DMS_SIMO%dgsave = .FALSE.
2109      ENDIF
2110      IF  (iom_use("DMS_ARAN")) THEN
2111          med_diag%DMS_ARAN%dgsave = .TRUE.
2112      ELSE
2113          med_diag%DMS_ARAN%dgsave = .FALSE.
2114      ENDIF
2115      IF  (iom_use("DMS_HALL")) THEN
2116          med_diag%DMS_HALL%dgsave = .TRUE.
2117      ELSE
2118          med_diag%DMS_HALL%dgsave = .FALSE.
2119      ENDIF
2120      IF  (iom_use("DMS_ANDM")) THEN
2121          med_diag%DMS_ANDM%dgsave = .TRUE.
2122      ELSE
2123          med_diag%DMS_ANDM%dgsave = .FALSE.
2124      ENDIF
2125      IF  (iom_use("ATM_XCO2")) THEN
2126          med_diag%ATM_XCO2%dgsave = .TRUE.
2127      ELSE
2128          med_diag%ATM_XCO2%dgsave = .FALSE.
2129      ENDIF
2130      IF  (iom_use("OCN_FCO2")) THEN
2131          med_diag%OCN_FCO2%dgsave = .TRUE.
2132      ELSE
2133          med_diag%OCN_FCO2%dgsave = .FALSE.
2134      ENDIF
2135      IF  (iom_use("ATM_FCO2")) THEN
2136          med_diag%ATM_FCO2%dgsave = .TRUE.
2137      ELSE
2138          med_diag%ATM_FCO2%dgsave = .FALSE.
2139      ENDIF
2140      IF  (iom_use("OCN_RHOSW")) THEN
2141          med_diag%OCN_RHOSW%dgsave = .TRUE.
2142      ELSE
2143          med_diag%OCN_RHOSW%dgsave = .FALSE.
2144      ENDIF
2145      IF  (iom_use("OCN_SCHCO2")) THEN
2146          med_diag%OCN_SCHCO2%dgsave = .TRUE.
2147      ELSE
2148          med_diag%OCN_SCHCO2%dgsave = .FALSE.
2149      ENDIF
2150      IF  (iom_use("OCN_KWCO2")) THEN
2151          med_diag%OCN_KWCO2%dgsave = .TRUE.
2152      ELSE
2153          med_diag%OCN_KWCO2%dgsave = .FALSE.
2154      ENDIF
2155      IF  (iom_use("OCN_K0")) THEN
2156          med_diag%OCN_K0%dgsave = .TRUE.
2157      ELSE
2158          med_diag%OCN_K0%dgsave = .FALSE.
2159      ENDIF
2160      IF  (iom_use("CO2STARAIR")) THEN
2161          med_diag%CO2STARAIR%dgsave = .TRUE.
2162      ELSE
2163          med_diag%CO2STARAIR%dgsave = .FALSE.
2164      ENDIF
2165      IF  (iom_use("OCN_DPCO2")) THEN
2166          med_diag%OCN_DPCO2%dgsave = .TRUE.
2167      ELSE
2168          med_diag%OCN_DPCO2%dgsave = .FALSE.
2169      ENDIF
2170      !!
2171      IF  (iom_use("TPP3")) THEN
2172          med_diag%TPP3%dgsave = .TRUE.
2173      ELSE
2174          med_diag%TPP3%dgsave = .FALSE.
2175      ENDIF
2176      IF  (iom_use("DETFLUX3")) THEN
2177          med_diag%DETFLUX3%dgsave = .TRUE.
2178      ELSE
2179          med_diag%DETFLUX3%dgsave = .FALSE.
2180      ENDIF
2181      IF  (iom_use("REMIN3N")) THEN
2182          med_diag%REMIN3N%dgsave = .TRUE.
2183      ELSE
2184          med_diag%REMIN3N%dgsave = .FALSE.
2185      ENDIF
2186      IF  (iom_use("PH3")) THEN
2187          med_diag%PH3%dgsave = .TRUE.
2188      ELSE
2189          med_diag%PH3%dgsave = .FALSE.
2190      ENDIF
2191      IF  (iom_use("OM_CAL3")) THEN
2192          med_diag%OM_CAL3%dgsave = .TRUE.
2193      ELSE
2194          med_diag%OM_CAL3%dgsave = .FALSE.
2195      ENDIF
2196      !!
2197      !!----------------------------------------------------------------------
2198      !! AXY (03/11/16): add in additional CMIP6 diagnostics
2199      !!----------------------------------------------------------------------
2200      !!
2201      !! 2D fields
2202      IF  (iom_use("epC100")) THEN
2203          med_diag%epC100%dgsave = .TRUE.
2204      ELSE
2205          med_diag%epC100%dgsave = .FALSE.
2206      ENDIF
2207      IF  (iom_use("epCALC100")) THEN
2208          med_diag%epCALC100%dgsave = .TRUE.
2209      ELSE
2210          med_diag%epCALC100%dgsave = .FALSE.
2211      ENDIF
2212      IF  (iom_use("epN100")) THEN
2213          med_diag%epN100%dgsave = .TRUE.
2214      ELSE
2215          med_diag%epN100%dgsave = .FALSE.
2216      ENDIF
2217      IF  (iom_use("epSI100")) THEN
2218          med_diag%epSI100%dgsave = .TRUE.
2219      ELSE
2220          med_diag%epSI100%dgsave = .FALSE.
2221      ENDIF
2222      IF  (iom_use("FGCO2")) THEN
2223          med_diag%FGCO2%dgsave = .TRUE.
2224      ELSE
2225          med_diag%FGCO2%dgsave = .FALSE.
2226      ENDIF
2227      IF  (iom_use("INTDISSIC")) THEN
2228          med_diag%INTDISSIC%dgsave = .TRUE.
2229      ELSE
2230          med_diag%INTDISSIC%dgsave = .FALSE.
2231      ENDIF
2232      IF  (iom_use("INTDISSIN")) THEN
2233          med_diag%INTDISSIN%dgsave = .TRUE.
2234      ELSE
2235          med_diag%INTDISSIN%dgsave = .FALSE.
2236      ENDIF
2237      IF  (iom_use("INTDISSISI")) THEN
2238          med_diag%INTDISSISI%dgsave = .TRUE.
2239      ELSE
2240          med_diag%INTDISSISI%dgsave = .FALSE.
2241      ENDIF
2242      IF  (iom_use("INTTALK")) THEN
2243          med_diag%INTTALK%dgsave = .TRUE.
2244      ELSE
2245          med_diag%INTTALK%dgsave = .FALSE.
2246      ENDIF
2247      IF  (iom_use("O2min")) THEN
2248          med_diag%O2min%dgsave = .TRUE.
2249      ELSE
2250          med_diag%O2min%dgsave = .FALSE.
2251      ENDIF
2252      IF  (iom_use("ZO2min")) THEN
2253          med_diag%ZO2min%dgsave = .TRUE.
2254      ELSE
2255          med_diag%ZO2min%dgsave = .FALSE.
2256      ENDIF
2257      IF  (iom_use("FBDDTALK")) THEN
2258          med_diag%FBDDTALK%dgsave = .TRUE.
2259      ELSE
2260          med_diag%FBDDTALK%dgsave = .FALSE.
2261      ENDIF
2262      IF  (iom_use("FBDDTDIC")) THEN
2263          med_diag%FBDDTDIC%dgsave = .TRUE.
2264      ELSE
2265          med_diag%FBDDTDIC%dgsave = .FALSE.
2266      ENDIF
2267      IF  (iom_use("FBDDTDIFE")) THEN
2268          med_diag%FBDDTDIFE%dgsave = .TRUE.
2269      ELSE
2270          med_diag%FBDDTDIFE%dgsave = .FALSE.
2271      ENDIF
2272      IF  (iom_use("FBDDTDIN")) THEN
2273          med_diag%FBDDTDIN%dgsave = .TRUE.
2274      ELSE
2275          med_diag%FBDDTDIN%dgsave = .FALSE.
2276      ENDIF
2277      IF  (iom_use("FBDDTDISI")) THEN
2278          med_diag%FBDDTDISI%dgsave = .TRUE.
2279      ELSE
2280          med_diag%FBDDTDISI%dgsave = .FALSE.
2281      ENDIF
2282      !!
2283      !! 3D
2284      IF  (iom_use("TPPD3")) THEN
2285          med_diag%TPPD3%dgsave = .TRUE.
2286      ELSE
2287          med_diag%TPPD3%dgsave = .FALSE.
2288      ENDIF
2289      IF  (iom_use("BDDTALK3")) THEN
2290          med_diag%BDDTALK3%dgsave = .TRUE.
2291      ELSE
2292          med_diag%BDDTALK3%dgsave = .FALSE.
2293      ENDIF
2294      IF  (iom_use("BDDTDIC3")) THEN
2295          med_diag%BDDTDIC3%dgsave = .TRUE.
2296      ELSE
2297          med_diag%BDDTDIC3%dgsave = .FALSE.
2298      ENDIF
2299      IF  (iom_use("BDDTDIFE3")) THEN
2300          med_diag%BDDTDIFE3%dgsave = .TRUE.
2301      ELSE
2302          med_diag%BDDTDIFE3%dgsave = .FALSE.
2303      ENDIF
2304      IF  (iom_use("BDDTDIN3")) THEN
2305          med_diag%BDDTDIN3%dgsave = .TRUE.
2306      ELSE
2307          med_diag%BDDTDIN3%dgsave = .FALSE.
2308      ENDIF
2309      IF  (iom_use("BDDTDISI3")) THEN
2310          med_diag%BDDTDISI3%dgsave = .TRUE.
2311      ELSE
2312          med_diag%BDDTDISI3%dgsave = .FALSE.
2313      ENDIF
2314      IF  (iom_use("FD_NIT3")) THEN
2315          med_diag%FD_NIT3%dgsave = .TRUE.
2316      ELSE
2317          med_diag%FD_NIT3%dgsave = .FALSE.
2318      ENDIF
2319      IF  (iom_use("FD_SIL3")) THEN
2320          med_diag%FD_SIL3%dgsave = .TRUE.
2321      ELSE
2322          med_diag%FD_SIL3%dgsave = .FALSE.
2323      ENDIF
2324      IF  (iom_use("FD_CAR3")) THEN
2325          med_diag%FD_CAR3%dgsave = .TRUE.
2326      ELSE
2327          med_diag%FD_CAR3%dgsave = .FALSE.
2328      ENDIF
2329      IF  (iom_use("FD_CAL3")) THEN
2330          med_diag%FD_CAL3%dgsave = .TRUE.
2331      ELSE
2332          med_diag%FD_CAL3%dgsave = .FALSE.
2333      ENDIF
2334      IF  (iom_use("CO33")) THEN
2335          med_diag%CO33%dgsave = .TRUE.
2336      ELSE
2337          med_diag%CO33%dgsave = .FALSE.
2338      ENDIF
2339      IF  (iom_use("CO3SATARAG3")) THEN
2340          med_diag%CO3SATARAG3%dgsave = .TRUE.
2341      ELSE
2342          med_diag%CO3SATARAG3%dgsave = .FALSE.
2343      ENDIF
2344      IF  (iom_use("CO3SATCALC3")) THEN
2345          med_diag%CO3SATCALC3%dgsave = .TRUE.
2346      ELSE
2347          med_diag%CO3SATCALC3%dgsave = .FALSE.
2348      ENDIF
2349      IF  (iom_use("DCALC3")) THEN
2350          med_diag%DCALC3%dgsave = .TRUE.
2351      ELSE
2352          med_diag%DCALC3%dgsave = .FALSE.
2353      ENDIF
2354      IF  (iom_use("EXPC3")) THEN
2355          med_diag%EXPC3%dgsave = .TRUE.
2356      ELSE
2357          med_diag%EXPC3%dgsave = .FALSE.
2358      ENDIF
2359      IF  (iom_use("EXPN3")) THEN
2360          med_diag%EXPN3%dgsave = .TRUE.
2361      ELSE
2362          med_diag%EXPN3%dgsave = .FALSE.
2363      ENDIF
2364      IF  (iom_use("FEDISS3")) THEN
2365          med_diag%FEDISS3%dgsave = .TRUE.
2366      ELSE
2367          med_diag%FEDISS3%dgsave = .FALSE.
2368      ENDIF
2369      IF  (iom_use("FESCAV3")) THEN
2370          med_diag%FESCAV3%dgsave = .TRUE.
2371      ELSE
2372          med_diag%FESCAV3%dgsave = .FALSE.
2373      ENDIF
2374      IF  (iom_use("MIGRAZP3")) THEN
2375          med_diag%MIGRAZP3%dgsave = .TRUE.
2376      ELSE
2377          med_diag%MIGRAZP3%dgsave = .FALSE.
2378      ENDIF
2379      IF  (iom_use("MIGRAZD3")) THEN
2380          med_diag%MIGRAZD3%dgsave = .TRUE.
2381      ELSE
2382          med_diag%MIGRAZD3%dgsave = .FALSE.
2383      ENDIF
2384      IF  (iom_use("MEGRAZP3")) THEN
2385          med_diag%MEGRAZP3%dgsave = .TRUE.
2386      ELSE
2387          med_diag%MEGRAZP3%dgsave = .FALSE.
2388      ENDIF
2389      IF  (iom_use("MEGRAZD3")) THEN
2390          med_diag%MEGRAZD3%dgsave = .TRUE.
2391      ELSE
2392          med_diag%MEGRAZD3%dgsave = .FALSE.
2393      ENDIF
2394      IF  (iom_use("MEGRAZZ3")) THEN
2395          med_diag%MEGRAZZ3%dgsave = .TRUE.
2396      ELSE
2397          med_diag%MEGRAZZ3%dgsave = .FALSE.
2398      ENDIF
2399      IF  (iom_use("O2SAT3")) THEN
2400          med_diag%O2SAT3%dgsave = .TRUE.
2401      ELSE
2402          med_diag%O2SAT3%dgsave = .FALSE.
2403      ENDIF
2404      IF  (iom_use("PBSI3")) THEN
2405          med_diag%PBSI3%dgsave = .TRUE.
2406      ELSE
2407          med_diag%PBSI3%dgsave = .FALSE.
2408      ENDIF
2409      IF  (iom_use("PCAL3")) THEN
2410          med_diag%PCAL3%dgsave = .TRUE.
2411      ELSE
2412          med_diag%PCAL3%dgsave = .FALSE.
2413      ENDIF
2414      IF  (iom_use("REMOC3")) THEN
2415          med_diag%REMOC3%dgsave = .TRUE.
2416      ELSE
2417          med_diag%REMOC3%dgsave = .FALSE.
2418      ENDIF
2419      IF  (iom_use("PNLIMJ3")) THEN
2420          med_diag%PNLIMJ3%dgsave = .TRUE.
2421      ELSE
2422          med_diag%PNLIMJ3%dgsave = .FALSE.
2423      ENDIF
2424      IF  (iom_use("PNLIMN3")) THEN
2425          med_diag%PNLIMN3%dgsave = .TRUE.
2426      ELSE
2427          med_diag%PNLIMN3%dgsave = .FALSE.
2428      ENDIF
2429      IF  (iom_use("PNLIMFE3")) THEN
2430          med_diag%PNLIMFE3%dgsave = .TRUE.
2431      ELSE
2432          med_diag%PNLIMFE3%dgsave = .FALSE.
2433      ENDIF
2434      IF  (iom_use("PDLIMJ3")) THEN
2435          med_diag%PDLIMJ3%dgsave = .TRUE.
2436      ELSE
2437          med_diag%PDLIMJ3%dgsave = .FALSE.
2438      ENDIF
2439      IF  (iom_use("PDLIMN3")) THEN
2440          med_diag%PDLIMN3%dgsave = .TRUE.
2441      ELSE
2442          med_diag%PDLIMN3%dgsave = .FALSE.
2443      ENDIF
2444      IF  (iom_use("PDLIMFE3")) THEN
2445          med_diag%PDLIMFE3%dgsave = .TRUE.
2446      ELSE
2447          med_diag%PDLIMFE3%dgsave = .FALSE.
2448      ENDIF
2449      IF  (iom_use("PDLIMSI3")) THEN
2450          med_diag%PDLIMSI3%dgsave = .TRUE.
2451      ELSE
2452          med_diag%PDLIMSI3%dgsave = .FALSE.
2453      ENDIF
2454
2455   END SUBROUTINE   trc_nam_iom_medusa
2456   
2457#else
2458   !!----------------------------------------------------------------------
2459   !!  Dummy module :                                             No MEDUSA
2460   !!----------------------------------------------------------------------
2461CONTAINS
2462   SUBROUTINE trc_nam_medusa                      ! Empty routine
2463   END  SUBROUTINE  trc_nam_medusa
2464#endif 
2465
2466   !!======================================================================
2467END MODULE trcnam_medusa
Note: See TracBrowser for help on using the repository browser.