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

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.

bdydta.F90 in branches/2014/dev_r4621_NOC4_BDY_VERT_INTERP/NEMOGCM/NEMO/OPA_SRC/BDY – NEMO

Context Navigation

source: branches/2014/dev_r4621_NOC4_BDY_VERT_INTERP/NEMOGCM/NEMO/OPA_SRC/BDY/bdydta.F90 @ 5901

Last change on this file since 5901 was 5620, checked in by jamesharle, 9 years ago
Merge with r5619 of trunk, update to unstructured BDY interpolation in fldread.F90. Structured BDY interpolation incomplete.
Property svn:keywords set to `Id`
File size: 40.5 KB

Line
1	MODULE bdydta
2	!!======================================================================
3	!! * MODULE bdydta *
4	!! Open boundary data : read the data for the unstructured open boundaries.
5	!!======================================================================
6	!! History : 1.0 ! 2005-01 (J. Chanut, A. Sellar) Original code
7	!! - ! 2007-01 (D. Storkey) Update to use IOM module
8	!! - ! 2007-07 (D. Storkey) add bdy_dta_fla
9	!! 3.0 ! 2008-04 (NEMO team) add in the reference version
10	!! 3.3 ! 2010-09 (E.O'Dea) modifications for Shelf configurations
11	!! 3.3 ! 2010-09 (D.Storkey) add ice boundary conditions
12	!! 3.4 ! 2011 (D. Storkey) rewrite in preparation for OBC-BDY merge
13	!! 3.6 ! 2012-01 (C. Rousset) add ice boundary conditions for lim3
14	!!----------------------------------------------------------------------
15	#if defined key_bdy
16	!!----------------------------------------------------------------------
17	!! 'key_bdy' Open Boundary Conditions
18	!!----------------------------------------------------------------------
19	!! bdy_dta : read external data along open boundaries from file
20	!! bdy_dta_init : initialise arrays etc for reading of external data
21	!!----------------------------------------------------------------------
22	USE timing ! Timing
23	USE oce ! ocean dynamics and tracers
24	USE dom_oce ! ocean space and time domain
25	USE phycst ! physical constants
26	USE bdy_oce ! ocean open boundary conditions
27	USE bdytides ! tidal forcing at boundaries
28	USE fldread ! read input fields
29	USE iom ! IOM library
30	USE in_out_manager ! I/O logical units
31	USE dynspg_oce, ONLY: lk_dynspg_ts ! Split-explicit free surface flag
32	#if defined key_lim2
33	USE ice_2
34	#elif defined key_lim3
35	USE ice
36	USE limvar ! redistribute ice input into categories
37	#endif
38	USE sbcapr
39
40	IMPLICIT NONE
41	PRIVATE
42
43	PUBLIC bdy_dta ! routine called by step.F90 and dynspg_ts.F90
44	PUBLIC bdy_dta_init ! routine called by nemogcm.F90
45
46	INTEGER, ALLOCATABLE, DIMENSION(:) :: nb_bdy_fld ! Number of fields to update for each boundary set.
47	INTEGER :: nb_bdy_fld_sum ! Total number of fields to update for all boundary sets.
48
49	LOGICAL, DIMENSION(jp_bdy) :: ln_full_vel_array ! =T => full velocities in 3D boundary conditions
50	! =F => baroclinic velocities in 3D boundary conditions
51	!$AGRIF_DO_NOT_TREAT
52	TYPE(FLD), PUBLIC, ALLOCATABLE, DIMENSION(:), TARGET :: bf ! structure of input fields (file informations, fields read)
53	!$AGRIF_END_DO_NOT_TREAT
54	TYPE(MAP_POINTER), ALLOCATABLE, DIMENSION(:) :: nbmap_ptr ! array of pointers to nbmap
55
56	#if defined key_lim3
57	LOGICAL :: ll_bdylim3 ! determine whether ice input is lim2 (F) or lim3 (T) type
58	INTEGER :: jfld_hti, jfld_hts, jfld_ai ! indices of ice thickness, snow thickness and concentration in bf structure
59	#endif
60
61	# include "domzgr_substitute.h90"
62	!!----------------------------------------------------------------------
63	!! NEMO/OPA 3.3 , NEMO Consortium (2010)
64	!! $Id$
65	!! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
66	!!----------------------------------------------------------------------
67	CONTAINS
68
69	SUBROUTINE bdy_dta( kt, jit, time_offset )
70	!!----------------------------------------------------------------------
71	!! * SUBROUTINE bdy_dta *
72	!!
73	!! ** Purpose : Update external data for open boundary conditions
74	!!
75	!! ** Method : Use fldread.F90
76	!!
77	!!----------------------------------------------------------------------
78	!!
79	INTEGER, INTENT( in ) :: kt ! ocean time-step index
80	INTEGER, INTENT( in ), OPTIONAL :: jit ! subcycle time-step index (for timesplitting option)
81	INTEGER, INTENT( in ), OPTIONAL :: time_offset ! time offset in units of timesteps. NB. if jit
82	! is present then units = subcycle timesteps.
83	! time_offset = 0 => get data at "now" time level
84	! time_offset = -1 => get data at "before" time level
85	! time_offset = +1 => get data at "after" time level
86	! etc.
87	!!
88	INTEGER :: ib_bdy, jfld, jstart, jend, ib, ii, ij, ik, igrd, jl ! local indices
89	INTEGER, DIMENSION(jpbgrd) :: ilen1
90	INTEGER, POINTER, DIMENSION(:) :: nblen, nblenrim ! short cuts
91	TYPE(OBC_DATA), POINTER :: dta ! short cut
92	!!
93	!!---------------------------------------------------------------------------
94	!!
95	IF( nn_timing == 1 ) CALL timing_start('bdy_dta')
96
97	! Initialise data arrays once for all from initial conditions where required
98	!---------------------------------------------------------------------------
99	IF( kt .eq. nit000 .and. .not. PRESENT(jit) ) THEN
100
101	! Calculate depth-mean currents
102	!-----------------------------
103
104	DO ib_bdy = 1, nb_bdy
105
106	nblen => idx_bdy(ib_bdy)%nblen
107	nblenrim => idx_bdy(ib_bdy)%nblenrim
108	dta => dta_bdy(ib_bdy)
109
110	IF( nn_dyn2d_dta(ib_bdy) .eq. 0 ) THEN
111	ilen1(:) = nblen(:)
112	IF( dta%ll_ssh ) THEN
113	igrd = 1
114	DO ib = 1, ilen1(igrd)
115	ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
116	ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
117	dta_bdy(ib_bdy)%ssh(ib) = sshn(ii,ij) * tmask(ii,ij,1)
118	END DO
119	END IF
120	IF( dta%ll_u2d ) THEN
121	igrd = 2
122	DO ib = 1, ilen1(igrd)
123	ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
124	ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
125	dta_bdy(ib_bdy)%u2d(ib) = un_b(ii,ij) * umask(ii,ij,1)
126	END DO
127	END IF
128	IF( dta%ll_v2d ) THEN
129	igrd = 3
130	DO ib = 1, ilen1(igrd)
131	ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
132	ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
133	dta_bdy(ib_bdy)%v2d(ib) = vn_b(ii,ij) * vmask(ii,ij,1)
134	END DO
135	END IF
136	ENDIF
137
138	IF( nn_dyn3d_dta(ib_bdy) .eq. 0 ) THEN
139	ilen1(:) = nblen(:)
140	IF( dta%ll_u3d ) THEN
141	igrd = 2
142	DO ib = 1, ilen1(igrd)
143	DO ik = 1, jpkm1
144	ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
145	ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
146	dta_bdy(ib_bdy)%u3d(ib,ik) = ( un(ii,ij,ik) - un_b(ii,ij) ) * umask(ii,ij,ik)
147	END DO
148	END DO
149	END IF
150	IF( dta%ll_v3d ) THEN
151	igrd = 3
152	DO ib = 1, ilen1(igrd)
153	DO ik = 1, jpkm1
154	ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
155	ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
156	dta_bdy(ib_bdy)%v3d(ib,ik) = ( vn(ii,ij,ik) - vn_b(ii,ij) ) * vmask(ii,ij,ik)
157	END DO
158	END DO
159	END IF
160	ENDIF
161
162	IF( nn_tra_dta(ib_bdy) .eq. 0 ) THEN
163	ilen1(:) = nblen(:)
164	IF( dta%ll_tem ) THEN
165	igrd = 1
166	DO ib = 1, ilen1(igrd)
167	DO ik = 1, jpkm1
168	ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
169	ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
170	dta_bdy(ib_bdy)%tem(ib,ik) = tsn(ii,ij,ik,jp_tem) * tmask(ii,ij,ik)
171	END DO
172	END DO
173	END IF
174	IF( dta%ll_sal ) THEN
175	igrd = 1
176	DO ib = 1, ilen1(igrd)
177	DO ik = 1, jpkm1
178	ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
179	ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
180	dta_bdy(ib_bdy)%sal(ib,ik) = tsn(ii,ij,ik,jp_sal) * tmask(ii,ij,ik)
181	END DO
182	END DO
183	END IF
184	ENDIF
185
186	#if defined key_lim2
187	IF( nn_ice_lim_dta(ib_bdy) .eq. 0 ) THEN
188	ilen1(:) = nblen(:)
189	IF( dta%ll_frld ) THEN
190	igrd = 1
191	DO ib = 1, ilen1(igrd)
192	ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
193	ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
194	dta_bdy(ib_bdy)%frld(ib) = frld(ii,ij) * tmask(ii,ij,1)
195	END DO
196	END IF
197	IF( dta%ll_hicif ) THEN
198	igrd = 1
199	DO ib = 1, ilen1(igrd)
200	ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
201	ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
202	dta_bdy(ib_bdy)%hicif(ib) = hicif(ii,ij) * tmask(ii,ij,1)
203	END DO
204	END IF
205	IF( dta%ll_hsnif ) THEN
206	igrd = 1
207	DO ib = 1, ilen1(igrd)
208	ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
209	ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
210	dta_bdy(ib_bdy)%hsnif(ib) = hsnif(ii,ij) * tmask(ii,ij,1)
211	END DO
212	END IF
213	ENDIF
214	#elif defined key_lim3
215	IF( nn_ice_lim_dta(ib_bdy) .eq. 0 ) THEN
216	ilen1(:) = nblen(:)
217	IF( dta%ll_a_i ) THEN
218	igrd = 1
219	DO jl = 1, jpl
220	DO ib = 1, ilen1(igrd)
221	ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
222	ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
223	dta_bdy(ib_bdy)%a_i (ib,jl) = a_i(ii,ij,jl) * tmask(ii,ij,1)
224	END DO
225	END DO
226	ENDIF
227	IF( dta%ll_ht_i ) THEN
228	igrd = 1
229	DO jl = 1, jpl
230	DO ib = 1, ilen1(igrd)
231	ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
232	ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
233	dta_bdy(ib_bdy)%ht_i (ib,jl) = ht_i(ii,ij,jl) * tmask(ii,ij,1)
234	END DO
235	END DO
236	ENDIF
237	IF( dta%ll_ht_s ) THEN
238	igrd = 1
239	DO jl = 1, jpl
240	DO ib = 1, ilen1(igrd)
241	ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
242	ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
243	dta_bdy(ib_bdy)%ht_s (ib,jl) = ht_s(ii,ij,jl) * tmask(ii,ij,1)
244	END DO
245	END DO
246	ENDIF
247	ENDIF
248	#endif
249
250	ENDDO ! ib_bdy
251
252
253	ENDIF ! kt .eq. nit000
254
255	! update external data from files
256	!--------------------------------
257
258	jstart = 1
259	DO ib_bdy = 1, nb_bdy
260	dta => dta_bdy(ib_bdy)
261	IF( nn_dta(ib_bdy) .eq. 1 ) THEN ! skip this bit if no external data required
262
263	IF( PRESENT(jit) ) THEN
264	! Update barotropic boundary conditions only
265	! jit is optional argument for fld_read and bdytide_update
266	IF( cn_dyn2d(ib_bdy) /= 'none' ) THEN
267	IF( nn_dyn2d_dta(ib_bdy) .eq. 2 ) THEN ! tidal harmonic forcing ONLY: initialise arrays
268	IF( dta%ll_ssh ) dta%ssh(:) = 0.0
269	IF( dta%ll_u2d ) dta%u2d(:) = 0.0
270	IF( dta%ll_u3d ) dta%v2d(:) = 0.0
271	ENDIF
272	IF (cn_tra(ib_bdy) /= 'runoff') THEN
273	IF( nn_dyn2d_dta(ib_bdy) .EQ. 1 .OR. nn_dyn2d_dta(ib_bdy) .EQ. 3 ) THEN
274
275	jend = jstart + dta%nread(2) - 1
276	IF( ln_full_vel_array(ib_bdy) ) THEN
277	CALL fld_read( kt=kt, kn_fsbc=1, sd=bf(jstart:jend), map=nbmap_ptr(jstart:jend), &
278	& kit=jit, kt_offset=time_offset , jpk_bdy=nb_jpk_bdy, fvl=ln_full_vel_array(ib_bdy) )
279	ELSE
280	CALL fld_read( kt=kt, kn_fsbc=1, sd=bf(jstart:jend), map=nbmap_ptr(jstart:jend), &
281	& kit=jit, kt_offset=time_offset )
282	ENDIF
283
284	! If full velocities in boundary data then extract barotropic velocities from 3D fields
285	IF( ln_full_vel_array(ib_bdy) .AND. &
286	& ( nn_dyn2d_dta(ib_bdy) .EQ. 1 .OR. nn_dyn2d_dta(ib_bdy) .EQ. 3 .OR. &
287	& nn_dyn3d_dta(ib_bdy) .EQ. 1 ) )THEN
288
289	igrd = 2 ! zonal velocity
290	dta%u2d(:) = 0.0
291	DO ib = 1, idx_bdy(ib_bdy)%nblen(igrd)
292	ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
293	ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
294	DO ik = 1, jpkm1
295	dta%u2d(ib) = dta%u2d(ib) &
296	& + fse3u(ii,ij,ik) * umask(ii,ij,ik) * dta%u3d(ib,ik)
297	END DO
298	dta%u2d(ib) = dta%u2d(ib) * hur(ii,ij)
299	END DO
300	igrd = 3 ! meridional velocity
301	dta%v2d(:) = 0.0
302	DO ib = 1, idx_bdy(ib_bdy)%nblen(igrd)
303	ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
304	ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
305	DO ik = 1, jpkm1
306	dta%v2d(ib) = dta%v2d(ib) &
307	& + fse3v(ii,ij,ik) * vmask(ii,ij,ik) * dta%v3d(ib,ik)
308	END DO
309	dta%v2d(ib) = dta%v2d(ib) * hvr(ii,ij)
310	END DO
311	ENDIF
312	ENDIF
313	IF( nn_dyn2d_dta(ib_bdy) .ge. 2 ) THEN ! update tidal harmonic forcing
314	CALL bdytide_update( kt=kt, idx=idx_bdy(ib_bdy), dta=dta, td=tides(ib_bdy), &
315	& jit=jit, time_offset=time_offset )
316	ENDIF
317	ENDIF
318	ENDIF
319	ELSE
320	IF (cn_tra(ib_bdy) == 'runoff') then ! runoff condition
321	jend = nb_bdy_fld(ib_bdy)
322	CALL fld_read( kt=kt, kn_fsbc=1, sd=bf(jstart:jend), &
323	& map=nbmap_ptr(jstart:jend), kt_offset=time_offset )
324	!
325	igrd = 2 ! zonal velocity
326	DO ib = 1, idx_bdy(ib_bdy)%nblen(igrd)
327	ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
328	ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
329	dta%u2d(ib) = dta%u2d(ib) / ( e2u(ii,ij) * hu_0(ii,ij) )
330	END DO
331	!
332	igrd = 3 ! meridional velocity
333	DO ib = 1, idx_bdy(ib_bdy)%nblen(igrd)
334	ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
335	ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
336	dta%v2d(ib) = dta%v2d(ib) / ( e1v(ii,ij) * hv_0(ii,ij) )
337	END DO
338	ELSE
339	IF( nn_dyn2d_dta(ib_bdy) .eq. 2 ) THEN ! tidal harmonic forcing ONLY: initialise arrays
340	IF( dta%ll_ssh ) dta%ssh(:) = 0.0
341	IF( dta%ll_u2d ) dta%u2d(:) = 0.0
342	IF( dta%ll_v2d ) dta%v2d(:) = 0.0
343	ENDIF
344	IF( dta%nread(1) .gt. 0 ) THEN ! update external data
345	jend = jstart + dta%nread(1) - 1
346	CALL fld_read( kt=kt, kn_fsbc=1, sd=bf(jstart:jend), &
347	& map=nbmap_ptr(jstart:jend), kt_offset=time_offset, jpk_bdy=nb_jpk_bdy, fvl=ln_full_vel_array(ib_bdy) )
348	ENDIF
349	! If full velocities in boundary data then split into barotropic and baroclinic data
350	IF( ln_full_vel_array(ib_bdy) .and. &
351	& ( nn_dyn2d_dta(ib_bdy) .EQ. 1 .OR. nn_dyn2d_dta(ib_bdy) .EQ. 3 .OR. &
352	& nn_dyn3d_dta(ib_bdy) .EQ. 1 ) ) THEN
353	igrd = 2 ! zonal velocity
354	dta%u2d(:) = 0.0
355	DO ib = 1, idx_bdy(ib_bdy)%nblen(igrd)
356	ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
357	ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
358	DO ik = 1, jpkm1
359	dta%u2d(ib) = dta%u2d(ib) &
360	& + fse3u(ii,ij,ik) * umask(ii,ij,ik) * dta%u3d(ib,ik)
361	END DO
362	dta%u2d(ib) = dta%u2d(ib) * hur(ii,ij)
363	DO ik = 1, jpkm1
364	dta%u3d(ib,ik) = dta%u3d(ib,ik) - dta%u2d(ib)
365	END DO
366	END DO
367	igrd = 3 ! meridional velocity
368	dta%v2d(:) = 0.0
369	DO ib = 1, idx_bdy(ib_bdy)%nblen(igrd)
370	ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
371	ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
372	DO ik = 1, jpkm1
373	dta%v2d(ib) = dta%v2d(ib) &
374	& + fse3v(ii,ij,ik) * vmask(ii,ij,ik) * dta%v3d(ib,ik)
375	END DO
376	dta%v2d(ib) = dta%v2d(ib) * hvr(ii,ij)
377	DO ik = 1, jpkm1
378	dta%v3d(ib,ik) = dta%v3d(ib,ik) - dta%v2d(ib)
379	END DO
380	END DO
381	ENDIF
382
383	ENDIF
384	#if defined key_lim3
385	IF( .NOT. ll_bdylim3 .AND. cn_ice_lim(ib_bdy) /= 'none' .AND. nn_ice_lim_dta(ib_bdy) == 1 ) THEN ! bdy ice input (case input is lim2 type)
386	CALL lim_var_itd ( bf(jfld_hti)%fnow(:,1,1), bf(jfld_hts)%fnow(:,1,1), bf(jfld_ai)%fnow(:,1,1), &
387	& dta_bdy(ib_bdy)%ht_i, dta_bdy(ib_bdy)%ht_s, dta_bdy(ib_bdy)%a_i )
388	ENDIF
389	#endif
390	ENDIF
391	jstart = jstart + dta%nread(1)
392	END IF ! nn_dta(ib_bdy) = 1
393	END DO ! ib_bdy
394
395	! bg jchanut tschanges
396	#if defined key_tide
397	! Add tides if not split-explicit free surface else this is done in ts loop
398	IF (.NOT.lk_dynspg_ts) CALL bdy_dta_tides( kt=kt, time_offset=time_offset )
399	#endif
400	! end jchanut tschanges
401
402	IF ( ln_apr_obc ) THEN
403	DO ib_bdy = 1, nb_bdy
404	IF (cn_tra(ib_bdy) /= 'runoff')THEN
405	igrd = 1 ! meridional velocity
406	DO ib = 1, idx_bdy(ib_bdy)%nblenrim(igrd)
407	ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
408	ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
409	dta_bdy(ib_bdy)%ssh(ib) = dta_bdy(ib_bdy)%ssh(ib) + ssh_ib(ii,ij)
410	ENDDO
411	ENDIF
412	ENDDO
413	ENDIF
414
415	IF( nn_timing == 1 ) CALL timing_stop('bdy_dta')
416
417	END SUBROUTINE bdy_dta
418
419
420	SUBROUTINE bdy_dta_init
421	!!----------------------------------------------------------------------
422	!! * SUBROUTINE bdy_dta_init *
423	!!
424	!! ** Purpose : Initialise arrays for reading of external data
425	!! for open boundary conditions
426	!!
427	!! ** Method :
428	!!
429	!!----------------------------------------------------------------------
430	USE dynspg_oce, ONLY: lk_dynspg_ts
431	!!
432	INTEGER :: ib_bdy, jfld, jstart, jend, ierror ! local indices
433	INTEGER :: ios ! Local integer output status for namelist read
434	!!
435	CHARACTER(len=100) :: cn_dir ! Root directory for location of data files
436	CHARACTER(len=100), DIMENSION(nb_bdy) :: cn_dir_array ! Root directory for location of data files
437	LOGICAL :: ln_full_vel ! =T => full velocities in 3D boundary data
438	! =F => baroclinic velocities in 3D boundary data
439	INTEGER :: ilen_global ! Max length required for global bdy dta arrays
440	INTEGER, ALLOCATABLE, DIMENSION(:) :: ilen1, ilen3 ! size of 1st and 3rd dimensions of local arrays
441	INTEGER, ALLOCATABLE, DIMENSION(:) :: ibdy ! bdy set for a particular jfld
442	INTEGER, ALLOCATABLE, DIMENSION(:) :: igrid ! index for grid type (1,2,3 = T,U,V)
443	INTEGER, POINTER, DIMENSION(:) :: nblen, nblenrim ! short cuts
444	TYPE(OBC_DATA), POINTER :: dta ! short cut
445	#if defined key_lim3
446	INTEGER, DIMENSION(3) :: zdimsz ! number of elements in each of the 4 dimensions (i.e. i,j,t,ice-cat) for an array
447	INTEGER :: zndims ! number of dimensions in an array (i.e. 3 = wo ice cat; 4 = w ice cat)
448	INTEGER :: inum,id1 ! local integer
449	#endif
450	TYPE(FLD_N), ALLOCATABLE, DIMENSION(:) :: blf_i ! array of namelist information structures
451	TYPE(FLD_N) :: bn_tem, bn_sal, bn_u3d, bn_v3d !
452	TYPE(FLD_N) :: bn_ssh, bn_u2d, bn_v2d ! informations about the fields to be read
453	#if defined key_lim2
454	TYPE(FLD_N) :: bn_frld, bn_hicif, bn_hsnif !
455	#elif defined key_lim3
456	TYPE(FLD_N) :: bn_a_i, bn_ht_i, bn_ht_s
457	#endif
458	NAMELIST/nambdy_dta/ cn_dir, bn_tem, bn_sal, bn_u3d, bn_v3d, bn_ssh, bn_u2d, bn_v2d
459	#if defined key_lim2
460	NAMELIST/nambdy_dta/ bn_frld, bn_hicif, bn_hsnif
461	#elif defined key_lim3
462	NAMELIST/nambdy_dta/ bn_a_i, bn_ht_i, bn_ht_s
463	#endif
464	NAMELIST/nambdy_dta/ ln_full_vel
465	!!---------------------------------------------------------------------------
466
467	IF( nn_timing == 1 ) CALL timing_start('bdy_dta_init')
468
469	IF(lwp) WRITE(numout,*)
470	IF(lwp) WRITE(numout,*) 'bdy_dta_ini : initialization of data at the open boundaries'
471	IF(lwp) WRITE(numout,*) '~~~~~~~~~~'
472	IF(lwp) WRITE(numout,*) ''
473
474	! Set nn_dta
475	DO ib_bdy = 1, nb_bdy
476	nn_dta(ib_bdy) = MAX( nn_dyn2d_dta(ib_bdy) &
477	,nn_dyn3d_dta(ib_bdy) &
478	,nn_tra_dta(ib_bdy) &
479	#if ( defined key_lim2 \|\| defined key_lim3 )
480	,nn_ice_lim_dta(ib_bdy) &
481	#endif
482	)
483	IF(nn_dta(ib_bdy) .gt. 1) nn_dta(ib_bdy) = 1
484	END DO
485
486	! Work out upper bound of how many fields there are to read in and allocate arrays
487	! ---------------------------------------------------------------------------
488	ALLOCATE( nb_bdy_fld(nb_bdy) )
489	nb_bdy_fld(:) = 0
490	DO ib_bdy = 1, nb_bdy
491	IF( cn_dyn2d(ib_bdy) /= 'none' .and. ( nn_dyn2d_dta(ib_bdy) .eq. 1 .or. nn_dyn2d_dta(ib_bdy) .eq. 3 ) ) THEN
492	nb_bdy_fld(ib_bdy) = nb_bdy_fld(ib_bdy) + 3
493	ENDIF
494	IF( cn_dyn3d(ib_bdy) /= 'none' .and. nn_dyn3d_dta(ib_bdy) .eq. 1 ) THEN
495	nb_bdy_fld(ib_bdy) = nb_bdy_fld(ib_bdy) + 2
496	ENDIF
497	IF( cn_tra(ib_bdy) /= 'none' .and. nn_tra_dta(ib_bdy) .eq. 1 ) THEN
498	nb_bdy_fld(ib_bdy) = nb_bdy_fld(ib_bdy) + 2
499	ENDIF
500	#if ( defined key_lim2 \|\| defined key_lim3 )
501	IF( cn_ice_lim(ib_bdy) /= 'none' .and. nn_ice_lim_dta(ib_bdy) .eq. 1 ) THEN
502	nb_bdy_fld(ib_bdy) = nb_bdy_fld(ib_bdy) + 3
503	ENDIF
504	#endif
505	IF(lwp) WRITE(numout,*) 'Maximum number of files to open =',nb_bdy_fld(ib_bdy)
506	ENDDO
507
508	nb_bdy_fld_sum = SUM( nb_bdy_fld )
509
510	ALLOCATE( bf(nb_bdy_fld_sum), STAT=ierror )
511	IF( ierror > 0 ) THEN
512	CALL ctl_stop( 'bdy_dta: unable to allocate bf structure' ) ; RETURN
513	ENDIF
514	ALLOCATE( blf_i(nb_bdy_fld_sum), STAT=ierror )
515	IF( ierror > 0 ) THEN
516	CALL ctl_stop( 'bdy_dta: unable to allocate blf_i structure' ) ; RETURN
517	ENDIF
518	ALLOCATE( nbmap_ptr(nb_bdy_fld_sum), STAT=ierror )
519	IF( ierror > 0 ) THEN
520	CALL ctl_stop( 'bdy_dta: unable to allocate nbmap_ptr structure' ) ; RETURN
521	ENDIF
522	ALLOCATE( ilen1(nb_bdy_fld_sum), ilen3(nb_bdy_fld_sum) )
523	ALLOCATE( ibdy(nb_bdy_fld_sum) )
524	ALLOCATE( igrid(nb_bdy_fld_sum) )
525
526	! Read namelists
527	! --------------
528	REWIND(numnam_ref)
529	REWIND(numnam_cfg)
530	jfld = 0
531	DO ib_bdy = 1, nb_bdy
532	IF( nn_dta(ib_bdy) .eq. 1 ) THEN
533	READ ( numnam_ref, nambdy_dta, IOSTAT = ios, ERR = 901)
534	901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nambdy_dta in reference namelist', lwp )
535
536	READ ( numnam_cfg, nambdy_dta, IOSTAT = ios, ERR = 902 )
537	902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nambdy_dta in configuration namelist', lwp )
538	IF(lwm) WRITE ( numond, nambdy_dta )
539
540	cn_dir_array(ib_bdy) = cn_dir
541	ln_full_vel_array(ib_bdy) = ln_full_vel
542	!dta%ll_fvl = ln_full_vel ! jdha need this in fldread routine to work out what type of correction to apply to interpolated bdy data (maybe we replace all instances of ln_full_vel_array with this rather than duplicate)
543
544	nblen => idx_bdy(ib_bdy)%nblen
545	nblenrim => idx_bdy(ib_bdy)%nblenrim
546	dta => dta_bdy(ib_bdy)
547	dta%nread(2) = 0
548
549	! Only read in necessary fields for this set.
550	! Important that barotropic variables come first.
551	IF( nn_dyn2d_dta(ib_bdy) .eq. 1 .or. nn_dyn2d_dta(ib_bdy) .eq. 3 ) THEN
552
553	IF( dta%ll_ssh ) THEN
554	if(lwp) write(numout,*) '++++++ reading in ssh field'
555	jfld = jfld + 1
556	blf_i(jfld) = bn_ssh
557	ibdy(jfld) = ib_bdy
558	igrid(jfld) = 1
559	ilen1(jfld) = nblen(igrid(jfld))
560	ilen3(jfld) = 1
561	dta%nread(2) = dta%nread(2) + 1
562	ENDIF
563
564	IF( dta%ll_u2d .and. .not. ln_full_vel_array(ib_bdy) ) THEN
565	if(lwp) write(numout,*) '++++++ reading in u2d field'
566	jfld = jfld + 1
567	blf_i(jfld) = bn_u2d
568	ibdy(jfld) = ib_bdy
569	igrid(jfld) = 2
570	ilen1(jfld) = nblen(igrid(jfld))
571	ilen3(jfld) = 1
572	dta%nread(2) = dta%nread(2) + 1
573	ENDIF
574
575	IF( dta%ll_v2d .and. .not. ln_full_vel_array(ib_bdy) ) THEN
576	if(lwp) write(numout,*) '++++++ reading in v2d field'
577	jfld = jfld + 1
578	blf_i(jfld) = bn_v2d
579	ibdy(jfld) = ib_bdy
580	igrid(jfld) = 3
581	ilen1(jfld) = nblen(igrid(jfld))
582	ilen3(jfld) = 1
583	dta%nread(2) = dta%nread(2) + 1
584	ENDIF
585
586	ENDIF
587
588	! read 3D velocities if baroclinic velocities require OR if
589	! barotropic velocities required and ln_full_vel set to .true.
590	IF( nn_dyn3d_dta(ib_bdy) .eq. 1 .or. &
591	& ( ln_full_vel_array(ib_bdy) .and. ( nn_dyn2d_dta(ib_bdy) .eq. 1 .or. nn_dyn2d_dta(ib_bdy) .eq. 3 ) ) ) THEN
592
593	IF( dta%ll_u3d .or. ( ln_full_vel_array(ib_bdy) .and. dta%ll_u2d ) ) THEN
594	if(lwp) write(numout,*) '++++++ reading in u3d field'
595	jfld = jfld + 1
596	blf_i(jfld) = bn_u3d
597	ibdy(jfld) = ib_bdy
598	igrid(jfld) = 2
599	ilen1(jfld) = nblen(igrid(jfld))
600	ilen3(jfld) = jpk
601	IF( ln_full_vel_array(ib_bdy) .and. dta%ll_u2d ) dta%nread(2) = dta%nread(2) + 1
602	ENDIF
603
604	IF( dta%ll_v3d .or. ( ln_full_vel_array(ib_bdy) .and. dta%ll_v2d ) ) THEN
605	if(lwp) write(numout,*) '++++++ reading in v3d field'
606	jfld = jfld + 1
607	blf_i(jfld) = bn_v3d
608	ibdy(jfld) = ib_bdy
609	igrid(jfld) = 3
610	ilen1(jfld) = nblen(igrid(jfld))
611	ilen3(jfld) = jpk
612	IF( ln_full_vel_array(ib_bdy) .and. dta%ll_v2d ) dta%nread(2) = dta%nread(2) + 1
613	ENDIF
614
615	ENDIF
616
617	! temperature and salinity
618	IF( nn_tra_dta(ib_bdy) .eq. 1 ) THEN
619
620	IF( dta%ll_tem ) THEN
621	if(lwp) write(numout,*) '++++++ reading in tem field'
622	jfld = jfld + 1
623	blf_i(jfld) = bn_tem
624	ibdy(jfld) = ib_bdy
625	igrid(jfld) = 1
626	ilen1(jfld) = nblen(igrid(jfld))
627	ilen3(jfld) = jpk
628	ENDIF
629
630	IF( dta%ll_sal ) THEN
631	if(lwp) write(numout,*) '++++++ reading in sal field'
632	jfld = jfld + 1
633	blf_i(jfld) = bn_sal
634	ibdy(jfld) = ib_bdy
635	igrid(jfld) = 1
636	ilen1(jfld) = nblen(igrid(jfld))
637	ilen3(jfld) = jpk
638	ENDIF
639
640	ENDIF
641
642	#if defined key_lim2
643	! sea ice
644	IF( nn_ice_lim_dta(ib_bdy) .eq. 1 ) THEN
645
646	IF( dta%ll_frld ) THEN
647	jfld = jfld + 1
648	blf_i(jfld) = bn_frld
649	ibdy(jfld) = ib_bdy
650	igrid(jfld) = 1
651	ilen1(jfld) = nblen(igrid(jfld))
652	ilen3(jfld) = 1
653	ENDIF
654
655	IF( dta%ll_hicif ) THEN
656	jfld = jfld + 1
657	blf_i(jfld) = bn_hicif
658	ibdy(jfld) = ib_bdy
659	igrid(jfld) = 1
660	ilen1(jfld) = nblen(igrid(jfld))
661	ilen3(jfld) = 1
662	ENDIF
663
664	IF( dta%ll_hsnif ) THEN
665	jfld = jfld + 1
666	blf_i(jfld) = bn_hsnif
667	ibdy(jfld) = ib_bdy
668	igrid(jfld) = 1
669	ilen1(jfld) = nblen(igrid(jfld))
670	ilen3(jfld) = 1
671	ENDIF
672
673	ENDIF
674	#elif defined key_lim3
675	! sea ice
676	IF( nn_ice_lim_dta(ib_bdy) .eq. 1 ) THEN
677
678	! Test for types of ice input (lim2 or lim3)
679	CALL iom_open ( bn_a_i%clname, inum )
680	id1 = iom_varid ( inum, bn_a_i%clvar, kdimsz=zdimsz, kndims=zndims, ldstop = .FALSE. )
681	CALL iom_close ( inum )
682	!CALL fld_clopn ( bn_a_i, nyear, nmonth, nday, ldstop=.TRUE. )
683	!CALL iom_open ( bn_a_i%clname, inum )
684	!id1 = iom_varid ( bn_a_i%num, bn_a_i%clvar, kdimsz=zdimsz, kndims=zndims, ldstop = .FALSE. )
685	IF ( zndims == 4 ) THEN
686	ll_bdylim3 = .TRUE. ! lim3 input
687	ELSE
688	ll_bdylim3 = .FALSE. ! lim2 input
689	ENDIF
690	! End test
691
692	IF( dta%ll_a_i ) THEN
693	jfld = jfld + 1
694	blf_i(jfld) = bn_a_i
695	ibdy(jfld) = ib_bdy
696	igrid(jfld) = 1
697	ilen1(jfld) = nblen(igrid(jfld))
698	IF ( ll_bdylim3 ) THEN ; ilen3(jfld)=jpl ; ELSE ; ilen3(jfld)=1 ; ENDIF
699	ENDIF
700
701	IF( dta%ll_ht_i ) THEN
702	jfld = jfld + 1
703	blf_i(jfld) = bn_ht_i
704	ibdy(jfld) = ib_bdy
705	igrid(jfld) = 1
706	ilen1(jfld) = nblen(igrid(jfld))
707	IF ( ll_bdylim3 ) THEN ; ilen3(jfld)=jpl ; ELSE ; ilen3(jfld)=1 ; ENDIF
708	ENDIF
709
710	IF( dta%ll_ht_s ) THEN
711	jfld = jfld + 1
712	blf_i(jfld) = bn_ht_s
713	ibdy(jfld) = ib_bdy
714	igrid(jfld) = 1
715	ilen1(jfld) = nblen(igrid(jfld))
716	IF ( ll_bdylim3 ) THEN ; ilen3(jfld)=jpl ; ELSE ; ilen3(jfld)=1 ; ENDIF
717	ENDIF
718
719	ENDIF
720	#endif
721	! Recalculate field counts
722	!-------------------------
723	IF( ib_bdy .eq. 1 ) THEN
724	nb_bdy_fld_sum = 0
725	nb_bdy_fld(ib_bdy) = jfld
726	nb_bdy_fld_sum = jfld
727	ELSE
728	nb_bdy_fld(ib_bdy) = jfld - nb_bdy_fld_sum
729	nb_bdy_fld_sum = nb_bdy_fld_sum + nb_bdy_fld(ib_bdy)
730	ENDIF
731
732	dta%nread(1) = nb_bdy_fld(ib_bdy)
733
734	ENDIF ! nn_dta .eq. 1
735	ENDDO ! ib_bdy
736
737	DO jfld = 1, nb_bdy_fld_sum
738	ALLOCATE( bf(jfld)%fnow(ilen1(jfld),1,ilen3(jfld)) )
739	IF( blf_i(jfld)%ln_tint ) ALLOCATE( bf(jfld)%fdta(ilen1(jfld),1,ilen3(jfld),2) )
740	nbmap_ptr(jfld)%ptr => idx_bdy(ibdy(jfld))%nbmap(:,igrid(jfld))
741	nbmap_ptr(jfld)%ll_unstruc = ln_coords_file(ibdy(jfld))
742	ENDDO
743
744	! fill bf with blf_i and control print
745	!-------------------------------------
746	jstart = 1
747	DO ib_bdy = 1, nb_bdy
748	jend = jstart - 1 + nb_bdy_fld(ib_bdy)
749	CALL fld_fill( bf(jstart:jend), blf_i(jstart:jend), cn_dir_array(ib_bdy), 'bdy_dta', &
750	& 'open boundary conditions', 'nambdy_dta' )
751	jstart = jend + 1
752	ENDDO
753
754	DO jfld = 1, nb_bdy_fld_sum
755	bf(jfld)%igrd = igrid(jfld)
756	bf(jfld)%ibdy = ibdy(jfld)
757	ENDDO
758
759	! Initialise local boundary data arrays
760	! nn_xxx_dta=0 : allocate space - will be filled from initial conditions later
761	! nn_xxx_dta=1 : point to "fnow" arrays
762	!-------------------------------------
763
764	jfld = 0
765	DO ib_bdy=1, nb_bdy
766
767	nblen => idx_bdy(ib_bdy)%nblen
768	dta => dta_bdy(ib_bdy)
769
770	if(lwp) then
771	write(numout,*) '++++++ dta%ll_ssh = ',dta%ll_ssh
772	write(numout,*) '++++++ dta%ll_u2d = ',dta%ll_u2d
773	write(numout,*) '++++++ dta%ll_v2d = ',dta%ll_v2d
774	write(numout,*) '++++++ dta%ll_u3d = ',dta%ll_u3d
775	write(numout,*) '++++++ dta%ll_v3d = ',dta%ll_v3d
776	write(numout,*) '++++++ dta%ll_tem = ',dta%ll_tem
777	write(numout,*) '++++++ dta%ll_sal = ',dta%ll_sal
778	endif
779
780	IF ( nn_dyn2d_dta(ib_bdy) .eq. 0 .or. nn_dyn2d_dta(ib_bdy) .eq. 2 ) THEN
781	if(lwp) write(numout,*) '++++++ dta%ssh/u2d/u3d allocated space'
782	IF( dta%ll_ssh ) ALLOCATE( dta%ssh(nblen(1)) )
783	IF( dta%ll_u2d ) ALLOCATE( dta%u2d(nblen(2)) )
784	IF( dta%ll_v2d ) ALLOCATE( dta%v2d(nblen(3)) )
785	ENDIF
786	IF ( nn_dyn2d_dta(ib_bdy) .eq. 1 .or. nn_dyn2d_dta(ib_bdy) .eq. 3 ) THEN
787	IF( dta%ll_ssh ) THEN
788	if(lwp) write(numout,*) '++++++ dta%ssh pointing to fnow'
789	jfld = jfld + 1
790	dta%ssh => bf(jfld)%fnow(:,1,1)
791	ENDIF
792	IF ( dta%ll_u2d ) THEN
793	IF ( ln_full_vel_array(ib_bdy) ) THEN
794	if(lwp) write(numout,*) '++++++ dta%u2d allocated space'
795	ALLOCATE( dta%u2d(nblen(2)) )
796	ELSE
797	if(lwp) write(numout,*) '++++++ dta%u2d pointing to fnow'
798	jfld = jfld + 1
799	dta%u2d => bf(jfld)%fnow(:,1,1)
800	ENDIF
801	ENDIF
802	IF ( dta%ll_v2d ) THEN
803	IF ( ln_full_vel_array(ib_bdy) ) THEN
804	if(lwp) write(numout,*) '++++++ dta%v2d allocated space'
805	ALLOCATE( dta%v2d(nblen(3)) )
806	ELSE
807	if(lwp) write(numout,*) '++++++ dta%v2d pointing to fnow'
808	jfld = jfld + 1
809	dta%v2d => bf(jfld)%fnow(:,1,1)
810	ENDIF
811	ENDIF
812	ENDIF
813
814	IF ( nn_dyn3d_dta(ib_bdy) .eq. 0 ) THEN
815	if(lwp) write(numout,*) '++++++ dta%u3d/v3d allocated space'
816	IF( dta%ll_u3d ) ALLOCATE( dta_bdy(ib_bdy)%u3d(nblen(2),jpk) )
817	IF( dta%ll_v3d ) ALLOCATE( dta_bdy(ib_bdy)%v3d(nblen(3),jpk) )
818	ENDIF
819	IF ( nn_dyn3d_dta(ib_bdy) .eq. 1 .or. &
820	& ( ln_full_vel_array(ib_bdy) .and. ( nn_dyn2d_dta(ib_bdy) .eq. 1 .or. nn_dyn2d_dta(ib_bdy) .eq. 3 ) ) ) THEN
821	IF ( dta%ll_u3d .or. ( ln_full_vel_array(ib_bdy) .and. dta%ll_u2d ) ) THEN
822	if(lwp) write(numout,*) '++++++ dta%u3d pointing to fnow'
823	jfld = jfld + 1
824	dta_bdy(ib_bdy)%u3d => bf(jfld)%fnow(:,1,:)
825	ENDIF
826	IF ( dta%ll_v3d .or. ( ln_full_vel_array(ib_bdy) .and. dta%ll_v2d ) ) THEN
827	if(lwp) write(numout,*) '++++++ dta%v3d pointing to fnow'
828	jfld = jfld + 1
829	dta_bdy(ib_bdy)%v3d => bf(jfld)%fnow(:,1,:)
830	ENDIF
831	ENDIF
832
833	IF( nn_tra_dta(ib_bdy) .eq. 0 ) THEN
834	if(lwp) write(numout,*) '++++++ dta%tem/sal allocated space'
835	IF( dta%ll_tem ) ALLOCATE( dta_bdy(ib_bdy)%tem(nblen(1),jpk) )
836	IF( dta%ll_sal ) ALLOCATE( dta_bdy(ib_bdy)%sal(nblen(1),jpk) )
837	ELSE
838	IF( dta%ll_tem ) THEN
839	if(lwp) write(numout,*) '++++++ dta%tem pointing to fnow'
840	jfld = jfld + 1
841	dta_bdy(ib_bdy)%tem => bf(jfld)%fnow(:,1,:)
842	ENDIF
843	IF( dta%ll_sal ) THEN
844	if(lwp) write(numout,*) '++++++ dta%sal pointing to fnow'
845	jfld = jfld + 1
846	dta_bdy(ib_bdy)%sal => bf(jfld)%fnow(:,1,:)
847	ENDIF
848	ENDIF
849
850	#if defined key_lim2
851	IF (cn_ice_lim(ib_bdy) /= 'none') THEN
852	IF( nn_ice_lim_dta(ib_bdy) .eq. 0 ) THEN
853	ALLOCATE( dta_bdy(ib_bdy)%frld(nblen(1)) )
854	ALLOCATE( dta_bdy(ib_bdy)%hicif(nblen(1)) )
855	ALLOCATE( dta_bdy(ib_bdy)%hsnif(nblen(1)) )
856	ELSE
857	jfld = jfld + 1
858	dta_bdy(ib_bdy)%frld => bf(jfld)%fnow(:,1,1)
859	jfld = jfld + 1
860	dta_bdy(ib_bdy)%hicif => bf(jfld)%fnow(:,1,1)
861	jfld = jfld + 1
862	dta_bdy(ib_bdy)%hsnif => bf(jfld)%fnow(:,1,1)
863	ENDIF
864	ENDIF
865	#elif defined key_lim3
866	IF (cn_ice_lim(ib_bdy) /= 'none') THEN
867	IF( nn_ice_lim_dta(ib_bdy) .eq. 0 ) THEN
868	ALLOCATE( dta_bdy(ib_bdy)%a_i (nblen(1),jpl) )
869	ALLOCATE( dta_bdy(ib_bdy)%ht_i(nblen(1),jpl) )
870	ALLOCATE( dta_bdy(ib_bdy)%ht_s(nblen(1),jpl) )
871	ELSE
872	IF ( ll_bdylim3 ) THEN ! case input is lim3 type
873	jfld = jfld + 1
874	dta_bdy(ib_bdy)%a_i => bf(jfld)%fnow(:,1,:)
875	jfld = jfld + 1
876	dta_bdy(ib_bdy)%ht_i => bf(jfld)%fnow(:,1,:)
877	jfld = jfld + 1
878	dta_bdy(ib_bdy)%ht_s => bf(jfld)%fnow(:,1,:)
879	ELSE ! case input is lim2 type
880	jfld_ai = jfld + 1
881	jfld_hti = jfld + 2
882	jfld_hts = jfld + 3
883	jfld = jfld + 3
884	ALLOCATE( dta_bdy(ib_bdy)%a_i (nblen(1),jpl) )
885	ALLOCATE( dta_bdy(ib_bdy)%ht_i(nblen(1),jpl) )
886	ALLOCATE( dta_bdy(ib_bdy)%ht_s(nblen(1),jpl) )
887	dta_bdy(ib_bdy)%a_i (:,:) = 0.0
888	dta_bdy(ib_bdy)%ht_i(:,:) = 0.0
889	dta_bdy(ib_bdy)%ht_s(:,:) = 0.0
890	ENDIF
891
892	ENDIF
893	ENDIF
894	#endif
895
896	ENDDO ! ib_bdy
897
898	IF( nn_timing == 1 ) CALL timing_stop('bdy_dta_init')
899
900	END SUBROUTINE bdy_dta_init
901
902	#else
903	!!----------------------------------------------------------------------
904	!! Dummy module NO Open Boundary Conditions
905	!!----------------------------------------------------------------------
906	CONTAINS
907	SUBROUTINE bdy_dta( kt, jit, time_offset ) ! Empty routine
908	INTEGER, INTENT( in ) :: kt
909	INTEGER, INTENT( in ), OPTIONAL :: jit
910	INTEGER, INTENT( in ), OPTIONAL :: time_offset
911	WRITE(,) 'bdy_dta: You should not have seen this print! error?', kt
912	END SUBROUTINE bdy_dta
913	SUBROUTINE bdy_dta_init() ! Empty routine
914	WRITE(,) 'bdy_dta_init: You should not have seen this print! error?'
915	END SUBROUTINE bdy_dta_init
916	#endif
917
918	!!==============================================================================
919	END MODULE bdydta

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

Download in other formats: