model/src/packages_init_fixed.F

C $Header: /u/gcmpack/MITgcm/model/src/packages_init_fixed.F,v 1.66 2009/01/04 01:00:19 jmc Exp $
C $Name:  $

#include "PACKAGES_CONFIG.h"
#include "CPP_OPTIONS.h"

CBOP
C     !ROUTINE: PACKAGES_INIT_FIXED
C     !INTERFACE:
      SUBROUTINE PACKAGES_INIT_FIXED( myThid )

C     !DESCRIPTION: \bv
C     *==========================================================*
C     | SUBROUTINE PACKAGES_INIT_FIXED
C     | o Does initialisation of package-related fixed fields
C     *==========================================================*
C     \ev

C     !CALLING SEQUENCE:
C     PACKAGES_INIT_FIXED
C       |
C       |-- RUNCLOCK_INIT
C       |
C       |-- DIAGNOSTICS_INIT_EARLY
C       |-- DIAGNOSTICS_MAIN_INIT
C       |
C       |-- GAD_INIT
C       |
C       |-- MOM_INIT_FIXED
C       |
C       |-- CD_CODE_INIT_FIXED
C       |
C       |-- EXF_INIT_FIXED
C       |
C       |-- GMREDI_INIT_FIXED
C       |
C       |-- DWNSLP_INIT_FIXED
C       |
C       |-- KPP_INIT_FIXED
C       |
C       |-- MY82_INIT_FIXED
C       |
C       |-- GGL90_INIT_FIXED
C       |
C       |-- OPPS_INIT
C       |
C       |-- OBCS_INIT_FIXED
C       |
C       |-- TIMEAVE_INIT_FIXED
C       |
C       |-- SHAP_FILT_INIT_FIZED
C       |
C       |-- ZONAL_FILT_INIT
C       |
C       |-- AIM_INITIALISE
C       |
C       |-- GRIDALT_INITIALISE
C       |
C       |-- FIZHI_INIT_FIXED
C       |
C       |-- LAND_INITIALISE
C       |
C       |-- SEAICE_COST_INIT_FIXED
C       |
C       |-- ECCO_COST_INIT_FIXED
C       |
C       |-- PROFILES_INIT_FIXED
C       |
C       |-- CTRL_INIT
C       |
C       |-- FLT_INIT_FIXED
C       |
C       |-- PTRACERS_INIT_FIXED
C       |
C       |-- GCHEM_INIT_FIXED
C       |
C       |-- THSICE_INIT_FIXED
C       |
C       |-- SHELFICE_INIT_FIXED
C       |
C       |-- CPL_INIT_FIXED
C       |
C       |-- ATM2D_INIT_FIXED
C       |
C       |-- REGRID_INIT_FIXED
C       |
C       |-- SALT_PLUME_INIT_FIXED
C       |
C       |-- MYPACKAGE_INIT_FIXED
C       |
C       |-- DIAGNOSTICS_INIT_FIXED

C     !USES:
      IMPLICIT NONE
C     === Global variables ===
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"

C     !INPUT/OUTPUT PARAMETERS:
C     === Routine arguments ===
C     myThid -  Number of this instances
      INTEGER myThid
CEOP

#ifdef ALLOW_RUNCLOCK
      IF ( useRunClock ) THEN
        CALL RUNCLOCK_INIT( myThid )
      ENDIF
#endif

#ifdef ALLOW_DIAGNOSTICS
      IF ( useDiagnostics ) THEN
C-    needs to call DIAGNOSTICS_INIT_EARLY before all package-diag-init calls
        CALL DIAGNOSTICS_INIT_EARLY( myThid )
        CALL DIAGNOSTICS_MAIN_INIT( myThid )
      ENDIF
#endif

#ifdef ALLOW_GENERIC_ADVDIFF
C--   Initialize fixed params for GAD
      CALL GAD_INIT_FIXED( myThid )
#endif

#ifdef ALLOW_MOM_COMMON
C--   Initialize fixed params for Momentum pkgs (common, fluxform, vecinv)
      IF ( momStepping ) THEN
       CALL MOM_INIT_FIXED( myThid )
      ENDIF
#endif

#ifdef ALLOW_CD_CODE
      IF (useCDscheme) THEN
       CALL CD_CODE_INIT_FIXED(myThid)
      ENDIF
#endif

#ifdef ALLOW_EXF
C--   Initialize fixed arrays for EXF
      IF ( useEXF ) THEN
       CALL EXF_INIT_FIXED( myThid )
      ENDIF
#endif

#ifdef ALLOW_GMREDI
C--   Initialize fixed arrays for GM-Redi
      IF ( useGMRedi ) THEN
       CALL GMREDI_INIT_FIXED( myThid )
      ENDIF
#endif

#ifdef ALLOW_DOWN_SLOPE
C--   Initialize fixed arrays for Down-Slope pkg
      IF ( useDOWN_SLOPE ) THEN
       CALL DWNSLP_INIT_FIXED( myThid )
      ENDIF
#endif

#ifdef ALLOW_KPP
C--   Initialize fixed arrays for KPP
      IF ( useKPP ) THEN
       CALL KPP_INIT_FIXED( myThid )
      ENDIF
#endif

#ifdef ALLOW_MY82
C--   Initialize fixed arrays for MY82
      IF ( useMY82 ) THEN
       CALL MY82_INIT_FIXED( myThid )
      ENDIF
#endif

#ifdef ALLOW_GGL90
C--   Initialize fixed arrays for GGL90
      IF ( useGGL90 ) THEN
       CALL GGL90_INIT_FIXED( myThid )
      ENDIF
#endif

#ifdef ALLOW_OPPS
      IF (useOPPS) CALL OPPS_INIT( myThid )
#endif /* ALLOW_OPPS */

C--   Initialize fixed arrays for OBCS
#ifdef ALLOW_OBCS
      IF (useOBCS) THEN
       CALL OBCS_INIT_FIXED( myThid )
      ENDIF
#endif

#ifdef ALLOW_TIMEAVE
C     IF (useTIMEAVE) THEN
      IF ( taveFreq.GT.0. ) THEN
       CALL TIMEAVE_INIT_FIXED( myThid )
      ENDIF
#endif

#ifdef ALLOW_SHAP_FILT
C--   Shapiro filter initialisation
      IF (useSHAP_FILT) THEN
       CALL SHAP_FILT_INIT_FIXED( myThid )
      ENDIF
#endif

#ifdef ALLOW_ZONAL_FILT
C--   Latitude circle filter initialisation
      IF (useZONAL_FILT) THEN
       CALL ZONAL_FILT_INIT(myThid)
      ENDIF
#endif

#ifdef ALLOW_AIM
C--   Initialise & Read AIM physical parameters
      IF (useAIM) CALL AIM_INITIALISE( myThid )
#endif

C AMM
#ifdef ALLOW_GRIDALT
C--   Initialise GRIDALT parameters - the alternative grid
      IF (useGRIDALT) then
       CALL TIMER_START('GRIDALT_INITIALISE  [PACKAGES_INIT_F]',myThid)
       CALL GRIDALT_INITIALISE( myThid )
       CALL TIMER_STOP ('GRIDALT_INITIALISE  [PACKAGES_INIT_F]',myThid)
      ENDIF

#endif

#ifdef ALLOW_FIZHI
C--   Initialise & Read FIZHI physical parameters
      IF (useFIZHI) CALL FIZHI_INIT_FIXED( myThid )
#endif
C AMM

#ifdef ALLOW_LAND
C--   Initialise & Read Land package parameters
      IF (useLand) CALL LAND_INITIALISE( myThid )
#endif

#if (defined (ALLOW_SEAICE) && defined (ALLOW_COST))
C--   Initialise ecco-specific cost function.
C--   This needs to preceed the call ctrl_init
C--   in order to provide the weight files
      IF (useSEAICE) CALL SEAICE_COST_INIT_FIXED( myThid )
#endif

#if (defined (ALLOW_ECCO) && defined (ALLOW_COST))
C--   Initialise ecco-specific cost function.
C--   This needs to preceed the call ctrl_init
C--   in order to provide the weight files
      CALL ECCO_COST_INIT_FIXED( myThid )
#endif

#ifdef ALLOW_PROFILES
      CALL PROFILES_INIT_FIXED( myThid )
#endif

#ifdef ALLOW_AUTODIFF
C--   Initialise the control variables
      CALL CTRL_INIT( myThid )
#endif

#ifdef ALLOW_ADMTLM
C--   Initialise ADMTLM
      CALL CTRL_ADMTLM( myThid )
#endif

#ifdef ALLOW_FLT
C--   Initialise fixed array for Float pkg
      IF ( useFLT ) CALL FLT_INIT_FIXED( myThid )
#endif

#ifdef ALLOW_PTRACERS
      IF (usePTRACERS) CALL PTRACERS_INIT_FIXED(myThid)
#endif

#ifdef ALLOW_GCHEM
      IF (useGCHEM) CALL GCHEM_INIT_FIXED(myThid)
#endif

#ifdef ALLOW_RBCS
      IF (useRBCS) CALL RBCS_INIT_FIXED(myThid)
#endif

#ifdef ALLOW_THSICE
      IF (useThSIce) CALL THSICE_INIT_FIXED(myThid)
#endif

#ifdef ALLOW_SEAICE
      IF (useSEAICE) CALL SEAICE_INIT_FIXED(myThid)
#endif

#ifdef ALLOW_SHELFICE
      IF (useShelfIce) CALL SHELFICE_INIT_FIXED( myThid )
#endif /* ALLOW_SHELFICE */

#ifdef ALLOW_EMBED_FILES
      IF ( useEMBED_FILES ) CALL EMBED_FILES_INIT( myThid )
#endif

#ifdef ALLOW_REGRID
      IF ( useREGRID ) CALL REGRID_INIT_FIXED( myThid )
#endif

#ifdef COMPONENT_MODULE
      IF ( useCoupler ) CALL CPL_INIT_FIXED( myThid )
#endif

#ifdef ALLOW_ATM2D
      IF ( useAtm2d ) CALL ATM2D_INIT_FIXED( myThid )
#endif

#ifdef ALLOW_SALT_PLUME
      IF (useSALT_PLUME) CALL SALT_PLUME_INIT_FIXED(myThid)
#endif

#ifdef ALLOW_MYPACKAGE
      IF (useMYPACKAGE) CALL MYPACKAGE_INIT_FIXED(myThid)
#endif

#ifdef ALLOW_DIAGNOSTICS
C-    needs to call DIAGNOSTICS_INIT_FIXED after all package-diag-init calls
      IF ( useDiagnostics ) CALL DIAGNOSTICS_INIT_FIXED( myThid )
#endif

      RETURN
      END
1	jmc	1.67	C $Header: /u/gcmpack/MITgcm/model/src/packages_init_fixed.F,v 1.66 2009/01/04 01:00:19 jmc Exp $
2	adcroft	1.3	C $Name: $
3	adcroft	1.2
4	edhill	1.11	#include "PACKAGES_CONFIG.h"
5	adcroft	1.2	#include "CPP_OPTIONS.h"
6
7	cnh	1.4	CBOP
8			C !ROUTINE: PACKAGES_INIT_FIXED
9			C !INTERFACE:
10	adcroft	1.2	SUBROUTINE PACKAGES_INIT_FIXED( myThid )
11	cnh	1.4
12			C !DESCRIPTION: \bv
13			C ==========================================================
14	jmc	1.51	C \| SUBROUTINE PACKAGES_INIT_FIXED
15			C \| o Does initialisation of package-related fixed fields
16	cnh	1.4	C ==========================================================
17			C \ev
18
19			C !CALLING SEQUENCE:
20			C PACKAGES_INIT_FIXED
21			C \|
22	adcroft	1.39	C \|-- RUNCLOCK_INIT
23			C \|
24	jmc	1.34	C \|-- DIAGNOSTICS_INIT_EARLY
25			C \|-- DIAGNOSTICS_MAIN_INIT
26			C \|
27	jmc	1.5	C \|-- GAD_INIT
28			C \|
29	jmc	1.44	C \|-- MOM_INIT_FIXED
30			C \|
31	jmc	1.35	C \|-- CD_CODE_INIT_FIXED
32			C \|
33	dimitri	1.38	C \|-- EXF_INIT_FIXED
34			C \|
35	jmc	1.35	C \|-- GMREDI_INIT_FIXED
36			C \|
37	jmc	1.64	C \|-- DWNSLP_INIT_FIXED
38			C \|
39	dimitri	1.37	C \|-- KPP_INIT_FIXED
40			C \|
41	mlosch	1.63	C \|-- MY82_INIT_FIXED
42			C \|
43	dfer	1.65	C \|-- GGL90_INIT_FIXED
44			C \|
45	jmc	1.35	C \|-- OPPS_INIT
46			C \|
47	cnh	1.4	C \|-- OBCS_INIT_FIXED
48			C \|
49	jmc	1.34	C \|-- TIMEAVE_INIT_FIXED
50			C \|
51	jmc	1.36	C \|-- SHAP_FILT_INIT_FIZED
52			C \|
53	jmc	1.6	C \|-- ZONAL_FILT_INIT
54	jmc	1.7	C \|
55			C \|-- AIM_INITIALISE
56			C \|
57	jmc	1.25	C \|-- GRIDALT_INITIALISE
58			C \|
59			C \|-- FIZHI_INIT_FIXED
60			C \|
61	jmc	1.7	C \|-- LAND_INITIALISE
62	jmc	1.25	C \|
63	jmc	1.44	C \|-- SEAICE_COST_INIT_FIXED
64			C \|
65	jmc	1.25	C \|-- ECCO_COST_INIT_FIXED
66			C \|
67	jmc	1.67	C \|-- PROFILES_INIT_FIXED
68			C \|
69	jmc	1.25	C \|-- CTRL_INIT
70			C \|
71	jmc	1.67	C \|-- FLT_INIT_FIXED
72			C \|
73	jmc	1.49	C \|-- PTRACERS_INIT_FIXED
74	jmc	1.34	C \|
75	jmc	1.25	C \|-- GCHEM_INIT_FIXED
76			C \|
77	jmc	1.35	C \|-- THSICE_INIT_FIXED
78	mlosch	1.27	C \|
79	mlosch	1.47	C \|-- SHELFICE_INIT_FIXED
80			C \|
81	jmc	1.34	C \|-- CPL_INIT_FIXED
82			C \|
83	jscott	1.54	C \|-- ATM2D_INIT_FIXED
84			C \|
85	dimitri	1.60	C \|-- REGRID_INIT_FIXED
86			C \|
87			C \|-- SALT_PLUME_INIT_FIXED
88			C \|
89	jmc	1.64	C \|-- MYPACKAGE_INIT_FIXED
90			C \|
91	jmc	1.34	C \|-- DIAGNOSTICS_INIT_FIXED
92	cnh	1.4
93			C !USES:
94	adcroft	1.2	IMPLICIT NONE
95			C === Global variables ===
96			#include "SIZE.h"
97			#include "EEPARAMS.h"
98			#include "PARAMS.h"
99
100	cnh	1.4	C !INPUT/OUTPUT PARAMETERS:
101	adcroft	1.2	C === Routine arguments ===
102			C myThid - Number of this instances
103			INTEGER myThid
104	cnh	1.4	CEOP
105	edhill	1.21
106	adcroft	1.39	#ifdef ALLOW_RUNCLOCK
107			IF ( useRunClock ) THEN
108			CALL RUNCLOCK_INIT( myThid )
109			ENDIF
110			#endif
111
112	jmc	1.33	#ifdef ALLOW_DIAGNOSTICS
113			IF ( useDiagnostics ) THEN
114	jmc	1.51	C- needs to call DIAGNOSTICS_INIT_EARLY before all package-diag-init calls
115	jmc	1.33	CALL DIAGNOSTICS_INIT_EARLY( myThid )
116			CALL DIAGNOSTICS_MAIN_INIT( myThid )
117			ENDIF
118			#endif
119
120	jmc	1.32	#ifdef ALLOW_GENERIC_ADVDIFF
121	edhill	1.21	C-- Initialize fixed params for GAD
122	jmc	1.57	CALL GAD_INIT_FIXED( myThid )
123	jmc	1.32	#endif
124	adcroft	1.2
125	jmc	1.44	#ifdef ALLOW_MOM_COMMON
126			C-- Initialize fixed params for Momentum pkgs (common, fluxform, vecinv)
127			IF ( momStepping ) THEN
128			CALL MOM_INIT_FIXED( myThid )
129			ENDIF
130			#endif
131
132	jmc	1.35	#ifdef ALLOW_CD_CODE
133			IF (useCDscheme) THEN
134	jmc	1.62	CALL CD_CODE_INIT_FIXED(myThid)
135	jmc	1.35	ENDIF
136			#endif
137
138	dimitri	1.38	#ifdef ALLOW_EXF
139			C-- Initialize fixed arrays for EXF
140	heimbach	1.58	IF ( useEXF ) THEN
141	dimitri	1.38	CALL EXF_INIT_FIXED( myThid )
142	heimbach	1.58	ENDIF
143	dimitri	1.38	#endif
144
145	jmc	1.35	#ifdef ALLOW_GMREDI
146			C-- Initialize fixed arrays for GM-Redi
147			IF ( useGMRedi ) THEN
148			CALL GMREDI_INIT_FIXED( myThid )
149			ENDIF
150			#endif
151
152	jmc	1.64	#ifdef ALLOW_DOWN_SLOPE
153			C-- Initialize fixed arrays for Down-Slope pkg
154			IF ( useDOWN_SLOPE ) THEN
155			CALL DWNSLP_INIT_FIXED( myThid )
156			ENDIF
157			#endif
158
159	dimitri	1.37	#ifdef ALLOW_KPP
160			C-- Initialize fixed arrays for KPP
161			IF ( useKPP ) THEN
162			CALL KPP_INIT_FIXED( myThid )
163			ENDIF
164			#endif
165
166	mlosch	1.63	#ifdef ALLOW_MY82
167			C-- Initialize fixed arrays for MY82
168			IF ( useMY82 ) THEN
169			CALL MY82_INIT_FIXED( myThid )
170			ENDIF
171			#endif
172
173	dfer	1.65	#ifdef ALLOW_GGL90
174			C-- Initialize fixed arrays for GGL90
175			IF ( useGGL90 ) THEN
176			CALL GGL90_INIT_FIXED( myThid )
177			ENDIF
178			#endif
179
180	jmc	1.35	#ifdef ALLOW_OPPS
181			IF (useOPPS) CALL OPPS_INIT( myThid )
182			#endif /* ALLOW_OPPS */
183
184	adcroft	1.2	C-- Initialize fixed arrays for OBCS
185			#ifdef ALLOW_OBCS
186			IF (useOBCS) THEN
187			CALL OBCS_INIT_FIXED( myThid )
188	adcroft	1.3	ENDIF
189			#endif
190	edhill	1.26
191			#ifdef ALLOW_TIMEAVE
192			C IF (useTIMEAVE) THEN
193	jmc	1.36	IF ( taveFreq.GT.0. ) THEN
194			CALL TIMEAVE_INIT_FIXED( myThid )
195			ENDIF
196			#endif
197
198			#ifdef ALLOW_SHAP_FILT
199			C-- Shapiro filter initialisation
200			IF (useSHAP_FILT) THEN
201			CALL SHAP_FILT_INIT_FIXED( myThid )
202			ENDIF
203	edhill	1.26	#endif
204
205	jmc	1.6	#ifdef ALLOW_ZONAL_FILT
206			C-- Latitude circle filter initialisation
207			IF (useZONAL_FILT) THEN
208			CALL ZONAL_FILT_INIT(myThid)
209	adcroft	1.2	ENDIF
210	jmc	1.7	#endif
211
212			#ifdef ALLOW_AIM
213			C-- Initialise & Read AIM physical parameters
214			IF (useAIM) CALL AIM_INITIALISE( myThid )
215			#endif
216	molod	1.15
217			C AMM
218			#ifdef ALLOW_GRIDALT
219			C-- Initialise GRIDALT parameters - the alternative grid
220	molod	1.24	IF (useGRIDALT) then
221	jmc	1.62	CALL TIMER_START('GRIDALT_INITIALISE [PACKAGES_INIT_F]',myThid)
222	molod	1.24	CALL GRIDALT_INITIALISE( myThid )
223	jmc	1.62	CALL TIMER_STOP ('GRIDALT_INITIALISE [PACKAGES_INIT_F]',myThid)
224	molod	1.24	ENDIF
225
226	molod	1.15	#endif
227
228			#ifdef ALLOW_FIZHI
229			C-- Initialise & Read FIZHI physical parameters
230			IF (useFIZHI) CALL FIZHI_INIT_FIXED( myThid )
231	molod	1.16	#endif
232			C AMM
233
234	jmc	1.7	#ifdef ALLOW_LAND
235			C-- Initialise & Read Land package parameters
236			IF (useLand) CALL LAND_INITIALISE( myThid )
237	heimbach	1.8	#endif
238
239	heimbach	1.43	#if (defined (ALLOW_SEAICE) && defined (ALLOW_COST))
240			C-- Initialise ecco-specific cost function.
241			C-- This needs to preceed the call ctrl_init
242			C-- in order to provide the weight files
243	heimbach	1.56	IF (useSEAICE) CALL SEAICE_COST_INIT_FIXED( myThid )
244	heimbach	1.43	#endif
245
246	heimbach	1.14	#if (defined (ALLOW_ECCO) && defined (ALLOW_COST))
247			C-- Initialise ecco-specific cost function.
248			C-- This needs to preceed the call ctrl_init
249			C-- in order to provide the weight files
250	heimbach	1.40	CALL ECCO_COST_INIT_FIXED( myThid )
251	heimbach	1.14	#endif
252
253	heimbach	1.48	#ifdef ALLOW_PROFILES
254			CALL PROFILES_INIT_FIXED( myThid )
255			#endif
256
257	heimbach	1.13	#ifdef ALLOW_AUTODIFF
258	heimbach	1.8	C-- Initialise the control variables
259			CALL CTRL_INIT( myThid )
260	adcroft	1.2	#endif
261	stephd	1.10
262	heimbach	1.59	#ifdef ALLOW_ADMTLM
263			C-- Initialise ADMTLM
264			CALL CTRL_ADMTLM( myThid )
265	jmc	1.67	#endif
266
267			#ifdef ALLOW_FLT
268			C-- Initialise fixed array for Float pkg
269			IF ( useFLT ) CALL FLT_INIT_FIXED( myThid )
270	heimbach	1.59	#endif
271
272	stephd	1.10	#ifdef ALLOW_PTRACERS
273	jmc	1.49	IF (usePTRACERS) CALL PTRACERS_INIT_FIXED(myThid)
274	jmc	1.25	#endif
275
276	stephd	1.10	#ifdef ALLOW_GCHEM
277	jmc	1.62	IF (useGCHEM) CALL GCHEM_INIT_FIXED(myThid)
278	edhill	1.42	#endif
279
280	stephd	1.45	#ifdef ALLOW_RBCS
281	jmc	1.62	IF (useRBCS) CALL RBCS_INIT_FIXED(myThid)
282	stephd	1.45	#endif
283
284	jmc	1.35	#ifdef ALLOW_THSICE
285	jmc	1.62	IF (useThSIce) CALL THSICE_INIT_FIXED(myThid)
286	heimbach	1.28	#endif
287
288	heimbach	1.55	#ifdef ALLOW_SEAICE
289	jmc	1.62	IF (useSEAICE) CALL SEAICE_INIT_FIXED(myThid)
290	heimbach	1.55	#endif
291
292	mlosch	1.47	#ifdef ALLOW_SHELFICE
293			IF (useShelfIce) CALL SHELFICE_INIT_FIXED( myThid )
294			#endif /* ALLOW_SHELFICE */
295
296	jmc	1.53	#ifdef ALLOW_EMBED_FILES
297	jmc	1.62	IF ( useEMBED_FILES ) CALL EMBED_FILES_INIT( myThid )
298	jmc	1.53	#endif
299
300			#ifdef ALLOW_REGRID
301	jmc	1.62	IF ( useREGRID ) CALL REGRID_INIT_FIXED( myThid )
302	jmc	1.53	#endif
303
304	jmc	1.33	#ifdef COMPONENT_MODULE
305	jmc	1.62	IF ( useCoupler ) CALL CPL_INIT_FIXED( myThid )
306	jmc	1.33	#endif
307
308	jscott	1.54	#ifdef ALLOW_ATM2D
309			IF ( useAtm2d ) CALL ATM2D_INIT_FIXED( myThid )
310			#endif
311
312	dimitri	1.60	#ifdef ALLOW_SALT_PLUME
313	jmc	1.62	IF (useSALT_PLUME) CALL SALT_PLUME_INIT_FIXED(myThid)
314	dimitri	1.60	#endif
315
316	jmc	1.64	#ifdef ALLOW_MYPACKAGE
317			IF (useMYPACKAGE) CALL MYPACKAGE_INIT_FIXED(myThid)
318			#endif
319
320	molod	1.18	#ifdef ALLOW_DIAGNOSTICS
321	jmc	1.51	C- needs to call DIAGNOSTICS_INIT_FIXED after all package-diag-init calls
322	heimbach	1.50	IF ( useDiagnostics ) CALL DIAGNOSTICS_INIT_FIXED( myThid )
323	molod	1.18	#endif
324	edhill	1.30
325	adcroft	1.2	RETURN
326			END