Context Navigation

← Previous Revision
Latest Revision
Next Revision →
Blame
Revision Log

bdydta.F90 @ 7029

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

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

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

Context Navigation

source: branches/NERC/dev_r6998_ORCHESTRA/NEMOGCM/NEMO/OPA_SRC/BDY/bdydta.F90 @ 7029

Download in other formats: