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

Context Navigation

source: branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcnam_medusa.F90 @ 6715

Last change on this file since 6715 was 6715, checked in by jpalmier, 8 years ago

JPALM -- 16-06-2016 -- MEDUSA branch update :

-- pass co2 flux and dms_surf through restart for atm coupling.
-- introduce CFC cycle for dynamic evolution comparison
-- add Tim Graham Age tracer
-- include MEDUSA Q10 modif
-- svn-key removed
-- still need debug stage

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

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

Download in other formats: