Context Navigation

bdydta.F90 @ 15473

Last change on this file since 15473 was 7567, checked in by hadjt, 7 years ago
CO6 version adapted for shelf seas climate projections, including added diagnostics
File size: 41.7 KB

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

Download in other formats: