Context Navigation

bdydta.F90 @ 6623

Last change on this file since 6623 was 6140, checked in by timgraham, 8 years ago

Merge of branches/2015/dev_merge_2015 back into trunk. Merge excludes NEMOGCM/TOOLS/OBSTOOLS/ for now due to issues with the change of file type. Will sort these manually with further commits.

Branch merged as follows:
In the working copy of branch ran:
svn merge svn+ssh://forge.ipsl.jussieu.fr/ipsl/forge/projets/nemo/svn/trunk@HEAD
Small conflicts due to bug fixes applied to trunk since the dev_merge_2015 was copied. Bug fixes were applied to the branch as well so these were easy to resolve.
Branch committed at this stage

In working copy run:
svn switch svn+ssh://forge.ipsl.jussieu.fr/ipsl/forge/projets/nemo/svn/trunk
to switch working copy

Run:
svn merge --reintegrate svn+ssh://forge.ipsl.jussieu.fr/ipsl/forge/projets/nemo/svn/branches/2015/dev_merge_2015
to merge the branch into the trunk and then commit - no conflicts at this stage.

Property svn:keywords set to Id

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

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

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

Context Navigation

source: trunk/NEMOGCM/NEMO/OPA_SRC/BDY/bdydta.F90 @ 6623

Download in other formats: