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

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.

bdydta.F90 in branches/2015/nemo_v3_6_STABLE/NEMOGCM/NEMO/OPA_SRC/BDY – NEMO

Context Navigation

source: branches/2015/nemo_v3_6_STABLE/NEMOGCM/NEMO/OPA_SRC/BDY/bdydta.F90 @ 8523

Last change on this file since 8523 was 8523, checked in by cbricaud, 7 years ago
bugfix in nemo_v3_6_STABLE for ticket #1922
Property svn:keywords set to `Id`
File size: 40.3 KB

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

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

Download in other formats: