Context Navigation

Back to Ticket #805

Ticket #805: checkdefs.F

File checkdefs.F, 109.1 KB (added by m.hadfield, 7 years ago)

Line
1	#include "cppdefs.h"
2	SUBROUTINE checkdefs
3	!
4	!svn $Id$
5	!================================================== Hernan G. Arango ===
6	! Copyright (c) 2002-2019 The ROMS/TOMS Group !
7	! Licensed under a MIT/X style license !
8	! See License_ROMS.txt !
9	!=======================================================================
10	! !
11	! This subroutine checks activated C-preprocessing options for !
12	! consistency. !
13	! !
14	!=======================================================================
15	!
16	USE mod_param
17	USE mod_parallel
18	USE mod_iounits
19	USE mod_scalars
20	USE mod_strings
21	!
22	USE strings_mod, ONLY : uppercase
23	!
24	implicit none
25	!
26	! Local variable declarations.
27	!
28	integer :: iatms = 0
29	integer :: ibbl = 0
30	integer :: ibiology = 0
31	integer :: idriver = 0
32	integer :: itrcHadv = 0
33	integer :: itrcVadv = 0
34	integer :: itrcHadvtl = 0
35	integer :: itrcVadvtl = 0
36	integer :: ivelHadv = 0
37	integer :: ivelVadv = 0
38	integer :: ivmix = 0
39	integer :: nearshore = 0
40
41	integer :: is, lstr, ng
42	!
43	!-----------------------------------------------------------------------
44	! Report activated C-preprocessing options.
45	!-----------------------------------------------------------------------
46	!
47	Coptions=' '
48	IF (Master) WRITE (stdout,10)
49	10 FORMAT (/,' Activated C-preprocessing Options:',/)
50	20 FORMAT (1x,a,t26,a)
51	!
52	IF (Master) THEN
53	WRITE (stdout,20) TRIM(ADJUSTL(MyAppCPP)), TRIM(ADJUSTL(title))
54	END IF
55	is=LEN_TRIM(Coptions)+1
56	lstr=LEN_TRIM(MyAppCPP)
57	Coptions(is:is+lstr)=TRIM(ADJUSTL(MyAppCPP))
58	is=LEN_TRIM(Coptions)+1
59	Coptions(is:is)=','
60
61	#if defined AD_AVERAGES && defined ADJOINT
62	!
63	IF (Master) WRITE (stdout,20) 'AD_AVERAGES', &
64	& 'Writing out time-averaged adjoint model fields'
65	is=LEN_TRIM(Coptions)+1
66	Coptions(is:is+13)=' AD_AVERAGES,'
67	#endif
68	#if defined ADD_FSOBC && defined SSH_TIDES
69	!
70	IF (Master) WRITE (stdout,20) 'ADD_FSOBC', &
71	& 'Adding tidal elevation to processed OBC data'
72	is=LEN_TRIM(Coptions)+1
73	Coptions(is:is+11)=' ADD_FSOBC,'
74	#endif
75	#if defined ADD_M2OBC && defined UV_TIDES
76	!
77	IF (Master) WRITE (stdout,20) 'ADD_M2OBC', &
78	& 'Adding tidal currents to processed OBC data'
79	is=LEN_TRIM(Coptions)+1
80	Coptions(is:is+11)=' ADD_M2OBC,'
81	#endif
82	#ifdef ADJOINT
83	!
84	IF (Master) WRITE (stdout,20) 'ADJOINT', &
85	& 'Adjoint Model'
86	is=LEN_TRIM(Coptions)+1
87	Coptions(is:is+9)=' ADJOINT,'
88	#endif
89	#if defined ADJUST_BOUNDARY && \
90	(defined CORRELATION \|\| defined IS4DVAR \|\| \
91	defined IS4DVAR_SENSITIVITY \|\| defined SENSITIVITY_4DVAR \|\| \
92	defined W4DPSAS \|\| defined W4DVAR)
93	!
94	IF (Master) WRITE (stdout,20) 'ADJUST_BOUNDARY', &
95	& 'Including boundary conditions in 4DVar state estimation'
96	is=LEN_TRIM(Coptions)+1
97	Coptions(is:is+17)=' ADJUST_BOUNDARY,'
98	#endif
99	#if defined ADJUST_STFLUX && \
100	(defined CORRELATION \|\| defined IS4DVAR \|\| \
101	defined IS4DVAR_SENSITIVITY \|\| defined SENSITIVITY_4DVAR \|\| \
102	defined W4DPSAS \|\| defined W4DVAR)
103	!
104	IF (Master) WRITE (stdout,20) 'ADJUST_STFLUX', &
105	& 'Including surface tracer flux in 4DVar state estimation'
106	is=LEN_TRIM(Coptions)+1
107	Coptions(is:is+15)=' ADJUST_STFLUX,'
108	#endif
109	#if defined ADJUST_WSTRESS && \
110	(defined CORRELATION \|\| defined IS4DVAR \|\| \
111	defined IS4DVAR_SENSITIVITY \|\| defined SENSITIVITY_4DVAR \|\| \
112	defined W4DPSAS \|\| defined W4DVAR)
113	!
114	IF (Master) WRITE (stdout,20) 'ADJUST_WSTRESS', &
115	& 'Including surface wind stress in 4DVar state estimation'
116	is=LEN_TRIM(Coptions)+1
117	Coptions(is:is+16)=' ADJUST_WSTRESS,'
118	#endif
119	#if defined AD_IMPULSE
120	!
121	IF (Master) WRITE (stdout,20) 'AD_IMPULSE', &
122	& 'Force adjoint model with intermittent impulses'
123	is=LEN_TRIM(Coptions)+1
124	Coptions(is:is+12)=' AD_IMPULSE,'
125	#endif
126	#ifdef ADM_DRIVER
127	!
128	IF (Master) WRITE (stdout,20) 'ADM_DRIVER', &
129	& 'Generic adjoint model driver'
130	is=LEN_TRIM(Coptions)+1
131	Coptions(is:is+12)=' ADM_DRIVER,'
132	idriver=idriver+1
133	#endif
134	#ifdef AD_SENSITIVITY
135	!
136	IF (Master) WRITE (stdout,20) 'AD_SENSITIVITY', &
137	& 'Adjoint Sensitivity Analysis'
138	is=LEN_TRIM(Coptions)+1
139	Coptions(is:is+16)=' AD_SENSITIVITY,'
140	idriver=idriver+1
141	#endif
142	#ifdef AFT_EIGENMODES
143	!
144	IF (Master) WRITE (stdout,20) 'AFT_EIGENMODES', &
145	& 'Adjoint Finite Time Eigenvalues'
146	is=LEN_TRIM(Coptions)+1
147	Coptions(is:is+16)=' AFT_EIGENMODES,'
148	idriver=idriver+1
149	#endif
150	#if defined AGE_MEAN && defined T_PASSIVE && defined SOLVE3D
151	!
152	IF (Master) WRITE (stdout,20) 'AGE_MEAN', &
153	& 'Computing Mean Age of inert passive tracer'
154	is=LEN_TRIM(Coptions)+1
155	Coptions(is:is+10)=' AGE_MEAN,'
156	idriver=idriver+1
157	#endif
158	#if defined ALBEDO && defined ANA_SRFLUX
159	!
160	IF (Master) WRITE (stdout,20) 'ALBEDO', &
161	& 'Shortwave radiation from albedo equation'
162	is=LEN_TRIM(Coptions)+1
163	Coptions(is:is+8)=' ALBEDO,'
164	#endif
165	#if defined BIOLOGY && defined ANA_BIOLOGY
166	!
167	IF (Master) WRITE (stdout,20) 'ANA_BIOLOGY', &
168	& 'Analytical biology initial conditions'
169	is=LEN_TRIM(Coptions)+1
170	Coptions(is:is+13)=' ANA_BIOLOGY,'
171	#endif
172	#if defined BIOLOGY \|\| defined SEDIMENT \|\| defined T_PASSIVE
173	# ifdef ANA_BPFLUX
174	!
175	IF (Master) WRITE (stdout,20) 'ANA_BPFLUX', &
176	& 'Analytical bottom passive tracers fluxes'
177	is=LEN_TRIM(Coptions)+1
178	Coptions(is:is+12)=' ANA_BPFLUX,'
179	# endif
180	#endif
181	#ifdef ANA_BSFLUX
182	!
183	IF (Master) WRITE (stdout,20) 'ANA_BSFLUX', &
184	& 'Analytical kinematic bottom salinity flux'
185	is=LEN_TRIM(Coptions)+1
186	Coptions(is:is+12)=' ANA_BSFLUX,'
187	#endif
188	#ifdef ANA_BTFLUX
189	!
190	IF (Master) WRITE (stdout,20) 'ANA_BTFLUX', &
191	& 'Analytical kinematic bottom temperature flux'
192	is=LEN_TRIM(Coptions)+1
193	Coptions(is:is+12)=' ANA_BTFLUX,'
194	#endif
195	#if defined ANA_CLOUD && defined BULK_FLUXES
196	!
197	IF (Master) WRITE (stdout,20) 'ANA_CLOUD', &
198	& 'Analytical cloud fraction'
199	is=LEN_TRIM(Coptions)+1
200	Coptions(is:is+11)=' ANA_CLOUD,'
201	#endif
202	#if defined ANA_DRAG && defined UV_DRAG_GRID
203	!
204	IF (Master) WRITE (stdout,20) 'ANA_DRAG_GRID', &
205	# if defined UV_LOGDRAG
206	& 'Analytical spatially varying bottom roughness length'
207	# elif defined UV_LDRAG
208	& 'Analytical spatially varying linear drag coefficient'
209	# elif defined UV_QDRAG
210	& 'Analytical spatially varying quadratic drag coefficient'
211	# endif
212	is=LEN_TRIM(Coptions)+1
213	Coptions(is:is+10)=' ANA_DRAG,'
214	#endif
215	#ifdef ANA_DIAG
216	!
217	IF (Master) WRITE (stdout,20) 'ANA_DIAG', &
218	& 'Customized diagnostics'
219	is=LEN_TRIM(Coptions)+1
220	Coptions(is:is+10)=' ANA_DIAG,'
221	#endif
222	#if defined ANA_DQDSST && defined QCORRECTION
223	!
224	IF (Master) WRITE (stdout,20) 'ANA_DQDSST', &
225	& 'Analytical surface heat flux sensitivity to SST'
226	is=LEN_TRIM(Coptions)+1
227	Coptions(is:is+12)=' ANA_DQDSST,'
228	#endif
229	#ifdef ANA_FSOBC
230	!
231	IF (Master) WRITE (stdout,20) 'ANA_FSOBC', &
232	& 'Analytical free-surface boundary conditions'
233	is=LEN_TRIM(Coptions)+1
234	Coptions(is:is+11)=' ANA_FSOBC,'
235	#endif
236	#ifdef ANA_GRID
237	!
238	IF (Master) WRITE (stdout,20) 'ANA_GRID', &
239	& 'Analytical grid set-up'
240	is=LEN_TRIM(Coptions)+1
241	Coptions(is:is+10)=' ANA_GRID,'
242	#endif
243	#if defined ANA_HUMIDITY && defined BULK_FLUXES
244	!
245	IF (Master) WRITE (stdout,20) 'ANA_HUMIDITY', &
246	& 'Analytical surface air humidity'
247	is=LEN_TRIM(Coptions)+1
248	Coptions(is:is+14)=' ANA_HUMIDITY,'
249	#endif
250	#ifdef ANA_INITIAL
251	!
252	IF (Master) WRITE (stdout,20) 'ANA_INITIAL', &
253	& 'Analytical initial conditions'
254	is=LEN_TRIM(Coptions)+1
255	Coptions(is:is+13)=' ANA_INITIAL,'
256	#endif
257	#ifdef ANA_M2CLIMA
258	!
259	IF (Master) WRITE (stdout,20) 'ANA_M2CLIMA', &
260	& 'Analytical 2D momentum climatology'
261	is=LEN_TRIM(Coptions)+1
262	Coptions(is:is+13)=' ANA_M2CLIMA,'
263	#endif
264	#ifdef ANA_M2OBC
265	!
266	IF (Master) WRITE (stdout,20) 'ANA_M2OBC', &
267	& 'Analytical 2D momentum boundary conditions'
268	is=LEN_TRIM(Coptions)+1
269	Coptions(is:is+11)=' ANA_M2OBC,'
270	#endif
271	#if defined ANA_M3CLIMA && defined SOLVE3D
272	!
273	IF (Master) WRITE (stdout,20) 'ANA_M3CLIMA', &
274	& 'Analytical 3D momentum climatology'
275	is=LEN_TRIM(Coptions)+1
276	Coptions(is:is+13)=' ANA_M3CLIMA,'
277	#endif
278	#if defined ANA_M3OBC && defined SOLVE3D
279	!
280	IF (Master) WRITE (stdout,20) 'ANA_M3OBC', &
281	& 'Analytical 3D momentum boundary conditions'
282	is=LEN_TRIM(Coptions)+1
283	Coptions(is:is+11)=' ANA_M3OBC,'
284	#endif
285	#ifdef ANA_MASK
286	!
287	IF (Master) WRITE (stdout,20) 'ANA_MASK', &
288	& 'Analytical Land/Sea Masking'
289	is=LEN_TRIM(Coptions)+1
290	Coptions(is:is+10)=' ANA_MASK,'
291	#endif
292	#ifdef ANA_NUDGCOEF
293	!
294	IF (ANY(Lnudging(:))) THEN
295	IF (Master) WRITE (stdout,20) 'ANA_NUDGCOEF', &
296	& 'Analytical spatially varying nudging time-scales'
297	is=LEN_TRIM(Coptions)+1
298	Coptions(is:is+14)=' ANA_NUDGCOEF,'
299	END IF
300	#endif
301	#if defined ANA_PAIR && defined BULK_FLUXES
302	!
303	IF (Master) WRITE (stdout,20) 'ANA_PAIR', &
304	& 'Analytical surface air pressure'
305	is=LEN_TRIM(Coptions)+1
306	Coptions(is:is+10)=' ANA_PAIR,'
307	#endif
308	#if defined ANA_PASSIVE && defined T_PASSIVE && defined SOLVE3D
309	!
310	IF (Master) WRITE (stdout,20) 'ANA_PASSIVE', &
311	& 'Analytical initial conditions for inert tracers'
312	is=LEN_TRIM(Coptions)+1
313	Coptions(is:is+13)=' ANA_PASSIVE,'
314	#endif
315	#if defined ANA_PERTURB
316	!
317	IF (Master) WRITE (stdout,20) 'ANA_PERTURB', &
318	& 'Perturb initial conditions'
319	is=LEN_TRIM(Coptions)+1
320	Coptions(is:is+13)=' ANA_PERTURB,'
321	#endif
322	#ifdef ANA_PSOURCE
323	!
324	IF (Master) WRITE (stdout,20) 'ANA_PSOURCE', &
325	& 'Analytical point sources and sinks'
326	is=LEN_TRIM(Coptions)+1
327	Coptions(is:is+13)=' ANA_PSOURCE,'
328	#endif
329	#if defined ANA_RAIN && defined BULK_FLUXES
330	!
331	IF (Master) WRITE (stdout,20) 'ANA_RAIN', &
332	& 'Analytical rain fall rate'
333	is=LEN_TRIM(Coptions)+1
334	Coptions(is:is+10)=' ANA_RAIN,'
335	#endif
336	#if defined ANA_RESPIRATION && defined HYPOXIA_SRM
337	!
338	IF (Master) WRITE (stdout,20) 'ANA_RESPIRATION', &
339	& 'Analytical respiration rate'
340	is=LEN_TRIM(Coptions)+1
341	Coptions(is:is+17)=' ANA_RESPIRATION,'
342	#endif
343	#if defined ANA_SCOPE && \
344	(defined AD_SENSITIVITY \|\| defined IS4DVAR_SENSITIVITY \|\| \
345	defined OPT_OBSERVATIONS \|\| defined SENSITIVITY_4DVAR \|\| \
346	defined SO_SEMI)
347	!
348	IF (Master) WRITE (stdout,20) 'ANA_SCOPE', &
349	& 'Analytical spatial scope masking arrays'
350	is=LEN_TRIM(Coptions)+1
351	Coptions(is:is+12)=' ANA_SCOPE,'
352	#endif
353	#if defined ANA_SEDIMENT && (defined SEDIMENT \|\| defined BBL_MODEL)
354	!
355	IF (Master) WRITE (stdout,20) 'ANA_SEDIMENT', &
356	& 'Analytical sediment initial conditions'
357	is=LEN_TRIM(Coptions)+1
358	Coptions(is:is+14)=' ANA_SEDIMENT,'
359	#endif
360	#ifdef ANA_SMFLUX
361	!
362	IF (Master) WRITE (stdout,20) 'ANA_SMFLUX', &
363	& 'Analytical kinematic surface momentum flux'
364	is=LEN_TRIM(Coptions)+1
365	Coptions(is:is+12)=' ANA_SMFLUX,'
366	#endif
367	#if defined BIOLOGY \|\| defined SEDIMENT \|\| defined T_PASSIVE
368	# ifdef ANA_SPFLUX
369	!
370	IF (Master) WRITE (stdout,20) 'ANA_SPFLUX', &
371	& 'Analytical surface passive tracer fluxes'
372	is=LEN_TRIM(Coptions)+1
373	Coptions(is:is+12)=' ANA_SPFLUX,'
374	# endif
375	#endif
376	#ifdef ANA_SPINNING
377	!
378	IF (Master) WRITE (stdout,20) 'ANA_SPINNING', &
379	& 'Analytical time-varying rotation force'
380	is=LEN_TRIM(Coptions)+1
381	Coptions(is:is+14)=' ANA_SPINNING,'
382	#endif
383	#ifdef ANA_SPONGE
384	!
385	IF (ANY(Lsponge)) THEN
386	IF (Master) WRITE (stdout,20) 'ANA_SPONGE', &
387	& 'Analytical enhanced viscosity/diffusion sponge'
388	is=LEN_TRIM(Coptions)+1
389	Coptions(is:is+12)=' ANA_SPONGE,'
390	END IF
391	#endif
392	#ifdef ANA_SRFLUX
393	!
394	IF (Master) WRITE (stdout,20) 'ANA_SRFLUX', &
395	& 'Analytical kinematic shortwave radiation flux'
396	is=LEN_TRIM(Coptions)+1
397	Coptions(is:is+12)=' ANA_SRFLUX,'
398	#endif
399	#ifdef ANA_SSH
400	!
401	IF (Master) WRITE (stdout,20) 'ANA_SSH', &
402	& 'Analytical sea surface height climatology'
403	is=LEN_TRIM(Coptions)+1
404	Coptions(is:is+9)=' ANA_SSH,'
405	#endif
406	#ifdef ANA_SSS
407	!
408	IF (Master) WRITE (stdout,20) 'ANA_SSS', &
409	& 'Analytical sea surface salinity'
410	is=LEN_TRIM(Coptions)+1
411	Coptions(is:is+9)=' ANA_SSS,'
412	#endif
413	#if defined ANA_SST && defined QCORRECTION
414	!
415	IF (Master) WRITE (stdout,20) 'ANA_SST', &
416	& 'Analytical sea surface temperature, SST'
417	is=LEN_TRIM(Coptions)+1
418	Coptions(is:is+9)=' ANA_SST,'
419	#endif
420	#ifdef ANA_SSFLUX
421	!
422	IF (Master) WRITE (stdout,20) 'ANA_SSFLUX', &
423	& 'Analytical kinematic surface salinity flux'
424	is=LEN_TRIM(Coptions)+1
425	Coptions(is:is+12)=' ANA_SSFLUX,'
426	#endif
427	#ifdef ANA_STFLUX
428	!
429	IF (Master) WRITE (stdout,20) 'ANA_STFLUX', &
430	& 'Analytical kinematic surface temperature flux'
431	is=LEN_TRIM(Coptions)+1
432	Coptions(is:is+12)=' ANA_STFLUX,'
433	#endif
434	#if defined ANA_TCLIMA && defined SOLVE3D
435	!
436	IF (Master) WRITE (stdout,20) 'ANA_TCLIMA', &
437	& 'Analytical tracer climatology'
438	is=LEN_TRIM(Coptions)+1
439	Coptions(is:is+12)=' ANA_TCLIMA,'
440	#endif
441	#if defined ANA_TOBC && defined SOLVE3D
442	!
443	IF (Master) WRITE (stdout,20) 'ANA_TOBC', &
444	& 'Analytical tracers boundary conditions'
445	is=LEN_TRIM(Coptions)+1
446	Coptions(is:is+10)=' ANA_TOBC,'
447	#endif
448	#if defined ANA_VMIX && defined SOLVE3D
449	!
450	IF (Master) WRITE (stdout,20) 'ANA_VMIX', &
451	& 'Analytical vertical mixing coefficients'
452	is=LEN_TRIM(Coptions)+1
453	Coptions(is:is+10)=' ANA_VMIX,'
454	ivmix=ivmix+1
455	#endif
456	#if defined ANA_WINDS && defined BULK_FLUXES
457	!
458	IF (Master) WRITE (stdout,20) 'ANA_WINDS', &
459	& 'Analytical surface wind components'
460	is=LEN_TRIM(Coptions)+1
461	Coptions(is:is+11)=' ANA_WINDS,'
462	#endif
463	#if defined ANA_WTYPE && defined WTYPE_GRID && \
464	(defined LMD_SKPP \|\| defined SOLAR_SOURCE)
465	!
466	IF (Master) WRITE (stdout,20) 'ANA_WTYPE', &
467	& 'Analytical spatially varying Jerlov water type index'
468	is=LEN_TRIM(Coptions)+1
469	Coptions(is:is+11)=' ANA_WTYPE,'
470	#endif
471	#ifdef ANA_WWAVE
472	!
473	IF (Master) WRITE (stdout,20) 'ANA_WWAVE', &
474	& 'Analytical wind induced waves'
475	is=LEN_TRIM(Coptions)+1
476	Coptions(is:is+11)=' ANA_WWAVE,'
477	#endif
478	#ifdef ARRAY_MODES
479	!
480	IF (Master) WRITE (stdout,20) 'ARRAY_MODES', &
481	& 'Stabilized Representer Matrix Array Modes'
482	is=LEN_TRIM(Coptions)+1
483	Coptions(is:is+13)=' ARRAY_MODES,'
484	idriver=idriver+1
485	# ifdef ARRAY_MODES_SPLIT
486	!
487	IF (Master) WRITE (stdout,20) 'ARRAY_MODES_SPLIT', &
488	& 'Separate analysis due to IC, surface forcing, and OBC'
489	is=LEN_TRIM(Coptions)+1
490	Coptions(is:is+19)=' ARRAY_MODES_SPLIT,'
491	# endif
492	#endif
493	#if defined DISTRIBUTE && defined NESTING
494	# if defined ASSEMBLE_ALLGATHER
495	!
496	IF (Master) WRITE (stdout,20) 'ASSEMBLE_ALLGATHER', &
497	& 'Using mpi_allgather in mp_assemble routine'
498	is=LEN_TRIM(Coptions)+1
499	Coptions(is:is+20)=' ASSEMBLE_ALLGATHER,'
500	# elif defined ASSEMBLE_ALLREDUCE
501	!
502	IF (Master) WRITE (stdout,20) 'ASSEMBLE_ALLREDUCE', &
503	& 'Using mpi_allreduce in mp_assemble routine'
504	is=LEN_TRIM(Coptions)+1
505	Coptions(is:is+20)=' ASSEMBLE_ALLGATHER,'
506	# else
507	!
508	IF (Master) WRITE (stdout,20) '!ASSEMBLE_ALL...', &
509	& 'Using mpi_isend/mpi_recv in mp_assemble routine'
510	is=LEN_TRIM(Coptions)+1
511	Coptions(is:is+18)=' !ASSEMBLE_ALL...,'
512	# endif
513	#endif
514	#ifdef ASSUMED_SHAPE
515	!
516	IF (Master) WRITE (stdout,20) 'ASSUMED_SHAPE', &
517	& 'Using assumed-shape arrays'
518	is=LEN_TRIM(Coptions)+1
519	Coptions(is:is+15)=' ASSUMED_SHAPE,'
520	#endif
521	#ifdef ATM_PRESS
522	!
523	IF (Master) WRITE (stdout,20) 'ATM_PRESS', &
524	& 'Impose atmospheric pressure onto sea surface'
525	is=LEN_TRIM(Coptions)+1
526	Coptions(is:is+11)=' ATM_PRESS,'
527	# ifdef PRESS_COMPENSATE
528	IF (Master) WRITE (stdout,20) 'PRESS_COMPENSATE', &
529	& 'Compensate for boundary data without atm. pressure effect'
530	is=LEN_TRIM(Coptions)+1
531	Coptions(is:is+18)=' PRESS_COMPENSATE,'
532	# endif
533	#endif
534	#ifdef AVERAGES
535	!
536	IF (Master) WRITE (stdout,20) 'AVERAGES', &
537	& 'Writing out time-averaged nonlinear model fields'
538	is=LEN_TRIM(Coptions)+1
539	Coptions(is:is+10)=' AVERAGES,'
540	# if defined AVERAGES_DETIDE && (defined SSH_TIDES \|\| defined UV_TIDES)
541	!
542	IF (Master) WRITE (stdout,20) 'AVERAGES_DETIDE', &
543	& 'Writing out time-averaged nonlinear model detided fields'
544	is=LEN_TRIM(Coptions)+1
545	Coptions(is:is+17)=' AVERAGES_DETIDE,'
546	# endif
547	#endif
548	#if defined BACKGROUND && defined FOUR_DVAR
549	!
550	IF (Master) WRITE (stdout,20) 'BACKGROUND', &
551	& 'Include background cost function'
552	is=LEN_TRIM(Coptions)+1
553	Coptions(is:is+12)=' BACKGROUND,'
554	#endif
555	#ifdef BALANCE_OPERATOR
556	!
557	IF (Master) WRITE (stdout,20) 'BALANCE_OPERATOR', &
558	& 'Error Covariance Multivariate Balance Operator'
559	is=LEN_TRIM(Coptions)+1
560	Coptions(is:is+18)=' BALANCE_OPERATOR,'
561	#endif
562	#if defined SEDIMENT && defined BEDLOAD_MPM
563	!
564	IF (Master) WRITE (stdout,20) 'BEDLOAD_MPM', &
565	& 'Activate bed load sediment transport Meyer-Peter-Mueller'
566	is=LEN_TRIM(Coptions)+1
567	Coptions(is:is+13)=' BEDLOAD_MPM,'
568	#endif
569	#if defined SEDIMENT && defined BEDLOAD_SOULSBY
570	!
571	IF (Master) WRITE (stdout,20) 'BEDLOAD_SOULSBY', &
572	& 'Activate bed load sediment transport Soulsby formula'
573	is=LEN_TRIM(Coptions)+1
574	Coptions(is:is+17)=' BEDLOAD_SOULSBY,'
575	#endif
576	#if defined BGQC && defined FOUR_DVAR
577	!
578	IF (Master) WRITE (stdout,20) 'BGQC', &
579	& 'Background quality control of observations'
580	is=LEN_TRIM(Coptions)+1
581	Coptions(is:is+6)=' BGQC,'
582	#endif
583	#if defined BEOFS_ONLY && defined WEAK_CONSTRAINT
584	!
585	IF (Master) WRITE (stdout,20) 'BEOFS_ONLY', &
586	& 'Compute EOFs of the background error covariance matrix'
587	is=LEN_TRIM(Coptions)+1
588	Coptions(is:is+12)=' BEOFS_ONLY,'
589	#endif
590	#ifdef BIO_FENNEL
591	!
592	IF (Master) WRITE (stdout,20) 'BIO_FENNEL', &
593	& 'Fennel et al. (2006) nitrogen-based model'
594	is=LEN_TRIM(Coptions)+1
595	Coptions(is:is+12)=' BIO_FENNEL,'
596	ibiology=ibiology+1
597	#endif
598	#if defined BNORM && defined HESSIAN_SV
599	!
600	IF (Master) WRITE (stdout,20) 'BNORM', &
601	& 'Background normalization of Hessian singular vectors'
602	is=LEN_TRIM(Coptions)+1
603	Coptions(is:is+7)=' BNORM,'
604	#endif
605	#if defined BIO_SEDIMENT && \
606	(defined BIO_FENNEL \|\| defined NEMURO \|\| defined ECOSIM)
607	!
608	IF (Master) WRITE (stdout,20) 'BIO_SEDIMENT', &
609	& 'Restore fallen particulate material to the nutrient pool'
610	is=LEN_TRIM(Coptions)+1
611	Coptions(is:is+14)=' BIO_SEDIMENT,'
612	#endif
613	#ifdef BODYFORCE
614	!
615	IF (Master) WRITE (stdout,20) 'BODYFORCE', &
616	& 'Momentum stresses as body-forces'
617	is=LEN_TRIM(Coptions)+1
618	Coptions(is:is+11)=' BODYFORCE,'
619	#endif
620	#ifdef DISTRIBUTE
621	# if defined BOUNDARY_ALLGATHER
622	!
623	IF (Master) WRITE (stdout,20) 'BOUNDARY_ALLGATHER', &
624	& 'Using mpi_allgather in mp_boundary routine'
625	is=LEN_TRIM(Coptions)+1
626	Coptions(is:is+20)=' BOUNDARY_ALLGATHER,'
627	# else
628	IF (Master) WRITE (stdout,20) '!BOUNDARY_ALLGATHER', &
629	& 'Using mpi_allreduce in mp_boundary routine'
630	is=LEN_TRIM(Coptions)+1
631	Coptions(is:is+11)=' !BOUNDARY_ALLGATHER,'
632	# endif
633	#endif
634	#ifdef BULK_FLUXES
635	!
636	IF (Master) WRITE (stdout,20) 'BULK_FLUXES', &
637	& 'Surface bulk fluxes parameterization'
638	is=LEN_TRIM(Coptions)+1
639	Coptions(is:is+13)=' BULK_FLUXES,'
640	#endif
641	#ifdef BVF_MIXING
642	!
643	IF (Master) WRITE (stdout,20) 'BVF_MIXING', &
644	& 'Brunt-Vaisala frequency based vertical mixing'
645	is=LEN_TRIM(Coptions)+1
646	Coptions(is:is+12)=' BVF_MIXING,'
647	ivmix=ivmix+1
648	#endif
649	#if defined CANUTO_A && defined GLS_MIXING
650	!
651	IF (Master) WRITE (stdout,20) 'CANUTO_A', &
652	& 'Canuto A-stability function formulation'
653	is=LEN_TRIM(Coptions)+1
654	Coptions(is:is+8)=' CANUTO_A,'
655	#endif
656	#if defined CANUTO_B && defined GLS_MIXING
657	!
658	IF (Master) WRITE (stdout,20) 'CANUTO_B', &
659	& 'Canuto B-stability function formulation'
660	is=LEN_TRIM(Coptions)+1
661	Coptions(is:is+8)=' CANUTO_B,'
662	#endif
663	#if defined CARBON && defined BIO_FENNEL
664	!
665	IF (Master) WRITE (stdout,20) 'CARBON', &
666	& 'Add Carbon constituents'
667	is=LEN_TRIM(Coptions)+1
668	Coptions(is:is+8)=' CARBON,'
669	#endif
670	#if defined CELERITY_READ && defined FORWARD_READ
671	!
672	IF (Master) WRITE (stdout,20) 'CELERITY_READ', &
673	& 'Read in radiation celerity from forward file'
674	is=LEN_TRIM(Coptions)+1
675	Coptions(is:is+15)=' CELERETY_READ,'
676	#endif
677	#if defined CELERITY_WRITE && defined FORWARD_WRITE
678	!
679	IF (Master) WRITE (stdout,20) 'CELERITY_WRITE', &
680	& 'Write out radiation celerity in forward file'
681	is=LEN_TRIM(Coptions)+1
682	Coptions(is:is+16)=' CELERITY_WRITE,'
683	#endif
684	#ifdef CLIPPING
685	!
686	IF (Master) WRITE (stdout,20) 'CLIPPING', &
687	& 'Stabilized Representer Matrix clipping spectrum analysis'
688	is=LEN_TRIM(Coptions)+1
689	Coptions(is:is+10)=' CLIPPING,'
690	idriver=idriver+1
691	# ifdef CLIPPING_SPLIT
692	!
693	IF (Master) WRITE (stdout,20) 'CLIPPING_SPLIT', &
694	& 'Separate analysis due to IC, surface forcing, and OBC'
695	is=LEN_TRIM(Coptions)+1
696	Coptions(is:is+16)=' CLIPPING_SPLIT,'
697	# endif
698	#endif
699	#ifdef COAMPS_COUPLING
700	!
701	IF (Master) WRITE (stdout,20) 'COAMPS_COUPLING', &
702	& 'Atmosphere coupling with COAMPS'
703	is=LEN_TRIM(Coptions)+1
704	Coptions(is:is+17)=' COAMPS_COUPLING,'
705	iatms=iatms+1
706	#endif
707	#if defined DISTRIBUTE
708	# if defined COLLECT_ALLGATHER
709	!
710	IF (Master) WRITE (stdout,20) 'COLLECT_ALLGATHER', &
711	& 'Using mpi_allgather in mp_collect routine'
712	is=LEN_TRIM(Coptions)+1
713	Coptions(is:is+19)=' COLLECT_ALLGATHER,'
714	# elif defined COLLECT_ALLREDUCE
715	!
716	IF (Master) WRITE (stdout,20) 'COLLECT_ALLREDUCE', &
717	& 'Using mpi_allreduce in mp_collect routine'
718	is=LEN_TRIM(Coptions)+1
719	Coptions(is:is+19)=' COLLECT_ALLGATHER,'
720	# else
721	!
722	IF (Master) WRITE (stdout,20) '!COLLECT_ALL...', &
723	& 'Using mpi_isend/mpi_recv in mp_collect routine'
724	is=LEN_TRIM(Coptions)+1
725	Coptions(is:is+17)=' !COLLECT_ALL...,'
726	# endif
727	#endif
728	#ifdef CORRELATION
729	!
730	IF (Master) WRITE (stdout,20) 'CORRELATION', &
731	& 'Background-error Correlation Model'
732	is=LEN_TRIM(Coptions)+1
733	Coptions(is:is+13)=' CORRELATION,'
734	idriver=idriver+1
735	#endif
736	#ifdef DATA_COUPLING
737	!
738	IF (Master) WRITE (stdout,20) 'DATA_COUPLING', &
739	& 'DATA model included in coupling'
740	is=LEN_TRIM(Coptions)+1
741	Coptions(is:is+15)=' DATA_COUPLING,'
742	#endif
743	#if defined DATALESS_LOOPS && defined W4DVAR
744	!
745	IF (Master) WRITE (stdout,20) 'DATALESS_LOOPS', &
746	& 'Testing convergence of RPM Picard iterations'
747	is=LEN_TRIM(Coptions)+1
748	Coptions(is:is+16)=' DATALESS_LOOPS,'
749	#endif
750	#if defined GLS_MIXING && defined CHARNOK
751	!
752	IF (Master) WRITE (stdout,20) 'CHARNOK', &
753	& 'Charnok surface roughness from wind stress'
754	is=LEN_TRIM(Coptions)+1
755	Coptions(is:is+9)=' CHARNOK,'
756	#endif
757	#if defined PROPAGATOR && defined CHECKPOINTING
758	!
759	IF (Master) WRITE (stdout,20) 'CHECKPOINTING', &
760	& 'Processing checkpointing NetCDF for GST analysis'
761	is=LEN_TRIM(Coptions)+1
762	Coptions(is:is+15)=' CHECKPOINTING,'
763	#endif
764	#if defined GLS_MIXING && defined CRAIG_BANNER
765	!
766	IF (Master) WRITE (stdout,20) 'CRAIG_BANNER', &
767	& 'Craig and Banner wave breaking surface flux'
768	is=LEN_TRIM(Coptions)+1
769	Coptions(is:is+14)=' CRAIG_BANNER,'
770	#endif
771	#if defined COARE_OOST
772	!
773	IF (Master) WRITE (stdout,20) 'COARE_OOST', &
774	& 'Oost et al (2002) relation for ZoW in bulk fluxes'
775	is=LEN_TRIM(Coptions)+1
776	Coptions(is:is+12)=' COARE_OOST,'
777	#endif
778	#if defined COARE_TAYLOR_YELLAND
779	!
780	IF (Master) WRITE (stdout,20) 'COARE_TAYLOR_YELLAND', &
781	& 'Taylor and Yelland (2001) relation for ZoW in bulk fluxes'
782	is=LEN_TRIM(Coptions)+1
783	Coptions(is:is+22)=' COARE_TAYLOR_YELLAND,'
784	#endif
785	#if defined COOL_SKIN && defined BULK_FLUXES
786	!
787	IF (Master) WRITE (stdout,20) 'COOL_SKIN', &
788	& 'Surface cool skin correction'
789	is=LEN_TRIM(Coptions)+1
790	Coptions(is:is+11)=' COOL_SKIN,'
791	#endif
792	#if defined COSINE2 && defined SOLVE3D
793	!
794	IF (Master) WRITE (stdout,20) 'COSINE2', &
795	& 'Cosine-squared shape time-averaging barotropic filter'
796	is=LEN_TRIM(Coptions)+1
797	Coptions(is:is+9)=' COSINE2,'
798	#endif
799	#ifdef CURVGRID
800	!
801	IF (Master) WRITE (stdout,20) 'CURVGRID', &
802	& 'Orthogonal curvilinear grid'
803	is=LEN_TRIM(Coptions)+1
804	Coptions(is:is+10)=' CURVGRID,'
805	#endif
806	#if defined DAILY_SHORTWAVE && defined RED_TIDE
807	!
808	IF (Master) WRITE (stdout,20) 'DAILY_SHORTWAVE', &
809	& 'Daily-averaged surface solar shortwave radiation flux'
810	is=LEN_TRIM(Coptions)+1
811	Coptions(is:is+17)=' DAILY_SHORTWAVE,'
812	#endif
813	#ifdef DEBUGGING
814	!
815	IF (Master) WRITE (stdout,20) 'DEBUGGING', &
816	& 'Internal debugging switch activated'
817	is=LEN_TRIM(Coptions)+1
818	Coptions(is:is+11)=' DEBUGGING,'
819	#endif
820	#ifdef DEEPWATER_WAVES
821	!
822	IF (Master) WRITE (stdout,20) 'DEEPWATER_WAVES', &
823	& 'Deep water waves approx in bulk fluxes'
824	is=LEN_TRIM(Coptions)+1
825	Coptions(is:is+17)=' DEEPWATER_WAVES,'
826	#endif
827	#if defined DEFLATE && defined HDF5
828	!
829	IF (Master) WRITE (stdout,20) 'DEFLATE', &
830	& 'Setting compression in output NetCDF-4/HDF5 files'
831	is=LEN_TRIM(Coptions)+1
832	Coptions(is:is+9)=' DEFLATE,'
833	#endif
834	#if defined DENITRIFICATION && defined BIO_FENNEL
835	!
836	IF (Master) WRITE (stdout,20) 'DENITRIFICATION', &
837	& 'Add denitrification processes'
838	is=LEN_TRIM(Coptions)+1
839	Coptions(is:is+17)=' DENITRIFICATION,'
840	#endif
841	#if defined DIAGNOSTICS_BIO && defined BIO_FENNEL
842	!
843	IF (Master) WRITE (stdout,20) 'DIAGNOSTICS_BIO', &
844	& 'Computing and writing biological diagnostic terms'
845	is=LEN_TRIM(Coptions)+1
846	Coptions(is:is+17)=' DIAGNOSTICS_BIO,'
847	#endif
848	#ifdef DIAGNOSTICS_TS
849	!
850	IF (Master) WRITE (stdout,20) 'DIAGNOSTICS_TS', &
851	& 'Computing and writing tracer diagnostic terms'
852	is=LEN_TRIM(Coptions)+1
853	Coptions(is:is+16)=' DIAGNOSTICS_TS,'
854	#endif
855	#ifdef DIAGNOSTICS_UV
856	!
857	IF (Master) WRITE (stdout,20) 'DIAGNOSTICS_UV', &
858	& 'Computing and writing momentum diagnostic terms'
859	is=LEN_TRIM(Coptions)+1
860	Coptions(is:is+16)=' DIAGNOSTICS_UV,'
861	#endif
862	#if defined DIFF_3DCOEF
863	!
864	IF (Master) WRITE (stdout,20) 'DIFF_3DCOEF', &
865	& 'Horizontal, time-dependent 3D diffusion coefficient'
866	is=LEN_TRIM(Coptions)+1
867	Coptions(is:is+13)=' DIFF_3DCOEF,'
868	#endif
869	#if defined TS_DIF2 \|\| defined TS_DIF4
870	# ifdef DIFF_GRID
871	!
872	IF (Master) WRITE (stdout,20) 'DIFF_GRID', &
873	& 'Horizontal diffusion coefficient scaled by grid size'
874	is=LEN_TRIM(Coptions)+1
875	Coptions(is:is+11)=' DIFF_GRID,'
876	# endif
877	#endif
878	#ifdef DIURNAL_SRFLUX
879	!
880	IF (Master) WRITE (stdout,20) 'DIURNAL_SRFLUX', &
881	& 'Modulate shortwave radiation by the local diurnal cycle'
882	is=LEN_TRIM(Coptions)+1
883	Coptions(is:is+16)=' DIURNAL_SRFLUX,'
884	#endif
885	#ifdef DJ_GRADPS
886	!
887	IF (Master) WRITE (stdout,20) 'DJ_GRADPS', &
888	& 'Parabolic Splines density Jacobian (Shchepetkin, 2002)'
889	is=LEN_TRIM(Coptions)+1
890	Coptions(is:is+11)=' DJ_GRADPS,'
891	#endif
892	#ifdef DOUBLE_PRECISION
893	!
894	IF (Master) WRITE (stdout,20) 'DOUBLE_PRECISION', &
895	& 'Double precision arithmetic numerical kernel.'
896	is=LEN_TRIM(Coptions)+1
897	Coptions(is:is+18)=' DOUBLE_PRECISION,'
898	#endif
899	#ifdef ECOSIM
900	!
901	IF (Master) WRITE (stdout,20) 'ECOSIM', &
902	& 'Bio-Optical EcoSim model'
903	is=LEN_TRIM(Coptions)+1
904	Coptions(is:is+8)=' ECOSIM,'
905	ibiology=ibiology+1
906	#endif
907	#if defined EMINUSP && defined BULK_FLUXES
908	!
909	IF (Master) WRITE (stdout,20) 'EMINUSP', &
910	& 'Compute Salt Flux using E-P'
911	is=LEN_TRIM(Coptions)+1
912	Coptions(is:is+9)=' EMINUSP,'
913	#endif
914	#if defined ESMF_LIB && defined MODEL_COUPLING
915	!
916	IF (Master) WRITE (stdout,20) 'ESMF_LIB', &
917	& 'Using Earth System Modeling Framework library'
918	is=LEN_TRIM(Coptions)+1
919	Coptions(is:is+10)=' ESMF_LIB,'
920	#endif
921	#ifdef ENKF_RESTART
922	!
923	IF (Master) WRITE (stdout,20) 'ENKF_RESTART', &
924	& 'Writting fields for offline Ensemble Kalman Filter'
925	is=LEN_TRIM(Coptions)+1
926	Coptions(is:is+14)=' ENKF_RESTART,'
927	#endif
928	#ifdef ENSEMBLE
929	!
930	IF (Master) WRITE (stdout,20) 'ENSEMBLE', &
931	& 'Ensemble Forecasting Propagator'
932	is=LEN_TRIM(Coptions)+1
933	Coptions(is:is+10)=' ENSEMBLE,'
934	#endif
935	#if defined EVOLVED_LCZ && defined FOUR_DVAR
936	!
937	IF (Master) WRITE (stdout,20) 'EVOLVED_LCZ', &
938	& 'Compute 4D-Var evolved Hessian singular vectors'
939	is=LEN_TRIM(Coptions)+1
940	Coptions(is:is+13)=' EVOLVED_LCZ,'
941	#endif
942	# if defined EXCLUDE_SPONGE && defined MODEL_COUPLING && \
943	(defined DATA_COUPLING && !defined ANA_SPONGE)
944	!
945	IF (Master) WRITE (stdout,20) 'EXCLUDE_SPONGE', &
946	& 'Exclude sponge points in export fields during coupling'
947	is=LEN_TRIM(Coptions)+1
948	Coptions(is:is+16)=' EXCLUDE_SPONGE,'
949	#endif
950	#ifdef FLOATS
951	!
952	IF (Master) WRITE (stdout,20) 'FLOATS', &
953	& 'Simulated Lagrangian drifters'
954	is=LEN_TRIM(Coptions)+1
955	Coptions(is:is+8)=' FLOATS,'
956	#endif
957	#if defined FLOAT_BIOLOGY && defined FLOATS
958	!
959	IF (Master) WRITE (stdout,20) 'FLOAT_BIOLOGY', &
960	& 'Biological behavior on Lagrangian drifters'
961	is=LEN_TRIM(Coptions)+1
962	Coptions(is:is+15)=' FLOAT_BIOLOGY,'
963	#endif
964	#if defined FLOAT_OYSTER && defined FLOATS
965	!
966	IF (Master) WRITE (stdout,20) 'FLOAT_OYSTER', &
967	& 'Oyster model behavior on Lagrangian drifters'
968	is=LEN_TRIM(Coptions)+1
969	Coptions(is:is+14)=' FLOAT_OYSTER,'
970	#endif
971	#if defined FLOAT_STICKY && defined FLOATS && defined SOLVE3D
972	!
973	IF (Master) WRITE (stdout,20) 'FLOAT_STICKY', &
974	& 'Reflect (stick) floats that hit the surface (bottom)'
975	is=LEN_TRIM(Coptions)+1
976	Coptions(is:is+14)=' FLOAT_STICKY,'
977	#endif
978	#if defined FLOATS && defined FLOAT_VWALK && defined SOLVE3D
979	!
980	IF (Master) WRITE (stdout,20) 'FLOAT_VWALK', &
981	& 'Lagrangian drifters with vertical diffusion, random walk'
982	is=LEN_TRIM(Coptions)+1
983	Coptions(is:is+13)=' FLOAT_VWALK,'
984	#endif
985	#ifdef FORCING_SV
986	!
987	IF (Master) WRITE (stdout,20) 'FORCING_SV', &
988	& 'Forcing Singular Vectors Propagator'
989	is=LEN_TRIM(Coptions)+1
990	Coptions(is:is+12)=' FORCING_SV,'
991	idriver=idriver+1
992	#endif
993	#if defined FORWARD_MIXING && defined SOLVE3D
994	!
995	IF (Master) WRITE (stdout,20) 'FORWARD_MIXING', &
996	& 'Read in Forward vertical mixing for Tangent/Adjoint'
997	is=LEN_TRIM(Coptions)+1
998	Coptions(is:is+16)=' FORWARD_MIXING,'
999	#endif
1000	#ifdef FORWARD_READ
1001	!
1002	IF (Master) WRITE (stdout,20) 'FORWARD_READ', &
1003	& 'Read in Forward solution for Tangent/Adjoint'
1004	is=LEN_TRIM(Coptions)+1
1005	Coptions(is:is+14)=' FORWARD_READ,'
1006	#endif
1007	#ifdef FORWARD_RHS
1008	!
1009	IF (Master) WRITE (stdout,20) 'FORWARD_RHS', &
1010	& 'Process Forward RHS terms for Tangent/Adjoint'
1011	is=LEN_TRIM(Coptions)+1
1012	Coptions(is:is+13)=' FORWARD_RHS,'
1013	#endif
1014	#ifdef FORWARD_WRITE
1015	!
1016	IF (Master) WRITE (stdout,20) 'FORWARD_WRITE', &
1017	& 'Write out Forward solution for Tangent/Adjoint'
1018	is=LEN_TRIM(Coptions)+1
1019	Coptions(is:is+15)=' FORWARD_WRITE,'
1020	#endif
1021	#if defined FRC_COUPLING && defined MODEL_COUPLING
1022	!
1023	IF (Master) WRITE (stdout,20) 'FRC_COUPLING', &
1024	& 'Atmospheric forcing fields are from ESM coupling'
1025	is=LEN_TRIM(Coptions)+1
1026	Coptions(is:is+14)=' FRC_COUPLING,'
1027	#endif
1028	#ifdef FSOBC_REDUCED
1029	!
1030	IF (Master) WRITE (stdout,20) 'FSOBC_REDUCED', &
1031	& 'Using free-surface data in reduced physics conditions'
1032	is=LEN_TRIM(Coptions)+1
1033	Coptions(is:is+15)=' FSOBC_REDUCED,'
1034	#endif
1035	#ifdef FT_EIGENMODES
1036	!
1037	IF (Master) WRITE (stdout,20) 'FT_EIGENMODES', &
1038	& 'Finite Time Eigenmodes: Normal Modes'
1039	is=LEN_TRIM(Coptions)+1
1040	Coptions(is:is+15)=' FT_EIGENMODES,'
1041	idriver=idriver+1
1042	#endif
1043	#if defined FOUR_DVAR \|\| defined PROPAGATOR \|\| defined VERIFICATION
1044	# ifdef FULL_GRID
1045	!
1046	IF (Master) WRITE (stdout,20) 'FULL_GRID', &
1047	# if (defined FOUR_DVAR \|\| defined VERIFICATION)
1048	& 'Considering observations at interior and boundary points'
1049	# elif defined PROPAGATOR
1050	& 'State vector includes interior and boundary points'
1051	# endif
1052	is=LEN_TRIM(Coptions)+1
1053	Coptions(is:is+11)=' FULL_GRID,'
1054	# else
1055	!
1056	IF (Master) WRITE (stdout,20) '!FULL_GRID', &
1057	# if (defined FOUR_DVAR \|\| defined VERIFICATION)
1058	& 'Considering observations at interior points only'
1059	# elif defined PROPAGATOR
1060	& 'State vector includes interior points only'
1061	# endif
1062	is=LEN_TRIM(Coptions)+1
1063	Coptions(is:is+12)=' !FULL_GRID,'
1064	# endif
1065	#endif
1066	#if defined GEOPOTENTIAL_HCONV && defined FOUR_DVAR && defined SOLVE3D
1067	!
1068	IF (Master) WRITE (stdout,20) 'GEOPOTENTIAL_HCONV', &
1069	& 'Horizontal convolutions along geopetential surfaces'
1070	is=LEN_TRIM(Coptions)+1
1071	Coptions(is:is+20)=' GEOPOTENTIAL_HCONV,'
1072	#endif
1073	#ifdef GLS_MIXING
1074	!
1075	IF (Master) WRITE (stdout,20) 'GLS_MIXING', &
1076	& 'Generic Length-Scale turbulence closure'
1077	is=LEN_TRIM(Coptions)+1
1078	Coptions(is:is+12)=' GLS_MIXING,'
1079	ivmix=ivmix+1
1080	#endif
1081	#ifdef HDF5
1082	!
1083	IF (Master) WRITE (stdout,20) 'HDF5', &
1084	& 'Creating NetCDF-4/HDF5 format files'
1085	is=LEN_TRIM(Coptions)+1
1086	Coptions(is:is+6)=' HDF5,'
1087	#endif
1088	#ifdef HESSIAN_FSV
1089	!
1090	IF (Master) WRITE (stdout,20) 'HESSIAN_FSV', &
1091	& 'Hessian (4D-Var) Forcing Singular Vectors Propagator'
1092	is=LEN_TRIM(Coptions)+1
1093	Coptions(is:is+13)=' HESSIAN_FSV,'
1094	#endif
1095	#ifdef HESSIAN_SO
1096	!
1097	IF (Master) WRITE (stdout,20) 'HESSIAN_SO', &
1098	& 'Hessian (4D-Var) Stochastic Optimals Propagator'
1099	is=LEN_TRIM(Coptions)+1
1100	Coptions(is:is+12)=' HESSIAN_SO,'
1101	#endif
1102	#ifdef HESSIAN_SV
1103	!
1104	IF (Master) WRITE (stdout,20) 'HESSIAN_SV', &
1105	& 'Hessian (4D-Var) Singular Vectors Propagator'
1106	is=LEN_TRIM(Coptions)+1
1107	Coptions(is:is+12)=' HESSIAN_SV,'
1108	#endif
1109	#if defined HOLLING_GRAZING && defined NEMURO
1110	!
1111	IF (Master) WRITE (stdout,20) 'HOLLING_GRAZING', &
1112	& 'Holling-type s-shaped grazing algorithm, Implicit'
1113	is=LEN_TRIM(Coptions)+1
1114	Coptions(is:is+17)=' HOLLING_GRAZING,'
1115	#endif
1116	#ifdef HYPOXIA_SRM
1117	!
1118	IF (Master) WRITE (stdout,20) 'HYPOXIA_SRM', &
1119	& 'Hypoxia Simple Respiration Model'
1120	is=LEN_TRIM(Coptions)+1
1121	Coptions(is:is+13)=' HYPOXIA_SRM,'
1122	ibiology=ibiology+1
1123	#endif
1124	#ifdef ICESHELF
1125	!
1126	IF (Master) WRITE (stdout,20) 'ICESHELF', &
1127	& 'Include Ice Shelf Cavities'
1128	is=LEN_TRIM(Coptions)+1
1129	Coptions(is:is+10)=' ICESHELF,'
1130	#endif
1131	#if defined IMPACT_INNER && defined OBS_IMPACT && \
1132	(defined W4DVAR_SENSITIVITY \|\| defined W4DPSAS_SENSITIVITY \|\| \
1133	defined IS4DVAR_SENSITIVITY)
1134	!
1135	IF (Master) WRITE (stdout,20) 'IMPACT_INNER', &
1136	& 'Writing observations impacts for each inner loop'
1137	is=LEN_TRIM(Coptions)+1
1138	Coptions(is:is+14)=' IMPACT_INNER,'
1139	#endif
1140	#ifdef IMPLICIT_NUDGING
1141	!
1142	IF (Master) WRITE (stdout,20) 'IMPLICIT_NUDGING', &
1143	& 'Implicit nudging term in momentum radiation conditions'
1144	is=LEN_TRIM(Coptions)+1
1145	Coptions(is:is+18)=' IMPLICIT_NUDGING,'
1146	#endif
1147	#if defined IMPLICIT_VCONV && defined FOUR_DVAR && defined VCONVOLUTION
1148	!
1149	IF (Master) WRITE (stdout,20) 'IMPLICIT_VCONV', &
1150	& 'Implicit Vertical Convolution Algorithm'
1151	is=LEN_TRIM(Coptions)+1
1152	Coptions(is:is+16)=' IMPLICIT_VCONV,'
1153	#endif
1154	#if defined LIMIT_BSTRESS && defined SOLVE3D && !defined BBL_MODEL
1155	!
1156	IF (Master) WRITE (stdout,20) 'LIMIT_BSTRESS', &
1157	& 'Limit bottom stress to maintain bottom velocity direction'
1158	is=LEN_TRIM(Coptions)+1
1159	Coptions(is:is+15)=' LIMIT_BSTRESS,'
1160	#endif
1161	#if defined LIMIT_STFLX_COOLING && defined SOLVE3D
1162	!
1163	IF (Master) WRITE (stdout,20) 'LIMIT_STFLX_COOLING', &
1164	& 'Suppress further cooling if SST is at freezing point'
1165	is=LEN_TRIM(Coptions)+1
1166	Coptions(is:is+21)=' LIMIT_STFLX_COOLING,'
1167	#endif
1168	#if defined NPZD_IRON && defined IRON_LIMIT
1169	!
1170	IF (Master) WRITE (stdout,20) 'IRON_LIMIT', &
1171	& 'Include dissolved iron limitation on phytoplankton growth'
1172	is=LEN_TRIM(Coptions)+1
1173	Coptions(is:is+12)=' IRON_LIMIT,'
1174	#endif
1175	#if defined NPZD_IRON && defined IRON_LIMIT && defined IRON_RELAX
1176	!
1177	IF (Master) WRITE (stdout,20) 'IRON_RELAX', &
1178	& 'Nudge dissolved iron over the shelf, h <= FeHmin'
1179	is=LEN_TRIM(Coptions)+1
1180	Coptions(is:is+12)=' IRON_RELAX,'
1181	#endif
1182	#ifdef IMPULSE
1183	!
1184	IF (Master) WRITE (stdout,20) 'IMPULSE', &
1185	& 'Processing Adjoint Impulse forcing'
1186	is=LEN_TRIM(Coptions)+1
1187	Coptions(is:is+9)=' IMPULSE,'
1188	#endif
1189	#ifdef INNER_PRODUCT
1190	!
1191	IF (Master) WRITE (stdout,20) 'INNER_PRODUCT', &
1192	& 'Tangent/Adjoint State Matrices Inner Product Test'
1193	is=LEN_TRIM(Coptions)+1
1194	Coptions(is:is+15)=' INNER_PRODUCT,'
1195	idriver=idriver+1
1196	#endif
1197	#ifdef IS4DVAR
1198	!
1199	IF (Master) WRITE (stdout,20) 'IS4DVAR', &
1200	& 'Incremental strong constraint 4D-Var data assimilation'
1201	is=LEN_TRIM(Coptions)+1
1202	Coptions(is:is+9)=' IS4DVAR,'
1203	idriver=idriver+1
1204	#endif
1205	#ifdef IS4DVAR_SENSITIVITY
1206	!
1207	IF (Master) WRITE (stdout,20) 'IS4DVAR_SENSITIVITY', &
1208	& 'I4D-Var Observations Sensitivity Analysis'
1209	is=LEN_TRIM(Coptions)+1
1210	Coptions(is:is+21)=' IS4DVAR_SENSITIVITY,'
1211	idriver=idriver+1
1212	#endif
1213	#if defined INLINE_2DIO && defined DISTRIBUTE
1214	!
1215	IF (Master) WRITE (stdout,20) 'INLINE_2DIO', &
1216	& 'Processing 3D IO level by level to reduce memory needs'
1217	is=LEN_TRIM(Coptions)+1
1218	Coptions(is:is+13)=' INLINE_2DIO,'
1219	#endif
1220	#if defined IVLEV_EXPLICIT && defined NEMURO
1221	!
1222	IF (Master) WRITE (stdout,20) 'IVLEV_EXPLICIT', &
1223	& 'Ivlev Grazing algorithm, Explicit'
1224	is=LEN_TRIM(Coptions)+1
1225	Coptions(is:is+16)=' IVLEV_EXPLICIT,'
1226	#endif
1227	#if defined KANTHA_CLAYSON && (defined GLS_MIXING \|\| defined MY25_MIXING)
1228	!
1229	IF (Master) WRITE (stdout,20) 'KANTHA_CLAYSON', &
1230	& 'Kantha and Clayson stability function formulation'
1231	is=LEN_TRIM(Coptions)+1
1232	Coptions(is:is+16)=' KANTHA_CLAYSON,'
1233	#endif
1234	#if defined LCZ_FINAL && defined FOUR_DVAR
1235	!
1236	IF (Master) WRITE (stdout,20) 'LCZ_FINAL', &
1237	& 'Compute 4D-Var Hessian singular vectors'
1238	is=LEN_TRIM(Coptions)+1
1239	Coptions(is:is+11)=' LCZ_FINAL,'
1240	#endif
1241	#if defined LIMIT_VDIFF && \
1242	(defined GLS_MIXING \|\| defined LMD_MIXING \|\| defined MY25_MIXING)
1243	!
1244	IF (Master) WRITE (stdout,20) 'LIMIT_VDIFF', &
1245	& 'Impose an upper limit on vertical diffusion coefficient'
1246	is=LEN_TRIM(Coptions)+1
1247	Coptions(is:is+13)=' LIMIT_VDIFF,'
1248	#endif
1249	#if defined LIMIT_VVISC && \
1250	(defined GLS_MIXING \|\| defined LMD_MIXING \|\| defined MY25_MIXING)
1251	!
1252	IF (Master) WRITE (stdout,20) 'LIMIT_VVISC', &
1253	& 'Impose an upper limit on vertical viscosity coefficient'
1254	is=LEN_TRIM(Coptions)+1
1255	Coptions(is:is+13)=' LIMIT_VVISC,'
1256	#endif
1257	#ifdef LMD_BKPP
1258	!
1259	IF (Master) WRITE (stdout,20) 'LMD_BKPP', &
1260	& 'KPP bottom boundary layer mixing'
1261	is=LEN_TRIM(Coptions)+1
1262	Coptions(is:is+10)=' LMD_BKPP,'
1263	#endif
1264	#ifdef LMD_CONVEC
1265	!
1266	IF (Master) WRITE (stdout,20) 'LMD_CONVEC', &
1267	& 'LMD convective mixing due to shear instability'
1268	is=LEN_TRIM(Coptions)+1
1269	Coptions(is:is+12)=' LMD_CONVEC,'
1270	#endif
1271	#ifdef LMD_DDMIX
1272	!
1273	IF (Master) WRITE (stdout,20) 'LMD_DDMIX', &
1274	& 'LMD double-diffusive mixing'
1275	is=LEN_TRIM(Coptions)+1
1276	Coptions(is:is+11)=' LMD_DDMIX,'
1277	#endif
1278	#ifdef LMD_MIXING
1279	!
1280	IF (Master) WRITE (stdout,20) 'LMD_MIXING', &
1281	& 'Large/McWilliams/Doney interior mixing'
1282	is=LEN_TRIM(Coptions)+1
1283	Coptions(is:is+12)=' LMD_MIXING,'
1284	ivmix=ivmix+1
1285	#endif
1286	#ifdef LMD_NONLOCAL
1287	!
1288	IF (Master) WRITE (stdout,20) 'LMD_NONLOCAL', &
1289	& 'LMD convective nonlocal transport'
1290	is=LEN_TRIM(Coptions)+1
1291	Coptions(is:is+14)=' LMD_NONLOCAL,'
1292	#endif
1293	#ifdef LMD_RIMIX
1294	!
1295	IF (Master) WRITE (stdout,20) 'LMD_RIMIX', &
1296	& 'LMD diffusivity due to shear instability'
1297	is=LEN_TRIM(Coptions)+1
1298	Coptions(is:is+11)=' LMD_RIMIX,'
1299	#endif
1300	#ifdef LMD_SHAPIRO
1301	!
1302	IF (Master) WRITE (stdout,20) 'LMD_SHAPIRO', &
1303	& 'Shapiro filtering boundary layer depth'
1304	is=LEN_TRIM(Coptions)+1
1305	Coptions(is:is+13)=' LMD_SHAPIRO,'
1306	#endif
1307	#ifdef LMD_SKPP
1308	!
1309	IF (Master) WRITE (stdout,20) 'LMD_SKPP', &
1310	& 'KPP surface boundary layer mixing'
1311	is=LEN_TRIM(Coptions)+1
1312	Coptions(is:is+10)=' LMD_SKPP,'
1313	#endif
1314	#if defined LONGWAVE && defined BULK_FLUXES
1315	!
1316	IF (Master) WRITE (stdout,20) 'LONGWAVE', &
1317	& 'Compute net longwave radiation internally'
1318	is=LEN_TRIM(Coptions)+1
1319	Coptions(is:is+10)=' LONGWAVE,'
1320	#endif
1321	#if defined LONGWAVE_OUT && defined BULK_FLUXES
1322	!
1323	IF (Master) WRITE (stdout,20) 'LONGWAVE_OUT', &
1324	& 'Compute outgoing longwave radiation internally'
1325	is=LEN_TRIM(Coptions)+1
1326	Coptions(is:is+14)=' LONGWAVE_OUT,'
1327	#endif
1328	#ifdef MASKING
1329	!
1330	IF (Master) WRITE (stdout,20) 'MASKING', &
1331	& 'Land/Sea masking'
1332	is=LEN_TRIM(Coptions)+1
1333	Coptions(is:is+9)=' MASKING,'
1334	#endif
1335	#ifdef MB_BBL
1336	!
1337	IF (Master) WRITE (stdout,20) 'MB_BBL', &
1338	& 'Blaas Bottom Boundary Layer'
1339	is=LEN_TRIM(Coptions)+1
1340	Coptions(is:is+8)=' MB_BBL,'
1341	ibbl=ibbl+1
1342	# ifdef MB_CALC_UB
1343	!
1344	IF (Master) WRITE (stdout,20) 'MB_CALC_UB', &
1345	& 'Internal computation of bottom orbital velocity'
1346	is=LEN_TRIM(Coptions)+1
1347	Coptions(is:is+12)=' MB_CALC_UB,'
1348	# endif
1349	# ifdef MB_CALC_ZNOT
1350	!
1351	IF (Master) WRITE (stdout,20) 'MB_CALC_ZNOT', &
1352	& 'Internal computation of bottom roughness'
1353	is=LEN_TRIM(Coptions)+1
1354	Coptions(is:is+14)=' MB_CALC_ZNOT,'
1355	# endif
1356	# ifdef MB_Z0BIO
1357	!
1358	IF (Master) WRITE (stdout,20) 'MB_Z0BIO', &
1359	& 'Biogenic bedform roughness ripple height and length'
1360	is=LEN_TRIM(Coptions)+1
1361	Coptions(is:is+10)=' MB_Z0BIO,'
1362	# endif
1363	# ifdef MB_Z0BL
1364	!
1365	IF (Master) WRITE (stdout,20) 'MB_Z0BL', &
1366	& 'Bedload roughness for ripples'
1367	is=LEN_TRIM(Coptions)+1
1368	Coptions(is:is+9)=' MB_Z0BL,'
1369	# endif
1370	# ifdef MB_Z0RIP
1371	!
1372	IF (Master) WRITE (stdout,20) 'MB_Z0RIP', &
1373	& 'Bedform roughness ripple height and length'
1374	is=LEN_TRIM(Coptions)+1
1375	Coptions(is:is+10)=' MB_Z0RIP,'
1376	# endif
1377	#endif
1378	#if defined MCT_LIB && defined MODEL_COUPLING
1379	!
1380	IF (Master) WRITE (stdout,20) 'MCT_LIB', &
1381	& 'Using Model Coupling Toolkit library'
1382	is=LEN_TRIM(Coptions)+1
1383	Coptions(is:is+9)=' MCT_LIB,'
1384	#endif
1385	#if defined MINRES && defined FOUR_DVAR
1386	!
1387	IF (Master) WRITE (stdout,20) 'MINRES', &
1388	& 'Minimal Residual Method for 4D-Var minimization'
1389	is=LEN_TRIM(Coptions)+1
1390	Coptions(is:is+8)=' MINRES,'
1391	#endif
1392	#if (defined TS_DIF2 \|\| defined TS_DIF4) && defined SOLVE3D
1393	# ifdef MIX_GEO_TS
1394	!
1395	IF (Master) WRITE (stdout,20) 'MIX_GEO_TS', &
1396	& 'Mixing of tracers along geopotential surfaces'
1397	is=LEN_TRIM(Coptions)+1
1398	Coptions(is:is+12)=' MIX_GEO_TS,'
1399	# endif
1400	# ifdef MIX_ISO_TS
1401	!
1402	IF (Master) WRITE (stdout,20) 'MIX_ISO_TS', &
1403	& 'Mixing of tracers along isopycnal surfaces'
1404	is=LEN_TRIM(Coptions)+1
1405	Coptions(is:is+12)=' MIX_ISO_TS,'
1406	# endif
1407	# ifdef MIX_S_TS
1408	!
1409	IF (Master) WRITE (stdout,20) 'MIX_S_TS', &
1410	& 'Mixing of tracers along constant S-surfaces'
1411	is=LEN_TRIM(Coptions)+1
1412	Coptions(is:is+10)=' MIX_S_TS,'
1413	# endif
1414	#endif
1415	#if (defined UV_VIS2 \|\| defined UV_VIS4) && defined SOLVE3D
1416	# ifdef MIX_GEO_UV
1417	!
1418	IF (Master) WRITE (stdout,20) 'MIX_GEO_UV', &
1419	& 'Mixing of momentum along geopotential surfaces'
1420	is=LEN_TRIM(Coptions)+1
1421	Coptions(is:is+11)=' MIX_GEO_UV,'
1422	# endif
1423	# ifdef MIX_S_UV
1424	!
1425	IF (Master) WRITE (stdout,20) 'MIX_S_UV', &
1426	& 'Mixing of momentum along constant S-surfaces'
1427	is=LEN_TRIM(Coptions)+1
1428	Coptions(is:is+10)=' MIX_S_UV,'
1429	# endif
1430	#endif
1431	#ifdef MPI
1432	!
1433	IF (Master) WRITE (stdout,20) 'MPI', &
1434	& 'MPI distributed-memory configuration'
1435	is=LEN_TRIM(Coptions)+1
1436	Coptions(is:is+4)=' MPI,'
1437	#endif
1438	#if defined MULTIPLE_TLM && defined TANGENT
1439	!
1440	IF (Master) WRITE (stdout,20) 'MULTIPLE_TLM', &
1441	& 'Create multiple TLM history files'
1442	is=LEN_TRIM(Coptions)+1
1443	Coptions(is:is+14)=' MULTIPLE_TLM,'
1444	#endif
1445	#ifdef MY25_MIXING
1446	!
1447	IF (Master) WRITE (stdout,20) 'MY25_MIXING', &
1448	& 'Mellor/Yamada Level-2.5 mixing closure'
1449	is=LEN_TRIM(Coptions)+1
1450	Coptions(is:is+13)=' MY25_MIXING,'
1451	ivmix=ivmix+1
1452	#endif
1453	#ifdef NEARSHORE_MELLOR05
1454	!
1455	IF (Master) WRITE (stdout,20) 'NEARSHORE_MELLOR05', &
1456	& 'Nearshore Radiation Stress Terms (Mellor 2005)'
1457	is=LEN_TRIM(Coptions)+1
1458	Coptions(is:is+20)=' NEARSHORE_MELLOR05,'
1459	nearshore=nearshore+1
1460	#endif
1461	#ifdef NEARSHORE_MELLOR08
1462	!
1463	IF (Master) WRITE (stdout,20) 'NEARSHORE_MELLOR08', &
1464	& 'Nearshore Radiation Stress Terms (Mellor 2008)'
1465	is=LEN_TRIM(Coptions)+1
1466	Coptions(is:is+20)=' NEARSHORE_MELLOR08,'
1467	nearshore=nearshore+1
1468	#endif
1469	#ifdef NEMURO
1470	!
1471	IF (Master) WRITE (stdout,20) 'NEMURO', &
1472	& 'Nemuro Lower Trophic Level Ecosystem Model'
1473	is=LEN_TRIM(Coptions)+1
1474	Coptions(is:is+8)=' NEMURO,'
1475	ibiology=ibiology+1
1476	#endif
1477	#ifdef NESTING
1478	!
1479	IF (Master) WRITE (stdout,20) 'NESTING', &
1480	& 'Nesting grids: Composite and Refinement'
1481	is=LEN_TRIM(Coptions)+1
1482	Coptions(is:is+9)=' NESTING,'
1483	# ifdef NESTING_DEBUG
1484	!
1485	IF (Master) WRITE (stdout,20) 'NESTING_DEBUG', &
1486	& 'Checking mass fluxes conservation in refinement.'
1487	is=LEN_TRIM(Coptions)+1
1488	Coptions(is:is+15)=' NESTING_DEBUG,'
1489	# endif
1490	#endif
1491	#if defined NL_BULK_FLUXES
1492	!
1493	IF (Master) WRITE (stdout,20) 'NL_BULK_FLUXES', &
1494	& 'Using bulk fluxes computed by nonlinear model'
1495	is=LEN_TRIM(Coptions)+1
1496	Coptions(is:is+16)=' NL_BULK_FLUXES,'
1497	#endif
1498	#if defined NLM_OUTER && defined WEAK_CONSTRAINT
1499	!
1500	IF (Master) WRITE (stdout,20) 'NLM_OUTER', &
1501	& 'Using the Nonlinear model as Basic State in Outer Loop'
1502	is=LEN_TRIM(Coptions)+1
1503	Coptions(is:is+11)=' NLM_OUTER,'
1504	#endif
1505	#ifdef NONLINEAR
1506	!
1507	IF (Master) WRITE (stdout,20) 'NONLINEAR', &
1508	& 'Nonlinear Model'
1509	is=LEN_TRIM(Coptions)+1
1510	Coptions(is:is+10)=' NONLINEAR,'
1511	#endif
1512	#ifdef SOLVE3D
1513	# if defined NONLIN_EOS
1514	!
1515	IF (Master) WRITE (stdout,20) 'NONLIN_EOS', &
1516	& 'Nonlinear Equation of State for seawater'
1517	is=LEN_TRIM(Coptions)+1
1518	Coptions(is:is+12)=' NONLIN_EOS,'
1519	# else
1520	!
1521	IF (Master) WRITE (stdout,20) '!NONLIN_EOS', &
1522	& 'Linear Equation of State for seawater'
1523	is=LEN_TRIM(Coptions)+1
1524	Coptions(is:is+13)=' !NONLIN_EOS,'
1525	# endif
1526	#endif
1527	#if defined NO_CORRECT_TRACER && defined NESTING && defined SOLVE3D
1528	!
1529	IF (Master) WRITE (stdout,20) 'NO_CORRECT_TRACER', &
1530	& 'Avoid two-way correction of coaser grid boundary tracer'
1531	is=LEN_TRIM(Coptions)+1
1532	Coptions(is:is+19)=' NO_CORRECT_TRACER,'
1533	#endif
1534	#ifdef NO_LBC_ATT
1535	!
1536	IF (Master) WRITE (stdout,20) 'NO_LBC_ATT', &
1537	& 'Not checking NetCDF global attribute NLM_LBC during restart'
1538	is=LEN_TRIM(Coptions)+1
1539	Coptions(is:is+12)=' NO_LBC_ATT,'
1540	#endif
1541	#if defined NO_READ_GHOST && defined DISTRIBUTE
1542	!
1543	IF (Master) WRITE (stdout,20) 'NO_READ_GHOST', &
1544	& 'Not including ghost points during scattering read data'
1545	is=LEN_TRIM(Coptions)+1
1546	Coptions(is:is+15)=' NO_READ_GHOST,'
1547	#endif
1548	#ifdef NO_WRITE_GRID
1549	!
1550	IF (Master) WRITE (stdout,20) 'NO_WRITE_GRID', &
1551	& 'Not Writing grid arrays into NetCDF ouput files'
1552	is=LEN_TRIM(Coptions)+1
1553	Coptions(is:is+15)=' NO_WRITE_GRID,'
1554	#endif
1555	#ifdef NPZD_FRANKS
1556	!
1557	IF (Master) WRITE (stdout,20) 'NPZD_FRANKS', &
1558	& 'NPZD Biological Model, Franks et al. fomulation'
1559	is=LEN_TRIM(Coptions)+1
1560	Coptions(is:is+13)=' NPZD_FRANKS,'
1561	ibiology=ibiology+1
1562	#endif
1563	#ifdef NPZD_IRON
1564	!
1565	IF (Master) WRITE (stdout,20) 'NPZD_IRON', &
1566	& 'NPZD Biological Model (Powell et al.) with iron limitation'
1567	is=LEN_TRIM(Coptions)+1
1568	Coptions(is:is+11)=' NPZD_IRON,'
1569	ibiology=ibiology+1
1570	#endif
1571	#ifdef NPZD_POWELL
1572	!
1573	IF (Master) WRITE (stdout,20) 'NPZD_POWELL', &
1574	& 'NPZD Biological Model, Powell et al. fomulation'
1575	is=LEN_TRIM(Coptions)+1
1576	Coptions(is:is+13)=' NPZD_POWELL,'
1577	ibiology=ibiology+1
1578	#endif
1579	#if defined N2S2_HORAVG && (defined GLS_MIXING \|\| defined MY25_MIXING)
1580	!
1581	IF (Master) WRITE (stdout,20) 'N2S2_HORAVG', &
1582	& 'Horizontal smoothing of buoyancy and shear'
1583	is=LEN_TRIM(Coptions)+1
1584	Coptions(is:is+13)=' N2S2_HORAVG,'
1585	#endif
1586	#ifdef OBSERVATIONS
1587	!
1588	IF (Master) WRITE (stdout,20) 'OBSERVATIONS', &
1589	& 'Processing 4DVar observations'
1590	is=LEN_TRIM(Coptions)+1
1591	Coptions(is:is+13)=' OBSERVATIONS,'
1592	#endif
1593	#ifdef OBS_IMPACT
1594	!
1595	IF (Master) WRITE (stdout,20) 'OBS_IMPACT', &
1596	& 'Compute observation impact to the 4DVAR system'
1597	is=LEN_TRIM(Coptions)+1
1598	Coptions(is:is+12)=' OBS_IMPACT,'
1599	#endif
1600	#if defined OBS_IMPACT && defined OBS_IMPACT_SPLIT
1601	!
1602	IF (Master) WRITE (stdout,20) 'OBS_IMPACT_SPLIT', &
1603	& 'Separate impact due to IC, surface forcing, and OBC'
1604	is=LEN_TRIM(Coptions)+1
1605	Coptions(is:is+18)=' OBS_IMPACT_SPLIT,'
1606	#endif
1607	#ifdef NESTING
1608	# ifdef ONE_WAY
1609	!
1610	IF (Master) WRITE (stdout,20) 'ONE_WAY', &
1611	& 'One-way nesting in refinement grids'
1612	is=LEN_TRIM(Coptions)+1
1613	Coptions(is:is+9)=' ONE_WAY,'
1614	# else
1615	!
1616	IF (Master) WRITE (stdout,20) '!ONE_WAY', &
1617	& 'Two-way nesting in refinement grids'
1618	is=LEN_TRIM(Coptions)+1
1619	Coptions(is:is+10)=' !ONE_WAY,'
1620	# endif
1621	#endif
1622	#ifdef OUT_DOUBLE
1623	!
1624	IF (Master) WRITE (stdout,20) 'OUT_DOUBLE', &
1625	& 'Double precision output fields in NetCDF files'
1626	is=LEN_TRIM(Coptions)+1
1627	Coptions(is:is+12)=' OUT_DOUBLE,'
1628	#endif
1629	#ifdef _OPENMP
1630	!
1631	IF (Master) WRITE (stdout,20) '_OPENMP', &
1632	& 'OpenMP parallel shared-memory directives'
1633	is=LEN_TRIM(Coptions)+1
1634	Coptions(is:is+9)=' _OPENMP,'
1635	#endif
1636	#ifdef OPT_OBSERVATIONS
1637	!
1638	IF (Master) WRITE (stdout,20) 'OPT_OBSERVATIONS', &
1639	& 'Optimal Observations Driver'
1640	is=LEN_TRIM(Coptions)+1
1641	Coptions(is:is+18)=' OPT_OBSERVATIONS,'
1642	idriver=idriver+1
1643	#endif
1644	#ifdef OPT_PERTURBATION
1645	!
1646	IF (Master) WRITE (stdout,20) 'OPT_PERTURBATION', &
1647	& 'Optimal Perturbation Propagator'
1648	is=LEN_TRIM(Coptions)+1
1649	Coptions(is:is+18)=' OPT_PERTURBATION,'
1650	idriver=idriver+1
1651	#endif
1652	#if defined OXYGEN && defined BIO_FENNEL
1653	# ifdef OCMIP_OXYGEN_SC
1654	!
1655	IF (Master) WRITE (stdout,20) 'OCMIP_OXYGEN_SC', &
1656	& 'Schmidt number from Keeling et al. (1998)'
1657	is=LEN_TRIM(Coptions)+1
1658	Coptions(is:is+17)=' OCMIP_OXYGEN_SC,'
1659	# else
1660	!
1661	IF (Master) WRITE (stdout,20) '!OCMIP_OXYGEN_SC', &
1662	& 'Schmidt number from Wanninkhof (1992)'
1663	is=LEN_TRIM(Coptions)+1
1664	Coptions(is:is+18)=' !OCMIP_OXYGEN_SC,'
1665	# endif
1666	!
1667	IF (Master) WRITE (stdout,20) 'OXYGEN', &
1668	& 'Add oxygen dynamics'
1669	is=LEN_TRIM(Coptions)+1
1670	Coptions(is:is+8)=' OXYGEN,'
1671	#elif !defined SURFACE_DO_SATURATION && defined HYPOXIA_SRM
1672	# ifdef OCMIP_OXYGEN_SC
1673	!
1674	IF (Master) WRITE (stdout,20) 'OCMIP_OXYGEN_SC', &
1675	& 'Schmidt number from Keeling et al. (1998)'
1676	is=LEN_TRIM(Coptions)+1
1677	Coptions(is:is+17)=' OCMIP_OXYGEN_SC,'
1678	# else
1679	!
1680	IF (Master) WRITE (stdout,20) '!OCMIP_OXYGEN_SC', &
1681	& 'Schmidt number from Wanninkhof (1992)'
1682	is=LEN_TRIM(Coptions)+1
1683	Coptions(is:is+18)=' !OCMIP_OXYGEN_SC,'
1684	# endif
1685	#endif
1686	#if defined PARALLEL_IO && defined DISTRIBUTE
1687	!
1688	IF (Master) WRITE (stdout,20) 'PARALLEL_IO', &
1689	& 'Parallel I/O processing'
1690	is=LEN_TRIM(Coptions)+1
1691	Coptions(is:is+13)=' PARALLEL_IO,'
1692	#endif
1693	#ifdef PERFECT_RESTART
1694	!
1695	IF (Master) WRITE (stdout,20) 'PERFECT_RESTART', &
1696	& 'Processing perfect restart variables'
1697	is=LEN_TRIM(Coptions)+1
1698	Coptions(is:is+17)=' PERFECT_RESTART,'
1699	#endif
1700	#ifdef PICARD_TEST
1701	!
1702	IF (Master) WRITE (stdout,20) 'PICARD_TEST', &
1703	& 'Representers tangent linear model Picard iterations test'
1704	is=LEN_TRIM(Coptions)+1
1705	Coptions(is:is+13)=' PICARD TEST,'
1706	idriver=idriver+1
1707	#endif
1708	#ifdef PJ_GRADP
1709	!
1710	IF (Master) WRITE (stdout,20) 'PJ_GRADP', &
1711	& 'Finite volume Pressure Jacobian (Lin, 1997)'
1712	is=LEN_TRIM(Coptions)+1
1713	Coptions(is:is+10)=' PJ_GRADP,'
1714	#endif
1715	#ifdef PJ_GRADPQ2
1716	!
1717	IF (Master) WRITE (stdout,20) 'PJ_GRADPQ2', &
1718	& 'Quartic 2 polynomial Pressure Jacobian (Shchepetkin, 2002)'
1719	is=LEN_TRIM(Coptions)+1
1720	Coptions(is:is+12)=' PJ_GRADPQ2,'
1721	#endif
1722	#ifdef PJ_GRADPQ4
1723	!
1724	IF (Master) WRITE (stdout,20) 'PJ_GRADPQ4', &
1725	& 'Quartic 4 polynomial Pressure Jacobian (Shchepetkin, 2002)'
1726	is=LEN_TRIM(Coptions)+1
1727	Coptions(is:is+12)=' PJ_GRADPQ4,'
1728	#endif
1729	#ifdef POSITIVE_ZERO
1730	!
1731	IF (Master) WRITE (stdout,20) 'POSITIVE_ZERO', &
1732	& 'Impose positive zero in output data'
1733	is=LEN_TRIM(Coptions)+1
1734	Coptions(is:is+15)=' POSITIVE_ZERO,'
1735	#endif
1736	#if defined POSTERIOR_EOFS && defined WEAK_CONSTRAINT
1737	!
1738	IF (Master) WRITE (stdout,20) 'POSTERIOR_EOFS', &
1739	& 'Estimating posterior analysis error covariace matrix EOFs'
1740	is=LEN_TRIM(Coptions)+1
1741	Coptions(is:is+16)=' POSTERIOR_EOFS,'
1742	#endif
1743	#if defined POSTERIOR_ERROR_F && defined WEAK_CONSTRAINT
1744	!
1745	IF (Master) WRITE (stdout,20) 'POSTERIOR_ERROR_F', &
1746	& 'Estimating final posterior analysis error covariace matrix'
1747	is=LEN_TRIM(Coptions)+1
1748	Coptions(is:is+16)=' POSTERIOR_ERROR_F,'
1749	#endif
1750	#if defined POSTERIOR_ERROR_I && defined WEAK_CONSTRAINT
1751	!
1752	IF (Master) WRITE (stdout,20) 'POSTERIOR_ERROR_I', &
1753	& 'Estimating initial posterior analysis error covariace matrix'
1754	is=LEN_TRIM(Coptions)+1
1755	Coptions(is:is+16)=' POSTERIOR_ERROR_F,'
1756	#endif
1757	#if defined POWER_LAW && defined SOLVE3D
1758	!
1759	IF (Master) WRITE (stdout,20) 'POWER_LAW', &
1760	& 'Power-law shape time-averaging barotropic filter'
1761	is=LEN_TRIM(Coptions)+1
1762	Coptions(is:is+11)=' POWER_LAW,'
1763	#endif
1764	#if !(defined PJ_GRADPQ4 \|\| defined PJ_GRADPQ2 \|\| defined PJ_GRADP \|\| \
1765	defined DJ_GRADPS) && defined SOLVE3D
1766	!
1767	IF (Master) WRITE (stdout,20) 'PRSGRD31', &
1768	& 'Standard density Jacobian formulation (Song, 1998)'
1769	is=LEN_TRIM(Coptions)+1
1770	Coptions(is:is+10)=' PRSGRD31,'
1771	#endif
1772	#ifdef PROFILE
1773	!
1774	IF (Master) WRITE (stdout,20) 'PROFILE', &
1775	& 'Time profiling activated'
1776	is=LEN_TRIM(Coptions)+1
1777	Coptions(is:is+9)=' PROFILE,'
1778	#endif
1779	#ifdef PSEUDOSPECTRA
1780	!
1781	IF (Master) WRITE (stdout,20) 'PSEUDOSPECTRA', &
1782	& 'Pseudospectra Propagator'
1783	is=LEN_TRIM(Coptions)+1
1784	Coptions(is:is+15)=' PSEUDOSPECTRA,'
1785	#endif
1786	#ifdef QCORRECTION
1787	!
1788	IF (Master) WRITE (stdout,20) 'QCORRECTION', &
1789	& 'Surface net heat flux correction'
1790	is=LEN_TRIM(Coptions)+1
1791	Coptions(is:is+13)=' QCORRECTION,'
1792	#endif
1793	#ifdef OBSERVATIONS
1794	# ifdef RADIAL_ANGLE_CCW_EAST
1795	!
1796	IF (Master) WRITE (stdout,20) 'RADIAL_ANGLE_CCW_EAST', &
1797	& 'Radial observations as an azimuth counterclockwise from EAST'
1798	is=LEN_TRIM(Coptions)+1
1799	Coptions(is:is+23)=' RADIAL_ANGLE_CCW_EAST,'
1800	# else
1801	IF (Master) WRITE (stdout,20) '!RADIAL_ANGLE_CCW_EAST', &
1802	& 'Radial observations as an azimuth clockwise from NORTH'
1803	is=LEN_TRIM(Coptions)+1
1804	Coptions(is:is+24)=' !RADIAL_ANGLE_CCW_EAST,'
1805	# endif
1806	#endif
1807	#if defined GLS_MIXING \|\| defined MY25_MIXING
1808	# if defined K_C2ADVECTION
1809	!
1810	IF (Master) WRITE (stdout,20) 'K_C2ADVECTION', &
1811	& 'Second-order centered differences advection of TKE fields'
1812	is=LEN_TRIM(Coptions)+1
1813	Coptions(is:is+15)=' K_C2ADVECTION,'
1814	# elif defined K_C4ADVECTION
1815	!
1816	IF (Master) WRITE (stdout,20) 'K_C4ADVECTION', &
1817	& 'Fourth-order centered differences advection of TKE fields'
1818	is=LEN_TRIM(Coptions)+1
1819	Coptions(is:is+15)=' K_C4ADVECTION,'
1820	# else
1821	!
1822	IF (Master) WRITE (stdout,20) 'K_GSCHEME', &
1823	& 'Third-order upstream advection of TKE fields'
1824	is=LEN_TRIM(Coptions)+1
1825	Coptions(is:is+11)=' K_GSCHEME,'
1826	# endif
1827	# ifdef TKE_DIF2
1828	!
1829	IF (Master) WRITE (stdout,20) 'TKE_DIF2', &
1830	& 'Harmonic mixing of TKE fields'
1831	is=LEN_TRIM(Coptions)+1
1832	Coptions(is:is+10)=' TKE_DIF2,'
1833	# endif
1834	# ifdef TKE_DIF4
1835	!
1836	IF (Master) WRITE (stdout,20) 'TKE_DIF4', &
1837	& 'Biharmonic mixing of TKE fields'
1838	is=LEN_TRIM(Coptions)+1
1839	Coptions(is:is+10)=' TKE_DIF4,'
1840	# endif
1841	#endif
1842	#ifdef RADIATION_2D
1843	!
1844	IF (Master) WRITE (stdout,20) 'RADIATION_2D', &
1845	& 'Use tangential phase speed in radiation conditions'
1846	is=LEN_TRIM(Coptions)+1
1847	Coptions(is:is+14)=' RADIATION_2D,'
1848	#endif
1849	#if defined RAMP_TIDES && (defined SSH_TIDES \|\| defined UV_TIDES)
1850	!
1851	IF (Master) WRITE (stdout,20) 'RAMP_TIDES', &
1852	& 'Ramping tidal forcing for one day'
1853	is=LEN_TRIM(Coptions)+1
1854	Coptions(is:is+14)=' RAMP_TIDES,'
1855	#endif
1856	#if defined READ_WATER && defined MASKING
1857	!
1858	IF (Master) WRITE (stdout,20) 'READ_WATER', &
1859	& 'Reading data at water points only'
1860	is=LEN_TRIM(Coptions)+1
1861	Coptions(is:is+12)=' READ_WATER,'
1862	#endif
1863	#if defined RECOMPUTE_4DVAR && defined FOUR_DVAR
1864	!
1865	IF (Master) WRITE (stdout,20) 'RECOMPUTE_4DVAR', &
1866	& 'Recomputing 4DVar during data assimilation analysis'
1867	is=LEN_TRIM(Coptions)+1
1868	Coptions(is:is+17)=' RECOMPUTE_4DVAR,'
1869	#endif
1870	#if defined DISTRIBUTE
1871	# if defined REDUCE_ALLGATHER
1872	!
1873	IF (Master) WRITE (stdout,20) 'REDUCE_ALLGATHER', &
1874	& 'Using mpi_allgather in mp_reduce routine'
1875	is=LEN_TRIM(Coptions)+1
1876	Coptions(is:is+18)=' REDUCE_ALLGATHER,'
1877	# elif defined REDUCE_ALLREDUCE
1878	!
1879	IF (Master) WRITE (stdout,20) 'REDUCE_ALLREDUCE', &
1880	& 'Using mpi_allreduce in mp_reduce routine'
1881	is=LEN_TRIM(Coptions)+1
1882	Coptions(is:is+18)=' REDUCE_ALLREDUCE,'
1883	# else
1884	!
1885	IF (Master) WRITE (stdout,20) '!REDUCE_ALL...', &
1886	& 'Using mpi_isend/mpi_recv in mp_reduce routine'
1887	is=LEN_TRIM(Coptions)+1
1888	Coptions(is:is+16)=' !REDUCE_ALL...,'
1889	# endif
1890	#endif
1891	#ifdef REGCM_COUPLING
1892	!
1893	IF (Master) WRITE (stdout,20) 'COAMPS_COUPLING', &
1894	& 'Atmosphere coupling with RegCM'
1895	is=LEN_TRIM(Coptions)+1
1896	Coptions(is:is+16)=' REGCM_COUPLING,'
1897	iatms=iatms+1
1898	#endif
1899	#ifdef RED_TIDE
1900	!
1901	IF (Master) WRITE (stdout,20) 'RED_TIDE', &
1902	& 'Red tide biological model'
1903	is=LEN_TRIM(Coptions)+1
1904	Coptions(is:is+10)=' RED_TIDE,'
1905	ibiology=ibiology+1
1906	#endif
1907	#if defined REGRESS_STARTCLOCK && defined MODEL_COUPLING && \
1908	defined ESMF_LIB
1909	!
1910	IF (Master) WRITE (stdout,20) 'REGRESS_STARTCLOCK', &
1911	& 'Regressing start clock by one coupling interval'
1912	is=LEN_TRIM(Coptions)+1
1913	Coptions(is:is+20)=' REGRESS_STARTCLOCK,'
1914	#endif
1915	#ifdef ROMS_STDOUT
1916	!
1917	IF (Master) WRITE (stdout,20) 'ROMS_STDOUT', &
1918	& 'Red tide biological model'
1919	is=LEN_TRIM(Coptions)+1
1920	Coptions(is:is+13)=' ROMS_STDOUT,'
1921	#endif
1922	#ifdef R_SYMMETRY
1923	!
1924	IF (Master) WRITE (stdout,20) 'REP_MATRIX', &
1925	& 'Representers Matrix Symmetry Test'
1926	is=LEN_TRIM(Coptions)+1
1927	Coptions(is:is+12)=' REP_MATRIX,'
1928	idriver=idriver+1
1929	#endif
1930	#if defined LMD_MIXING
1931	# ifdef RI_HORAVG
1932	!
1933	IF (Master) WRITE (stdout,20) 'RI_HORAVG', &
1934	& 'Smooth Richardson number horizontally'
1935	is=LEN_TRIM(Coptions)+1
1936	Coptions(is:is+11)=' RI_HORAVG,'
1937	# endif
1938	# ifdef RI_VERAVG
1939	!
1940	IF (Master) WRITE (stdout,20) 'RI_VERAVG', &
1941	& 'Smooth Richardson number vertically'
1942	is=LEN_TRIM(Coptions)+1
1943	Coptions(is:is+11)=' RI_VERAVG,'
1944	# endif
1945	#endif
1946	#if defined SOLVE3D && \
1947	(defined GLS_MIXING \|\| defined LMD_MIXING \|\| defined MY25_MIXING)
1948	# ifdef RI_SPLINES
1949	!
1950	IF (Master) WRITE (stdout,20) 'RI_SPLINES', &
1951	& 'Parabolic Spline Reconstruction for Richardson Number'
1952	is=LEN_TRIM(Coptions)+1
1953	Coptions(is:is+12)=' RI_SPLINES,'
1954	# endif
1955	#endif
1956	#if !(defined PJ_GRADPQ4 \|\| defined PJ_GRADPQ2 \|\| defined PJ_GRADP \|\| \
1957	defined DJ_GRADPS) && defined RHO_SURF && defined SOLVE3D
1958	!
1959	IF (Master) WRITE (stdout,20) 'RHO_SURF', &
1960	& 'Include difference between rho0 and surface density'
1961	is=LEN_TRIM(Coptions)+1
1962	Coptions(is:is+10)=' RHO_SURF,'
1963	#endif
1964	#if defined RPCG && defined WEAK_CONSTRAINT
1965	!
1966	IF (Master) WRITE (stdout,20) 'RPCG', &
1967	& 'Restricted B-preconditioned Lanczos minimization'
1968	is=LEN_TRIM(Coptions)+1
1969	Coptions(is:is+6)=' RPCG,'
1970	#endif
1971	#if defined RP_AVERAGES && defined TL_IOMS
1972	!
1973	IF (Master) WRITE (stdout,20) 'RP_AVERAGES', &
1974	& 'Writing out time-averaged representer model fields'
1975	is=LEN_TRIM(Coptions)+1
1976	Coptions(is:is+13)=' RP_AVERAGES,'
1977	#endif
1978	#ifdef RPM_DRIVER
1979	!
1980	IF (Master) WRITE (stdout,20) 'RPM_DRIVER', &
1981	& 'Generic representers tangent linear model driver'
1982	is=LEN_TRIM(Coptions)+1
1983	Coptions(is:is+12)=' RPM_DRIVER,'
1984	idriver=idriver+1
1985	#endif
1986	#ifdef RPM_RELAXATION
1987	!
1988	IF (Master) WRITE (stdout,20) 'RPM_RELAXATION', &
1989	& 'Diffusive Relaxation of RPM in Picard Iterations'
1990	is=LEN_TRIM(Coptions)+1
1991	Coptions(is:is+16)=' RPM_RELAXATION,'
1992	#endif
1993	#ifdef RST_SINGLE
1994	!
1995	IF (Master) WRITE (stdout,20) 'RST_SINGLE', &
1996	& 'Single precision fields in restart NetCDF file'
1997	is=LEN_TRIM(Coptions)+1
1998	Coptions(is:is+12)=' RST_SINGLE,'
1999	#else
2000	!
2001	IF (Master) WRITE (stdout,20) '!RST_SINGLE', &
2002	& 'Double precision fields in restart NetCDF file'
2003	is=LEN_TRIM(Coptions)+1
2004	Coptions(is:is+13)=' !RST_SINGLE,'
2005	#endif
2006	#ifdef SALINITY
2007	!
2008	IF (Master) WRITE (stdout,20) 'SALINITY', &
2009	& 'Using salinity'
2010	is=LEN_TRIM(Coptions)+1
2011	Coptions(is:is+10)=' SALINITY,'
2012	#endif
2013	#ifdef SANITY_CHECK
2014	!
2015	IF (Master) WRITE (stdout,20) 'SANITY_CHECK', &
2016	& 'Tangent Linear and Adjoint Models Correctness Check'
2017	is=LEN_TRIM(Coptions)+1
2018	Coptions(is:is+14)=' SANITY_CHECK,'
2019	idriver=idriver+1
2020	#endif
2021	#ifdef SCORRECTION
2022	!
2023	IF (Master) WRITE (stdout,20) 'SCORRECTION', &
2024	& 'Surface salinity flux correction'
2025	is=LEN_TRIM(Coptions)+1
2026	Coptions(is:is+13)=' SCORRECTION,'
2027	#endif
2028	#ifdef SEDIMENT
2029	!
2030	IF (Master) WRITE (stdout,20) 'SEDIMENT', &
2031	& 'Cohesive and noncohesive sediments'
2032	is=LEN_TRIM(Coptions)+1
2033	Coptions(is:is+10)=' SEDIMENT,'
2034	# ifdef SED_DENS
2035	!
2036	IF (Master) WRITE (stdout,20) 'SED_DENS', &
2037	& 'Include density of suspended sediment in equation of state'
2038	is=LEN_TRIM(Coptions)+1
2039	Coptions(is:is+10)=' SED_DENS,'
2040	# endif
2041	# ifdef SED_MORPH
2042	!
2043	IF (Master) WRITE (stdout,20) 'SED_MORPH', &
2044	& 'Allow bottom model elevation to evolve'
2045	is=LEN_TRIM(Coptions)+1
2046	Coptions(is:is+11)=' SED_MORPH,'
2047	# endif
2048	# ifdef SUSPLOAD
2049	!
2050	IF (Master) WRITE (stdout,20) 'SUSPLOAD', &
2051	& 'Activate suspended sediment transport'
2052	is=LEN_TRIM(Coptions)+1
2053	Coptions(is:is+10)=' SUSPLOAD,'
2054	# endif
2055	#endif
2056	#ifdef SG_BBL
2057	!
2058	IF (Master) WRITE (stdout,20) 'SG_BBL', &
2059	& 'Styles and Glenn Bottom Boundary Layer'
2060	is=LEN_TRIM(Coptions)+1
2061	Coptions(is:is+8)=' SG_BBL,'
2062	ibbl=ibbl+1
2063	# ifdef SG_CALC_UB
2064	!
2065	IF (Master) WRITE (stdout,20) 'SG_CALC_UB', &
2066	& 'Internal computation of bottom orbital velocity'
2067	is=LEN_TRIM(Coptions)+1
2068	Coptions(is:is+12)=' SG_CALC_UB,'
2069	# endif
2070	# ifdef SG_CALC_ZNOT
2071	!
2072	IF (Master) WRITE (stdout,20) 'SG_CALC_ZNOT', &
2073	& 'Internal computation of bottom roughness'
2074	is=LEN_TRIM(Coptions)+1
2075	Coptions(is:is+14)=' SG_CALC_ZNOT,'
2076	# endif
2077	# ifdef SG_LOGINT
2078	!
2079	IF (Master) WRITE (stdout,20) 'SG_LOGINT', &
2080	& 'Bottom currents logarithmic interpolation'
2081	is=LEN_TRIM(Coptions)+1
2082	Coptions(is:is+11)=' SG_LOGINT,'
2083	# endif
2084	#endif
2085	#ifdef SINGLE_PRECISION
2086	!
2087	IF (Master) WRITE (stdout,20) 'SINGLE_PRECISION', &
2088	& 'Single precision arithmetic numerical kernel.'
2089	is=LEN_TRIM(Coptions)+1
2090	Coptions(is:is+18)=' SINGLE_PRECISION,'
2091	#endif
2092	#if defined SKIP_NLM && defined SENSITIVITY_4DVAR
2093	!
2094	IF (Master) WRITE (stdout,20) 'SKIP_NLM', &
2095	& 'Skipping running NLM, reading NLM state trajectory'
2096	is=LEN_TRIM(Coptions)+1
2097	Coptions(is:is+10)=' SKIP_NLM,'
2098	#endif
2099	#ifdef SRELAXATION
2100	!
2101	IF (Master) WRITE (stdout,20) 'SRELAXATION', &
2102	& 'Surface salinity relaxation as surface flux'
2103	is=LEN_TRIM(Coptions)+1
2104	Coptions(is:is+13)=' SRELAXATION,'
2105	#endif
2106	#ifdef SOLAR_SOURCE
2107	!
2108	IF (Master) WRITE (stdout,20) 'SOLAR_SOURCE', &
2109	& 'Solar Radiation Source Term'
2110	is=LEN_TRIM(Coptions)+1
2111	Coptions(is:is+14)=' SOLAR_SOURCE,'
2112	#endif
2113	#ifdef SOLVE3D
2114	!
2115	IF (Master) WRITE (stdout,20) 'SOLVE3D', &
2116	& 'Solving 3D Primitive Equations'
2117	is=LEN_TRIM(Coptions)+1
2118	Coptions(is:is+9)=' SOLVE3D,'
2119	#endif
2120	#ifdef SO_SEMI
2121	!
2122	IF (Master) WRITE (stdout,20) 'SO_SEMI', &
2123	& 'Stochastic Optimals Propagator, Semi-norm Estimation'
2124	is=LEN_TRIM(Coptions)+1
2125	Coptions(is:is+9)=' SO_SEMI,'
2126	idriver=idriver+1
2127	#endif
2128	#ifdef SO_TRACE
2129	!
2130	IF (Master) WRITE (stdout,20) 'SO_TRACE', &
2131	& 'Stochastic Optimals Propagator, Randomized Trace Estimation'
2132	is=LEN_TRIM(Coptions)+1
2133	Coptions(is:is+10)=' SO_TRACE,'
2134	idriver=idriver+1
2135	#endif
2136	#ifdef SO_SEMI
2137	# ifdef SO_SEMI_WHITE
2138	!
2139	IF (Master) WRITE (stdout,20) 'SO_SEMI_WHITE', &
2140	& 'Stochastic Optimals, white noise processes'
2141	is=LEN_TRIM(Coptions)+1
2142	Coptions(is:is+15)=' SO_SEMI_WHITE,'
2143	# else
2144	!
2145	IF (Master) WRITE (stdout,20) '!SO_WHITE', &
2146	& 'Stochastic Optimals, red noise processes'
2147	is=LEN_TRIM(Coptions)+1
2148	Coptions(is:is+16)=' !SO_SEMI_WHITE,'
2149	# endif
2150	#endif
2151	#if defined SPLINES_VCONV && defined FOUR_DVAR && defined VCONVOLUTION
2152	!
2153	IF (Master) WRITE (stdout,20) 'IMPLICIT_VCONV', &
2154	& 'Implicit, Parabolic Splines Vertical Convolution Algorithm'
2155	is=LEN_TRIM(Coptions)+1
2156	Coptions(is:is+16)=' SPLINES_VCONV,'
2157	#endif
2158	#if defined SPLINES_VDIFF && defined SOLVE3D
2159	!
2160	IF (Master) WRITE (stdout,20) 'SPLINES_VDIFF', &
2161	& 'Parabolic Spline Reconstruction for Vertical Diffusion'
2162	is=LEN_TRIM(Coptions)+1
2163	Coptions(is:is+15)=' SPLINES_VDIFF,'
2164	#endif
2165	#if defined SPLINES_VVISC && defined SOLVE3D
2166	!
2167	IF (Master) WRITE (stdout,20) 'SPLINES_VVISC', &
2168	& 'Parabolic Spline Reconstruction for Vertical Viscosity'
2169	is=LEN_TRIM(Coptions)+1
2170	Coptions(is:is+15)=' SPLINES_VVISC,'
2171	#endif
2172	#ifdef SPHERICAL
2173	!
2174	IF (Master) WRITE (stdout,20) 'SPHERICAL', &
2175	& 'Spherical grid configuration'
2176	is=LEN_TRIM(Coptions)+1
2177	Coptions(is:is+11)=' SPHERICAL,'
2178	#endif
2179	#if defined SSH_TIDES
2180	!
2181	IF (Master) WRITE (stdout,20) 'SSH_TIDES', &
2182	& 'Add tidal elevation to SSH climatology'
2183	is=LEN_TRIM(Coptions)+1
2184	Coptions(is:is+11)=' SSH_TIDES,'
2185	#endif
2186	#ifdef SSW_BBL
2187	!
2188	IF (Master) WRITE (stdout,20) 'SSW_BBL', &
2189	& 'Styles and Glenn Bottom Boundary Layer - modified'
2190	is=LEN_TRIM(Coptions)+1
2191	Coptions(is:is+8)=' SSW_BBL,'
2192	ibbl=ibbl+1
2193	# ifdef SSW_CALC_UB
2194	!
2195	IF (Master) WRITE (stdout,20) 'SSW_CALC_UB', &
2196	& 'Internal computation of bottom orbital velocity'
2197	is=LEN_TRIM(Coptions)+1
2198	Coptions(is:is+14)=' SSW_CALC_UB,'
2199	# endif
2200	# ifdef SSW_CALC_ZNOT
2201	!
2202	IF (Master) WRITE (stdout,20) 'SSW_CALC_ZNOT', &
2203	& 'Internal computation of bottom roughness'
2204	is=LEN_TRIM(Coptions)+1
2205	Coptions(is:is+14)=' SSW_CALC_ZNOT,'
2206	# endif
2207	# ifdef SSW_FORM_DRAG_COR
2208	!
2209	IF (Master) WRITE (stdout,20) 'SSW_FORM_DRAG_COR', &
2210	& 'Form drag coefficient'
2211	is=LEN_TRIM(Coptions)+1
2212	Coptions(is:is+19)=' SSW_FORM_DRAG_COR,'
2213	# endif
2214	# ifdef SSW_ZOBIO
2215	!
2216	IF (Master) WRITE (stdout,20) 'SSW_Z0BIO', &
2217	& 'Biogenic bedfrom roughness from ripples'
2218	is=LEN_TRIM(Coptions)+1
2219	Coptions(is:is+11)=' SSW_Z0BIO,'
2220	# endif
2221	# ifdef SSW_ZOBL
2222	!
2223	IF (Master) WRITE (stdout,20) 'SSW_Z0BL', &
2224	& 'Bedload roughness from ripples'
2225	is=LEN_TRIM(Coptions)+1
2226	Coptions(is:is+10)=' SSW_Z0BL,'
2227	# endif
2228	# ifdef SSW_ZORIP
2229	!
2230	IF (Master) WRITE (stdout,20) 'SSW_Z0RIP', &
2231	& 'Bedfrom roughness from ripples'
2232	is=LEN_TRIM(Coptions)+1
2233	Coptions(is:is+11)=' SSW_Z0RIP,'
2234	# endif
2235	# ifdef SSW_LOGINT
2236	!
2237	IF (Master) WRITE (stdout,20) 'SSW_LOGINT', &
2238	& 'Bottom currents logarithmic interpolation'
2239	is=LEN_TRIM(Coptions)+1
2240	Coptions(is:is+11)=' SSW_LOGINT,'
2241	# endif
2242	#endif
2243	#ifdef STATIONS
2244	!
2245	IF (Master) WRITE (stdout,20) 'STATIONS', &
2246	& 'Writing out station data'
2247	is=LEN_TRIM(Coptions)+1
2248	Coptions(is:is+10)=' STATIONS,'
2249	#endif
2250	#if defined STATIONS_CGRID && defined STATIONS
2251	!
2252	IF (Master) WRITE (stdout,20) 'STATIONS_CGRID', &
2253	& 'Extracting station data at native C-grid locations'
2254	is=LEN_TRIM(Coptions)+1
2255	Coptions(is:is+16)=' STATIONS_CGRID,'
2256	#endif
2257	#ifdef STOCHASTIC_OPT
2258	!
2259	IF (Master) WRITE (stdout,20) 'STOCHASTIC_OPT', &
2260	& 'Stochastic Optimals Propagator'
2261	is=LEN_TRIM(Coptions)+1
2262	Coptions(is:is+16)=' STOCHASTIC_OPT,'
2263	idriver=idriver+1
2264	# ifdef STOCH_OPT_WHITE
2265	!
2266	IF (Master) WRITE (stdout,20) 'STOCH_OPT_WHITE', &
2267	& 'Stochastic Optimals, white noise processes'
2268	is=LEN_TRIM(Coptions)+1
2269	Coptions(is:is+17)=' STOCH_OPT_WHITE,'
2270	# else
2271	!
2272	IF (Master) WRITE (stdout,20) '!STOCH_OPT_WHITE', &
2273	& 'Stochastic Optimals, red noise processes'
2274	is=LEN_TRIM(Coptions)+1
2275	Coptions(is:is+18)=' !STOCH_OPT_WHITE,'
2276	# endif
2277	#endif
2278	#if defined SURFACE_DO_SATURATION && defined HYPOXIA_SRM
2279	!
2280	IF (Master) WRITE (stdout,20) 'SURFACE_DO_SATURATION', &
2281	& 'Use surface dissolved oxygen saturation'
2282	is=LEN_TRIM(Coptions)+1
2283	Coptions(is:is+23)=' SURFACE_DO_SATURATION,'
2284	#endif
2285	#ifdef SWAN_COUPLING
2286	!
2287	IF (Master) WRITE (stdout,20) 'SWAN_COUPLING', &
2288	& 'Two-way SWAN/ROMS coupling'
2289	is=LEN_TRIM(Coptions)+1
2290	Coptions(is:is+15)=' SWAN_COUPLING,'
2291	#endif
2292	#if defined CARBON && defined BIO_FENNEL
2293	# ifdef TALK_NONCONSERV
2294	!
2295	IF (Master) WRITE (stdout,20) 'TALK_NONCONSERV', &
2296	& 'Alkalinity is affected by changes in nitrate or ammonium'
2297	is=LEN_TRIM(Coptions)+1
2298	Coptions(is:is+17)=' TALK_NONCONSERV,'
2299	# else
2300	!
2301	IF (Master) WRITE (stdout,20) '!TALK_NONCONSERV', &
2302	& 'Alkalinity is passive and unaffected by nitrate or ammonium'
2303	is=LEN_TRIM(Coptions)+1
2304	Coptions(is:is+18)=' !TALK_NONCONSERV,'
2305	# endif
2306	#endif
2307	#ifdef TANGENT
2308	!
2309	IF (Master) WRITE (stdout,20) 'TANGENT', &
2310	& 'Tangent Linear Model'
2311	is=LEN_TRIM(Coptions)+1
2312	Coptions(is:is+9)=' TANGENT,'
2313	#endif
2314	#if defined TIME_CONV && defined WEAK_CONSTRAINT
2315	!
2316	IF (Master) WRITE (stdout,20) 'TIME_CONV', &
2317	& 'Weak-constraint 4D-Var time convolution'
2318	is=LEN_TRIM(Coptions)+1
2319	Coptions(is:is+11)=' TIME_CONV,'
2320	#endif
2321	#if defined TIME_INTERP && defined MODEL_COUPLING
2322	!
2323	IF (Master) WRITE (stdout,20) 'TIME_INTERP', &
2324	& 'Time interpolation from coupling snapshots arrays'
2325	is=LEN_TRIM(Coptions)+1
2326	Coptions(is:is+13)=' TIME_INTERP,'
2327	#endif
2328	#if defined TIME_INTERP_FLUX && defined NESTING && defined SOLVE3D
2329	!
2330	IF (Master) WRITE (stdout,20) 'TIME_INTERP_FLUX', &
2331	& 'Time interpolate coaser grid mass flux in refinement'
2332	is=LEN_TRIM(Coptions)+1
2333	Coptions(is:is+18)=' TIME_INTERP_FLUX,'
2334	#endif
2335	#if defined TL_AVERAGES && defined TANGENT
2336	!
2337	IF (Master) WRITE (stdout,20) 'TL_AVERAGES', &
2338	& 'Writing out time-averaged tangent linear model fields'
2339	is=LEN_TRIM(Coptions)+1
2340	Coptions(is:is+13)=' TL_AVERAGES,'
2341	#endif
2342	#ifdef TLM_DRIVER
2343	!
2344	IF (Master) WRITE (stdout,20) 'TLM_DRIVER', &
2345	& 'Generic tangent linear model driver'
2346	is=LEN_TRIM(Coptions)+1
2347	Coptions(is:is+12)=' TLM_DRIVER,'
2348	idriver=idriver+1
2349	#endif
2350	#if defined GLS_MIXING && defined TKE_WAVEDISS
2351	!
2352	IF (Master) WRITE (stdout,20) 'TKE_WAVEDISS', &
2353	& 'Wave breaking surface tke flux based on wave amplitude'
2354	is=LEN_TRIM(Coptions)+1
2355	Coptions(is:is+14)=' TKE_WAVEDISS,'
2356	#endif
2357	#ifdef TL_IOMS
2358	!
2359	IF (Master) WRITE (stdout,20) 'TL_IOMS', &
2360	& 'Representers Tangent Linear Model'
2361	is=LEN_TRIM(Coptions)+1
2362	Coptions(is:is+9)=' TL_IOMS,'
2363	#endif
2364	#ifdef TLM_CHECK
2365	!
2366	IF (Master) WRITE (stdout,20) 'TLM_CHECK', &
2367	& 'Tangent Linear Model Linearization Check'
2368	is=LEN_TRIM(Coptions)+1
2369	Coptions(is:is+11)=' TLM_CHECK,'
2370	#endif
2371	#ifdef SOLVE3D
2372	# if defined T_PASSIVE
2373	!
2374	IF (Master) WRITE (stdout,20) 'T_PASSIVE', &
2375	& 'Advecting and diffusing inert passive tracer'
2376	is=LEN_TRIM(Coptions)+1
2377	Coptions(is:is+11)=' T_PASSIVE,'
2378	# endif
2379	# if defined TS_MIX_CLIMA && \
2380	(defined TS_DIF2 \|\| defined TS_DIF4)
2381	!
2382	IF (Master) WRITE (stdout,20) 'TS_MIX_CLIMA', &
2383	& 'Horizontal diffusion of tracer perturbation (T-Tclm)'
2384	is=LEN_TRIM(Coptions)+1
2385	Coptions(is:is+14)=' TS_MIX_CLIMA,'
2386	# endif
2387	# if defined TS_MIX_MAX_SLOPE && defined MIX_ISO_TS && \
2388	(defined TS_DIF2 \|\| defined TS_DIF4)
2389	!
2390	IF (Master) WRITE (stdout,20) 'TS_MIX_MAX_SLOPE', &
2391	& 'Maximum density slope in isopycnal diffusion'
2392	is=LEN_TRIM(Coptions)+1
2393	Coptions(is:is+18)=' TS_MIX_MAX_SLOPE,'
2394	# endif
2395	# if defined TS_MIX_MAX_SLOPE && defined MIX_ISO_TS && \
2396	(defined TS_DIF2 \|\| defined TS_DIF4)
2397	!
2398	IF (Master) WRITE (stdout,20) 'TS_MIX_MAX_STRAT', &
2399	& 'Minimum density stratification in isopycnal diffusion'
2400	is=LEN_TRIM(Coptions)+1
2401	Coptions(is:is+18)=' TS_MIX_MAX_STRAT,'
2402	# endif
2403	# if defined TS_MIX_STABILITY && \
2404	(defined TS_DIF2 \|\| defined TS_DIF4)
2405	!
2406	IF (Master) WRITE (stdout,20) 'TS_MIX_STABILITY', &
2407	& 'Weighting diffusion between two time levels'
2408	is=LEN_TRIM(Coptions)+1
2409	Coptions(is:is+18)=' TS_MIX_STABILITY,'
2410	# endif
2411	# ifdef TS_MPDATA
2412	!
2413	IF (Master) WRITE (stdout,20) 'TS_MPDATA', &
2414	& 'Recursive flux corrected MPDATA 3D advection of tracers'
2415	is=LEN_TRIM(Coptions)+1
2416	Coptions(is:is+11)=' TS_MPDATA,'
2417	itrcHadv=itrcHadv+1
2418	itrcVadv=itrcVadv+1
2419	# ifdef TS_MPDATA_LIMIT
2420	!
2421	IF (Master) WRITE (stdout,20) 'TS_MPDATA_LIMIT', &
2422	& 'Further limiter in upwind corrector fluxes for stability'
2423	is=LEN_TRIM(Coptions)+1
2424	Coptions(is:is+17)=' TS_MPDATA_LIMIT,'
2425	# endif
2426	# endif
2427	# ifdef TS_A4HADVECTION
2428	!
2429	IF (Master) WRITE (stdout,20) 'TS_A4HADVECTION', &
2430	& 'Fourth-order Akima horizontal advection of tracers'
2431	is=LEN_TRIM(Coptions)+1
2432	Coptions(is:is+17)=' TS_A4HADVECTION,'
2433	itrcHadv=itrcHadv+1
2434	# endif
2435	# if defined TS_A4HADVECTION_TL && \
2436	(defined TANGENT \|\| defined ADJOINT \|\| defined TL_IOMS)
2437	!
2438	IF (Master) WRITE (stdout,20) 'TS_A4HADVECTION_TL', &
2439	& 'TL/AD fourth-order Akima horizontal tracer advection'
2440	is=LEN_TRIM(Coptions)+1
2441	Coptions(is:is+20)=' TS_A4HADVECTION_TL,'
2442	itrcHadvtl=itrcHadvtl+1
2443	# endif
2444	# ifdef TS_C2HADVECTION
2445	!
2446	IF (Master) WRITE (stdout,20) 'TS_C2HADVECTION', &
2447	& 'Second-order centered horizontal advection of tracers'
2448	is=LEN_TRIM(Coptions)+1
2449	Coptions(is:is+17)=' TS_C2HADVECTION,'
2450	itrcHadv=itrcHadv+1
2451	# endif
2452	# if defined TS_C2HADVECTION_TL && \
2453	(defined TANGENT \|\| defined ADJOINT \|\| defined TL_IOMS)
2454	!
2455	IF (Master) WRITE (stdout,20) 'TS_C2HADVECTION_TL', &
2456	& 'TL/AD second-order centered horizontal tracer advection'
2457	is=LEN_TRIM(Coptions)+1
2458	Coptions(is:is+20)=' TS_C2HADVECTION_TL,'
2459	itrcHadvtl=itrcHadvtl+1
2460	# endif
2461	# if defined TS_C4HADVECTION \|\| defined TS_U3ADV_SPLIT
2462	!
2463	IF (Master) WRITE (stdout,20) 'TS_C4HADVECTION', &
2464	& 'Fourth-order centered horizontal advection of tracers'
2465	is=LEN_TRIM(Coptions)+1
2466	Coptions(is:is+17)=' TS_C4HADVECTION,'
2467	itrcHadv=itrcHadv+1
2468	# endif
2469	# if defined TS_C4HADVECTION_TL && \
2470	(defined TANGENT \|\| defined ADJOINT \|\| defined TL_IOMS)
2471	!
2472	IF (Master) WRITE (stdout,20) 'TS_C4HADVECTION_TL', &
2473	& 'TL/AD fourth-order centered horizontal tracer advection'
2474	is=LEN_TRIM(Coptions)+1
2475	Coptions(is:is+20)=' TS_C4HADVECTION_TL,'
2476	itrcHadvtl=itrcHadvtl+1
2477	# endif
2478	# ifdef TS_U3HADVECTION
2479	!
2480	IF (Master) WRITE (stdout,20) 'TS_U3HADVECTION', &
2481	& 'Third-order upstream horizontal advection of tracers'
2482	is=LEN_TRIM(Coptions)+1
2483	Coptions(is:is+17)=' TS_U3HADVECTION,'
2484	itrcHadv=itrcHadv+1
2485	# endif
2486	# if defined TS_U3HADVECTION_TL && \
2487	(defined TANGENT \|\| defined ADJOINT \|\| defined TL_IOMS)
2488	!
2489	IF (Master) WRITE (stdout,20) 'TS_U3HADVECTION_TL', &
2490	& 'TL/AD third-order upstream horizontal tracer advection'
2491	is=LEN_TRIM(Coptions)+1
2492	Coptions(is:is+20)=' TS_U3HADVECTION_TL,'
2493	itrcHadvtl=itrcHadvtl+1
2494	# endif
2495	# if !(defined TS_MPDATA \|\| defined TS_A4HADVECTION \|\| \
2496	defined TS_C2HADVECTION \|\| defined TS_C4HADVECTION \|\| \
2497	defined TS_U3HADVECTION)
2498	!
2499	IF (Master) WRITE (stdout,20) 'TS_C4HADVECTION', &
2500	& 'Fourth-order centered horizontal advection of tracers'
2501	is=LEN_TRIM(Coptions)+1
2502	Coptions(is:is+17)=' TS_C4HADVECTION,'
2503	itrcHadv=itrcHadv+1
2504	# endif
2505	# ifdef TS_A4VADVECTION
2506	!
2507	IF (Master) WRITE (stdout,20) 'TS_A4VADVECTION', &
2508	& 'Fourth-order Akima vertical advection of tracers'
2509	is=LEN_TRIM(Coptions)+1
2510	Coptions(is:is+17)=' TS_A4VADVECTION,'
2511	itrcVadv=itrcVadv+1
2512	# endif
2513	# if defined TS_A4VADVECTION_TL && \
2514	(defined TANGENT \|\| defined ADJOINT \|\| defined TL_IOMS)
2515	!
2516	IF (Master) WRITE (stdout,20) 'TS_A4VADVECTION_TL', &
2517	& 'TL/AD fourth-order Akima vertical tracer advection'
2518	is=LEN_TRIM(Coptions)+1
2519	Coptions(is:is+20)=' TS_A4VADVECTION_TL,'
2520	itrcVadvtl=itrcVadvtl+1
2521	# endif
2522	# ifdef TS_C2VADVECTION
2523	!
2524	IF (Master) WRITE (stdout,20) 'TS_C2VADVECTION', &
2525	& 'Second-order centered vertical advection of tracers'
2526	is=LEN_TRIM(Coptions)+1
2527	Coptions(is:is+17)=' TS_C2VADVECTION,'
2528	itrcVadv=itrcVadv+1
2529	# endif
2530	# if defined TS_C2VADVECTION_TL && \
2531	(defined TANGENT \|\| defined ADJOINT \|\| defined TL_IOMS)
2532	!
2533	IF (Master) WRITE (stdout,20) 'TS_C2VADVECTION_TL', &
2534	& 'TL/AD second-order centered vertical tracer advection'
2535	is=LEN_TRIM(Coptions)+1
2536	Coptions(is:is+20)=' TS_C2VADVECTION_TL,'
2537	itrcVadvtl=itrcVadvtl+1
2538	# endif
2539	# if defined TS_C4VADVECTION \|\| defined TS_U3ADV_SPLIT
2540	!
2541	IF (Master) WRITE (stdout,20) 'TS_C4VADVECTION', &
2542	& 'Fourth-order centered vertical advection of tracers'
2543	is=LEN_TRIM(Coptions)+1
2544	Coptions(is:is+17)=' TS_C4VADVECTION,'
2545	itrcVadv=itrcVadv+1
2546	# endif
2547	# if defined TS_C4VADVECTION_TL && \
2548	(defined TANGENT \|\| defined ADJOINT \|\| defined TL_IOMS)
2549	!
2550	IF (Master) WRITE (stdout,20) 'TS_C4VADVECTION_TL', &
2551	& 'TL/AD fourth-order centered vertical tracer advection'
2552	is=LEN_TRIM(Coptions)+1
2553	Coptions(is:is+20)=' TS_C4VADVECTION_TL,'
2554	itrcVadvtl=itrcVadvtl+1
2555	# endif
2556	# ifdef TS_SVADVECTION
2557	!
2558	IF (Master) WRITE (stdout,20) 'TS_SVADVECTION', &
2559	& 'Parabolic splines vertical advection of tracers'
2560	is=LEN_TRIM(Coptions)+1
2561	Coptions(is:is+16)=' TS_SVADVECTION,'
2562	itrcVadv=itrcVadv+1
2563	# endif
2564	# if defined TS_SVADVECTION_TL && \
2565	(defined TANGENT \|\| defined ADJOINT \|\| defined TL_IOMS)
2566	!
2567	IF (Master) WRITE (stdout,20) 'TS_SVADVECTION_TL', &
2568	& 'TL/AD parabolic splines vertical tracer advection'
2569	is=LEN_TRIM(Coptions)+1
2570	Coptions(is:is+19)=' TS_SVADVECTION_TL,'
2571	itrcVadvtl=itrcVadvtl+1
2572	# endif
2573	# if !(defined TS_MPDATA \|\| defined TS_A4VADVECTION \|\| \
2574	defined TS_C2VADVECTION \|\| defined TS_C4VADVECTION \|\| \
2575	defined TS_SVADVECTION)
2576	!
2577	IF (Master) WRITE (stdout,20) 'TS_C4VADVECTION', &
2578	& 'Fourth-order centered vertical advection of tracers'
2579	is=LEN_TRIM(Coptions)+1
2580	Coptions(is:is+17)=' TS_C4VADVECTION,'
2581	itrcVadv=itrcVadv+1
2582	# endif
2583	# ifdef TS_U3ADV_SPLIT
2584	!
2585	IF (Master) WRITE (stdout,20) 'TS_U3ADV_SPLIT', &
2586	& 'Split third-order upstream advection of tracers'
2587	is=LEN_TRIM(Coptions)+1
2588	Coptions(is:is+16)=' TS_U3ADV_SPLIT,'
2589	# endif
2590	#endif
2591	#if defined TS_DIF2 && defined SOLVE3D
2592	!
2593	IF (Master) WRITE (stdout,20) 'TS_DIF2', &
2594	& 'Harmonic mixing of tracers'
2595	is=LEN_TRIM(Coptions)+1
2596	Coptions(is:is+9)=' TS_DIF2,'
2597	#endif
2598	#if defined TS_DIF4 && defined SOLVE3D
2599	!
2600	IF (Master) WRITE (stdout,20) 'TS_DIF4', &
2601	& 'Biharmonic mixing of tracers'
2602	is=LEN_TRIM(Coptions)+1
2603	Coptions(is:is+9)=' TS_DIF4,'
2604	#endif
2605	#ifdef TS_FIXED
2606	!
2607	IF (Master) WRITE (stdout,20) 'TS_FIXED', &
2608	& 'Diagnostic configuration, no evolution of tracer'
2609	is=LEN_TRIM(Coptions)+1
2610	Coptions(is:is+10)=' TS_FIXED,'
2611	#endif
2612	#if defined TS_SMAGORINSKY && (defined TS_DIF2 \|\| defined TS_DIF4)
2613	!
2614	IF (Master) WRITE (stdout,20) 'TS_SMAGORINSKY', &
2615	& 'Smagorinksy-like time-dependent diffusion coefficients'
2616	is=LEN_TRIM(Coptions)+1
2617	Coptions(is:is+16)=' TS_SMAGORINSKY,'
2618	#endif
2619	#ifdef UV_ADV
2620	!
2621	IF (Master) WRITE (stdout,20) 'UV_ADV', &
2622	& 'Advection of momentum'
2623	is=LEN_TRIM(Coptions)+1
2624	Coptions(is:is+8)=' UV_ADV,'
2625	#endif
2626	#ifdef UV_COR
2627	!
2628	IF (Master) WRITE (stdout,20) 'UV_COR', &
2629	& 'Coriolis term'
2630	is=LEN_TRIM(Coptions)+1
2631	Coptions(is:is+8)=' UV_COR,'
2632	#endif
2633	#ifdef UV_ADV
2634	# ifdef UV_C2ADVECTION
2635	!
2636	IF (Master) WRITE (stdout,20) 'UV_C2ADVECTION', &
2637	& 'Second-order centered differences advection of momentum'
2638	is=LEN_TRIM(Coptions)+1
2639	Coptions(is:is+16)=' UV_C2ADVECTION,'
2640	ivelHadv=ivelHadv+1
2641	ivelVadv=ivelVadv+1
2642	# endif
2643	# if defined UV_C4ADVECTION \|\| defined UV_U3ADV_SPLIT
2644	!
2645	IF (Master) WRITE (stdout,20) 'UV_C4ADVECTION', &
2646	& 'Fourth-order centered differences advection of momentum'
2647	is=LEN_TRIM(Coptions)+1
2648	Coptions(is:is+16)=' UV_C4ADVECTION,'
2649	ivelHadv=ivelHadv+1
2650	ivelVadv=ivelVadv+1
2651	# endif
2652	# ifdef SOLVE3D
2653	# if !(defined UV_C2ADVECTION \|\| defined UV_C4ADVECTION)
2654	!
2655	IF (Master) WRITE (stdout,20) 'UV_U3HADVECTION', &
2656	& 'Third-order upstream horizontal advection of 3D momentum'
2657	is=LEN_TRIM(Coptions)+1
2658	Coptions(is:is+17)=' UV_U3HADVECTION,'
2659	ivelHadv=ivelHadv+1
2660	# endif
2661	# if !(defined UV_C2ADVECTION \|\| defined UV_C4ADVECTION \|\| \
2662	defined UV_SADVECTION)
2663	!
2664	IF (Master) WRITE (stdout,20) 'UV_C4VADVECTION', &
2665	& 'Fourth-order centered vertical advection of momentum'
2666	is=LEN_TRIM(Coptions)+1
2667	Coptions(is:is+16)=' UV_C4VADVECTION,'
2668	ivelVadv=ivelVadv+1
2669	# endif
2670	# ifdef UV_SADVECTION
2671	!
2672	IF (Master) WRITE (stdout,20) 'UV_SADVECTION', &
2673	& 'Parabolic splines vertical advection of momentum'
2674	is=LEN_TRIM(Coptions)+1
2675	Coptions(is:is+15)=' UV_SADVECTION,'
2676	ivelVadv=ivelVadv+1
2677	# endif
2678	# ifdef UV_U3ADV_SPLIT
2679	!
2680	IF (Master) WRITE (stdout,20) 'UV_U3ADV_SPLIT', &
2681	& 'Split third-order upstream advection of momentum'
2682	is=LEN_TRIM(Coptions)+1
2683	Coptions(is:is+16)=' UV_U3ADV_SPLIT,'
2684	# endif
2685	# endif
2686	#endif
2687	#ifdef UV_DRAG_GRID
2688	!
2689	IF (Master) WRITE (stdout,20) 'UV_DRAG_GRID', &
2690	# if defined UV_LOGDRAG
2691	& 'Spatially varying bottom roughness length'
2692	# elif defined UV_LDRAG
2693	& 'Spatially varying linear drag coefficient'
2694	# elif defined UV_QDRAG
2695	& 'Spatially varying quadratic drag coefficient'
2696	# endif
2697	is=LEN_TRIM(Coptions)+1
2698	Coptions(is:is+14)=' UV_DRAG_GRID,'
2699	#endif
2700	#ifdef UV_LDRAG
2701	!
2702	IF (Master) WRITE (stdout,20) 'UV_LDRAG', &
2703	& 'Linear bottom stress'
2704	is=LEN_TRIM(Coptions)+1
2705	Coptions(is:is+10)=' UV_LDRAG,'
2706	ibbl=ibbl+1
2707	#endif
2708	#ifdef UV_LOGDRAG
2709	!
2710	IF (Master) WRITE (stdout,20) 'UV_LOGDRAG', &
2711	& 'Logarithmic bottom stress'
2712	is=LEN_TRIM(Coptions)+1
2713	Coptions(is:is+12)=' UV_LOGDRAG,'
2714	ibbl=ibbl+1
2715	#endif
2716	#ifdef UV_QDRAG
2717	!
2718	IF (Master) WRITE (stdout,20) 'UV_QDRAG', &
2719	& 'Quadratic bottom stress'
2720	is=LEN_TRIM(Coptions)+1
2721	Coptions(is:is+10)=' UV_QDRAG,'
2722	ibbl=ibbl+1
2723	#endif
2724	#if defined UV_TIDES
2725	!
2726	IF (Master) WRITE (stdout,20) 'UV_TIDES', &
2727	& 'Add tidal currents to 2D momentum climatologies'
2728	is=LEN_TRIM(Coptions)+1
2729	Coptions(is:is+10)=' UV_TIDES,'
2730	#endif
2731	#ifdef UV_VIS2
2732	!
2733	IF (Master) WRITE (stdout,20) 'UV_VIS2', &
2734	& 'Harmonic mixing of momentum'
2735	is=LEN_TRIM(Coptions)+1
2736	Coptions(is:is+9)=' UV_VIS2,'
2737	#endif
2738	#ifdef UV_VIS4
2739	!
2740	IF (Master) WRITE (stdout,20) 'UV_VIS4', &
2741	& 'Biharmonic mixing of momentum'
2742	is=LEN_TRIM(Coptions)+1
2743	Coptions(is:is+9)=' UV_VIS4,'
2744	#endif
2745	#if defined UV_SMAGORINSKY && (defined UV_VIS2 \|\| defined UV_VIS4)
2746	!
2747	IF (Master) WRITE (stdout,20) 'UV_SMAGORINSKY', &
2748	& 'Smagorinksy-like time-dependent viscosity coefficients'
2749	is=LEN_TRIM(Coptions)+1
2750	Coptions(is:is+16)=' UV_SMAGORINSKY,'
2751	#endif
2752	#if defined VAR_RHO_2D && defined SOLVE3D
2753	!
2754	IF (Master) WRITE (stdout,20) 'VAR_RHO_2D', &
2755	& 'Variable density barotropic mode'
2756	is=LEN_TRIM(Coptions)+1
2757	Coptions(is:is+12)=' VAR_RHO_2D,'
2758	#endif
2759	#if defined UV_VIS2 \|\| defined UV_VIS4
2760	# ifdef VISC_GRID
2761	!
2762	IF (Master) WRITE (stdout,20) 'VISC_GRID', &
2763	& 'Horizontal viscosity coefficient scaled by grid size'
2764	is=LEN_TRIM(Coptions)+1
2765	Coptions(is:is+11)=' VISC_GRID,'
2766	# endif
2767	#endif
2768	#if defined VCONVOLUTION && defined FOUR_DVAR && defined SOLVE3D
2769	!
2770	IF (Master) WRITE (stdout,20) 'VCONVOLUTION', &
2771	& 'Include vertical correlations in convolutions'
2772	is=LEN_TRIM(Coptions)+1
2773	Coptions(is:is+14)=' VCONVOLUTION,'
2774	#endif
2775	#ifdef VERIFICATION
2776	!
2777	IF (Master) WRITE (stdout,20) 'VERIFICATION', &
2778	& 'Proccess model solution at observation locations'
2779	is=LEN_TRIM(Coptions)+1
2780	Coptions(is:is+14)=' VERIFICATION,'
2781	#endif
2782	#if defined FLOATS && defined FLOAT_VWALK && defined SOLVE3D
2783	# ifdef VWALK_FORWARD
2784	!
2785	IF (Master) WRITE (stdout,20) 'VWALK_FORWARD', &
2786	& 'Lagrangian drifters, forward stepping random walk'
2787	is=LEN_TRIM(Coptions)+1
2788	Coptions(is:is+15)=' VWALK_FORWARD,'
2789	# else
2790	!
2791	IF (Master) WRITE (stdout,20) '!VWALK_FORWARD', &
2792	& 'Lagrangian drifters, predictor/corrector random walk'
2793	is=LEN_TRIM(Coptions)+1
2794	Coptions(is:is+16)=' !VWALK_FORWARD,'
2795	# endif
2796	#endif
2797	#if defined VISC_3DCOEF
2798	!
2799	IF (Master) WRITE (stdout,20) 'VISC_3DCOEF', &
2800	& 'Horizontal, time-dependent 3D viscosity coefficient'
2801	is=LEN_TRIM(Coptions)+1
2802	Coptions(is:is+13)=' VISC_3DCOEF,'
2803	#endif
2804	#ifdef W4DPSAS
2805	!
2806	IF (Master) WRITE (stdout,20) 'W4DPSAS', &
2807	& 'Weak constraint 4D-PSAS data assimilation'
2808	is=LEN_TRIM(Coptions)+1
2809	Coptions(is:is+9)=' W4DPSAS,'
2810	idriver=idriver+1
2811	#endif
2812	#ifdef W4DPSAS_SENSITIVITY
2813	!
2814	IF (Master) WRITE (stdout,20) 'W4DPSAS_SENSITIVITY', &
2815	& 'Weak constraint 4D-PSAS observation sensitivity'
2816	is=LEN_TRIM(Coptions)+1
2817	Coptions(is:is+21)=' W4DPSAS_SENSITIVITY,'
2818	idriver=idriver+1
2819	#endif
2820	#ifdef W4DVAR
2821	!
2822	IF (Master) WRITE (stdout,20) 'W4DVAR', &
2823	& 'Weak constraint 4DVAR data assimilation'
2824	is=LEN_TRIM(Coptions)+1
2825	Coptions(is:is+8)=' W4DVAR,'
2826	idriver=idriver+1
2827	#endif
2828	#ifdef W4DVAR_SENSITIVITY
2829	!
2830	IF (Master) WRITE (stdout,20) 'W4DVAR_SENSITIVITY', &
2831	& 'Weak constraint 4DVAR observation sensitivity'
2832	is=LEN_TRIM(Coptions)+1
2833	Coptions(is:is+20)=' W4DVAR_SENSITIVITY,'
2834	idriver=idriver+1
2835	#endif
2836	#ifdef WAV_COUPLING
2837	!
2838	IF (Master) WRITE (stdout,20) 'WAV_COUPLING', &
2839	& 'Two-way wave-ocean models coupling'
2840	is=LEN_TRIM(Coptions)+1
2841	Coptions(is:is+13)=' WAV_COUPLING,'
2842	#endif
2843	#ifdef WEAK_CONSTRAINT
2844	!
2845	IF (Master) WRITE (stdout,20) 'WEAK_CONSTRAINT', &
2846	& 'Activated weak constraint assimilation set-up'
2847	is=LEN_TRIM(Coptions)+1
2848	Coptions(is:is+17)=' WEAK_CONSTRAINT,'
2849	#endif
2850	#if defined WRITE_WATER && defined MASKING
2851	!
2852	IF (Master) WRITE (stdout,20) 'WRITE_WATER', &
2853	& 'Writing data at water points only'
2854	is=LEN_TRIM(Coptions)+1
2855	Coptions(is:is+13)=' WRITE_WATER,'
2856	#endif
2857	#ifdef WET_DRY
2858	!
2859	IF (Master) WRITE (stdout,20) 'WET_DRY', &
2860	& 'Wetting and drying activated'
2861	is=LEN_TRIM(Coptions)+1
2862	Coptions(is:is+9)=' WET_DRY,'
2863	#endif
2864	#ifdef WIND_MINUS_CURRENT && defined BULK_FLUXES && defined SOLVE3D
2865	!
2866	IF (Master) WRITE (stdout,20) 'WIND_MINUS_CURRENT', &
2867	& 'Compute effective wind by removing surface ocean current'
2868	is=LEN_TRIM(Coptions)+1
2869	Coptions(is:is+20)=' WIND_MINUS_CURRENT,'
2870	#endif
2871	#ifdef WIND_WAVES
2872	!
2873	IF (Master) WRITE (stdout,20) 'WIND_WAVES', &
2874	& 'Wind-induced wave data available'
2875	is=LEN_TRIM(Coptions)+1
2876	Coptions(is:is+11)=' WIND_WAVES,'
2877	#endif
2878	#if !(defined PJ_GRADPQ4 \|\| defined PJ_GRADPQ2 \|\| defined PJ_GRADP \|\| \
2879	defined DJ_GRADPS) && defined WJ_GRADP
2880	!
2881	IF (Master) WRITE (stdout,20) 'WJ_GRADP', &
2882	& 'Weighted Jacobians pressure gradient adjustment'
2883	is=LEN_TRIM(Coptions)+1
2884	Coptions(is:is+10)=' WJ_GRADP,'
2885	#endif
2886	#ifdef WRF_COUPLING
2887	!
2888	IF (Master) WRITE (stdout,20) 'WRF_COUPLING', &
2889	& 'Atmosphere coupling with WRF'
2890	is=LEN_TRIM(Coptions)+1
2891	Coptions(is:is+14)=' WRF_COUPLING,'
2892	iatms=iatms+1
2893	# ifdef WRF_TIMEAVG
2894	!
2895	IF (Master) WRITE (stdout,20) 'WRF_TIMEAVG', &
2896	& 'WRF Exports time-averaged fields over the coupling interval'
2897	is=LEN_TRIM(Coptions)+1
2898	Coptions(is:is+13)=' WRF_TIMEAVG,'
2899	# endif
2900	#endif
2901	#if defined WTYPE_GRID && \
2902	(defined LMD_SKPP \|\| defined SOLAR_SOURCE) && \
2903	!defined ANA_WTYPE
2904	!
2905	IF (Master) WRITE (stdout,20) 'WTYPE_GRID', &
2906	& 'Spatially varying Jerlov water type index'
2907	is=LEN_TRIM(Coptions)+1
2908	Coptions(is:is+12)=' WTYPE_GRID,'
2909	#endif
2910	#if defined ZETA_ELLIPTIC && defined BALANCE_OPERATOR && \
2911	defined FOUR_DVAR
2912	!
2913	IF (Master) WRITE (stdout,20) 'ZETA_ELLIPTIC', &
2914	& 'Solving SSH elliptic equation in balance operator'
2915	is=LEN_TRIM(Coptions)+1
2916	Coptions(is:is+15)=' ZETA_ELLIPTIC,'
2917	#endif
2918	#ifdef ZOS_HSIG
2919	!
2920	IF (Master) WRITE (stdout,20) 'ZOS_HSIG', &
2921	& 'Wave amplitude used for Zos calculation'
2922	is=LEN_TRIM(Coptions)+1
2923	Coptions(is:is+14)=' ZOS_HSIG,'
2924	#endif
2925	!
2926	!-----------------------------------------------------------------------
2927	! Stop if unsupported C-preprocessing options or report issues with
2928	! particular options.
2929	!-----------------------------------------------------------------------
2930	!
2931	CALL checkadj
2932	!
2933	!-----------------------------------------------------------------------
2934	! Check C-preprocessing options.
2935	!-----------------------------------------------------------------------
2936	!
2937	! Stop if more than one vertical closure scheme is selected.
2938	!
2939	IF (Master.and.(ivmix.gt.1)) THEN
2940	WRITE (stdout,30)
2941	30 FORMAT (/,' CHECKDEFS - only one vertical closure scheme', &
2942	& ' is allowed.')
2943	exit_flag=5
2944	END IF
2945	!
2946	! Stop if more that one bottom stress formulation is selected.
2947	!
2948	IF (Master.and.(ibbl.gt.1)) THEN
2949	WRITE (stdout,40)
2950	40 FORMAT (/,' CHECKDEFS - only one bottom stress formulation is', &
2951	& ' allowed.')
2952	exit_flag=5
2953	END IF
2954	!
2955	! Stop if no bottom stress formulation is selected.
2956	!
2957	IF (Master.and.(ibbl.eq.0)) THEN
2958	WRITE (stdout,50)
2959	50 FORMAT (/,' CHECKDEFS - no bottom stress formulation is', &
2960	& ' selected.')
2961	exit_flag=5
2962	END IF
2963	!
2964	! Stop if more than one biological module is selected.
2965	!
2966	IF (Master.and.(ibiology.gt.1)) THEN
2967	WRITE (stdout,60)
2968	60 FORMAT (/,' CHECKDEFS - only one biology MODULE is allowed.')
2969	exit_flag=5
2970	END IF
2971	!
2972	! Stop if more that one model driver is selected.
2973	!
2974	IF (Master.and.(idriver.gt.1)) THEN
2975	WRITE (stdout,70)
2976	70 FORMAT (/,' CHECKDEFS - only one model example is allowed.')
2977	exit_flag=5
2978	END IF
2979
2980	#ifndef SOLVE3D
2981	!
2982	! Stop it explicit time-step splitting on shallow water set-up.
2983	!
2984	DO ng=1,Ngrids
2985	IF (Master.and.(ndtfast(ng).gt.1)) THEN
2986	WRITE (stdout,80)
2987	80 FORMAT (/,' CHECKDEFS - explicit time-step splitting is ', &
2988	& ' inconsistent.', &
2989	& /,13x,'Change parameter NDTFAST to unity.')
2990	exit_flag=5
2991	END IF
2992	END DO
2993	#endif
2994	#if defined ADJOINT && defined _OPENMP
2995	IF (Master) THEN
2996	WRITE (stdout,90)
2997	90 FORMAT (/,' CHECKDEFS - cannot run the adjoint model in', &
2998	& ' shared-memory mode.')
2999	exit_flag=5
3000	END IF
3001	#endif
3002	#if defined NEMURO && defined HOLLING_GRAZING && defined IVLEV_EXPLICIT
3003	!
3004	! Stop if using more than one grazing algorithm.
3005	!
3006	IF (Master) THEN
3007	WRITE (stdout,100)
3008	100 FORMAT (/,' CHECKDEFS - More than one grazing algorithm', &
3009	& ' selected for Nemuro model.')
3010	exit_flag=5
3011	END IF
3012	#endif
3013	#if defined FOUR_DVAR && defined VCONVOLUTION
3014	# if defined IMPLICIT_VCONV && defined SPLINES_VCONV
3015	!
3016	! Stop if using more than vertical convolution algorithm.
3017	!
3018	IF (Master) THEN
3019	WRITE (stdout,110)
3020	110 FORMAT (/,' CHECKDEFS - More than one vertical convolution', &
3021	& ' algorithm selected.')
3022	exit_flag=5
3023	END IF
3024	# endif
3025	#endif
3026	#if defined MODEL_COUPLING && (defined ESMF_LIB && defined MCT_LIB)
3027	IF (Master) THEN
3028	WRITE (stdout,120)
3029	120 FORMAT (/,' CHECKDEFS - More than one coupling library', &
3030	& ' selected.')
3031	exit_flag=5
3032	END IF
3033	#endif
3034	#if defined POSTERIOR_ERROR_F && defined POSTERIOR_ERROR_I && \
3035	defined WEAK_CONSTRAINT
3036	IF (Master) THEN
3037	WRITE (stdout,130)
3038	130 FORMAT (/,' CHECKDEFS - Both initial and final posterior', &
3039	& ' error covariance analysis selected, choose one.')
3040	exit_flag=5
3041	END IF
3042	#endif
3043	#ifdef UV_ADV
3044	!
3045	! Stop if more than one advection scheme has been activated.
3046	!
3047	IF (Master.and.(ivelHadv.gt.1)) THEN
3048	WRITE (stdout,140) 'horizontal','momentum','ivelHadv =',ivelHadv
3049	exit_flag=5
3050	END IF
3051	# ifdef SOLVE3D
3052	IF (Master.and.(ivelVadv.gt.1)) THEN
3053	WRITE (stdout,140) 'vertical','momentum','ivelVadv =',ivelVadv
3054	exit_flag=5
3055	END IF
3056	# endif
3057	#endif
3058	#ifdef SOLVE3D
3059	IF (Master.and.(itrcHadv.gt.1)) THEN
3060	WRITE (stdout,140) 'horizontal','tracers','itrcHadv =',itrcHadv
3061	exit_flag=5
3062	END IF
3063	IF (Master.and.(itrcVadv.gt.1)) THEN
3064	WRITE (stdout,140) 'vertical','tracers','itrcVadv =',itrcVadv
3065	exit_flag=5
3066	END IF
3067	# if defined TANGENT \|\| defined TL_IOMS \|\| defined ADJOINT
3068	IF (Master.and.(itrcHadvtl.gt.1)) THEN
3069	WRITE (stdout,140) 'TL/AD horizontal','tracers','itrcHadvtl =', &
3070	& itrcHadvtl
3071	exit_flag=5
3072	END IF
3073	IF (Master.and.(itrcVadvtl.gt.1)) THEN
3074	WRITE (stdout,140) 'TL/AD vertical','tracers','itrcVadvtl =', &
3075	& itrcVadvtl
3076	exit_flag=5
3077	END IF
3078	# endif
3079	#endif
3080	140 FORMAT (/,' CHECKDEFS - only one ',a,' advection scheme', &
3081	& /,13x,'is allowed for ',a,', ',a,1x,i1)
3082	#if defined SOLVE3D && defined NEMURO
3083	# if !(defined IVLEV_EXPLICIT \|\| defined HOLLING_GRAZING)
3084	IF (Master) THEN
3085	WRITE (stdout,150) uppercase('ivlev_explicit'), &
3086	& uppercase('holling_grazing')
3087	150 FORMAT (/,' CHECKDEFS - Need to choose a zooplankton grazing', &
3088	& ' option;',/,13x,'use either ',a,' or implicit ',a, &
3089	& '.',/,13x,'The default implicit IVLEV algorithm', &
3090	& ' does not work well yet.')
3091	exit_flag=5
3092	END IF
3093	# endif
3094	#endif
3095	#ifdef NEARSHORE
3096	!
3097	! Stop if more that one radiation stress formulation is activated.
3098	!
3099	IF (Master.and.(nearshore.gt.1)) THEN
3100	WRITE (stdout,160)
3101	160 FORMAT (/,' CHECKDEFS - only one radiation stress formulation' &
3102	& ' is allowed.')
3103	exit_flag=5
3104	END IF
3105	#endif
3106	#ifdef PARALLEL_IO
3107	# if !(defined HDF5 \|\| defined PNETCDF)
3108	IF (Master) THEN
3109	WRITE (stdout,170) uppercase('hdf5'), &
3110	& uppercase('pnetcdf')
3111	170 FORMAT (/,' CHECKDEFS - Need to activate either ',a,' or ',a, &
3112	& ' options',/,13x,'for parallel I/O processing.')
3113	exit_flag=5
3114	END IF
3115	# endif
3116	# if defined HDF5 && defined PNETCDF
3117	IF (Master) THEN
3118	WRITE (stdout,180) uppercase('hdf5'), &
3119	& uppercase('pnetcdf')
3120	180 FORMAT (/,' CHECKDEFS - Both ',a,' and ',a,' options', &
3121	& ' selected',/,13x,'for parallel I/O processing.', &
3122	& /,13x,'Choose only one.')
3123	exit_flag=5
3124	END IF
3125	# endif
3126	#endif
3127	#if defined WRITE_WATER && defined WET_DRY
3128	!
3129	! Stop if writting only water points when wetting and drying is
3130	! activated. This is not possible since the Land/Sea masking is
3131	! time dependent.
3132	!
3133	IF (Master) THEN
3134	WRITE (stdout,190) uppercase('write_water'), &
3135	& uppercase('wet_dry')
3136	190 FORMAT (/,' CHECKDEFS - cannot activate ',a,' and ',a, &
3137	& ' together',/,13x,'because of time dependent ', &
3138	& ' Land/Sea masking.')
3139	exit_flag=5
3140	END IF
3141	#endif
3142	#if ((defined AD_AVERAGES && defined ADJOINT) && defined AVERAGES) \|\| \
3143	((defined RP_AVERAGES && defined TL_IOMS) && defined AVERAGES) \|\| \
3144	((defined TL_AVERAGES && defined TANGENT) && defined AVERAGES) \|\| \
3145	(defined AD_AVERAGES && defined TL_AVERAGES) \|\| \
3146	(defined AD_AVERAGES && defined RP_AVERAGES) \|\| \
3147	(defined RP_AVERAGES && defined TL_AVERAGES)
3148	!
3149	! Stop if activating more that one time-averaged option. This is not
3150	! possible in the same executable because the same internal arrays
3151	! are used to store the cummulative sums.
3152	!
3153	IF (Master) THEN
3154	WRITE (stdout,200) uppercase('averages'), &
3155	& uppercase('ad_averages'), &
3156	& uppercase('rp_averages'), &
3157	& uppercase('tl_averages')
3158	200 FORMAT (/,' CHECKDEFS - cannot activate ',a,' ',a,' ',a, &
3159	& ' or',a,/,13x,'at the same time. Choose one!')
3160	exit_flag=5
3161	END IF
3162	#endif
3163	#if !defined DISTRIBUTE && defined TS_MPDATA
3164	!
3165	! Stop if activating MPDATA in serial with partitions or shared-
3166	! memory.
3167	!
3168	IF (Master) THEN
3169	DO ng=1,Ngrids
3170	IF (NtileX(ng)*NtileE(ng).gt.1) THEN
3171	WRITE (stdout,210) uppercase('ts_mpdata')
3172	exit_flag=5
3173	EXIT
3174	END IF
3175	END DO
3176	210 FORMAT (/,' CHECKDEFS - cannot activate option: ',a, &
3177	& /,13x,'in serial with partitions or shared-memory...', &
3178	& /,13x,'Use distributed-memory (MPI) in parallel runs.')
3179	END IF
3180	#endif
3181	#ifdef Q_PSOURCE
3182	!
3183	! Stop if activating obsolete mass point Sources/Sinks option. This
3184	! capability is now activated with standard input switch "LwSrc".
3185	!
3186	IF (Master) THEN
3187	WRITE (stdout,220) uppercase('q_psource'), 'LwSrc'
3188	220 FORMAT (/,' CHECKDEFS - cannot use obsolete option: ',a, &
3189	& /,13x,'Use instead standard input switch: ',a, &
3190	& /,13x,'Edit header file or build script!')
3191	exit_flag=5
3192	END IF
3193	#endif
3194	#if defined TS_PSOURCE && defined SOLVE3D
3195	!
3196	! Stop if activating obsolete tracers point Sources/Sinks option. This
3197	! capability is now activated with standard input switch "LtracerSrc".
3198	!
3199	IF (Master) THEN
3200	WRITE (stdout,230) uppercase('ts_psource'), 'LtracerSrc'
3201	230 FORMAT (/,' CHECKDEFS - cannot use obsolete option: ',a, &
3202	& /,13x,'Use instead standard input switch: ',a, &
3203	& /,13x,'Edit header file or build script!')
3204	exit_flag=5
3205	END IF
3206	#endif
3207	#ifdef UV_PSOURCE
3208	!
3209	! Stop if activating obsolete mass point Sources/Sinks option. This
3210	! capability is now activated with standard input switch "LwSrc".
3211	!
3212	IF (Master) THEN
3213	WRITE (stdout,240) uppercase('uv_psource'), 'LuvSrc'
3214	240 FORMAT (/,' CHECKDEFS - cannot use obsolete option: ',a, &
3215	& /,13x,'Use instead standard input switch: ',a, &
3216	& /,13x,'Edit header file or build script!')
3217	exit_flag=5
3218	END IF
3219	#endif
3220	#if defined SPLINES && defined SOLVE3D
3221	!
3222	! Stop if activating parabolic splines reconstruction with deprecated
3223	! option. See https://www.myroms.org/projects/src/ticket/681 for more
3224	! details.
3225	!
3226	IF (Master) THEN
3227	WRITE (stdout,250) uppercase('splines'), &
3228	& 'https://www.myroms.org/projects/src/ticket/681'
3229	250 FORMAT (/,' CHECKDEFS - cannot use obsolete option: ',a, &
3230	& /,13x,'See following trac ticket for details: ', &
3231	& /,13x,a,/,13x,'please change header file accordingly.')
3232	exit_flag=5
3233	END IF
3234	#endif
3235	#if !defined DISTRIBUTE && defined NESTING
3236	!
3237	! Stop if activating nesting in serial with partitions or shared-
3238	! memory.
3239	!
3240	IF (Master) THEN
3241	DO ng=1,Ngrids
3242	IF (NtileX(ng)*NtileE(ng).gt.1) THEN
3243	WRITE (stdout,260) uppercase('nesting')
3244	exit_flag=5
3245	EXIT
3246	END IF
3247	END DO
3248	260 FORMAT (/,' CHECKDEFS - cannot activate option: ',a, &
3249	& /,13x,'in serial with partitions or shared-memory...', &
3250	& /,13x,'There bugs in shared-memory that need fixing.', &
3251	& /,13x,'Use distributed-memory (MPI) in parallel runs.')
3252	END IF
3253	#endif
3254	#ifdef ATM_COUPLING
3255	!
3256	! Stop if more than one atmosphere model is selected for coupling.
3257	!
3258	IF (Master.and.(iatms.gt.1)) THEN
3259	WRITE (stdout,270)
3260	270 FORMAT (/,' CHECKDEFS - only one atmosphere model coupling', &
3261	& ' is allowed.')
3262	exit_flag=5
3263	END IF
3264	#endif
3265
3266	RETURN
3267	END SUBROUTINE checkdefs

Download in other formats:

Original Format