Context Navigation

trcnam_medusa.F90 @ 8280

Last change on this file since 8280 was 8131, checked in by jpalmier, 7 years ago
JPALM -- DMS add an updated version of ANDERSON sheme
File size: 89.7 KB

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

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

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

Context Navigation

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

Download in other formats: