- Timestamp:
- 2019-12-05T19:31:55+01:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r10742_ENHANCE-12_SimonM-Tides/src/OCE/BDY/bdydta.F90
r12065 r12078 75 75 CONTAINS 76 76 77 SUBROUTINE bdy_dta( kt, k it, kt_offset )77 SUBROUTINE bdy_dta( kt, kt_offset ) 78 78 !!---------------------------------------------------------------------- 79 79 !! *** SUBROUTINE bdy_dta *** … … 85 85 !!---------------------------------------------------------------------- 86 86 INTEGER, INTENT(in) :: kt ! ocean time-step index 87 INTEGER, INTENT(in), OPTIONAL :: kit ! subcycle time-step index (for timesplitting option) 88 INTEGER, INTENT(in), OPTIONAL :: kt_offset ! time offset in units of timesteps. NB. if kit 89 ! ! is present then units = subcycle timesteps. 87 INTEGER, INTENT(in), OPTIONAL :: kt_offset ! time offset in units of timesteps 90 88 ! ! kt_offset = 0 => get data at "now" time level 91 89 ! ! kt_offset = -1 => get data at "before" time level … … 105 103 ! Initialise data arrays once for all from initial conditions where required 106 104 !--------------------------------------------------------------------------- 107 IF( kt == nit000 .AND. .NOT.PRESENT(kit)) THEN105 IF( kt == nit000 ) THEN 108 106 109 107 ! Calculate depth-mean currents … … 216 214 ! read/update all bdy data 217 215 ! ------------------------ 218 CALL fld_read( kt, 1, bf_alias, k it = kit, kt_offset = kt_offset )216 CALL fld_read( kt, 1, bf_alias, kt_offset = kt_offset ) 219 217 220 218 ! apply some corrections in some specific cases... … … 275 273 END DO 276 274 ENDIF ! ltotvel 277 278 ! update tidal harmonic forcing279 IF( PRESENT(kit) .AND. nn_dyn2d_dta(jbdy) .GE. 2 ) THEN280 CALL bdytide_update( kt = kt, idx = idx_bdy(jbdy), dta = dta_alias, td = tides(jbdy), &281 & kit = kit, kt_offset = kt_offset )282 ENDIF283 275 284 276 ! atm surface pressure : add inverted barometer effect to ssh if it was read
Note: See TracChangeset
for help on using the changeset viewer.