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/NERC/dev_r5518_GO6_MEDUSA_DMS_TUNING/NEMOGCM/NEMO/TOP_SRC/MEDUSA – NEMO

Context Navigation

source: branches/NERC/dev_r5518_GO6_MEDUSA_DMS_TUNING/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcnam_medusa.F90 @ 9187

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

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

Download in other formats: