source: branches/UKMO/dev_r5518_GO6_package_FOAMv14_biophys/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcnam_medusa.F90 @ 13677

Last change on this file since 13677 was 13677, checked in by dford, 7 months ago

Add option for MEDUSA to use RGB light scheme.

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