Context Navigation

bdydta.F90 @ 9383

Last change on this file since 9383 was 9383, checked in by andmirek, 6 years ago
#2050 fixes and changes
File size: 42.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
40	IMPLICIT NONE
41	PRIVATE
42
43	PUBLIC bdy_dta ! routine called by step.F90 and dynspg_ts.F90
44	PUBLIC bdy_dta_init ! routine called by nemogcm.F90
45	PRIVATE dta_namelist
46
47	INTEGER, ALLOCATABLE, DIMENSION(:) :: nb_bdy_fld ! Number of fields to update for each boundary set.
48	INTEGER :: nb_bdy_fld_sum ! Total number of fields to update for all boundary sets.
49	LOGICAL :: ln_bdy_sio ! single processor read flag
50
51	LOGICAL, DIMENSION(jp_bdy) :: ln_full_vel_array ! =T => full velocities in 3D boundary conditions
52	! =F => baroclinic velocities in 3D boundary conditions
53	!$AGRIF_DO_NOT_TREAT
54	TYPE(FLD), PUBLIC, ALLOCATABLE, DIMENSION(:), TARGET :: bf ! structure of input fields (file informations, fields read)
55	!$AGRIF_END_DO_NOT_TREAT
56	TYPE(MAP_POINTER), ALLOCATABLE, DIMENSION(:) :: nbmap_ptr ! array of pointers to nbmap
57
58	#if defined key_lim3
59	LOGICAL :: ll_bdylim3 ! determine whether ice input is lim2 (F) or lim3 (T) type
60	INTEGER :: jfld_hti, jfld_hts, jfld_ai ! indices of ice thickness, snow thickness and concentration in bf structure
61	#endif
62
63	# include "domzgr_substitute.h90"
64	!!----------------------------------------------------------------------
65	!! NEMO/OPA 3.3 , NEMO Consortium (2010)
66	!! $Id$
67	!! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
68	!!----------------------------------------------------------------------
69	CONTAINS
70
71	SUBROUTINE bdy_dta( kt, jit, time_offset )
72	!!----------------------------------------------------------------------
73	!! * SUBROUTINE bdy_dta *
74	!!
75	!! ** Purpose : Update external data for open boundary conditions
76	!!
77	!! ** Method : Use fldread.F90
78	!!
79	!!----------------------------------------------------------------------
80	!!
81	INTEGER, INTENT( in ) :: kt ! ocean time-step index
82	INTEGER, INTENT( in ), OPTIONAL :: jit ! subcycle time-step index (for timesplitting option)
83	INTEGER, INTENT( in ), OPTIONAL :: time_offset ! time offset in units of timesteps. NB. if jit
84	! is present then units = subcycle timesteps.
85	! time_offset = 0 => get data at "now" time level
86	! time_offset = -1 => get data at "before" time level
87	! time_offset = +1 => get data at "after" time level
88	! etc.
89	!!
90	INTEGER :: ib_bdy, jfld, jstart, jend, ib, ii, ij, ik, igrd, jl ! local indices
91	INTEGER, DIMENSION(jpbgrd) :: ilen1
92	INTEGER, POINTER, DIMENSION(:) :: nblen, nblenrim ! short cuts
93	TYPE(OBC_DATA), POINTER :: dta ! short cut
94	!!
95	!!---------------------------------------------------------------------------
96	!!
97	IF( nn_timing == 1 ) CALL timing_start('bdy_dta')
98
99	! Initialise data arrays once for all from initial conditions where required
100	!---------------------------------------------------------------------------
101	IF( kt .eq. nit000 .and. .not. PRESENT(jit) ) THEN
102
103	! Calculate depth-mean currents
104	!-----------------------------
105
106	DO ib_bdy = 1, nb_bdy
107
108	nblen => idx_bdy(ib_bdy)%nblen
109	nblenrim => idx_bdy(ib_bdy)%nblenrim
110	dta => dta_bdy(ib_bdy)
111
112	IF( nn_dyn2d_dta(ib_bdy) .eq. 0 ) THEN
113	ilen1(:) = nblen(:)
114	IF( dta%ll_ssh ) THEN
115	igrd = 1
116	DO ib = 1, ilen1(igrd)
117	ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
118	ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
119	dta_bdy(ib_bdy)%ssh(ib) = sshn(ii,ij) * tmask(ii,ij,1)
120	END DO
121	END IF
122	IF( dta%ll_u2d ) THEN
123	igrd = 2
124	DO ib = 1, ilen1(igrd)
125	ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
126	ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
127	dta_bdy(ib_bdy)%u2d(ib) = un_b(ii,ij) * umask(ii,ij,1)
128	END DO
129	END IF
130	IF( dta%ll_v2d ) THEN
131	igrd = 3
132	DO ib = 1, ilen1(igrd)
133	ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
134	ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
135	dta_bdy(ib_bdy)%v2d(ib) = vn_b(ii,ij) * vmask(ii,ij,1)
136	END DO
137	END IF
138	ENDIF
139
140	IF( nn_dyn3d_dta(ib_bdy) .eq. 0 ) THEN
141	ilen1(:) = nblen(:)
142	IF( dta%ll_u3d ) THEN
143	igrd = 2
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)%u3d(ib,ik) = ( un(ii,ij,ik) - un_b(ii,ij) ) * umask(ii,ij,ik)
149	END DO
150	END DO
151	END IF
152	IF( dta%ll_v3d ) THEN
153	igrd = 3
154	DO ib = 1, ilen1(igrd)
155	DO ik = 1, jpkm1
156	ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
157	ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
158	dta_bdy(ib_bdy)%v3d(ib,ik) = ( vn(ii,ij,ik) - vn_b(ii,ij) ) * vmask(ii,ij,ik)
159	END DO
160	END DO
161	END IF
162	ENDIF
163
164	IF( nn_tra_dta(ib_bdy) .eq. 0 ) THEN
165	ilen1(:) = nblen(:)
166	IF( dta%ll_tem ) 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)%tem(ib,ik) = tsn(ii,ij,ik,jp_tem) * tmask(ii,ij,ik)
173	END DO
174	END DO
175	END IF
176	IF( dta%ll_sal ) THEN
177	igrd = 1
178	DO ib = 1, ilen1(igrd)
179	DO ik = 1, jpkm1
180	ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
181	ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
182	dta_bdy(ib_bdy)%sal(ib,ik) = tsn(ii,ij,ik,jp_sal) * tmask(ii,ij,ik)
183	END DO
184	END DO
185	END IF
186	ENDIF
187
188	#if defined key_lim2
189	IF( nn_ice_lim_dta(ib_bdy) .eq. 0 ) THEN
190	ilen1(:) = nblen(:)
191	IF( dta%ll_frld ) THEN
192	igrd = 1
193	DO ib = 1, ilen1(igrd)
194	ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
195	ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
196	dta_bdy(ib_bdy)%frld(ib) = frld(ii,ij) * tmask(ii,ij,1)
197	END DO
198	END IF
199	IF( dta%ll_hicif ) THEN
200	igrd = 1
201	DO ib = 1, ilen1(igrd)
202	ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
203	ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
204	dta_bdy(ib_bdy)%hicif(ib) = hicif(ii,ij) * tmask(ii,ij,1)
205	END DO
206	END IF
207	IF( dta%ll_hsnif ) THEN
208	igrd = 1
209	DO ib = 1, ilen1(igrd)
210	ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
211	ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
212	dta_bdy(ib_bdy)%hsnif(ib) = hsnif(ii,ij) * tmask(ii,ij,1)
213	END DO
214	END IF
215	ENDIF
216	#elif defined key_lim3
217	IF( nn_ice_lim_dta(ib_bdy) .eq. 0 ) THEN
218	ilen1(:) = nblen(:)
219	IF( dta%ll_a_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)%a_i (ib,jl) = a_i(ii,ij,jl) * tmask(ii,ij,1)
226	END DO
227	END DO
228	ENDIF
229	IF( dta%ll_ht_i ) 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_i (ib,jl) = ht_i(ii,ij,jl) * tmask(ii,ij,1)
236	END DO
237	END DO
238	ENDIF
239	IF( dta%ll_ht_s ) THEN
240	igrd = 1
241	DO jl = 1, jpl
242	DO ib = 1, ilen1(igrd)
243	ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
244	ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
245	dta_bdy(ib_bdy)%ht_s (ib,jl) = ht_s(ii,ij,jl) * tmask(ii,ij,1)
246	END DO
247	END DO
248	ENDIF
249	ENDIF
250	#endif
251
252	ENDDO ! ib_bdy
253
254
255	ENDIF ! kt .eq. nit000
256
257	! update external data from files
258	!--------------------------------
259	lspr = ln_bdy_sio
260	jstart = 1
261	DO ib_bdy = 1, nb_bdy
262	dta => dta_bdy(ib_bdy)
263	IF( nn_dta(ib_bdy) .eq. 1 ) THEN ! skip this bit if no external data required
264
265	IF( PRESENT(jit) ) THEN
266	! Update barotropic boundary conditions only
267	! jit is optional argument for fld_read and bdytide_update
268	IF( cn_dyn2d(ib_bdy) /= 'none' ) THEN
269	IF( nn_dyn2d_dta(ib_bdy) .eq. 2 ) THEN ! tidal harmonic forcing ONLY: initialise arrays
270	IF( dta%ll_ssh ) dta%ssh(:) = 0.0
271	IF( dta%ll_u2d ) dta%u2d(:) = 0.0
272	IF( dta%ll_u3d ) dta%v2d(:) = 0.0
273	ENDIF
274	IF (cn_tra(ib_bdy) /= 'runoff') THEN
275	IF( nn_dyn2d_dta(ib_bdy) .EQ. 1 .OR. nn_dyn2d_dta(ib_bdy) .EQ. 3 ) THEN
276
277	jend = jstart + dta%nread(2) - 1
278	CALL fld_read( kt=kt, kn_fsbc=1, sd=bf(jstart:jend), map=nbmap_ptr(jstart:jend), &
279	& kit=jit, kt_offset=time_offset )
280
281	! If full velocities in boundary data then extract barotropic velocities from 3D fields
282	IF( ln_full_vel_array(ib_bdy) .AND. &
283	& ( nn_dyn2d_dta(ib_bdy) .EQ. 1 .OR. nn_dyn2d_dta(ib_bdy) .EQ. 3 .OR. &
284	& nn_dyn3d_dta(ib_bdy) .EQ. 1 ) )THEN
285
286	igrd = 2 ! zonal velocity
287	dta%u2d(:) = 0.0
288	DO ib = 1, idx_bdy(ib_bdy)%nblen(igrd)
289	ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
290	ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
291	DO ik = 1, jpkm1
292	dta%u2d(ib) = dta%u2d(ib) &
293	& + fse3u(ii,ij,ik) * umask(ii,ij,ik) * dta%u3d(ib,ik)
294	END DO
295	dta%u2d(ib) = dta%u2d(ib) * hur(ii,ij)
296	END DO
297	igrd = 3 ! meridional velocity
298	dta%v2d(:) = 0.0
299	DO ib = 1, idx_bdy(ib_bdy)%nblen(igrd)
300	ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
301	ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
302	DO ik = 1, jpkm1
303	dta%v2d(ib) = dta%v2d(ib) &
304	& + fse3v(ii,ij,ik) * vmask(ii,ij,ik) * dta%v3d(ib,ik)
305	END DO
306	dta%v2d(ib) = dta%v2d(ib) * hvr(ii,ij)
307	END DO
308	ENDIF
309	ENDIF
310	IF( nn_dyn2d_dta(ib_bdy) .ge. 2 ) THEN ! update tidal harmonic forcing
311	CALL bdytide_update( kt=kt, idx=idx_bdy(ib_bdy), dta=dta, td=tides(ib_bdy), &
312	& jit=jit, time_offset=time_offset )
313	ENDIF
314	ENDIF
315	ENDIF
316	ELSE
317	IF (cn_tra(ib_bdy) == 'runoff') then ! runoff condition
318	jend = nb_bdy_fld(ib_bdy)
319	CALL fld_read( kt=kt, kn_fsbc=1, sd=bf(jstart:jend), &
320	& map=nbmap_ptr(jstart:jend), kt_offset=time_offset )
321	!
322	igrd = 2 ! zonal 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%u2d(ib) = dta%u2d(ib) / ( e2u(ii,ij) * hu_0(ii,ij) )
327	END DO
328	!
329	igrd = 3 ! meridional velocity
330	DO ib = 1, idx_bdy(ib_bdy)%nblen(igrd)
331	ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
332	ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
333	dta%v2d(ib) = dta%v2d(ib) / ( e1v(ii,ij) * hv_0(ii,ij) )
334	END DO
335	ELSE
336	IF( nn_dyn2d_dta(ib_bdy) .eq. 2 ) THEN ! tidal harmonic forcing ONLY: initialise arrays
337	IF( dta%ll_ssh ) dta%ssh(:) = 0.0
338	IF( dta%ll_u2d ) dta%u2d(:) = 0.0
339	IF( dta%ll_v2d ) dta%v2d(:) = 0.0
340	ENDIF
341	IF( dta%nread(1) .gt. 0 ) THEN ! update external data
342	jend = jstart + dta%nread(1) - 1
343	CALL fld_read( kt=kt, kn_fsbc=1, sd=bf(jstart:jend), &
344	& map=nbmap_ptr(jstart:jend), kt_offset=time_offset )
345	ENDIF
346	! If full velocities in boundary data then split into barotropic and baroclinic data
347	IF( ln_full_vel_array(ib_bdy) .and. &
348	& ( nn_dyn2d_dta(ib_bdy) .EQ. 1 .OR. nn_dyn2d_dta(ib_bdy) .EQ. 3 .OR. &
349	& nn_dyn3d_dta(ib_bdy) .EQ. 1 ) ) THEN
350	igrd = 2 ! zonal velocity
351	dta%u2d(:) = 0.0
352	DO ib = 1, idx_bdy(ib_bdy)%nblen(igrd)
353	ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
354	ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
355	DO ik = 1, jpkm1
356	dta%u2d(ib) = dta%u2d(ib) &
357	& + fse3u(ii,ij,ik) * umask(ii,ij,ik) * dta%u3d(ib,ik)
358	END DO
359	dta%u2d(ib) = dta%u2d(ib) * hur(ii,ij)
360	DO ik = 1, jpkm1
361	dta%u3d(ib,ik) = dta%u3d(ib,ik) - dta%u2d(ib)
362	END DO
363	END DO
364	igrd = 3 ! meridional velocity
365	dta%v2d(:) = 0.0
366	DO ib = 1, idx_bdy(ib_bdy)%nblen(igrd)
367	ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
368	ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
369	DO ik = 1, jpkm1
370	dta%v2d(ib) = dta%v2d(ib) &
371	& + fse3v(ii,ij,ik) * vmask(ii,ij,ik) * dta%v3d(ib,ik)
372	END DO
373	dta%v2d(ib) = dta%v2d(ib) * hvr(ii,ij)
374	DO ik = 1, jpkm1
375	dta%v3d(ib,ik) = dta%v3d(ib,ik) - dta%v2d(ib)
376	END DO
377	END DO
378	ENDIF
379
380	ENDIF
381	#if defined key_lim3
382	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)
383	CALL lim_var_itd ( bf(jfld_hti)%fnow(:,1,1), bf(jfld_hts)%fnow(:,1,1), bf(jfld_ai)%fnow(:,1,1), &
384	& dta_bdy(ib_bdy)%ht_i, dta_bdy(ib_bdy)%ht_s, dta_bdy(ib_bdy)%a_i )
385	ENDIF
386	#endif
387	ENDIF
388	jstart = jstart + dta%nread(1)
389	END IF ! nn_dta(ib_bdy) = 1
390	END DO ! ib_bdy
391	lspr = .false.
392	! bg jchanut tschanges
393	#if defined key_tide
394	! Add tides if not split-explicit free surface else this is done in ts loop
395	IF (.NOT.lk_dynspg_ts) CALL bdy_dta_tides( kt=kt, time_offset=time_offset )
396	#endif
397	! end jchanut tschanges
398
399	IF ( ln_apr_obc ) THEN
400	DO ib_bdy = 1, nb_bdy
401	IF (cn_tra(ib_bdy) /= 'runoff')THEN
402	igrd = 1 ! meridional velocity
403	DO ib = 1, idx_bdy(ib_bdy)%nblenrim(igrd)
404	ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
405	ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
406	dta_bdy(ib_bdy)%ssh(ib) = dta_bdy(ib_bdy)%ssh(ib) + ssh_ib(ii,ij)
407	ENDDO
408	ENDIF
409	ENDDO
410	ENDIF
411
412	IF( nn_timing == 1 ) CALL timing_stop('bdy_dta')
413
414	END SUBROUTINE bdy_dta
415
416
417	SUBROUTINE bdy_dta_init
418	!!----------------------------------------------------------------------
419	!! * SUBROUTINE bdy_dta_init *
420	!!
421	!! ** Purpose : Initialise arrays for reading of external data
422	!! for open boundary conditions
423	!!
424	!! ** Method :
425	!!
426	!!----------------------------------------------------------------------
427	USE dynspg_oce, ONLY: lk_dynspg_ts
428	!!
429	INTEGER :: ib_bdy, jfld, jstart, jend, ierror ! local indices
430	INTEGER :: ios ! Local integer output status for namelist read
431	LOGICAL :: ln_bdy_sio
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, ln_bdy_sio
458	#if defined key_lim2
459	NAMELIST/nambdy_dta/ bn_frld, bn_hicif, bn_hsnif, ln_bdy_sio
460	#elif defined key_lim3
461	NAMELIST/nambdy_dta/ bn_a_i, bn_ht_i, bn_ht_s, ln_bdy_sio
462	#endif
463	NAMELIST/nambdy_dta/ ln_full_vel, ln_bdy_sio
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) .gt. 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) .eq. 1 .or. nn_dyn2d_dta(ib_bdy) .eq. 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) .eq. 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) .eq. 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) .eq. 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	ENDDO
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	IF(lwm) THEN
528	REWIND(numnam_ref)
529	REWIND(numnam_cfg)
530	ENDIF
531	jfld = 0
532	DO ib_bdy = 1, nb_bdy
533	IF( nn_dta(ib_bdy) .eq. 1 ) THEN
534	IF(lwm) THEN
535	READ ( numnam_ref, nambdy_dta, IOSTAT = ios, ERR = 901)
536	901 CONTINUE
537	ENDIF
538	call mpp_bcast(ios)
539	IF( ios /= 0 ) CALL ctl_nam ( ios , 'nambdy_dta in reference namelist', lwp )
540	IF(lwm) THEN
541	READ ( numnam_cfg, nambdy_dta, IOSTAT = ios, ERR = 902 )
542	902 CONTINUE
543	ENDIF
544	call mpp_bcast(ios)
545	IF( ios /= 0 ) CALL ctl_nam ( ios , 'nambdy_dta in configuration namelist', lwp )
546
547	IF(lwm) WRITE ( numond, nambdy_dta )
548
549	CALL dta_namelist(cn_dir, bn_tem, bn_sal, bn_u3d, bn_v3d, bn_ssh, bn_u2d, bn_v2d, ln_full_vel)
550
551	#if defined key_lim2
552	CALL lim_namelist(bn_frld, bn_hicif, bn_hsnif)
553	#elif defined key_lim3
554	CALL lim_namelist(bn_a_i, bn_ht_i, bn_ht_s)
555	#endif
556
557	cn_dir_array(ib_bdy) = cn_dir
558	ln_full_vel_array(ib_bdy) = ln_full_vel
559
560	nblen => idx_bdy(ib_bdy)%nblen
561	nblenrim => idx_bdy(ib_bdy)%nblenrim
562	dta => dta_bdy(ib_bdy)
563	dta%nread(2) = 0
564
565	! Only read in necessary fields for this set.
566	! Important that barotropic variables come first.
567	IF( nn_dyn2d_dta(ib_bdy) .eq. 1 .or. nn_dyn2d_dta(ib_bdy) .eq. 3 ) THEN
568
569	IF( dta%ll_ssh ) THEN
570	if(lwp) write(numout,*) '++++++ reading in ssh field'
571	jfld = jfld + 1
572	blf_i(jfld) = bn_ssh
573	ibdy(jfld) = ib_bdy
574	igrid(jfld) = 1
575	ilen1(jfld) = nblen(igrid(jfld))
576	ilen3(jfld) = 1
577	dta%nread(2) = dta%nread(2) + 1
578	ENDIF
579
580	IF( dta%ll_u2d .and. .not. ln_full_vel_array(ib_bdy) ) THEN
581	if(lwp) write(numout,*) '++++++ reading in u2d field'
582	jfld = jfld + 1
583	blf_i(jfld) = bn_u2d
584	ibdy(jfld) = ib_bdy
585	igrid(jfld) = 2
586	ilen1(jfld) = nblen(igrid(jfld))
587	ilen3(jfld) = 1
588	dta%nread(2) = dta%nread(2) + 1
589	ENDIF
590
591	IF( dta%ll_v2d .and. .not. ln_full_vel_array(ib_bdy) ) THEN
592	if(lwp) write(numout,*) '++++++ reading in v2d field'
593	jfld = jfld + 1
594	blf_i(jfld) = bn_v2d
595	ibdy(jfld) = ib_bdy
596	igrid(jfld) = 3
597	ilen1(jfld) = nblen(igrid(jfld))
598	ilen3(jfld) = 1
599	dta%nread(2) = dta%nread(2) + 1
600	ENDIF
601
602	ENDIF
603
604	! read 3D velocities if baroclinic velocities require OR if
605	! barotropic velocities required and ln_full_vel set to .true.
606	IF( nn_dyn3d_dta(ib_bdy) .eq. 1 .or. &
607	& ( ln_full_vel_array(ib_bdy) .and. ( nn_dyn2d_dta(ib_bdy) .eq. 1 .or. nn_dyn2d_dta(ib_bdy) .eq. 3 ) ) ) THEN
608
609	IF( dta%ll_u3d .or. ( ln_full_vel_array(ib_bdy) .and. dta%ll_u2d ) ) THEN
610	if(lwp) write(numout,*) '++++++ reading in u3d field'
611	jfld = jfld + 1
612	blf_i(jfld) = bn_u3d
613	ibdy(jfld) = ib_bdy
614	igrid(jfld) = 2
615	ilen1(jfld) = nblen(igrid(jfld))
616	ilen3(jfld) = jpk
617	IF( ln_full_vel_array(ib_bdy) .and. dta%ll_u2d ) dta%nread(2) = dta%nread(2) + 1
618	ENDIF
619
620	IF( dta%ll_v3d .or. ( ln_full_vel_array(ib_bdy) .and. dta%ll_v2d ) ) THEN
621	if(lwp) write(numout,*) '++++++ reading in v3d field'
622	jfld = jfld + 1
623	blf_i(jfld) = bn_v3d
624	ibdy(jfld) = ib_bdy
625	igrid(jfld) = 3
626	ilen1(jfld) = nblen(igrid(jfld))
627	ilen3(jfld) = jpk
628	IF( ln_full_vel_array(ib_bdy) .and. dta%ll_v2d ) dta%nread(2) = dta%nread(2) + 1
629	ENDIF
630
631	ENDIF
632
633	! temperature and salinity
634	IF( nn_tra_dta(ib_bdy) .eq. 1 ) THEN
635
636	IF( dta%ll_tem ) THEN
637	if(lwp) write(numout,*) '++++++ reading in tem field'
638	jfld = jfld + 1
639	blf_i(jfld) = bn_tem
640	ibdy(jfld) = ib_bdy
641	igrid(jfld) = 1
642	ilen1(jfld) = nblen(igrid(jfld))
643	ilen3(jfld) = jpk
644	ENDIF
645
646	IF( dta%ll_sal ) THEN
647	if(lwp) write(numout,*) '++++++ reading in sal field'
648	jfld = jfld + 1
649	blf_i(jfld) = bn_sal
650	ibdy(jfld) = ib_bdy
651	igrid(jfld) = 1
652	ilen1(jfld) = nblen(igrid(jfld))
653	ilen3(jfld) = jpk
654	ENDIF
655
656	ENDIF
657
658	#if defined key_lim2
659	! sea ice
660	IF( nn_ice_lim_dta(ib_bdy) .eq. 1 ) THEN
661
662	IF( dta%ll_frld ) THEN
663	jfld = jfld + 1
664	blf_i(jfld) = bn_frld
665	ibdy(jfld) = ib_bdy
666	igrid(jfld) = 1
667	ilen1(jfld) = nblen(igrid(jfld))
668	ilen3(jfld) = 1
669	ENDIF
670
671	IF( dta%ll_hicif ) THEN
672	jfld = jfld + 1
673	blf_i(jfld) = bn_hicif
674	ibdy(jfld) = ib_bdy
675	igrid(jfld) = 1
676	ilen1(jfld) = nblen(igrid(jfld))
677	ilen3(jfld) = 1
678	ENDIF
679
680	IF( dta%ll_hsnif ) THEN
681	jfld = jfld + 1
682	blf_i(jfld) = bn_hsnif
683	ibdy(jfld) = ib_bdy
684	igrid(jfld) = 1
685	ilen1(jfld) = nblen(igrid(jfld))
686	ilen3(jfld) = 1
687	ENDIF
688
689	ENDIF
690	#elif defined key_lim3
691	! sea ice
692	IF( nn_ice_lim_dta(ib_bdy) .eq. 1 ) THEN
693	! Test for types of ice input (lim2 or lim3)
694	! Build file name to find dimensions
695	clname=TRIM(bn_a_i%clname)
696	IF( .NOT. bn_a_i%ln_clim ) THEN
697	WRITE(clname, '(a,"_y",i4.4)' ) TRIM( bn_a_i%clname ), nyear ! add year
698	IF( bn_a_i%cltype /= 'yearly' ) WRITE(clname, '(a,"m" ,i2.2)' ) TRIM( clname ), nmonth ! add month
699	ELSE
700	IF( bn_a_i%cltype /= 'yearly' ) WRITE(clname, '(a,"_m",i2.2)' ) TRIM( bn_a_i%clname ), nmonth ! add month
701	ENDIF
702	IF( bn_a_i%cltype == 'daily' .OR. bn_a_i%cltype(1:4) == 'week' ) &
703	& WRITE(clname, '(a,"d" ,i2.2)' ) TRIM( clname ), nday ! add day
704	!
705	CALL iom_open ( clname, inum )
706	id1 = iom_varid( inum, bn_a_i%clvar, kdimsz=zdimsz, kndims=zndims, ldstop = .FALSE. )
707	CALL iom_close ( inum )
708
709	IF ( zndims == 4 ) THEN
710	ll_bdylim3 = .TRUE. ! lim3 input
711	ELSE
712	ll_bdylim3 = .FALSE. ! lim2 input
713	ENDIF
714	! End test
715
716	IF( dta%ll_a_i ) THEN
717	jfld = jfld + 1
718	blf_i(jfld) = bn_a_i
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	IF( dta%ll_ht_i ) THEN
726	jfld = jfld + 1
727	blf_i(jfld) = bn_ht_i
728	ibdy(jfld) = ib_bdy
729	igrid(jfld) = 1
730	ilen1(jfld) = nblen(igrid(jfld))
731	IF ( ll_bdylim3 ) THEN ; ilen3(jfld)=jpl ; ELSE ; ilen3(jfld)=1 ; ENDIF
732	ENDIF
733
734	IF( dta%ll_ht_s ) THEN
735	jfld = jfld + 1
736	blf_i(jfld) = bn_ht_s
737	ibdy(jfld) = ib_bdy
738	igrid(jfld) = 1
739	ilen1(jfld) = nblen(igrid(jfld))
740	IF ( ll_bdylim3 ) THEN ; ilen3(jfld)=jpl ; ELSE ; ilen3(jfld)=1 ; ENDIF
741	ENDIF
742
743	ENDIF
744	#endif
745	! Recalculate field counts
746	!-------------------------
747	IF( ib_bdy .eq. 1 ) THEN
748	nb_bdy_fld_sum = 0
749	nb_bdy_fld(ib_bdy) = jfld
750	nb_bdy_fld_sum = jfld
751	ELSE
752	nb_bdy_fld(ib_bdy) = jfld - nb_bdy_fld_sum
753	nb_bdy_fld_sum = nb_bdy_fld_sum + nb_bdy_fld(ib_bdy)
754	ENDIF
755
756	dta%nread(1) = nb_bdy_fld(ib_bdy)
757
758	ENDIF ! nn_dta .eq. 1
759	ENDDO ! ib_bdy
760
761	DO jfld = 1, nb_bdy_fld_sum
762	ALLOCATE( bf(jfld)%fnow(ilen1(jfld),1,ilen3(jfld)) )
763	IF( blf_i(jfld)%ln_tint ) ALLOCATE( bf(jfld)%fdta(ilen1(jfld),1,ilen3(jfld),2) )
764	nbmap_ptr(jfld)%ptr => idx_bdy(ibdy(jfld))%nbmap(:,igrid(jfld))
765	nbmap_ptr(jfld)%ll_unstruc = ln_coords_file(ibdy(jfld))
766	ENDDO
767
768	! fill bf with blf_i and control print
769	!-------------------------------------
770	jstart = 1
771	DO ib_bdy = 1, nb_bdy
772	jend = jstart - 1 + nb_bdy_fld(ib_bdy)
773	CALL fld_fill( bf(jstart:jend), blf_i(jstart:jend), cn_dir_array(ib_bdy), 'bdy_dta', &
774	& 'open boundary conditions', 'nambdy_dta' )
775	jstart = jend + 1
776	ENDDO
777
778	! Initialise local boundary data arrays
779	! nn_xxx_dta=0 : allocate space - will be filled from initial conditions later
780	! nn_xxx_dta=1 : point to "fnow" arrays
781	!-------------------------------------
782
783	jfld = 0
784	DO ib_bdy=1, nb_bdy
785
786	nblen => idx_bdy(ib_bdy)%nblen
787	dta => dta_bdy(ib_bdy)
788
789	if(lwp) then
790	write(numout,*) '++++++ dta%ll_ssh = ',dta%ll_ssh
791	write(numout,*) '++++++ dta%ll_u2d = ',dta%ll_u2d
792	write(numout,*) '++++++ dta%ll_v2d = ',dta%ll_v2d
793	write(numout,*) '++++++ dta%ll_u3d = ',dta%ll_u3d
794	write(numout,*) '++++++ dta%ll_v3d = ',dta%ll_v3d
795	write(numout,*) '++++++ dta%ll_tem = ',dta%ll_tem
796	write(numout,*) '++++++ dta%ll_sal = ',dta%ll_sal
797	endif
798
799	IF ( nn_dyn2d_dta(ib_bdy) .eq. 0 .or. nn_dyn2d_dta(ib_bdy) .eq. 2 ) THEN
800	if(lwp) write(numout,*) '++++++ dta%ssh/u2d/u3d allocated space'
801	IF( dta%ll_ssh ) ALLOCATE( dta%ssh(nblen(1)) )
802	IF( dta%ll_u2d ) ALLOCATE( dta%u2d(nblen(2)) )
803	IF( dta%ll_v2d ) ALLOCATE( dta%v2d(nblen(3)) )
804	ENDIF
805	IF ( nn_dyn2d_dta(ib_bdy) .eq. 1 .or. nn_dyn2d_dta(ib_bdy) .eq. 3 ) THEN
806	IF( dta%ll_ssh ) THEN
807	if(lwp) write(numout,*) '++++++ dta%ssh pointing to fnow'
808	jfld = jfld + 1
809	dta%ssh => bf(jfld)%fnow(:,1,1)
810	ENDIF
811	IF ( dta%ll_u2d ) THEN
812	IF ( ln_full_vel_array(ib_bdy) ) THEN
813	if(lwp) write(numout,*) '++++++ dta%u2d allocated space'
814	ALLOCATE( dta%u2d(nblen(2)) )
815	ELSE
816	if(lwp) write(numout,*) '++++++ dta%u2d pointing to fnow'
817	jfld = jfld + 1
818	dta%u2d => bf(jfld)%fnow(:,1,1)
819	ENDIF
820	ENDIF
821	IF ( dta%ll_v2d ) THEN
822	IF ( ln_full_vel_array(ib_bdy) ) THEN
823	if(lwp) write(numout,*) '++++++ dta%v2d allocated space'
824	ALLOCATE( dta%v2d(nblen(3)) )
825	ELSE
826	if(lwp) write(numout,*) '++++++ dta%v2d pointing to fnow'
827	jfld = jfld + 1
828	dta%v2d => bf(jfld)%fnow(:,1,1)
829	ENDIF
830	ENDIF
831	ENDIF
832
833	IF ( nn_dyn3d_dta(ib_bdy) .eq. 0 ) THEN
834	if(lwp) write(numout,*) '++++++ dta%u3d/v3d allocated space'
835	IF( dta%ll_u3d ) ALLOCATE( dta_bdy(ib_bdy)%u3d(nblen(2),jpk) )
836	IF( dta%ll_v3d ) ALLOCATE( dta_bdy(ib_bdy)%v3d(nblen(3),jpk) )
837	ENDIF
838	IF ( nn_dyn3d_dta(ib_bdy) .eq. 1 .or. &
839	& ( ln_full_vel_array(ib_bdy) .and. ( nn_dyn2d_dta(ib_bdy) .eq. 1 .or. nn_dyn2d_dta(ib_bdy) .eq. 3 ) ) ) THEN
840	IF ( dta%ll_u3d .or. ( ln_full_vel_array(ib_bdy) .and. dta%ll_u2d ) ) THEN
841	if(lwp) write(numout,*) '++++++ dta%u3d pointing to fnow'
842	jfld = jfld + 1
843	dta_bdy(ib_bdy)%u3d => bf(jfld)%fnow(:,1,:)
844	ENDIF
845	IF ( dta%ll_v3d .or. ( ln_full_vel_array(ib_bdy) .and. dta%ll_v2d ) ) THEN
846	if(lwp) write(numout,*) '++++++ dta%v3d pointing to fnow'
847	jfld = jfld + 1
848	dta_bdy(ib_bdy)%v3d => bf(jfld)%fnow(:,1,:)
849	ENDIF
850	ENDIF
851
852	IF( nn_tra_dta(ib_bdy) .eq. 0 ) THEN
853	if(lwp) write(numout,*) '++++++ dta%tem/sal allocated space'
854	IF( dta%ll_tem ) ALLOCATE( dta_bdy(ib_bdy)%tem(nblen(1),jpk) )
855	IF( dta%ll_sal ) ALLOCATE( dta_bdy(ib_bdy)%sal(nblen(1),jpk) )
856	ELSE
857	IF( dta%ll_tem ) THEN
858	if(lwp) write(numout,*) '++++++ dta%tem pointing to fnow'
859	jfld = jfld + 1
860	dta_bdy(ib_bdy)%tem => bf(jfld)%fnow(:,1,:)
861	ENDIF
862	IF( dta%ll_sal ) THEN
863	if(lwp) write(numout,*) '++++++ dta%sal pointing to fnow'
864	jfld = jfld + 1
865	dta_bdy(ib_bdy)%sal => bf(jfld)%fnow(:,1,:)
866	ENDIF
867	ENDIF
868
869	#if defined key_lim2
870	IF (cn_ice_lim(ib_bdy) /= 'none') THEN
871	IF( nn_ice_lim_dta(ib_bdy) .eq. 0 ) THEN
872	ALLOCATE( dta_bdy(ib_bdy)%frld(nblen(1)) )
873	ALLOCATE( dta_bdy(ib_bdy)%hicif(nblen(1)) )
874	ALLOCATE( dta_bdy(ib_bdy)%hsnif(nblen(1)) )
875	ELSE
876	jfld = jfld + 1
877	dta_bdy(ib_bdy)%frld => bf(jfld)%fnow(:,1,1)
878	jfld = jfld + 1
879	dta_bdy(ib_bdy)%hicif => bf(jfld)%fnow(:,1,1)
880	jfld = jfld + 1
881	dta_bdy(ib_bdy)%hsnif => bf(jfld)%fnow(:,1,1)
882	ENDIF
883	ENDIF
884	#elif defined key_lim3
885	IF (cn_ice_lim(ib_bdy) /= 'none') THEN
886	IF( nn_ice_lim_dta(ib_bdy) .eq. 0 ) THEN
887	ALLOCATE( dta_bdy(ib_bdy)%a_i (nblen(1),jpl) )
888	ALLOCATE( dta_bdy(ib_bdy)%ht_i(nblen(1),jpl) )
889	ALLOCATE( dta_bdy(ib_bdy)%ht_s(nblen(1),jpl) )
890	ELSE
891	IF ( ll_bdylim3 ) THEN ! case input is lim3 type
892	jfld = jfld + 1
893	dta_bdy(ib_bdy)%a_i => bf(jfld)%fnow(:,1,:)
894	jfld = jfld + 1
895	dta_bdy(ib_bdy)%ht_i => bf(jfld)%fnow(:,1,:)
896	jfld = jfld + 1
897	dta_bdy(ib_bdy)%ht_s => bf(jfld)%fnow(:,1,:)
898	ELSE ! case input is lim2 type
899	jfld_ai = jfld + 1
900	jfld_hti = jfld + 2
901	jfld_hts = jfld + 3
902	jfld = jfld + 3
903	ALLOCATE( dta_bdy(ib_bdy)%a_i (nblen(1),jpl) )
904	ALLOCATE( dta_bdy(ib_bdy)%ht_i(nblen(1),jpl) )
905	ALLOCATE( dta_bdy(ib_bdy)%ht_s(nblen(1),jpl) )
906	dta_bdy(ib_bdy)%a_i (:,:) = 0.0
907	dta_bdy(ib_bdy)%ht_i(:,:) = 0.0
908	dta_bdy(ib_bdy)%ht_s(:,:) = 0.0
909	ENDIF
910
911	ENDIF
912	ENDIF
913	#endif
914
915	ENDDO ! ib_bdy
916
917	IF( nn_timing == 1 ) CALL timing_stop('bdy_dta_init')
918
919	END SUBROUTINE bdy_dta_init
920
921	SUBROUTINE dta_namelist(cd_dir, bd_tem, bd_sal, bd_u3d, bd_v3d, bd_ssh, &
922	& bd_u2d, bd_v2d, ln_full_vel)
923	!!---------------------------------------------------------------------
924	!! * ROUTINE dta_namelist *
925	!!
926	!! ** Purpose : Broadcast namelist variables read by procesor lwm
927	!!
928	!! ** Method : use lib_mpp
929	!!----------------------------------------------------------------------
930	CHARACTER(len=100) :: cd_dir
931	TYPE(FLD_N) :: bd_tem, bd_sal, bd_u3d, bd_v3d !
932	TYPE(FLD_N) :: bd_ssh, bd_u2d, bd_v2d
933	#if defined key_mpp_mpi
934	CALL mpp_bcast(cd_dir, 100)
935	CALL fld_n_bcast(bd_tem)
936	CALL fld_n_bcast(bd_sal)
937	CALL fld_n_bcast(bd_u3d)
938	CALL fld_n_bcast(bd_v3d)
939	CALL fld_n_bcast(bd_ssh)
940	CALL fld_n_bcast(bd_u2d)
941	CALL fld_n_bcast(bd_v2d)
942	CALL mpp_bcast(ln_bdy_sio)
943	CALL mpp_bcast(ln_full_vel)
944	#endif
945	END SUBROUTINE dta_namelist
946
947	SUBROUTINE lim_namelist(bd_1, bd_2, bd_3)
948	!!---------------------------------------------------------------------
949	!! * ROUTINE lim_namelist *
950	!!
951	!! ** Purpose : Broadcast namelist variables read by procesor lwm
952	!!
953	!! ** Method : use lib_mpp
954	!!----------------------------------------------------------------------
955	TYPE(FLD_N) :: bd_1, bd_2, bd_3
956	#if defined key_mpp_mpi
957	CALL fld_n_bcast(bd_1)
958	CALL fld_n_bcast(bd_2)
959	CALL fld_n_bcast(bd_3)
960	#endif
961	END SUBROUTINE lim_namelist
962
963	#else
964	!!----------------------------------------------------------------------
965	!! Dummy module NO Open Boundary Conditions
966	!!----------------------------------------------------------------------
967	CONTAINS
968	SUBROUTINE bdy_dta( kt, jit, time_offset ) ! Empty routine
969	INTEGER, INTENT( in ) :: kt
970	INTEGER, INTENT( in ), OPTIONAL :: jit
971	INTEGER, INTENT( in ), OPTIONAL :: time_offset
972	WRITE(,) 'bdy_dta: You should not have seen this print! error?', kt
973	END SUBROUTINE bdy_dta
974	SUBROUTINE bdy_dta_init() ! Empty routine
975	WRITE(,) 'bdy_dta_init: You should not have seen this print! error?'
976	END SUBROUTINE bdy_dta_init
977	#endif
978
979	!!==============================================================================
980	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/test_moci_test_suite_namelist_read/NEMOGCM/NEMO/OPA_SRC/BDY/bdydta.F90 @ 9383

Download in other formats: