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

Context Navigation

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

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

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

Download in other formats: