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 NEMO/branches/2020/ticket2487/src/OCE/BDY – NEMO

Context Navigation

source: NEMO/branches/2020/ticket2487/src/OCE/BDY/bdydta.F90 @ 15410

Last change on this file since 15410 was 15410, checked in by smueller, 3 years ago
Synchronizing with ^{/NEMO/releases/r4.0/r4.0-HEAD@15405 (ticket #2487)}
Property svn:keywords set to `Id`
File size: 42.2 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
[9656]	13	!! 3.6 ! 2012-01 (C. Rousset) add ice boundary conditions for sea ice
	14	!! 4.0 ! 2018 (C. Rousset) SI3 compatibility
[1125]	15	!!----------------------------------------------------------------------
[9019]	16
[1125]	17	!!----------------------------------------------------------------------
[9019]	18	!! bdy_dta : read external data along open boundaries from file
	19	!! bdy_dta_init : initialise arrays etc for reading of external data
	20	!!----------------------------------------------------------------------
	21	USE oce ! ocean dynamics and tracers
	22	USE dom_oce ! ocean space and time domain
	23	USE phycst ! physical constants
	24	USE sbcapr ! atmospheric pressure forcing
	25	USE sbctide ! Tidal forcing or not
	26	USE bdy_oce ! ocean open boundary conditions
	27	USE bdytides ! tidal forcing at boundaries
[9570]	28	#if defined key_si3
[9019]	29	USE ice ! sea-ice variables
	30	USE icevar ! redistribute ice input into categories
[2528]	31	#endif
[9019]	32	!
[10529]	33	USE lib_mpp, ONLY: ctl_stop, ctl_nam
[9019]	34	USE fldread ! read input fields
	35	USE iom ! IOM library
	36	USE in_out_manager ! I/O logical units
	37	USE timing ! Timing
[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
[15264]	45	INTEGER , PARAMETER :: jpbdyfld = 17 ! maximum number of files to read
[11536]	46	INTEGER , PARAMETER :: jp_bdyssh = 1 !
	47	INTEGER , PARAMETER :: jp_bdyu2d = 2 !
	48	INTEGER , PARAMETER :: jp_bdyv2d = 3 !
	49	INTEGER , PARAMETER :: jp_bdyu3d = 4 !
	50	INTEGER , PARAMETER :: jp_bdyv3d = 5 !
	51	INTEGER , PARAMETER :: jp_bdytem = 6 !
	52	INTEGER , PARAMETER :: jp_bdysal = 7 !
	53	INTEGER , PARAMETER :: jp_bdya_i = 8 !
	54	INTEGER , PARAMETER :: jp_bdyh_i = 9 !
	55	INTEGER , PARAMETER :: jp_bdyh_s = 10 !
	56	INTEGER , PARAMETER :: jp_bdyt_i = 11 !
	57	INTEGER , PARAMETER :: jp_bdyt_s = 12 !
	58	INTEGER , PARAMETER :: jp_bdytsu = 13 !
	59	INTEGER , PARAMETER :: jp_bdys_i = 14 !
	60	INTEGER , PARAMETER :: jp_bdyaip = 15 !
	61	INTEGER , PARAMETER :: jp_bdyhip = 16 !
[15264]	62	INTEGER , PARAMETER :: jp_bdyhil = 17 !
[11536]	63	#if ! defined key_si3
	64	INTEGER , PARAMETER :: jpl = 1
	65	#endif
	66
[4354]	67	!$AGRIF_DO_NOT_TREAT
[11536]	68	TYPE(FLD), PUBLIC, ALLOCATABLE, DIMENSION(:,:), TARGET :: bf ! structure of input fields (file informations, fields read)
[4354]	69	!$AGRIF_END_DO_NOT_TREAT
[911]	70
[1125]	71	!!----------------------------------------------------------------------
[9598]	72	!! NEMO/OCE 4.0 , NEMO Consortium (2018)
[1146]	73	!! $Id$
[10068]	74	!! Software governed by the CeCILL license (see ./LICENSE)
[1125]	75	!!----------------------------------------------------------------------
[911]	76	CONTAINS
	77
[11536]	78	SUBROUTINE bdy_dta( kt, kit, kt_offset )
[1125]	79	!!----------------------------------------------------------------------
[3294]	80	!! * SUBROUTINE bdy_dta *
[911]	81	!!
[3294]	82	!! ** Purpose : Update external data for open boundary conditions
[911]	83	!!
[3294]	84	!! ** Method : Use fldread.F90
	85	!!
[1125]	86	!!----------------------------------------------------------------------
[6140]	87	INTEGER, INTENT(in) :: kt ! ocean time-step index
[11536]	88	INTEGER, INTENT(in), OPTIONAL :: kit ! subcycle time-step index (for timesplitting option)
	89	INTEGER, INTENT(in), OPTIONAL :: kt_offset ! time offset in units of timesteps. NB. if kit
[6140]	90	! ! is present then units = subcycle timesteps.
[11536]	91	! ! kt_offset = 0 => get data at "now" time level
	92	! ! kt_offset = -1 => get data at "before" time level
	93	! ! kt_offset = +1 => get data at "after" time level
[6140]	94	! ! etc.
	95	!
[11536]	96	INTEGER :: jbdy, jfld, jstart, jend, ib, jl ! dummy loop indices
	97	INTEGER :: ii, ij, ik, igrd, ipl ! local integers
	98	TYPE(OBC_DATA) , POINTER :: dta_alias ! short cut
	99	TYPE(FLD), DIMENSION(:), POINTER :: bf_alias
[911]	100	!!---------------------------------------------------------------------------
[6140]	101	!
[9124]	102	IF( ln_timing ) CALL timing_start('bdy_dta')
[6140]	103	!
[3294]	104	! Initialise data arrays once for all from initial conditions where required
	105	!---------------------------------------------------------------------------
[11536]	106	IF( kt == nit000 .AND. .NOT.PRESENT(kit) ) THEN
[1125]	107
[3294]	108	! Calculate depth-mean currents
	109	!-----------------------------
[11536]	110
[9019]	111	DO jbdy = 1, nb_bdy
[6140]	112	!
[9019]	113	IF( nn_dyn2d_dta(jbdy) == 0 ) THEN
[11536]	114	IF( dta_bdy(jbdy)%lneed_ssh ) THEN
[4292]	115	igrd = 1
[12548]	116	DO ib = 1, idx_bdy(jbdy)%nblenrim(igrd) ! ssh is allocated and used only on the rim
[9019]	117	ii = idx_bdy(jbdy)%nbi(ib,igrd)
	118	ij = idx_bdy(jbdy)%nbj(ib,igrd)
	119	dta_bdy(jbdy)%ssh(ib) = sshn(ii,ij) * tmask(ii,ij,1)
[11536]	120	END DO
[9019]	121	ENDIF
[12910]	122	IF( ASSOCIATED(dta_bdy(jbdy)%u2d) ) THEN ! no SIZE with a unassociated pointer. v2d and u2d can differ on subdomain
[4292]	123	igrd = 2
[12910]	124	DO ib = 1, SIZE(dta_bdy(jbdy)%u2d) ! u2d is used either over the whole bdy or only on the rim
[9019]	125	ii = idx_bdy(jbdy)%nbi(ib,igrd)
	126	ij = idx_bdy(jbdy)%nbj(ib,igrd)
	127	dta_bdy(jbdy)%u2d(ib) = un_b(ii,ij) * umask(ii,ij,1)
[11536]	128	END DO
[12910]	129	ENDIF
	130	IF( ASSOCIATED(dta_bdy(jbdy)%v2d) ) THEN ! no SIZE with a unassociated pointer. v2d and u2d can differ on subdomain
[4292]	131	igrd = 3
[12910]	132	DO ib = 1, SIZE(dta_bdy(jbdy)%v2d) ! v2d is used either over the whole bdy or only on the rim
[9019]	133	ii = idx_bdy(jbdy)%nbi(ib,igrd)
	134	ij = idx_bdy(jbdy)%nbj(ib,igrd)
	135	dta_bdy(jbdy)%v2d(ib) = vn_b(ii,ij) * vmask(ii,ij,1)
[11536]	136	END DO
[9019]	137	ENDIF
[4292]	138	ENDIF
[9019]	139	!
	140	IF( nn_dyn3d_dta(jbdy) == 0 ) THEN
[11536]	141	IF( dta_bdy(jbdy)%lneed_dyn3d ) THEN
[4292]	142	igrd = 2
[12548]	143	DO ib = 1, idx_bdy(jbdy)%nblen(igrd)
[4292]	144	DO ik = 1, jpkm1
[9019]	145	ii = idx_bdy(jbdy)%nbi(ib,igrd)
	146	ij = idx_bdy(jbdy)%nbj(ib,igrd)
	147	dta_bdy(jbdy)%u3d(ib,ik) = ( un(ii,ij,ik) - un_b(ii,ij) ) * umask(ii,ij,ik)
[3294]	148	END DO
[11536]	149	END DO
[4292]	150	igrd = 3
[12548]	151	DO ib = 1, idx_bdy(jbdy)%nblen(igrd)
[4292]	152	DO ik = 1, jpkm1
[9019]	153	ii = idx_bdy(jbdy)%nbi(ib,igrd)
	154	ij = idx_bdy(jbdy)%nbj(ib,igrd)
	155	dta_bdy(jbdy)%v3d(ib,ik) = ( vn(ii,ij,ik) - vn_b(ii,ij) ) * vmask(ii,ij,ik)
[11536]	156	END DO
	157	END DO
[9019]	158	ENDIF
[3294]	159	ENDIF
[911]	160
[9019]	161	IF( nn_tra_dta(jbdy) == 0 ) THEN
[11536]	162	IF( dta_bdy(jbdy)%lneed_tra ) THEN
[4292]	163	igrd = 1
[12548]	164	DO ib = 1, idx_bdy(jbdy)%nblen(igrd)
[4292]	165	DO ik = 1, jpkm1
[9019]	166	ii = idx_bdy(jbdy)%nbi(ib,igrd)
	167	ij = idx_bdy(jbdy)%nbj(ib,igrd)
[12049]	168	dta_bdy(jbdy)%tem(ib,ik) = tsn(ii,ij,ik,jp_tem) * tmask(ii,ij,ik)
	169	dta_bdy(jbdy)%sal(ib,ik) = tsn(ii,ij,ik,jp_sal) * tmask(ii,ij,ik)
[4292]	170	END DO
[11536]	171	END DO
[9019]	172	ENDIF
[3294]	173	ENDIF
[911]	174
[9570]	175	#if defined key_si3
[9657]	176	IF( nn_ice_dta(jbdy) == 0 ) THEN ! set ice to initial values
[11536]	177	IF( dta_bdy(jbdy)%lneed_ice ) THEN
[4292]	178	igrd = 1
	179	DO jl = 1, jpl
[12548]	180	DO ib = 1, idx_bdy(jbdy)%nblen(igrd)
[9019]	181	ii = idx_bdy(jbdy)%nbi(ib,igrd)
	182	ij = idx_bdy(jbdy)%nbj(ib,igrd)
[11536]	183	dta_bdy(jbdy)%a_i(ib,jl) = a_i (ii,ij,jl) * tmask(ii,ij,1)
	184	dta_bdy(jbdy)%h_i(ib,jl) = h_i (ii,ij,jl) * tmask(ii,ij,1)
	185	dta_bdy(jbdy)%h_s(ib,jl) = h_s (ii,ij,jl) * tmask(ii,ij,1)
	186	dta_bdy(jbdy)%t_i(ib,jl) = SUM(t_i (ii,ij,:,jl)) * r1_nlay_i * tmask(ii,ij,1)
	187	dta_bdy(jbdy)%t_s(ib,jl) = SUM(t_s (ii,ij,:,jl)) * r1_nlay_s * tmask(ii,ij,1)
	188	dta_bdy(jbdy)%tsu(ib,jl) = t_su(ii,ij,jl) * tmask(ii,ij,1)
	189	dta_bdy(jbdy)%s_i(ib,jl) = s_i (ii,ij,jl) * tmask(ii,ij,1)
	190	! melt ponds
	191	dta_bdy(jbdy)%aip(ib,jl) = a_ip(ii,ij,jl) * tmask(ii,ij,1)
	192	dta_bdy(jbdy)%hip(ib,jl) = h_ip(ii,ij,jl) * tmask(ii,ij,1)
[15264]	193	dta_bdy(jbdy)%hil(ib,jl) = h_il(ii,ij,jl) * tmask(ii,ij,1)
[4292]	194	END DO
	195	END DO
	196	ENDIF
	197	ENDIF
[3294]	198	#endif
[9019]	199	END DO ! jbdy
[6140]	200	!
	201	ENDIF ! kt == nit000
[911]	202
[3294]	203	! update external data from files
	204	!--------------------------------
[3703]	205
[11536]	206	DO jbdy = 1, nb_bdy
[3703]	207
[11536]	208	dta_alias => dta_bdy(jbdy)
	209	bf_alias => bf(:,jbdy)
[3703]	210
[11536]	211	! read/update all bdy data
	212	! ------------------------
	213	CALL fld_read( kt, 1, bf_alias, kit = kit, kt_offset = kt_offset )
[4292]	214
[11536]	215	! apply some corrections in some specific cases...
	216	! --------------------------------------------------
	217	!
	218	! if runoff condition: change river flow we read (in m3/s) into barotropic velocity (m/s)
[12910]	219	IF( cn_tra(jbdy) == 'runoff' ) THEN ! runoff
[11536]	220	!
[12910]	221	IF( ASSOCIATED(dta_bdy(jbdy)%u2d) ) THEN ! no SIZE with a unassociated pointer. v2d and u2d can differ on subdomain
	222	igrd = 2 ! zonal flow (m3/s) to barotropic zonal velocity (m/s)
	223	DO ib = 1, SIZE(dta_alias%u2d) ! u2d is used either over the whole bdy or only on the rim
	224	ii = idx_bdy(jbdy)%nbi(ib,igrd)
	225	ij = idx_bdy(jbdy)%nbj(ib,igrd)
	226	dta_alias%u2d(ib) = dta_alias%u2d(ib) / ( e2u(ii,ij) * hu_0(ii,ij) )
	227	END DO
	228	ENDIF
	229	IF( ASSOCIATED(dta_bdy(jbdy)%v2d) ) THEN ! no SIZE with a unassociated pointer. v2d and u2d can differ on subdomain
	230	igrd = 3 ! meridional flow (m3/s) to barotropic meridional velocity (m/s)
	231	DO ib = 1, SIZE(dta_alias%v2d) ! v2d is used either over the whole bdy or only on the rim
	232	ii = idx_bdy(jbdy)%nbi(ib,igrd)
	233	ij = idx_bdy(jbdy)%nbj(ib,igrd)
	234	dta_alias%v2d(ib) = dta_alias%v2d(ib) / ( e1v(ii,ij) * hv_0(ii,ij) )
	235	END DO
	236	ENDIF
[11536]	237	ENDIF
	238
	239	! tidal harmonic forcing ONLY: initialise arrays
	240	IF( nn_dyn2d_dta(jbdy) == 2 ) THEN ! we did not read ssh, u/v2d
[12910]	241	IF( ASSOCIATED(dta_alias%ssh) ) dta_alias%ssh(:) = 0._wp
	242	IF( ASSOCIATED(dta_alias%u2d) ) dta_alias%u2d(:) = 0._wp
	243	IF( ASSOCIATED(dta_alias%v2d) ) dta_alias%v2d(:) = 0._wp
[11536]	244	ENDIF
	245
	246	! If full velocities in boundary data, then split it into barotropic and baroclinic component
	247	IF( bf_alias(jp_bdyu3d)%ltotvel ) THEN ! if we read 3D total velocity (can be true only if u3d was read)
	248	igrd = 2 ! zonal velocity
	249	DO ib = 1, idx_bdy(jbdy)%nblen(igrd)
	250	ii = idx_bdy(jbdy)%nbi(ib,igrd)
	251	ij = idx_bdy(jbdy)%nbj(ib,igrd)
[12639]	252	dta_alias%u2d(ib) = 0._wp ! compute barotrope zonal velocity and put it in u2d
[11536]	253	DO ik = 1, jpkm1
	254	dta_alias%u2d(ib) = dta_alias%u2d(ib) + e3u_n(ii,ij,ik) * umask(ii,ij,ik) * dta_alias%u3d(ib,ik)
	255	END DO
	256	dta_alias%u2d(ib) = dta_alias%u2d(ib) * r1_hu_n(ii,ij)
	257	DO ik = 1, jpkm1 ! compute barocline zonal velocity and put it in u3d
	258	dta_alias%u3d(ib,ik) = dta_alias%u3d(ib,ik) - dta_alias%u2d(ib)
	259	END DO
	260	END DO
[15410]	261	ENDIF ! ltotvel
	262	IF( bf_alias(jp_bdyv3d)%ltotvel ) THEN ! if we read 3D total velocity (can be true only if u3d was read)
[11536]	263	igrd = 3 ! meridional velocity
	264	DO ib = 1, idx_bdy(jbdy)%nblen(igrd)
	265	ii = idx_bdy(jbdy)%nbi(ib,igrd)
	266	ij = idx_bdy(jbdy)%nbj(ib,igrd)
[12639]	267	dta_alias%v2d(ib) = 0._wp ! compute barotrope meridional velocity and put it in v2d
[11536]	268	DO ik = 1, jpkm1
	269	dta_alias%v2d(ib) = dta_alias%v2d(ib) + e3v_n(ii,ij,ik) * vmask(ii,ij,ik) * dta_alias%v3d(ib,ik)
	270	END DO
	271	dta_alias%v2d(ib) = dta_alias%v2d(ib) * r1_hv_n(ii,ij)
	272	DO ik = 1, jpkm1 ! compute barocline meridional velocity and put it in v3d
	273	dta_alias%v3d(ib,ik) = dta_alias%v3d(ib,ik) - dta_alias%v2d(ib)
	274	END DO
	275	END DO
	276	ENDIF ! ltotvel
	277
	278	! update tidal harmonic forcing
	279	IF( PRESENT(kit) .AND. nn_dyn2d_dta(jbdy) .GE. 2 ) THEN
	280	CALL bdytide_update( kt = kt, idx = idx_bdy(jbdy), dta = dta_alias, td = tides(jbdy), &
	281	& kit = kit, kt_offset = kt_offset )
	282	ENDIF
	283
	284	! atm surface pressure : add inverted barometer effect to ssh if it was read
	285	IF ( ln_apr_obc .AND. TRIM(bf_alias(jp_bdyssh)%clrootname) /= 'NOT USED' ) THEN
	286	igrd = 1
	287	DO ib = 1, idx_bdy(jbdy)%nblenrim(igrd) ! ssh is used only on the rim
	288	ii = idx_bdy(jbdy)%nbi(ib,igrd)
	289	ij = idx_bdy(jbdy)%nbj(ib,igrd)
	290	dta_alias%ssh(ib) = dta_alias%ssh(ib) + ssh_ib(ii,ij)
	291	END DO
	292	ENDIF
	293
[9570]	294	#if defined key_si3
[12639]	295	IF( dta_alias%lneed_ice .AND. idx_bdy(jbdy)%nblen(1) > 0 ) THEN
[11536]	296	! fill temperature and salinity arrays
	297	IF( TRIM(bf_alias(jp_bdyt_i)%clrootname) == 'NOT USED' ) bf_alias(jp_bdyt_i)%fnow(:,1,:) = rice_tem (jbdy)
	298	IF( TRIM(bf_alias(jp_bdyt_s)%clrootname) == 'NOT USED' ) bf_alias(jp_bdyt_s)%fnow(:,1,:) = rice_tem (jbdy)
	299	IF( TRIM(bf_alias(jp_bdytsu)%clrootname) == 'NOT USED' ) bf_alias(jp_bdytsu)%fnow(:,1,:) = rice_tem (jbdy)
	300	IF( TRIM(bf_alias(jp_bdys_i)%clrootname) == 'NOT USED' ) bf_alias(jp_bdys_i)%fnow(:,1,:) = rice_sal (jbdy)
	301	IF( TRIM(bf_alias(jp_bdyaip)%clrootname) == 'NOT USED' ) bf_alias(jp_bdyaip)%fnow(:,1,:) = rice_apnd(jbdy) * & ! rice_apnd is the pond fraction
	302	& bf_alias(jp_bdya_i)%fnow(:,1,:) ! ( a_ip = rice_apnd * a_i )
	303	IF( TRIM(bf_alias(jp_bdyhip)%clrootname) == 'NOT USED' ) bf_alias(jp_bdyhip)%fnow(:,1,:) = rice_hpnd(jbdy)
[15264]	304	IF( TRIM(bf_alias(jp_bdyhil)%clrootname) == 'NOT USED' ) bf_alias(jp_bdyhil)%fnow(:,1,:) = rice_hlid(jbdy)
[12395]	305
	306	! if T_i is read and not T_su, set T_su = T_i
	307	IF( TRIM(bf_alias(jp_bdyt_i)%clrootname) /= 'NOT USED' .AND. TRIM(bf_alias(jp_bdytsu)%clrootname) == 'NOT USED' ) &
	308	& bf_alias(jp_bdytsu)%fnow(:,1,:) = bf_alias(jp_bdyt_i)%fnow(:,1,:)
	309	! if T_s is read and not T_su, set T_su = T_s
	310	IF( TRIM(bf_alias(jp_bdyt_s)%clrootname) /= 'NOT USED' .AND. TRIM(bf_alias(jp_bdytsu)%clrootname) == 'NOT USED' ) &
	311	& bf_alias(jp_bdytsu)%fnow(:,1,:) = bf_alias(jp_bdyt_s)%fnow(:,1,:)
	312	! if T_i is read and not T_s, set T_s = T_i
	313	IF( TRIM(bf_alias(jp_bdyt_i)%clrootname) /= 'NOT USED' .AND. TRIM(bf_alias(jp_bdyt_s)%clrootname) == 'NOT USED' ) &
	314	& bf_alias(jp_bdyt_s)%fnow(:,1,:) = bf_alias(jp_bdyt_i)%fnow(:,1,:)
	315	! if T_su is read and not T_s, set T_s = T_su
	316	IF( TRIM(bf_alias(jp_bdytsu)%clrootname) /= 'NOT USED' .AND. TRIM(bf_alias(jp_bdyt_s)%clrootname) == 'NOT USED' ) &
	317	& bf_alias(jp_bdyt_s)%fnow(:,1,:) = bf_alias(jp_bdytsu)%fnow(:,1,:)
[11536]	318	! if T_su is read and not T_i, set T_i = (T_su + T_freeze)/2
	319	IF( TRIM(bf_alias(jp_bdytsu)%clrootname) /= 'NOT USED' .AND. TRIM(bf_alias(jp_bdyt_i)%clrootname) == 'NOT USED' ) &
	320	& bf_alias(jp_bdyt_i)%fnow(:,1,:) = 0.5_wp * ( bf_alias(jp_bdytsu)%fnow(:,1,:) + 271.15 )
	321	! if T_s is read and not T_i, set T_i = (T_s + T_freeze)/2
	322	IF( TRIM(bf_alias(jp_bdyt_s)%clrootname) /= 'NOT USED' .AND. TRIM(bf_alias(jp_bdyt_i)%clrootname) == 'NOT USED' ) &
	323	& bf_alias(jp_bdyt_i)%fnow(:,1,:) = 0.5_wp * ( bf_alias(jp_bdyt_s)%fnow(:,1,:) + 271.15 )
	324
	325	! make sure ponds = 0 if no ponds scheme
	326	IF ( .NOT.ln_pnd ) THEN
	327	bf_alias(jp_bdyaip)%fnow(:,1,:) = 0._wp
	328	bf_alias(jp_bdyhip)%fnow(:,1,:) = 0._wp
[15264]	329	bf_alias(jp_bdyhil)%fnow(:,1,:) = 0._wp
[11536]	330	ENDIF
[15264]	331	IF ( .NOT.ln_pnd_lids ) THEN
	332	bf_alias(jp_bdyhil)%fnow(:,1,:) = 0._wp
	333	ENDIF
[11536]	334
	335	! convert N-cat fields (input) into jpl-cat (output)
	336	ipl = SIZE(bf_alias(jp_bdya_i)%fnow, 3)
	337	IF( ipl /= jpl ) THEN ! ice: convert N-cat fields (input) into jpl-cat (output)
[15264]	338	CALL ice_var_itd( bf_alias(jp_bdyh_i)%fnow(:,1,:), bf_alias(jp_bdyh_s)%fnow(:,1,:), bf_alias(jp_bdya_i)%fnow(:,1,:), & ! in
	339	& dta_alias%h_i , dta_alias%h_s , dta_alias%a_i , & ! out
	340	& bf_alias(jp_bdyt_i)%fnow(:,1,:), bf_alias(jp_bdyt_s)%fnow(:,1,:), & ! in (optional)
	341	& bf_alias(jp_bdytsu)%fnow(:,1,:), bf_alias(jp_bdys_i)%fnow(:,1,:), & ! in -
	342	& bf_alias(jp_bdyaip)%fnow(:,1,:), bf_alias(jp_bdyhip)%fnow(:,1,:), bf_alias(jp_bdyhil)%fnow(:,1,:), & ! in -
	343	& dta_alias%t_i , dta_alias%t_s , & ! out -
	344	& dta_alias%tsu , dta_alias%s_i , & ! out -
	345	& dta_alias%aip , dta_alias%hip , dta_alias%hil ) ! out -
[11536]	346	ENDIF
	347	ENDIF
[4292]	348	#endif
[9019]	349	END DO ! jbdy
[911]	350
[7646]	351	IF ( ln_tide ) THEN
	352	IF (ln_dynspg_ts) THEN ! Fill temporary arrays with slow-varying bdy data
[11536]	353	DO jbdy = 1, nb_bdy ! Tidal component added in ts loop
	354	IF ( nn_dyn2d_dta(jbdy) .GE. 2 ) THEN
[12910]	355	IF( ASSOCIATED(dta_bdy(jbdy)%ssh) ) dta_bdy_s(jbdy)%ssh(:) = dta_bdy(jbdy)%ssh(:)
	356	IF( ASSOCIATED(dta_bdy(jbdy)%u2d) ) dta_bdy_s(jbdy)%u2d(:) = dta_bdy(jbdy)%u2d(:)
	357	IF( ASSOCIATED(dta_bdy(jbdy)%v2d) ) dta_bdy_s(jbdy)%v2d(:) = dta_bdy(jbdy)%v2d(:)
	358	ENDIF
	359	END DO
	360	ELSE ! Add tides if not split-explicit free surface else this is done in ts loop
	361	!
	362	CALL bdy_dta_tides( kt=kt, kt_offset=kt_offset )
[7646]	363	ENDIF
[12910]	364	ENDIF
	365	!
	366	IF( ln_timing ) CALL timing_stop('bdy_dta')
	367	!
	368	END SUBROUTINE bdy_dta
	369
[4292]	370
[6140]	371	SUBROUTINE bdy_dta_init
[3294]	372	!!----------------------------------------------------------------------
	373	!! * SUBROUTINE bdy_dta_init *
	374	!!
	375	!! ** Purpose : Initialise arrays for reading of external data
	376	!! for open boundary conditions
	377	!!
[4292]	378	!! ** Method :
[3294]	379	!!
	380	!!----------------------------------------------------------------------
[11536]	381	INTEGER :: jbdy, jfld ! Local integers
	382	INTEGER :: ierror, ios !
[6140]	383	!
[11536]	384	CHARACTER(len=3) :: cl3 !
[3294]	385	CHARACTER(len=100) :: cn_dir ! Root directory for location of data files
	386	LOGICAL :: ln_full_vel ! =T => full velocities in 3D boundary data
[9019]	387	! ! =F => baroclinic velocities in 3D boundary data
[11536]	388	LOGICAL :: ln_zinterp ! =T => requires a vertical interpolation of the bdydta
[15264]	389	REAL(wp) :: rn_ice_tem, rn_ice_sal, rn_ice_age, rn_ice_apnd, rn_ice_hpnd, rn_ice_hlid
[11536]	390	INTEGER :: ipk,ipl !
	391	INTEGER :: idvar ! variable ID
	392	INTEGER :: indims ! number of dimensions of the variable
	393	INTEGER :: iszdim ! number of dimensions of the variable
	394	INTEGER, DIMENSION(4) :: i4dimsz ! size of variable dimensions
	395	INTEGER :: igrd ! index for grid type (1,2,3 = T,U,V)
	396	LOGICAL :: lluld ! is the variable using the unlimited dimension
	397	LOGICAL :: llneed !
	398	LOGICAL :: llread !
[12910]	399	LOGICAL :: llfullbdy !
[11536]	400	TYPE(FLD_N), DIMENSION(1), TARGET :: bn_tem, bn_sal, bn_u3d, bn_v3d ! must be an array to be used with fld_fill
	401	TYPE(FLD_N), DIMENSION(1), TARGET :: bn_ssh, bn_u2d, bn_v2d ! informations about the fields to be read
[15264]	402	TYPE(FLD_N), DIMENSION(1), TARGET :: bn_a_i, bn_h_i, bn_h_s, bn_t_i, bn_t_s, bn_tsu, bn_s_i, bn_aip, bn_hip, bn_hil
[11536]	403	TYPE(FLD_N), DIMENSION(:), POINTER :: bn_alias ! must be an array to be used with fld_fill
	404	TYPE(FLD ), DIMENSION(:), POINTER :: bf_alias
	405	!
[15264]	406	NAMELIST/nambdy_dta/ cn_dir, bn_tem, bn_sal, bn_u3d, bn_v3d, bn_ssh, bn_u2d, bn_v2d, &
	407	& bn_a_i, bn_h_i, bn_h_s, bn_t_i, bn_t_s, bn_tsu, bn_s_i, bn_aip, bn_hip, bn_hil, &
	408	& rn_ice_tem, rn_ice_sal, rn_ice_age, rn_ice_apnd, rn_ice_hpnd, rn_ice_hlid, &
	409	& ln_full_vel, ln_zinterp
[3294]	410	!!---------------------------------------------------------------------------
[6140]	411	!
[3651]	412	IF(lwp) WRITE(numout,*)
	413	IF(lwp) WRITE(numout,*) 'bdy_dta_ini : initialization of data at the open boundaries'
	414	IF(lwp) WRITE(numout,*) '~~~~~~~~~~'
	415	IF(lwp) WRITE(numout,*) ''
	416
[11536]	417	ALLOCATE( bf(jpbdyfld,nb_bdy), STAT=ierror )
[3294]	418	IF( ierror > 0 ) THEN
	419	CALL ctl_stop( 'bdy_dta: unable to allocate bf structure' ) ; RETURN
[911]	420	ENDIF
[11536]	421	bf(:,:)%clrootname = 'NOT USED' ! default definition used as a flag in fld_read to do nothing.
	422	bf(:,:)%lzint = .FALSE. ! default definition
	423	bf(:,:)%ltotvel = .FALSE. ! default definition
	424
[3294]	425	! Read namelists
	426	! --------------
[4147]	427	REWIND(numnam_cfg)
[11536]	428	DO jbdy = 1, nb_bdy
[4147]	429
[11536]	430	WRITE(ctmp1, '(a,i2)') 'BDY number ', jbdy
	431	WRITE(ctmp2, '(a,i2)') 'block nambdy_dta number ', jbdy
[911]	432
[11536]	433	! There is only one nambdy_dta block in namelist_ref -> use it for each bdy so we do a rewind
	434	REWIND(numnam_ref)
	435	READ ( numnam_ref, nambdy_dta, IOSTAT = ios, ERR = 901)
	436	901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nambdy_dta in reference namelist' )
[911]	437
[11536]	438	! by-pass nambdy_dta reading if no input data used in this bdy
	439	IF( ( dta_bdy(jbdy)%lneed_dyn2d .AND. MOD(nn_dyn2d_dta(jbdy),2) == 1 ) &
	440	& .OR. ( dta_bdy(jbdy)%lneed_dyn3d .AND. nn_dyn3d_dta(jbdy) == 1 ) &
	441	& .OR. ( dta_bdy(jbdy)%lneed_tra .AND. nn_tra_dta(jbdy) == 1 ) &
	442	& .OR. ( dta_bdy(jbdy)%lneed_ice .AND. nn_ice_dta(jbdy) == 1 ) ) THEN
	443	! WARNING: we don't do a rewind here, each bdy reads its own nambdy_dta block one after another
	444	READ ( numnam_cfg, nambdy_dta, IOSTAT = ios, ERR = 902 )
	445	902 IF( ios > 0 ) CALL ctl_nam ( ios , 'nambdy_dta in configuration namelist' )
	446	IF(lwm) WRITE( numond, nambdy_dta )
	447	ENDIF
[911]	448
[11536]	449	! get the number of ice categories in bdy data file (use a_i information to do this)
	450	ipl = jpl ! default definition
	451	IF( dta_bdy(jbdy)%lneed_ice ) THEN ! if we need ice bdy data
	452	IF( nn_ice_dta(jbdy) == 1 ) THEN ! if we get ice bdy data from netcdf file
	453	CALL fld_fill( bf(jp_bdya_i,jbdy:jbdy), bn_a_i, cn_dir, 'bdy_dta', 'a_i'//' '//ctmp1, ctmp2 ) ! use namelist info
	454	CALL fld_clopn( bf(jp_bdya_i,jbdy), nyear, nmonth, nday ) ! not a problem when we call it again after
	455	idvar = iom_varid( bf(jp_bdya_i,jbdy)%num, bf(jp_bdya_i,jbdy)%clvar, kndims=indims, kdimsz=i4dimsz, lduld=lluld )
	456	IF( indims == 4 .OR. ( indims == 3 .AND. .NOT. lluld ) ) THEN ; ipl = i4dimsz(3) ! xylt or xyl
	457	ELSE ; ipl = 1 ! xy or xyt
[3294]	458	ENDIF
[11891]	459	bf(jp_bdya_i,jbdy)%clrootname = 'NOT USED' ! reset to default value as this subdomain may not need to read this bdy
[11536]	460	ENDIF
	461	ENDIF
[911]	462
[11536]	463	#if defined key_si3
	464	IF( .NOT.ln_pnd ) THEN
[15264]	465	rn_ice_apnd = 0. ; rn_ice_hpnd = 0. ; rn_ice_hlid = 0.
	466	CALL ctl_warn( 'rn_ice_apnd & rn_ice_hpnd = 0 & rn_ice_hlid = 0 when no ponds' )
[11536]	467	ENDIF
[15264]	468	IF( .NOT.ln_pnd_lids ) THEN
	469	rn_ice_hlid = 0.
	470	ENDIF
[11536]	471	#endif
[911]	472
[11536]	473	! temp, salt, age and ponds of incoming ice
	474	rice_tem (jbdy) = rn_ice_tem
	475	rice_sal (jbdy) = rn_ice_sal
	476	rice_age (jbdy) = rn_ice_age
	477	rice_apnd(jbdy) = rn_ice_apnd
	478	rice_hpnd(jbdy) = rn_ice_hpnd
[15264]	479	rice_hlid(jbdy) = rn_ice_hlid
	480
[11536]	481
	482	DO jfld = 1, jpbdyfld
[911]	483
[11536]	484	! =====================
	485	! ssh
	486	! =====================
	487	IF( jfld == jp_bdyssh ) THEN
	488	cl3 = 'ssh'
	489	igrd = 1 ! T point
	490	ipk = 1 ! surface data
	491	llneed = dta_bdy(jbdy)%lneed_ssh ! dta_bdy(jbdy)%ssh will be needed
	492	llread = MOD(nn_dyn2d_dta(jbdy),2) == 1 ! get data from NetCDF file
	493	bf_alias => bf(jp_bdyssh,jbdy:jbdy) ! alias for ssh structure of bdy number jbdy
	494	bn_alias => bn_ssh ! alias for ssh structure of nambdy_dta
	495	iszdim = idx_bdy(jbdy)%nblenrim(igrd) ! length of this bdy on this MPI processus : used only on the rim
[3294]	496	ENDIF
[11536]	497	! =====================
	498	! dyn2d
	499	! =====================
	500	IF( jfld == jp_bdyu2d ) THEN
	501	cl3 = 'u2d'
	502	igrd = 2 ! U point
	503	ipk = 1 ! surface data
[12910]	504	llneed = dta_bdy(jbdy)%lneed_dyn2d ! dta_bdy(jbdy)%u2d will be needed
[11536]	505	llread = .NOT. ln_full_vel .AND. MOD(nn_dyn2d_dta(jbdy),2) == 1 ! don't get u2d from u3d and read NetCDF file
	506	bf_alias => bf(jp_bdyu2d,jbdy:jbdy) ! alias for u2d structure of bdy number jbdy
	507	bn_alias => bn_u2d ! alias for u2d structure of nambdy_dta
[12910]	508	llfullbdy = ln_full_vel .OR. cn_dyn2d(jbdy) == 'frs' ! need u2d over the whole bdy or only over the rim?
	509	IF( llfullbdy ) THEN ; iszdim = idx_bdy(jbdy)%nblen(igrd)
	510	ELSE ; iszdim = idx_bdy(jbdy)%nblenrim(igrd)
[4292]	511	ENDIF
[3294]	512	ENDIF
[11536]	513	IF( jfld == jp_bdyv2d ) THEN
	514	cl3 = 'v2d'
	515	igrd = 3 ! V point
	516	ipk = 1 ! surface data
[12910]	517	llneed = dta_bdy(jbdy)%lneed_dyn2d ! dta_bdy(jbdy)%v2d will be needed
[11536]	518	llread = .NOT. ln_full_vel .AND. MOD(nn_dyn2d_dta(jbdy),2) == 1 ! don't get v2d from v3d and read NetCDF file
	519	bf_alias => bf(jp_bdyv2d,jbdy:jbdy) ! alias for v2d structure of bdy number jbdy
	520	bn_alias => bn_v2d ! alias for v2d structure of nambdy_dta
[12910]	521	llfullbdy = ln_full_vel .OR. cn_dyn2d(jbdy) == 'frs' ! need v2d over the whole bdy or only over the rim?
	522	IF( llfullbdy ) THEN ; iszdim = idx_bdy(jbdy)%nblen(igrd)
	523	ELSE ; iszdim = idx_bdy(jbdy)%nblenrim(igrd)
[4292]	524	ENDIF
[911]	525	ENDIF
[11536]	526	! =====================
	527	! dyn3d
	528	! =====================
	529	IF( jfld == jp_bdyu3d ) THEN
	530	cl3 = 'u3d'
	531	igrd = 2 ! U point
	532	ipk = jpk ! 3d data
	533	llneed = dta_bdy(jbdy)%lneed_dyn3d .OR. & ! dta_bdy(jbdy)%u3d will be needed
	534	& ( dta_bdy(jbdy)%lneed_dyn2d .AND. ln_full_vel ) ! u3d needed to compute u2d
	535	llread = nn_dyn3d_dta(jbdy) == 1 ! get data from NetCDF file
	536	bf_alias => bf(jp_bdyu3d,jbdy:jbdy) ! alias for u3d structure of bdy number jbdy
	537	bn_alias => bn_u3d ! alias for u3d structure of nambdy_dta
	538	iszdim = idx_bdy(jbdy)%nblen(igrd) ! length of this bdy on this MPI processus
	539	ENDIF
	540	IF( jfld == jp_bdyv3d ) THEN
	541	cl3 = 'v3d'
	542	igrd = 3 ! V point
	543	ipk = jpk ! 3d data
	544	llneed = dta_bdy(jbdy)%lneed_dyn3d .OR. & ! dta_bdy(jbdy)%v3d will be needed
	545	& ( dta_bdy(jbdy)%lneed_dyn2d .AND. ln_full_vel ) ! v3d needed to compute v2d
	546	llread = nn_dyn3d_dta(jbdy) == 1 ! get data from NetCDF file
	547	bf_alias => bf(jp_bdyv3d,jbdy:jbdy) ! alias for v3d structure of bdy number jbdy
	548	bn_alias => bn_v3d ! alias for v3d structure of nambdy_dta
	549	iszdim = idx_bdy(jbdy)%nblen(igrd) ! length of this bdy on this MPI processus
	550	ENDIF
[911]	551
[11536]	552	! =====================
	553	! tra
	554	! =====================
	555	IF( jfld == jp_bdytem ) THEN
	556	cl3 = 'tem'
	557	igrd = 1 ! T point
	558	ipk = jpk ! 3d data
	559	llneed = dta_bdy(jbdy)%lneed_tra ! dta_bdy(jbdy)%tem will be needed
	560	llread = nn_tra_dta(jbdy) == 1 ! get data from NetCDF file
	561	bf_alias => bf(jp_bdytem,jbdy:jbdy) ! alias for ssh structure of bdy number jbdy
	562	bn_alias => bn_tem ! alias for ssh structure of nambdy_dta
	563	iszdim = idx_bdy(jbdy)%nblen(igrd) ! length of this bdy on this MPI processus
[4333]	564	ENDIF
[11536]	565	IF( jfld == jp_bdysal ) THEN
	566	cl3 = 'sal'
	567	igrd = 1 ! T point
	568	ipk = jpk ! 3d data
	569	llneed = dta_bdy(jbdy)%lneed_tra ! dta_bdy(jbdy)%sal will be needed
	570	llread = nn_tra_dta(jbdy) == 1 ! get data from NetCDF file
	571	bf_alias => bf(jp_bdysal,jbdy:jbdy) ! alias for ssh structure of bdy number jbdy
	572	bn_alias => bn_sal ! alias for ssh structure of nambdy_dta
	573	iszdim = idx_bdy(jbdy)%nblen(igrd) ! length of this bdy on this MPI processus
[3294]	574	ENDIF
[911]	575
[11536]	576	! =====================
	577	! ice
	578	! =====================
	579	IF( jfld == jp_bdya_i .OR. jfld == jp_bdyh_i .OR. jfld == jp_bdyh_s .OR. &
	580	& jfld == jp_bdyt_i .OR. jfld == jp_bdyt_s .OR. jfld == jp_bdytsu .OR. &
[15264]	581	& jfld == jp_bdys_i .OR. jfld == jp_bdyaip .OR. jfld == jp_bdyhip .OR. jfld == jp_bdyhil ) THEN
[11536]	582	igrd = 1 ! T point
	583	ipk = ipl ! jpl-cat data
	584	llneed = dta_bdy(jbdy)%lneed_ice ! ice will be needed
	585	llread = nn_ice_dta(jbdy) == 1 ! get data from NetCDF file
	586	iszdim = idx_bdy(jbdy)%nblen(igrd) ! length of this bdy on this MPI processus
[4292]	587	ENDIF
[11536]	588	IF( jfld == jp_bdya_i ) THEN
	589	cl3 = 'a_i'
	590	bf_alias => bf(jp_bdya_i,jbdy:jbdy) ! alias for a_i structure of bdy number jbdy
	591	bn_alias => bn_a_i ! alias for a_i structure of nambdy_dta
[4292]	592	ENDIF
[11536]	593	IF( jfld == jp_bdyh_i ) THEN
	594	cl3 = 'h_i'
	595	bf_alias => bf(jp_bdyh_i,jbdy:jbdy) ! alias for h_i structure of bdy number jbdy
	596	bn_alias => bn_h_i ! alias for h_i structure of nambdy_dta
[3294]	597	ENDIF
[11536]	598	IF( jfld == jp_bdyh_s ) THEN
	599	cl3 = 'h_s'
	600	bf_alias => bf(jp_bdyh_s,jbdy:jbdy) ! alias for h_s structure of bdy number jbdy
	601	bn_alias => bn_h_s ! alias for h_s structure of nambdy_dta
[4292]	602	ENDIF
[11536]	603	IF( jfld == jp_bdyt_i ) THEN
	604	cl3 = 't_i'
	605	bf_alias => bf(jp_bdyt_i,jbdy:jbdy) ! alias for t_i structure of bdy number jbdy
	606	bn_alias => bn_t_i ! alias for t_i structure of nambdy_dta
[4292]	607	ENDIF
[11536]	608	IF( jfld == jp_bdyt_s ) THEN
	609	cl3 = 't_s'
	610	bf_alias => bf(jp_bdyt_s,jbdy:jbdy) ! alias for t_s structure of bdy number jbdy
	611	bn_alias => bn_t_s ! alias for t_s structure of nambdy_dta
[4292]	612	ENDIF
[11536]	613	IF( jfld == jp_bdytsu ) THEN
	614	cl3 = 'tsu'
	615	bf_alias => bf(jp_bdytsu,jbdy:jbdy) ! alias for tsu structure of bdy number jbdy
	616	bn_alias => bn_tsu ! alias for tsu structure of nambdy_dta
[3294]	617	ENDIF
[11536]	618	IF( jfld == jp_bdys_i ) THEN
	619	cl3 = 's_i'
	620	bf_alias => bf(jp_bdys_i,jbdy:jbdy) ! alias for s_i structure of bdy number jbdy
	621	bn_alias => bn_s_i ! alias for s_i structure of nambdy_dta
	622	ENDIF
	623	IF( jfld == jp_bdyaip ) THEN
	624	cl3 = 'aip'
	625	bf_alias => bf(jp_bdyaip,jbdy:jbdy) ! alias for aip structure of bdy number jbdy
	626	bn_alias => bn_aip ! alias for aip structure of nambdy_dta
	627	ENDIF
	628	IF( jfld == jp_bdyhip ) THEN
	629	cl3 = 'hip'
	630	bf_alias => bf(jp_bdyhip,jbdy:jbdy) ! alias for hip structure of bdy number jbdy
	631	bn_alias => bn_hip ! alias for hip structure of nambdy_dta
	632	ENDIF
[15264]	633	IF( jfld == jp_bdyhil ) THEN
	634	cl3 = 'hil'
	635	bf_alias => bf(jp_bdyhil,jbdy:jbdy) ! alias for hil structure of bdy number jbdy
	636	bn_alias => bn_hil ! alias for hil structure of nambdy_dta
	637	ENDIF
[911]	638
[11879]	639	IF( llneed .AND. iszdim > 0 ) THEN ! dta_bdy(jbdy)%xxx will be needed
[11536]	640	! ! -> must be associated with an allocated target
	641	ALLOCATE( bf_alias(1)%fnow( iszdim, 1, ipk ) ) ! allocate the target
	642	!
	643	IF( llread ) THEN ! get data from NetCDF file
	644	CALL fld_fill( bf_alias, bn_alias, cn_dir, 'bdy_dta', cl3//' '//ctmp1, ctmp2 ) ! use namelist info
	645	IF( bf_alias(1)%ln_tint ) ALLOCATE( bf_alias(1)%fdta( iszdim, 1, ipk, 2 ) )
	646	bf_alias(1)%imap => idx_bdy(jbdy)%nbmap(1:iszdim,igrd) ! associate the mapping used for this bdy
	647	bf_alias(1)%igrd = igrd ! used only for vertical integration of 3D arrays
[11793]	648	bf_alias(1)%ibdy = jbdy ! " " " " " " " "
[11536]	649	bf_alias(1)%ltotvel = ln_full_vel ! T if u3d is full velocity
	650	bf_alias(1)%lzint = ln_zinterp ! T if it requires a vertical interpolation
[4292]	651	ENDIF
	652
[11536]	653	! associate the pointer and get rid of the dimensions with a size equal to 1
	654	IF( jfld == jp_bdyssh ) dta_bdy(jbdy)%ssh => bf_alias(1)%fnow(:,1,1)
	655	IF( jfld == jp_bdyu2d ) dta_bdy(jbdy)%u2d => bf_alias(1)%fnow(:,1,1)
	656	IF( jfld == jp_bdyv2d ) dta_bdy(jbdy)%v2d => bf_alias(1)%fnow(:,1,1)
	657	IF( jfld == jp_bdyu3d ) dta_bdy(jbdy)%u3d => bf_alias(1)%fnow(:,1,:)
	658	IF( jfld == jp_bdyv3d ) dta_bdy(jbdy)%v3d => bf_alias(1)%fnow(:,1,:)
	659	IF( jfld == jp_bdytem ) dta_bdy(jbdy)%tem => bf_alias(1)%fnow(:,1,:)
	660	IF( jfld == jp_bdysal ) dta_bdy(jbdy)%sal => bf_alias(1)%fnow(:,1,:)
	661	IF( jfld == jp_bdya_i ) THEN
	662	IF( ipk == jpl ) THEN ; dta_bdy(jbdy)%a_i => bf_alias(1)%fnow(:,1,:)
	663	ELSE ; ALLOCATE( dta_bdy(jbdy)%a_i(iszdim,jpl) )
	664	ENDIF
	665	ENDIF
	666	IF( jfld == jp_bdyh_i ) THEN
	667	IF( ipk == jpl ) THEN ; dta_bdy(jbdy)%h_i => bf_alias(1)%fnow(:,1,:)
	668	ELSE ; ALLOCATE( dta_bdy(jbdy)%h_i(iszdim,jpl) )
	669	ENDIF
	670	ENDIF
	671	IF( jfld == jp_bdyh_s ) THEN
	672	IF( ipk == jpl ) THEN ; dta_bdy(jbdy)%h_s => bf_alias(1)%fnow(:,1,:)
	673	ELSE ; ALLOCATE( dta_bdy(jbdy)%h_s(iszdim,jpl) )
	674	ENDIF
	675	ENDIF
	676	IF( jfld == jp_bdyt_i ) THEN
	677	IF( ipk == jpl ) THEN ; dta_bdy(jbdy)%t_i => bf_alias(1)%fnow(:,1,:)
	678	ELSE ; ALLOCATE( dta_bdy(jbdy)%t_i(iszdim,jpl) )
	679	ENDIF
	680	ENDIF
	681	IF( jfld == jp_bdyt_s ) THEN
	682	IF( ipk == jpl ) THEN ; dta_bdy(jbdy)%t_s => bf_alias(1)%fnow(:,1,:)
	683	ELSE ; ALLOCATE( dta_bdy(jbdy)%t_s(iszdim,jpl) )
	684	ENDIF
	685	ENDIF
	686	IF( jfld == jp_bdytsu ) THEN
	687	IF( ipk == jpl ) THEN ; dta_bdy(jbdy)%tsu => bf_alias(1)%fnow(:,1,:)
	688	ELSE ; ALLOCATE( dta_bdy(jbdy)%tsu(iszdim,jpl) )
	689	ENDIF
	690	ENDIF
	691	IF( jfld == jp_bdys_i ) THEN
	692	IF( ipk == jpl ) THEN ; dta_bdy(jbdy)%s_i => bf_alias(1)%fnow(:,1,:)
	693	ELSE ; ALLOCATE( dta_bdy(jbdy)%s_i(iszdim,jpl) )
	694	ENDIF
	695	ENDIF
	696	IF( jfld == jp_bdyaip ) THEN
	697	IF( ipk == jpl ) THEN ; dta_bdy(jbdy)%aip => bf_alias(1)%fnow(:,1,:)
	698	ELSE ; ALLOCATE( dta_bdy(jbdy)%aip(iszdim,jpl) )
	699	ENDIF
	700	ENDIF
	701	IF( jfld == jp_bdyhip ) THEN
	702	IF( ipk == jpl ) THEN ; dta_bdy(jbdy)%hip => bf_alias(1)%fnow(:,1,:)
	703	ELSE ; ALLOCATE( dta_bdy(jbdy)%hip(iszdim,jpl) )
	704	ENDIF
	705	ENDIF
[15264]	706	IF( jfld == jp_bdyhil ) THEN
	707	IF( ipk == jpl ) THEN ; dta_bdy(jbdy)%hil => bf_alias(1)%fnow(:,1,:)
	708	ELSE ; ALLOCATE( dta_bdy(jbdy)%hil(iszdim,jpl) )
	709	ENDIF
	710	ENDIF
[4292]	711	ENDIF
[11536]	712
	713	END DO ! jpbdyfld
[6140]	714	!
[9019]	715	END DO ! jbdy
[6140]	716	!
	717	END SUBROUTINE bdy_dta_init
[11536]	718
[911]	719	!!==============================================================================
	720	END MODULE bdydta

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

Download in other formats: