Context Navigation

bdydta.F90 @ 7753

Last change on this file since 7753 was 7646, checked in by timgraham, 7 years ago

Merge of dev_merge_2016 into trunk. UPDATE TO ARCHFILES NEEDED for XIOS2.
LIM_SRC_s/limrhg.F90 to follow in next commit due to change of kind (I'm unable to do it in this commit).
Merged using the following steps:

1) svn merge --reintegrate svn+ssh://forge.ipsl.jussieu.fr/ipsl/forge/projets/nemo/svn/trunk .
2) Resolve minor conflicts in sette.sh and namelist_cfg for ORCA2LIM3 (due to a change in trunk after branch was created)
3) svn commit
4) svn switch svn+ssh://forge.ipsl.jussieu.fr/ipsl/forge/projets/nemo/svn/trunk
5) svn merge svn+ssh://forge.ipsl.jussieu.fr/ipsl/forge/projets/nemo/svn/branches/2016/dev_merge_2016 .
6) At this stage I checked out a clean copy of the branch to compare against what is about to be committed to the trunk.
6) svn commit #Commit code to the trunk

In this commit I have also reverted a change to Fcheck_archfile.sh which was causing problems on the Paris machine.

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

Download in other formats: