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 trunk/NEMO/OPA_SRC/BDY – NEMO

Context Navigation

source: trunk/NEMO/OPA_SRC/BDY/bdydta.F90 @ 1473

Last change on this file since 1473 was 1241, checked in by rblod, 15 years ago
Fix a stupid bug for time splitting and ensure restartability for dynspg_ts in addition, see tickets #280 and #292
Property svn:executable set to ``* Property svn:keywords set to `Id`
File size: 45.3 KB

Rev	Line
[911]	1	MODULE bdydta
[1125]	2	!!======================================================================
	3	!! * MODULE bdydta *
[911]	4	!! Open boundary data : read the data for the unstructured open boundaries.
[1125]	5	!!======================================================================
	6	!! History : 1.0 ! 2005-01 (J. Chanut, A. Sellar) Original code
	7	!! - ! 2007-01 (D. Storkey) Update to use IOM module
	8	!! - ! 2007-07 (D. Storkey) add bdy_dta_bt
	9	!! 3.0 ! 2008-04 (NEMO team) add in the reference version
	10	!!----------------------------------------------------------------------
	11	#if defined key_bdy
	12	!!----------------------------------------------------------------------
	13	!! 'key_bdy' Unstructured Open Boundary Conditions
	14	!!----------------------------------------------------------------------
	15	!! bdy_dta : read u, v, t, s data along open boundaries
	16	!! bdy_dta_bt : read depth-mean velocities and elevation along open
	17	!! boundaries
	18	!!----------------------------------------------------------------------
[911]	19	USE oce ! ocean dynamics and tracers
	20	USE dom_oce ! ocean space and time domain
	21	USE daymod ! calendar
	22	USE phycst ! physical constants
	23	USE bdy_oce ! ocean open boundary conditions
	24	USE bdytides ! tidal forcing at boundaries
	25	USE iom
	26	USE ioipsl
	27	USE in_out_manager ! I/O logical units
	28
	29	IMPLICIT NONE
	30	PRIVATE
	31
[1125]	32	PUBLIC bdy_dta ! routines called by step.F90
	33	PUBLIC bdy_dta_bt
[911]	34
[1125]	35	INTEGER :: numbdyt, numbdyu, numbdyv !: logical units for T-, U-, & V-points data file, resp.
	36	INTEGER :: ntimes_bdy !: exact number of time dumps in data files
	37	INTEGER :: nbdy_b, nbdy_a !: record of bdy data file for before and after model time step
	38	INTEGER :: numbdyt_bt, numbdyu_bt, numbdyv_bt !: logical unit for T-, U- & V-points data file, resp.
	39	INTEGER :: ntimes_bdy_bt !: exact number of time dumps in data files
	40	INTEGER :: nbdy_b_bt, nbdy_a_bt !: record of bdy data file for before and after model time step
[911]	41
[1125]	42	INTEGER, DIMENSION (jpbtime) :: istep, istep_bt !: time array in seconds in each data file
[911]	43
[1125]	44	REAL(wp) :: zoffset !: time offset between time origin in file & start time of model run
[911]	45
[1125]	46	REAL(wp), DIMENSION(jpbdim,jpk,2) :: tbdydta, sbdydta !: time interpolated values of T and S bdy data
	47	REAL(wp), DIMENSION(jpbdim,jpk,2) :: ubdydta, vbdydta !: time interpolated values of U and V bdy data
	48	REAL(wp), DIMENSION(jpbdim,2) :: ubtbdydta, vbtbdydta !: Arrays used for time interpolation of bdy data
	49	REAL(wp), DIMENSION(jpbdim,2) :: sshbdydta !: bdy data of ssh
[911]	50
[1125]	51	!!----------------------------------------------------------------------
	52	!! NEMO/OPA 3.0 , LOCEAN-IPSL (2008)
[1146]	53	!! $Id$
[1125]	54	!! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
	55	!!----------------------------------------------------------------------
[911]	56
	57	CONTAINS
	58
	59	SUBROUTINE bdy_dta( kt )
[1125]	60	!!----------------------------------------------------------------------
	61	!! * SUBROUTINE bdy_dta *
[911]	62	!!
[1125]	63	!! ** Purpose : Read unstructured boundary data for FRS condition.
[911]	64	!!
[1125]	65	!! ** Method : At the first timestep, read in boundary data for two
	66	!! times from the file and time-interpolate. At other
	67	!! timesteps, check to see if we need another time from
	68	!! the file. If so read it in. Time interpolate.
	69	!!----------------------------------------------------------------------
	70	INTEGER, INTENT( in ) :: kt ! ocean time-step index (for timesplitting option, otherwise zero)
[911]	71	!!
[1125]	72	CHARACTER(LEN=80), DIMENSION(3) :: clfile ! names of input files
	73	CHARACTER(LEN=70 ) :: clunits ! units attribute of time coordinate
	74	LOGICAL :: lect ! flag for reading
	75	INTEGER :: it, ib, ik, igrd ! dummy loop indices
	76	INTEGER :: igrd_start, igrd_end ! start and end of loops on igrd
	77	INTEGER :: idvar ! netcdf var ID
	78	INTEGER :: iman, i15, imois ! Time variables for monthly clim forcing
	79	INTEGER :: ntimes_bdyt, ntimes_bdyu, ntimes_bdyv
[911]	80	INTEGER :: itimer, totime
[1125]	81	INTEGER :: ii, ij ! array addresses
	82	INTEGER :: ipi, ipj, ipk, inum ! temporary integers (NetCDF read)
[911]	83	INTEGER :: iyear0, imonth0, iday0
	84	INTEGER :: ihours0, iminutes0, isec0
[1125]	85	INTEGER :: iyear, imonth, iday, isecs
[911]	86	INTEGER, DIMENSION(jpbtime) :: istept, istepu, istepv ! time arrays from data files
[1125]	87	REAL(wp) :: dayfrac, zxy, zoffsett
	88	REAL(wp) :: zoffsetu, zoffsetv
[911]	89	REAL(wp) :: dayjul0, zdayjulini
[1125]	90	REAL(wp), DIMENSION(jpbtime) :: zstepr ! REAL time array from data files
	91	REAL(wp), DIMENSION(jpbdta,1,jpk) :: zdta ! temporary array for data fields
[911]	92	!!---------------------------------------------------------------------------
	93
[1125]	94	IF( ln_bdy_dyn_frs .OR. ln_bdy_tra_frs ) THEN ! If these are both false then this routine
	95	! does nothing.
	96
[911]	97	! -------------------- !
	98	! Initialization !
	99	! -------------------- !
	100
	101	lect = .false. ! If true, read a time record
	102
	103	! Some time variables for monthly climatological forcing:
	104	! *******************************************************
[1125]	105	!!gm here use directely daymod variables
[911]	106
[1125]	107	iman = INT( raamo ) ! Number of months in a year
[911]	108
	109	i15 = INT( 2*FLOAT( nday ) / ( FLOAT( nobis(nmonth) ) + 0.5 ) )
	110	! i15=0 if the current day is in the first half of the month, else i15=1
	111
	112	imois = nmonth + i15 - 1 ! imois is the first month record
[1125]	113	IF( imois == 0 ) imois = iman
[911]	114
	115	! Time variable for non-climatological forcing:
	116	! *********************************************
	117	itimer = (kt-nit000+1)*rdt ! current time in seconds for interpolation
	118
	119
[1125]	120	! !-------------------!
	121	IF( kt == nit000 ) THEN ! First call only !
	122	! !-------------------!
	123	istep(:) = 0
	124	nbdy_b = 0
	125	nbdy_a = 0
[911]	126
[1125]	127	! Get time information from bdy data file
	128	! ***************************************
[911]	129
[1125]	130	IF(lwp) WRITE(numout,*)
	131	IF(lwp) WRITE(numout,*) 'bdy_dta : Initialize unstructured boundary data'
	132	IF(lwp) WRITE(numout,*) '~~~~~~~'
[911]	133
[1125]	134	IF ( nbdy_dta == 0 ) THEN
	135	!
	136	IF(lwp) WRITE(numout,*) ' Bdy data are taken from initial conditions'
	137	!
	138	ELSEIF (nbdy_dta == 1) THEN
	139	!
	140	IF(lwp) WRITE(numout,*) ' Bdy data are read in netcdf files'
	141	!
	142	dayfrac = adatrj - FLOAT( itimer ) / 86400. ! day fraction at time step kt-1
	143	dayfrac = dayfrac - INT ( dayfrac ) !
	144	totime = ( nitend - nit000 + 1 ) * rdt ! Total time of the run to verify that all the
	145	! ! necessary time dumps in file are included
	146	!
	147	clfile(1) = filbdy_data_T
	148	clfile(2) = filbdy_data_U
	149	clfile(3) = filbdy_data_V
	150	!
	151	! how many files are we to read in?
	152	igrd_start = 1
	153	igrd_end = 3
	154	IF(.NOT. ln_bdy_tra_frs .AND. .NOT. ln_bdy_ice_frs) THEN
	155	! No T-grid file.
	156	igrd_start = 2
	157	ELSEIF ( .NOT. ln_bdy_dyn_frs ) THEN
	158	! No U-grid or V-grid file.
	159	igrd_end = 1
	160	ENDIF
[911]	161
[1125]	162	DO igrd = igrd_start, igrd_end ! loop over T, U & V grid !
	163	! !---------------------------!
	164	CALL iom_open( clfile(igrd), inum )
	165	CALL iom_gettime( inum, zstepr, kntime=ntimes_bdy, cdunits=clunits )
[911]	166
[1125]	167	SELECT CASE( igrd )
	168	CASE (1)
	169	numbdyt = inum
	170	CASE (2)
	171	numbdyu = inum
	172	CASE (3)
	173	numbdyv = inum
	174	END SELECT
[911]	175
[1125]	176	! Calculate time offset
	177	READ(clunits,7000) iyear0, imonth0, iday0, ihours0, iminutes0, isec0
	178	! Convert time origin in file to julian days
	179	isec0 = isec0 + ihours060.60. + iminutes0*60.
	180	CALL ymds2ju(iyear0, imonth0, iday0, real(isec0), dayjul0)
	181	! Compute model initialization time
	182	iyear = ndastp / 10000
	183	imonth = ( ndastp - iyear * 10000 ) / 100
	184	iday = ndastp - iyear * 10000 - imonth * 100
	185	isecs = dayfrac * 86400
	186	CALL ymds2ju(iyear, imonth, iday, real(isecs) , zdayjulini)
	187	! offset from initialization date:
	188	zoffset = (dayjul0-zdayjulini)*86400
	189	!
	190	7000 FORMAT('seconds since ', I4.4,'-',I2.2,'-',I2.2,' ',I2.2,':',I2.2,':',I2.2)
[911]	191
[1125]	192	!! TO BE DONE... Check consistency between calendar from file
	193	!! (available optionally from iom_gettime) and calendar in model
	194	!! when calendar in model available outside of IOIPSL.
[911]	195
[1125]	196	IF(lwp) WRITE(numout,*) 'number of times: ',ntimes_bdy
	197	IF(lwp) WRITE(numout,*) 'offset: ',zoffset
	198	IF(lwp) WRITE(numout,*) 'totime: ',totime
	199	IF(lwp) WRITE(numout,*) 'zstepr: ',zstepr
[911]	200
[1125]	201	! Check that there are not too many times in the file.
	202	IF( ntimes_bdy > jpbtime ) THEN
	203	WRITE(ctmp1,*) 'Check file: ', clfile(igrd), 'jpbtime= ', jpbtime, ' ntimes_bdy= ', ntimes_bdy
	204	CALL ctl_stop( 'Number of time dumps in files exceed jpbtime parameter', ctmp1 )
	205	ENDIF
[911]	206
[1125]	207	! Check that time array increases:
	208
	209	it = 1
	210	DO WHILE( zstepr(it+1) > zstepr(it) .AND. it /= ntimes_bdy - 1 )
	211	it = it + 1
	212	END DO
[911]	213
[1125]	214	IF( it.NE.ntimes_bdy-1 .AND. ntimes_bdy > 1 ) THEN
	215	WRITE(ctmp1,*) 'Check file: ', clfile(igrd)
	216	CALL ctl_stop( 'Time array in unstructured boundary data files', &
	217	& 'does not continuously increase.' , ctmp1 )
	218	ENDIF
	219	!
	220	! Check that times in file span model run time:
	221	IF( zstepr(1) + zoffset > 0 ) THEN
	222	WRITE(ctmp1,*) 'Check file: ', clfile(igrd)
	223	CALL ctl_stop( 'First time dump in bdy file is after model initial time', ctmp1 )
	224	END IF
	225	IF( zstepr(ntimes_bdy) + zoffset < totime ) THEN
	226	WRITE(ctmp1,*) 'Check file: ', clfile(igrd)
	227	CALL ctl_stop( 'Last time dump in bdy file is before model final time', ctmp1 )
	228	END IF
	229	!
	230	IF ( igrd == 1 ) THEN
	231	ntimes_bdyt = ntimes_bdy
	232	zoffsett = zoffset
	233	istept(:) = INT( zstepr(:) + zoffset )
	234	ELSEIF(igrd == 2 ) THEN
	235	ntimes_bdyu = ntimes_bdy
	236	zoffsetu = zoffset
	237	istepu(:) = INT( zstepr(:) + zoffset )
	238	ELSEIF(igrd == 3 ) THEN
	239	ntimes_bdyv = ntimes_bdy
	240	zoffsetv = zoffset
	241	istepv(:) = INT( zstepr(:) + zoffset )
	242	ENDIF
	243	!
	244	END DO ! end loop over T, U & V grid
[911]	245
[1125]	246	IF (igrd_start == 1 .and. igrd_end == 3) THEN
	247	! Only test differences if we are reading in 3 files
	248	! Verify time consistency between files
	249	IF( ntimes_bdyu /= ntimes_bdyt .OR. ntimes_bdyv /= ntimes_bdyt ) THEN
	250	CALL ctl_stop( 'Bdy data files must have the same number of time dumps', &
	251	& 'Multiple time frequencies not implemented yet' )
	252	ENDIF
	253	ntimes_bdy = ntimes_bdyt
	254	!
	255	IF( zoffsetu /= zoffsett .OR. zoffsetv /= zoffsett ) THEN
	256	CALL ctl_stop( 'Bdy data files must have the same time origin', &
	257	& 'Multiple time frequencies not implemented yet' )
	258	ENDIF
	259	zoffset = zoffsett
	260	ENDIF
[911]	261
[1125]	262	IF( igrd_start == 1 ) THEN
	263	istep(:) = istept(:)
	264	ELSE
	265	istep(:) = istepu(:)
	266	ENDIF
[911]	267
[1125]	268	! Check number of time dumps:
	269	IF( ntimes_bdy == 1 .AND. .NOT. ln_bdy_clim ) THEN
	270	CALL ctl_stop( 'There is only one time dump in data files', &
	271	& 'Choose ln_bdy_clim=.true. in namelist for constant bdy forcing.' )
	272	ENDIF
[911]	273
[1125]	274	IF( ln_bdy_clim ) THEN
	275	IF( ntimes_bdy /= 1 .AND. ntimes_bdy /= 12 ) THEN
	276	CALL ctl_stop( 'For climatological boundary forcing (ln_bdy_clim=.true.),', &
	277	& 'bdy data files must contain 1 or 12 time dumps.' )
	278	ELSEIF( ntimes_bdy == 1 ) THEN
	279	IF(lwp) WRITE(numout,*)
	280	IF(lwp) WRITE(numout,*) 'We assume constant boundary forcing from bdy data files'
	281	ELSEIF( ntimes_bdy == 12 ) THEN
	282	IF(lwp) WRITE(numout,*)
	283	IF(lwp) WRITE(numout,*) 'We assume monthly (and cyclic) boundary forcing from bdy data files'
	284	ENDIF
[911]	285	ENDIF
	286
[1125]	287	! Find index of first record to read (before first model time).
	288	it = 1
	289	DO WHILE( istep(it+1) <= 0 .AND. it <= ntimes_bdy - 1 )
	290	it = it + 1
[911]	291	END DO
[1125]	292	nbdy_b = it
	293	!
	294	WRITE(numout,*) 'Time offset is ',zoffset
	295	WRITE(numout,*) 'First record to read is ',nbdy_b
[911]	296
[1125]	297	ENDIF ! endif (nbdy_dta == 1)
[911]	298
	299
[1125]	300	! 1.2 Read first record in file if necessary (ie if nbdy_dta == 1)
	301	! *****************************************************************
[911]	302
[1125]	303	IF( nbdy_dta == 0) THEN ! boundary data arrays are filled with initial conditions
	304	!
	305	IF (ln_bdy_tra_frs) THEN
	306	igrd = 1 ! T-points data
	307	DO ib = 1, nblen(igrd)
	308	ii = nbi(ib,igrd)
	309	ij = nbj(ib,igrd)
	310	DO ik = 1, jpkm1
	311	tbdy(ib,ik) = tn(ii, ij, ik)
	312	sbdy(ib,ik) = sn(ii, ij, ik)
	313	ENDDO
	314	END DO
[911]	315	ENDIF
	316
[1125]	317	IF(ln_bdy_dyn_frs) THEN
	318	igrd = 2 ! U-points data
	319	DO ib = 1, nblen(igrd)
	320	ii = nbi(ib,igrd)
	321	ij = nbj(ib,igrd)
	322	DO ik = 1, jpkm1
	323	ubdy(ib,ik) = un(ii, ij, ik)
	324	ENDDO
	325	END DO
[911]	326
[1125]	327	igrd = 3 ! V-points data
	328	DO ib = 1, nblen(igrd)
	329	ii = nbi(ib,igrd)
	330	ij = nbj(ib,igrd)
	331	DO ik = 1, jpkm1
	332	vbdy(ib,ik) = vn(ii, ij, ik)
	333	ENDDO
	334	END DO
[911]	335	ENDIF
[1125]	336	!
	337	ELSEIF( nbdy_dta == 1 ) THEN ! Set first record in the climatological case:
	338	!
	339	IF( ln_bdy_clim .AND. ntimes_bdy == 1 ) THEN
	340	nbdy_a = 1
	341	ELSEIF( ln_bdy_clim .AND. ntimes_bdy == iman ) THEN
	342	nbdy_b = 0
	343	nbdy_a = imois
	344	ELSE
	345	nbdy_a = nbdy_b
	346	ENDIF
	347
	348	! Read first record:
	349	ipj = 1
	350	ipk = jpk
	351	igrd = 1
	352	ipi = nblendta(igrd)
[911]	353
[1125]	354	IF(ln_bdy_tra_frs) THEN
	355	igrd = 1 ! Temperature
	356	IF( nblendta(igrd) <= 0 ) THEN
	357	idvar = iom_varid( numbdyt, 'votemper' )
	358	nblendta(igrd) = iom_file(numbdyt)%dimsz(1,idvar)
	359	ENDIF
	360	WRITE(numout,*) 'Dim size for votemper is ', nblendta(igrd)
	361	ipi = nblendta(igrd)
	362	CALL iom_get ( numbdyt, jpdom_unknown, 'votemper', zdta(1:ipi,1:ipj,1:ipk), nbdy_a )
[911]	363
[1125]	364	DO ib = 1, nblen(igrd)
	365	DO ik = 1, jpkm1
	366	tbdydta(ib,ik,2) = zdta(nbmap(ib,igrd),1,ik)
	367	END DO
	368	END DO
	369	!
	370	igrd = 1 ! salinity
	371	IF( nblendta(igrd) .le. 0 ) THEN
	372	idvar = iom_varid( numbdyt, 'vosaline' )
	373	nblendta(igrd) = iom_file(numbdyt)%dimsz(1,idvar)
	374	ENDIF
	375	WRITE(numout,*) 'Dim size for vosaline is ', nblendta(igrd)
	376	ipi = nblendta(igrd)
	377	CALL iom_get ( numbdyt, jpdom_unknown, 'vosaline', zdta(1:ipi,1:ipj,1:ipk), nbdy_a )
[911]	378
[1125]	379	DO ib = 1, nblen(igrd)
	380	DO ik = 1, jpkm1
	381	sbdydta(ib,ik,2) = zdta(nbmap(ib,igrd),1,ik)
	382	END DO
	383	END DO
	384	ENDIF ! ln_bdy_tra_frs
	385
	386	IF(ln_bdy_dyn_frs) THEN
[911]	387
[1125]	388	igrd = 2 ! u-velocity
	389	IF ( nblendta(igrd) .le. 0 ) THEN
	390	idvar = iom_varid( numbdyu,'vozocrtx' )
	391	nblendta(igrd) = iom_file(numbdyu)%dimsz(1,idvar)
	392	ENDIF
	393	WRITE(numout,*) 'Dim size for vozocrtx is ', nblendta(igrd)
	394	ipi = nblendta(igrd)
	395	CALL iom_get ( numbdyu, jpdom_unknown,'vozocrtx',zdta(1:ipi,1:ipj,1:ipk),nbdy_a )
	396	DO ib = 1, nblen(igrd)
	397	DO ik = 1, jpkm1
	398	ubdydta(ib,ik,2) = zdta(nbmap(ib,igrd),1,ik)
	399	END DO
	400	END DO
	401	!
	402	igrd = 3 ! v-velocity
	403	IF ( nblendta(igrd) .le. 0 ) THEN
	404	idvar = iom_varid( numbdyv,'vomecrty' )
	405	nblendta(igrd) = iom_file(numbdyv)%dimsz(1,idvar)
	406	ENDIF
	407	WRITE(numout,*) 'Dim size for vomecrty is ', nblendta(igrd)
	408	ipi = nblendta(igrd)
	409	CALL iom_get ( numbdyv, jpdom_unknown,'vomecrty',zdta(1:ipi,1:ipj,1:ipk),nbdy_a )
	410	DO ib = 1, nblen(igrd)
	411	DO ik = 1, jpkm1
	412	vbdydta(ib,ik,2) = zdta(nbmap(ib,igrd),1,ik)
	413	END DO
	414	END DO
	415	ENDIF ! ln_bdy_dyn_frs
[911]	416
	417
[1125]	418	IF ((.NOT.ln_bdy_clim) .AND. (istep(1) > 0)) THEN
	419	! First data time is after start of run
	420	! Put first value in both time levels
	421	nbdy_b = nbdy_a
	422	IF(ln_bdy_tra_frs) THEN
	423	tbdydta(:,:,1) = tbdydta(:,:,2)
	424	sbdydta(:,:,1) = sbdydta(:,:,2)
	425	ENDIF
	426	IF(ln_bdy_dyn_frs) THEN
	427	ubdydta(:,:,1) = ubdydta(:,:,2)
	428	vbdydta(:,:,1) = vbdydta(:,:,2)
	429	ENDIF
	430	END IF
[911]	431
[1125]	432	END IF ! nbdy_dta == 0/1
[911]	433
[1125]	434	! In the case of constant boundary forcing fill bdy arrays once for all
	435	IF ((ln_bdy_clim).AND.(ntimes_bdy==1)) THEN
	436	IF(ln_bdy_tra_frs) THEN
	437	tbdy (:,:) = tbdydta (:,:,2)
	438	sbdy (:,:) = sbdydta (:,:,2)
	439	ENDIF
	440	IF(ln_bdy_dyn_frs) THEN
	441	ubdy (:,:) = ubdydta (:,:,2)
	442	vbdy (:,:) = vbdydta (:,:,2)
	443	ENDIF
[911]	444
[1125]	445	IF(ln_bdy_tra_frs .or. ln_bdy_ice_frs) CALL iom_close( numbdyt )
	446	IF(ln_bdy_dyn_frs) CALL iom_close( numbdyu )
	447	IF(ln_bdy_dyn_frs) CALL iom_close( numbdyv )
	448	END IF
[911]	449
[1125]	450	ENDIF ! End if nit000
[911]	451
	452
[1125]	453	! !---------------------!
	454	! ! at each time step !
	455	! !---------------------!
[911]	456
[1125]	457	IF( nbdy_dta == 1 .AND. ntimes_bdy > 1 ) THEN
	458	!
	459	! Read one more record if necessary
	460	!**********************************
[911]	461
[1125]	462	IF( ln_bdy_clim .AND. imois /= nbdy_b ) THEN ! remember that nbdy_b=0 for kt=nit000
	463	nbdy_b = imois
	464	nbdy_a = imois + 1
	465	nbdy_b = MOD( nbdy_b, iman ) ; IF( nbdy_b == 0 ) nbdy_b = iman
	466	nbdy_a = MOD( nbdy_a, iman ) ; IF( nbdy_a == 0 ) nbdy_a = iman
	467	lect=.true.
	468	ELSEIF( .NOT.ln_bdy_clim .AND. itimer >= istep(nbdy_a) ) THEN
[911]	469
[1125]	470	IF ( nbdy_a < ntimes_bdy ) THEN
	471	nbdy_b = nbdy_a
	472	nbdy_a = nbdy_a + 1
	473	lect =.true.
	474	ELSE
	475	! We have reached the end of the file
	476	! put the last data time into both time levels
	477	nbdy_b = nbdy_a
	478	IF(ln_bdy_tra_frs) THEN
	479	tbdydta(:,:,1) = tbdydta(:,:,2)
	480	sbdydta(:,:,1) = sbdydta(:,:,2)
	481	ENDIF
	482	IF(ln_bdy_dyn_frs) THEN
	483	ubdydta(:,:,1) = ubdydta(:,:,2)
	484	vbdydta(:,:,1) = vbdydta(:,:,2)
	485	ENDIF
	486	END IF ! nbdy_a < ntimes_bdy
[911]	487
	488	END IF
	489
[1125]	490	IF( lect ) THEN
	491	! Swap arrays
	492	IF(ln_bdy_tra_frs) THEN
	493	tbdydta(:,:,1) = tbdydta(:,:,2)
	494	sbdydta(:,:,1) = sbdydta(:,:,2)
	495	ENDIF
	496	IF(ln_bdy_dyn_frs) THEN
	497	ubdydta(:,:,1) = ubdydta(:,:,2)
	498	vbdydta(:,:,1) = vbdydta(:,:,2)
	499	ENDIF
[911]	500
[1125]	501	! read another set
	502	ipj = 1
	503	ipk = jpk
[911]	504
[1125]	505	IF(ln_bdy_tra_frs) THEN
	506	!
	507	igrd = 1 ! temperature
	508	ipi = nblendta(igrd)
	509	CALL iom_get ( numbdyt, jpdom_unknown, 'votemper', zdta(1:ipi,1:ipj,1:ipk), nbdy_a )
	510	DO ib = 1, nblen(igrd)
	511	DO ik = 1, jpkm1
	512	tbdydta(ib,ik,2) = zdta(nbmap(ib,igrd),1,ik)
	513	END DO
	514	END DO
	515	!
	516	igrd = 1 ! salinity
	517	ipi = nblendta(igrd)
	518	CALL iom_get ( numbdyt, jpdom_unknown, 'vosaline', zdta(1:ipi,1:ipj,1:ipk), nbdy_a )
	519	DO ib = 1, nblen(igrd)
	520	DO ik = 1, jpkm1
	521	sbdydta(ib,ik,2) = zdta(nbmap(ib,igrd),1,ik)
	522	END DO
	523	END DO
	524	ENDIF ! ln_bdy_tra_frs
[911]	525
[1125]	526	IF(ln_bdy_dyn_frs) THEN
	527	!
	528	igrd = 2 ! u-velocity
	529	ipi = nblendta(igrd)
	530	CALL iom_get ( numbdyu, jpdom_unknown,'vozocrtx',zdta(1:ipi,1:ipj,1:ipk),nbdy_a )
	531	DO ib = 1, nblen(igrd)
	532	DO ik = 1, jpkm1
	533	ubdydta(ib,ik,2) = zdta(nbmap(ib,igrd),1,ik)
	534	END DO
	535	END DO
	536	!
	537	igrd = 3 ! v-velocity
	538	ipi = nblendta(igrd)
	539	CALL iom_get ( numbdyv, jpdom_unknown,'vomecrty',zdta(1:ipi,1:ipj,1:ipk),nbdy_a )
	540	DO ib = 1, nblen(igrd)
	541	DO ik = 1, jpkm1
	542	vbdydta(ib,ik,2) = zdta(nbmap(ib,igrd),1,ik)
	543	END DO
	544	END DO
	545	ENDIF ! ln_bdy_dyn_frs
[911]	546
[1125]	547	!
	548	IF(lwp) WRITE(numout,*) 'bdy_dta : first record file used nbdy_b ',nbdy_b
	549	IF(lwp) WRITE(numout,*) '~~~~~~~~ last record file used nbdy_a ',nbdy_a
	550	IF (.NOT.ln_bdy_clim) THEN
	551	IF(lwp) WRITE(numout,*) 'first record time (s): ', istep(nbdy_b)
	552	IF(lwp) WRITE(numout,*) 'model time (s) : ', itimer
	553	IF(lwp) WRITE(numout,*) 'second record time (s): ', istep(nbdy_a)
	554	ENDIF
	555	!
	556	ENDIF ! end lect=.true.
[911]	557
	558
[1125]	559	! Interpolate linearly
	560	! ********************
	561	!
	562	IF( ln_bdy_clim ) THEN ; zxy = FLOAT( nday ) / FLOAT( nobis(nbdy_b) ) + 0.5 - i15
	563	ELSE ; zxy = FLOAT( istep(nbdy_b) - itimer ) / FLOAT( istep(nbdy_b) - istep(nbdy_a) )
	564	END IF
[911]	565
[1125]	566	IF(ln_bdy_tra_frs) THEN
	567	igrd = 1 ! temperature & salinity
	568	DO ib = 1, nblen(igrd)
	569	DO ik = 1, jpkm1
	570	tbdy(ib,ik) = zxy * tbdydta(ib,ik,2) + (1.-zxy) * tbdydta(ib,ik,1)
	571	sbdy(ib,ik) = zxy * sbdydta(ib,ik,2) + (1.-zxy) * sbdydta(ib,ik,1)
	572	END DO
	573	END DO
	574	ENDIF
[911]	575
[1125]	576	IF(ln_bdy_dyn_frs) THEN
	577	igrd = 2 ! u-velocity
	578	DO ib = 1, nblen(igrd)
	579	DO ik = 1, jpkm1
	580	ubdy(ib,ik) = zxy * ubdydta(ib,ik,2) + (1.-zxy) * ubdydta(ib,ik,1)
	581	END DO
	582	END DO
	583	!
	584	igrd = 3 ! v-velocity
	585	DO ib = 1, nblen(igrd)
	586	DO ik = 1, jpkm1
	587	vbdy(ib,ik) = zxy * vbdydta(ib,ik,2) + (1.-zxy) * vbdydta(ib,ik,1)
	588	END DO
	589	END DO
	590	ENDIF
[911]	591
[1125]	592	END IF !end if ((nbdy_dta==1).AND.(ntimes_bdy>1))
[911]	593
	594
[1125]	595	! !---------------------!
	596	! ! last call !
	597	! !---------------------!
[911]	598	IF( kt == nitend ) THEN
[1125]	599	IF(ln_bdy_tra_frs .or. ln_bdy_ice_frs) CALL iom_close( numbdyt ) ! Closing of the 3 files
	600	IF(ln_bdy_dyn_frs) CALL iom_close( numbdyu )
	601	IF(ln_bdy_dyn_frs) CALL iom_close( numbdyv )
[911]	602	ENDIF
[1125]	603	!
	604	ENDIF ! ln_bdy_dyn_frs .OR. ln_bdy_tra_frs
[911]	605
	606	END SUBROUTINE bdy_dta
	607
[1125]	608
[911]	609	SUBROUTINE bdy_dta_bt( kt, jit )
	610	!!---------------------------------------------------------------------------
	611	!! * SUBROUTINE bdy_dta_bt *
	612	!!
[1125]	613	!! ** Purpose : Read unstructured boundary data for Flather condition
[911]	614	!!
[1125]	615	!! ** Method : At the first timestep, read in boundary data for two
	616	!! times from the file and time-interpolate. At other
	617	!! timesteps, check to see if we need another time from
	618	!! the file. If so read it in. Time interpolate.
[911]	619	!!---------------------------------------------------------------------------
[1125]	620	!!gm DOCTOR names : argument integer : start with "k"
[911]	621	INTEGER, INTENT( in ) :: kt ! ocean time-step index
	622	INTEGER, INTENT( in ) :: jit ! barotropic time step index
[1125]	623	! ! (for timesplitting option, otherwise zero)
	624	!!
[911]	625	LOGICAL :: lect ! flag for reading
[1125]	626	INTEGER :: it, ib, igrd ! dummy loop indices
[911]	627	INTEGER :: idvar ! netcdf var ID
	628	INTEGER :: iman, i15, imois ! Time variables for monthly clim forcing
[1125]	629	INTEGER :: ntimes_bdyt, ntimes_bdyu, ntimes_bdyv
[911]	630	INTEGER :: itimer, totime
	631	INTEGER :: ipi, ipj, ipk, inum ! temporary integers (NetCDF read)
	632	INTEGER :: iyear0, imonth0, iday0
	633	INTEGER :: ihours0, iminutes0, isec0
[1125]	634	INTEGER :: iyear, imonth, iday, isecs
[911]	635	INTEGER, DIMENSION(jpbtime) :: istept, istepu, istepv ! time arrays from data files
[1125]	636	REAL(wp) :: dayfrac, zxy, zoffsett
	637	REAL(wp) :: zoffsetu, zoffsetv
[1126]	638	REAL(wp) :: dayjul0, zdayjulini
[1125]	639	REAL(wp) :: zinterval_s, zinterval_e ! First and last interval in time axis
	640	REAL(wp), DIMENSION(jpbtime) :: zstepr ! REAL time array from data files
	641	REAL(wp), DIMENSION(jpbdta,1) :: zdta ! temporary array for data fields
	642	CHARACTER(LEN=80), DIMENSION(3) :: clfile
[911]	643	CHARACTER(LEN=70 ) :: clunits ! units attribute of time coordinate
	644	!!---------------------------------------------------------------------------
	645
[1125]	646	!!gm add here the same style as in bdy_dta
	647	!!gm clearly bdy_dta_bt and bdy_dta can be combined...
	648	!!gm too many things duplicated in the read of data... simplification can be done
	649
[911]	650	! -------------------- !
	651	! Initialization !
	652	! -------------------- !
	653
	654	lect = .false. ! If true, read a time record
	655
	656	! Some time variables for monthly climatological forcing:
	657	! *******************************************************
[1125]	658	!!gm here use directely daymod variables
[911]	659
	660	iman = INT( raamo ) ! Number of months in a year
	661
	662	i15 = INT( 2*FLOAT( nday ) / ( FLOAT( nobis(nmonth) ) + 0.5 ) )
	663	! i15=0 if the current day is in the first half of the month, else i15=1
	664
	665	imois = nmonth + i15 - 1 ! imois is the first month record
	666	IF( imois == 0 ) imois = iman
	667
	668	! Time variable for non-climatological forcing:
	669	! *********************************************
	670
[1241]	671	itimer = ((kt-1)-nit000+1)*rdt ! current time in seconds for interpolation
	672	itimer = itimer + jit*rdt/REAL(nn_baro,wp) ! in non-climatological case
[911]	673
[1125]	674	IF ( ln_bdy_tides ) THEN
[911]	675
[1125]	676	! -------------------------------------!
	677	! Update BDY fields with tidal forcing !
	678	! -------------------------------------!
	679
	680	CALL tide_update( kt, jit )
[911]	681
[1125]	682	ENDIF
[911]	683
[1125]	684	IF ( ln_bdy_dyn_fla ) THEN
[911]	685
[1125]	686	! -------------------------------------!
	687	! Update BDY fields with model data !
	688	! -------------------------------------!
[911]	689
[1125]	690	! !-------------------!
	691	IF( kt == nit000 ) THEN ! First call only !
	692	! !-------------------!
	693	istep_bt(:) = 0
	694	nbdy_b_bt = 0
	695	nbdy_a_bt = 0
[911]	696
[1125]	697	! Get time information from bdy data file
	698	! ***************************************
[911]	699
	700	IF(lwp) WRITE(numout,*)
	701	IF(lwp) WRITE(numout,*) 'bdy_dta_bt :Initialize unstructured boundary data for barotropic variables.'
	702	IF(lwp) WRITE(numout,*) '~~~~~~~'
	703
[1125]	704	IF( nbdy_dta == 0 ) THEN
[911]	705	IF(lwp) WRITE(numout,*) 'Bdy data are taken from initial conditions'
	706
	707	ELSEIF (nbdy_dta == 1) THEN
	708	IF(lwp) WRITE(numout,*) 'Bdy data are read in netcdf files'
	709
	710	dayfrac = adatrj-float(itimer)/86400. ! day fraction at time step kt-1
	711	dayfrac = dayfrac - INT(dayfrac) !
	712	totime = (nitend-nit000+1)*rdt ! Total time of the run to verify that all the
	713	! necessary time dumps in file are included
	714
[1125]	715	clfile(1) = filbdy_data_bt_T
	716	clfile(2) = filbdy_data_bt_U
	717	clfile(3) = filbdy_data_bt_V
[911]	718
[1125]	719	DO igrd = 1,3
[911]	720
[1125]	721	CALL iom_open( clfile(igrd), inum )
	722	CALL iom_gettime( inum, zstepr, kntime=ntimes_bdy, cdunits=clunits )
[911]	723
[1125]	724	SELECT CASE( igrd )
	725	CASE (1)
	726	numbdyt = inum
	727	CASE (2)
	728	numbdyu = inum
	729	CASE (3)
	730	numbdyv = inum
	731	END SELECT
	732
[911]	733	! Calculate time offset
	734	READ(clunits,7000) iyear0, imonth0, iday0, ihours0, iminutes0, isec0
	735	! Convert time origin in file to julian days
	736	isec0 = isec0 + ihours060.60. + iminutes0*60.
	737	CALL ymds2ju(iyear0, imonth0, iday0, real(isec0), dayjul0)
	738	! Compute model initialization time
[1125]	739	iyear = ndastp / 10000
	740	imonth = ( ndastp - iyear * 10000 ) / 100
	741	iday = ndastp - iyear * 10000 - imonth * 100
	742	isecs = dayfrac * 86400
	743	CALL ymds2ju(iyear, imonth, iday, real(isecs) , zdayjulini)
	744	! zoffset from initialization date:
	745	zoffset = (dayjul0-zdayjulini)*86400
[911]	746	!
	747
	748	7000 FORMAT('seconds since ', I4.4,'-',I2.2,'-',I2.2,' ',I2.2,':',I2.2,':',I2.2)
	749
	750	!! TO BE DONE... Check consistency between calendar from file
	751	!! (available optionally from iom_gettime) and calendar in model
	752	!! when calendar in model available outside of IOIPSL.
	753
	754	! Check that there are not too many times in the file.
[1125]	755	IF (ntimes_bdy_bt > jpbtime) CALL ctl_stop( &
	756	'Number of time dumps in bdy file exceed jpbtime parameter', &
	757	'Check file:' // TRIM(clfile(igrd)) )
[911]	758
[1125]	759	! Check that time array increases (or interp will fail):
	760	DO it = 2, ntimes_bdy
	761	IF ( zstepr(it-1) >= zstepr(it) ) THEN
	762	CALL ctl_stop('Time array in unstructured boundary data file', &
	763	'does not continuously increase.', &
	764	'Check file:' // TRIM(clfile(igrd)) )
	765	EXIT
	766	END IF
[911]	767	END DO
	768
[1125]	769	IF ( .NOT. ln_bdy_clim ) THEN
	770	! Check that times in file span model run time:
[911]	771
[1125]	772	! Note: the fields may be time means, so we allow nit000 to be before
	773	! first time in the file, provided that it falls inside the meaning
	774	! period of the first field. Until we can get the meaning period
	775	! from the file, use the interval between fields as a proxy.
	776	! If nit000 is before the first time, use the value at first time
	777	! instead of extrapolating. This is done by putting time 1 into
	778	! both time levels.
	779	! The same applies to the last time level: see setting of lect below.
[911]	780
[1125]	781	IF ( ntimes_bdy == 1 ) CALL ctl_stop( &
	782	'There is only one time dump in data files', &
	783	'Set ln_bdy_clim=.true. in namelist for constant bdy forcing.' )
[911]	784
[1125]	785	zinterval_s = zstepr(2) - zstepr(1)
	786	zinterval_e = zstepr(ntimes_bdy) - zstepr(ntimes_bdy-1)
	787
	788	IF ( zstepr(1) - zinterval_s / 2.0 > 0 ) THEN
	789	IF(lwp) WRITE(numout,*) 'First bdy time relative to nit000:', zstepr(1)
	790	IF(lwp) WRITE(numout,*) 'Interval between first two times: ', zinterval_s
	791	CALL ctl_stop( 'First data time is after start of run', &
	792	'by more than half a meaning period', &
	793	'Check file: ' // TRIM(clfile(igrd)) )
	794	END IF
	795
[1126]	796	IF ( zstepr(ntimes_bdy) + zinterval_e / 2.0 < totime ) THEN
[1125]	797	IF(lwp) WRITE(numout,*) 'Last bdy time relative to nit000:', zstepr(ntimes_bdy)
	798	IF(lwp) WRITE(numout,*) 'Interval between last two times: ', zinterval_e
	799	CALL ctl_stop( 'Last data time is before end of run', &
	800	'by more than half a meaning period', &
	801	'Check file: ' // TRIM(clfile(igrd)) )
	802	END IF
	803
	804	END IF ! .NOT. ln_bdy_clim
	805
	806	IF ( igrd .EQ. 1) THEN
	807	ntimes_bdyt = ntimes_bdy_bt
	808	zoffsett = zoffset
	809	istept(:) = INT( zstepr(:) + zoffset )
	810	ELSE IF (igrd .EQ. 2) THEN
	811	ntimes_bdyu = ntimes_bdy_bt
	812	zoffsetu = zoffset
	813	istepu(:) = INT( zstepr(:) + zoffset )
	814	ELSE IF (igrd .EQ. 3) THEN
	815	ntimes_bdyv = ntimes_bdy_bt
	816	zoffsetv = zoffset
	817	istepv(:) = INT( zstepr(:) + zoffset )
[911]	818	ENDIF
	819
	820	ENDDO
	821
	822	! Verify time consistency between files
	823
[1125]	824	IF ( ntimes_bdyu /= ntimes_bdyt .OR. ntimes_bdyv /= ntimes_bdyt ) THEN
	825	CALL ctl_stop( &
	826	'Time axis lengths differ between bdy data files', &
	827	'Multiple time frequencies not implemented yet' )
	828	ELSE
	829	ntimes_bdy_bt = ntimes_bdyt
[911]	830	ENDIF
	831
[1125]	832	IF (zoffsetu.NE.zoffsett .OR. zoffsetv.NE.zoffsett) THEN
	833	CALL ctl_stop( &
	834	'Bdy data files must have the same time origin', &
	835	'Multiple time frequencies not implemented yet' )
[911]	836	ENDIF
[1125]	837	zoffset = zoffsett
[911]	838
	839	!! Check that times are the same in the three files... HERE.
	840	istep_bt(:) = istept(:)
	841
	842	! Check number of time dumps:
	843	IF (ln_bdy_clim) THEN
[1125]	844	SELECT CASE ( ntimes_bdy_bt )
	845	CASE( 1 )
[911]	846	IF(lwp) WRITE(numout,*)
	847	IF(lwp) WRITE(numout,*) 'We assume constant boundary forcing from bdy data files'
	848	IF(lwp) WRITE(numout,*)
[1125]	849	CASE( 12 )
[911]	850	IF(lwp) WRITE(numout,*)
	851	IF(lwp) WRITE(numout,*) 'We assume monthly (and cyclic) boundary forcing from bdy data files'
	852	IF(lwp) WRITE(numout,*)
[1125]	853	CASE DEFAULT
	854	CALL ctl_stop( &
	855	'For climatological boundary forcing (ln_bdy_clim=.true.),',&
	856	'bdy data files must contain 1 or 12 time dumps.' )
	857	END SELECT
[911]	858	ENDIF
	859
	860	! Find index of first record to read (before first model time).
	861
[1125]	862	it=1
	863	DO WHILE ( ((istep_bt(it+1)) <= 0 ).AND.(it.LE.(ntimes_bdy_bt-1)))
	864	it=it+1
[911]	865	END DO
[1125]	866	nbdy_b_bt = it
[911]	867
[1125]	868	WRITE(numout,*) 'Time offset is ',zoffset
	869	WRITE(numout,*) 'First record to read is ',nbdy_b_bt
[911]	870
	871	ENDIF ! endif (nbdy_dta == 1)
	872
	873	! 1.2 Read first record in file if necessary (ie if nbdy_dta == 1)
	874	! *****************************************************************
	875
	876	IF ( nbdy_dta == 0) THEN
	877	! boundary data arrays are filled with initial conditions
[1125]	878	igrd = 2 ! U-points data
	879	DO ib = 1, nblen(igrd)
	880	ubtbdy(ib) = un(nbi(ib,igrd), nbj(ib,igrd), 1)
[911]	881	END DO
	882
[1125]	883	igrd = 3 ! V-points data
	884	DO ib = 1, nblen(igrd)
	885	vbtbdy(ib) = vn(nbi(ib,igrd), nbj(ib,igrd), 1)
[911]	886	END DO
	887
[1125]	888	igrd = 1 ! T-points data
	889	DO ib = 1, nblen(igrd)
	890	sshbdy(ib) = sshn(nbi(ib,igrd), nbj(ib,igrd))
[911]	891	END DO
	892
	893	ELSEIF (nbdy_dta == 1) THEN
	894
	895	! Set first record in the climatological case:
[1125]	896	IF ((ln_bdy_clim).AND.(ntimes_bdy_bt==1)) THEN
	897	nbdy_a_bt = 1
	898	ELSEIF ((ln_bdy_clim).AND.(ntimes_bdy_bt==iman)) THEN
	899	nbdy_b_bt = 0
	900	nbdy_a_bt = imois
[911]	901	ELSE
[1125]	902	nbdy_a_bt = nbdy_b_bt
[911]	903	END IF
	904
	905	! Open Netcdf files:
	906
[1125]	907	CALL iom_open ( filbdy_data_bt_T, numbdyt_bt )
	908	CALL iom_open ( filbdy_data_bt_U, numbdyu_bt )
	909	CALL iom_open ( filbdy_data_bt_V, numbdyv_bt )
[911]	910
	911	! Read first record:
	912	ipj=1
[1125]	913	igrd=1
	914	ipi=nblendta(igrd)
[911]	915
	916	! ssh
[1125]	917	igrd=1
	918	IF ( nblendta(igrd) .le. 0 ) THEN
	919	idvar = iom_varid( numbdyt_bt,'sossheig' )
	920	nblendta(igrd) = iom_file(numbdyt_bt)%dimsz(1,idvar)
[911]	921	ENDIF
[1125]	922	WRITE(numout,*) 'Dim size for sossheig is ',nblendta(igrd)
	923	ipi=nblendta(igrd)
[911]	924
[1125]	925	CALL iom_get ( numbdyt_bt, jpdom_unknown,'sossheig',zdta(1:ipi,1:ipj),nbdy_a_bt )
[911]	926
[1125]	927	DO ib=1, nblen(igrd)
	928	sshbdydta(ib,2) = zdta(nbmap(ib,igrd),1)
[911]	929	END DO
	930
	931	! u-velocity
[1125]	932	igrd=2
	933	IF ( nblendta(igrd) .le. 0 ) THEN
	934	idvar = iom_varid( numbdyu_bt,'vobtcrtx' )
	935	nblendta(igrd) = iom_file(numbdyu_bt)%dimsz(1,idvar)
[911]	936	ENDIF
[1125]	937	WRITE(numout,*) 'Dim size for vobtcrtx is ',nblendta(igrd)
	938	ipi=nblendta(igrd)
[911]	939
[1125]	940	CALL iom_get ( numbdyu_bt, jpdom_unknown,'vobtcrtx',zdta(1:ipi,1:ipj),nbdy_a_bt )
[911]	941
[1125]	942	DO ib=1, nblen(igrd)
	943	ubtbdydta(ib,2) = zdta(nbmap(ib,igrd),1)
[911]	944	END DO
	945
	946	! v-velocity
[1125]	947	igrd=3
	948	IF ( nblendta(igrd) .le. 0 ) THEN
	949	idvar = iom_varid( numbdyv_bt,'vobtcrty' )
	950	nblendta(igrd) = iom_file(numbdyv_bt)%dimsz(1,idvar)
[911]	951	ENDIF
[1125]	952	WRITE(numout,*) 'Dim size for vobtcrty is ',nblendta(igrd)
	953	ipi=nblendta(igrd)
[911]	954
[1125]	955	CALL iom_get ( numbdyv_bt, jpdom_unknown,'vobtcrty',zdta(1:ipi,1:ipj),nbdy_a_bt )
[911]	956
[1125]	957	DO ib=1, nblen(igrd)
	958	vbtbdydta(ib,2) = zdta(nbmap(ib,igrd),1)
[911]	959	END DO
	960
	961	END IF
	962
	963	! In the case of constant boundary forcing fill bdy arrays once for all
[1125]	964	IF ((ln_bdy_clim).AND.(ntimes_bdy_bt==1)) THEN
[911]	965
	966	ubtbdy (:) = ubtbdydta (:,2)
	967	vbtbdy (:) = vbtbdydta (:,2)
	968	sshbdy (:) = sshbdydta (:,2)
	969
[1125]	970	CALL iom_close( numbdyt_bt )
	971	CALL iom_close( numbdyu_bt )
	972	CALL iom_close( numbdyv_bt )
[911]	973
	974	END IF
	975
	976	ENDIF ! End if nit000
	977
	978	! -------------------- !
	979	! 2. At each time step !
	980	! -------------------- !
	981
[1125]	982	IF ((nbdy_dta==1).AND.(ntimes_bdy_bt>1)) THEN
[911]	983
	984	! 2.1 Read one more record if necessary
	985	!**************************************
	986
[1125]	987	IF ( (ln_bdy_clim).AND.(imois/=nbdy_b_bt) ) THEN ! remember that nbdy_b_bt=0 for kt=nit000
	988	nbdy_b_bt = imois
	989	nbdy_a_bt = imois+1
	990	nbdy_b_bt = MOD( nbdy_b_bt, iman )
	991	IF( nbdy_b_bt == 0 ) nbdy_b_bt = iman
	992	nbdy_a_bt = MOD( nbdy_a_bt, iman )
	993	IF( nbdy_a_bt == 0 ) nbdy_a_bt = iman
[911]	994	lect=.true.
	995
[1125]	996	ELSEIF ((.NOT.ln_bdy_clim).AND.(itimer >= istep_bt(nbdy_a_bt))) THEN
	997	nbdy_b_bt=nbdy_a_bt
	998	nbdy_a_bt=nbdy_a_bt+1
[911]	999	lect=.true.
	1000	END IF
	1001
	1002	IF (lect) THEN
	1003
	1004	! Swap arrays
	1005	sshbdydta(:,1) = sshbdydta(:,2)
	1006	ubtbdydta(:,1) = ubtbdydta(:,2)
	1007	vbtbdydta(:,1) = vbtbdydta(:,2)
	1008
	1009	! read another set
	1010
	1011	ipj=1
	1012	ipk=jpk
[1125]	1013	igrd=1
	1014	ipi=nblendta(igrd)
[911]	1015
	1016
	1017	! ssh
[1125]	1018	igrd=1
	1019	ipi=nblendta(igrd)
[911]	1020
[1125]	1021	CALL iom_get ( numbdyt_bt, jpdom_unknown,'sossheig',zdta(1:ipi,1:ipj),nbdy_a_bt )
[911]	1022
[1125]	1023	DO ib=1, nblen(igrd)
	1024	sshbdydta(ib,2) = zdta(nbmap(ib,igrd),1)
[911]	1025	END DO
	1026
	1027	! u-velocity
[1125]	1028	igrd=2
	1029	ipi=nblendta(igrd)
[911]	1030
[1125]	1031	CALL iom_get ( numbdyu_bt, jpdom_unknown,'vobtcrtx',zdta(1:ipi,1:ipj),nbdy_a_bt )
[911]	1032
[1125]	1033	DO ib=1, nblen(igrd)
	1034	ubtbdydta(ib,2) = zdta(nbmap(ib,igrd),1)
[911]	1035	END DO
	1036
	1037	! v-velocity
[1125]	1038	igrd=3
	1039	ipi=nblendta(igrd)
[911]	1040
[1125]	1041	CALL iom_get ( numbdyv_bt, jpdom_unknown,'vobtcrty',zdta(1:ipi,1:ipj),nbdy_a_bt )
[911]	1042
[1125]	1043	DO ib=1, nblen(igrd)
	1044	vbtbdydta(ib,2) = zdta(nbmap(ib,igrd),1)
[911]	1045	END DO
	1046
	1047
[1125]	1048	IF(lwp) WRITE(numout,*) 'bdy_dta : first record file used nbdy_b_bt ',nbdy_b_bt
	1049	IF(lwp) WRITE(numout,*) '~~~~~~~~ last record file used nbdy_a_bt ',nbdy_a_bt
[911]	1050	IF (.NOT.ln_bdy_clim) THEN
[1125]	1051	IF(lwp) WRITE(numout,*) 'first record time (s): ', istep_bt(nbdy_b_bt)
[911]	1052	IF(lwp) WRITE(numout,*) 'model time (s) : ', itimer
[1125]	1053	IF(lwp) WRITE(numout,*) 'second record time (s): ', istep_bt(nbdy_a_bt)
[911]	1054	ENDIF
	1055	END IF ! end lect=.true.
	1056
	1057
	1058	! 2.2 Interpolate linearly:
	1059	! ***************************
	1060
	1061	IF (ln_bdy_clim) THEN
[1125]	1062	zxy = FLOAT( nday ) / FLOAT( nobis(nbdy_b_bt) ) + 0.5 - i15
[911]	1063	ELSE
[1125]	1064	zxy = FLOAT(istep_bt(nbdy_b_bt)-itimer) / FLOAT(istep_bt(nbdy_b_bt)-istep_bt(nbdy_a_bt))
[911]	1065	END IF
	1066
[1125]	1067	igrd=1
	1068	DO ib=1, nblen(igrd)
	1069	sshbdy(ib) = zxy * sshbdydta(ib,2) + &
	1070	(1.-zxy) * sshbdydta(ib,1)
[911]	1071	END DO
	1072
[1125]	1073	igrd=2
	1074	DO ib=1, nblen(igrd)
	1075	ubtbdy(ib) = zxy * ubtbdydta(ib,2) + &
	1076	(1.-zxy) * ubtbdydta(ib,1)
[911]	1077	END DO
	1078
[1125]	1079	igrd=3
	1080	DO ib=1, nblen(igrd)
	1081	vbtbdy(ib) = zxy * vbtbdydta(ib,2) + &
	1082	(1.-zxy) * vbtbdydta(ib,1)
[911]	1083	END DO
	1084
	1085
[1125]	1086	END IF !end if ((nbdy_dta==1).AND.(ntimes_bdy_bt>1))
[911]	1087
	1088	! ------------------- !
	1089	! Last call kt=nitend !
	1090	! ------------------- !
	1091
	1092	! Closing of the 3 files
	1093	IF( kt == nitend ) THEN
[1125]	1094	CALL iom_close( numbdyt_bt )
	1095	CALL iom_close( numbdyu_bt )
	1096	CALL iom_close( numbdyv_bt )
[911]	1097	ENDIF
	1098
[1125]	1099	ENDIF ! ln_bdy_dyn_frs
[911]	1100
	1101	END SUBROUTINE bdy_dta_bt
	1102
	1103
	1104	#else
[1125]	1105	!!----------------------------------------------------------------------
	1106	!! Dummy module NO Unstruct Open Boundary Conditions
	1107	!!----------------------------------------------------------------------
[911]	1108	CONTAINS
[1125]	1109	SUBROUTINE bdy_dta( kt ) ! Empty routine
[911]	1110	WRITE(,) 'bdy_dta: You should not have seen this print! error?', kt
	1111	END SUBROUTINE bdy_dta
[1125]	1112	SUBROUTINE bdy_dta_bt( kt, kit ) ! Empty routine
	1113	WRITE(,) 'bdy_dta: You should not have seen this print! error?', kt, kit
[911]	1114	END SUBROUTINE bdy_dta_bt
	1115	#endif
	1116
	1117	!!==============================================================================
	1118	END MODULE bdydta

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

Download in other formats: