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/2019/dev_r12072_MERGE_OPTION2_2019/src/OCE/BDY – NEMO

Context Navigation

source: NEMO/branches/2019/dev_r12072_MERGE_OPTION2_2019/src/OCE/BDY/bdydta.F90 @ 12256

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

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

Download in other formats: