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 @ 6079

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

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

Download in other formats: