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

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.

trcnam_medusa.F90 in branches/UKMO/dev_r5518_GO6_package_FOAMv14_biophys/NEMOGCM/NEMO/TOP_SRC/MEDUSA – NEMO

Context Navigation

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

Last change on this file since 13695 was 13695, checked in by dford, 3 years ago
Add 3D PFT-specific grazing diagnostics
File size: 89.3 KB

Line
1	MODULE 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
46	CONTAINS
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)
268	903 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 )
272	904 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)
1050	905 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 )
1054	906 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)
1106	907 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 )
1110	908 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("MEGRAZPN3")) THEN
2348	med_diag%MEGRAZPN3%dgsave = .TRUE.
2349	ELSE
2350	med_diag%MEGRAZPN3%dgsave = .FALSE.
2351	ENDIF
2352	IF (iom_use("MEGRAZPD3")) THEN
2353	med_diag%MEGRAZPD3%dgsave = .TRUE.
2354	ELSE
2355	med_diag%MEGRAZPD3%dgsave = .FALSE.
2356	ENDIF
2357	IF (iom_use("MEGRAZD3")) THEN
2358	med_diag%MEGRAZD3%dgsave = .TRUE.
2359	ELSE
2360	med_diag%MEGRAZD3%dgsave = .FALSE.
2361	ENDIF
2362	IF (iom_use("MEGRAZZ3")) THEN
2363	med_diag%MEGRAZZ3%dgsave = .TRUE.
2364	ELSE
2365	med_diag%MEGRAZZ3%dgsave = .FALSE.
2366	ENDIF
2367	IF (iom_use("O2SAT3")) THEN
2368	med_diag%O2SAT3%dgsave = .TRUE.
2369	ELSE
2370	med_diag%O2SAT3%dgsave = .FALSE.
2371	ENDIF
2372	IF (iom_use("PBSI3")) THEN
2373	med_diag%PBSI3%dgsave = .TRUE.
2374	ELSE
2375	med_diag%PBSI3%dgsave = .FALSE.
2376	ENDIF
2377	IF (iom_use("PCAL3")) THEN
2378	med_diag%PCAL3%dgsave = .TRUE.
2379	ELSE
2380	med_diag%PCAL3%dgsave = .FALSE.
2381	ENDIF
2382	IF (iom_use("REMOC3")) THEN
2383	med_diag%REMOC3%dgsave = .TRUE.
2384	ELSE
2385	med_diag%REMOC3%dgsave = .FALSE.
2386	ENDIF
2387	IF (iom_use("PNLIMJ3")) THEN
2388	med_diag%PNLIMJ3%dgsave = .TRUE.
2389	ELSE
2390	med_diag%PNLIMJ3%dgsave = .FALSE.
2391	ENDIF
2392	IF (iom_use("PNLIMN3")) THEN
2393	med_diag%PNLIMN3%dgsave = .TRUE.
2394	ELSE
2395	med_diag%PNLIMN3%dgsave = .FALSE.
2396	ENDIF
2397	IF (iom_use("PNLIMFE3")) THEN
2398	med_diag%PNLIMFE3%dgsave = .TRUE.
2399	ELSE
2400	med_diag%PNLIMFE3%dgsave = .FALSE.
2401	ENDIF
2402	IF (iom_use("PDLIMJ3")) THEN
2403	med_diag%PDLIMJ3%dgsave = .TRUE.
2404	ELSE
2405	med_diag%PDLIMJ3%dgsave = .FALSE.
2406	ENDIF
2407	IF (iom_use("PDLIMN3")) THEN
2408	med_diag%PDLIMN3%dgsave = .TRUE.
2409	ELSE
2410	med_diag%PDLIMN3%dgsave = .FALSE.
2411	ENDIF
2412	IF (iom_use("PDLIMFE3")) THEN
2413	med_diag%PDLIMFE3%dgsave = .TRUE.
2414	ELSE
2415	med_diag%PDLIMFE3%dgsave = .FALSE.
2416	ENDIF
2417	IF (iom_use("PDLIMSI3")) THEN
2418	med_diag%PDLIMSI3%dgsave = .TRUE.
2419	ELSE
2420	med_diag%PDLIMSI3%dgsave = .FALSE.
2421	ENDIF
2422
2423	END SUBROUTINE trc_nam_iom_medusa
2424
2425	#else
2426	!!----------------------------------------------------------------------
2427	!! Dummy module : No MEDUSA
2428	!!----------------------------------------------------------------------
2429	CONTAINS
2430	SUBROUTINE trc_nam_medusa ! Empty routine
2431	END SUBROUTINE trc_nam_medusa
2432	#endif
2433
2434	!!======================================================================
2435	END MODULE trcnam_medusa

Note: See TracBrowser for help on using the repository browser.

Download in other formats: