Changeset 12205 for NEMO/branches/2019/dev_r11943_MERGE_2019/src/OCE/BDY/bdydta.F90

Timestamp:

2019-12-12T11:52:50+01:00 (4 years ago)

Author:

acc

Message:

2019/dev_r11943_MERGE_2019: Merge in dev_r11879_ENHANCE-05_SimonM-Harmonic_Analysis. SETTE tested

File:

• NEMO/branches/2019/dev_r11943_MERGE_2019/src/OCE/BDY/bdydta.F90 (modified) (6 diffs)

NEMO/branches/2019/dev_r11943_MERGE_2019/src/OCE/BDY/bdydta.F90

@@ -23 +23 @@
    USE phycst         ! physical constants
    USE sbcapr         ! atmospheric pressure forcing
-   USE sbctide        ! Tidal forcing or not
+   USE tide_mod, ONLY: ln_tide ! tidal forcing
    USE bdy_oce        ! ocean open boundary conditions  
    USE bdytides       ! tidal forcing at boundaries
@@ -75 +75 @@
 CONTAINS
 
-   SUBROUTINE bdy_dta( kt, Kmm, kit, kt_offset )
+   SUBROUTINE bdy_dta( kt, Kmm, kt_offset )
       !!----------------------------------------------------------------------
       !!                   ***  SUBROUTINE bdy_dta  ***
@@ -86 +86 @@
       INTEGER, INTENT(in)           ::   kt           ! ocean time-step index 
       INTEGER, INTENT(in)           ::   Kmm          ! ocean time level index
-      INTEGER, INTENT(in), OPTIONAL ::   kit          ! subcycle time-step index (for timesplitting option)
-      INTEGER, INTENT(in), OPTIONAL ::   kt_offset    ! time offset in units of timesteps. NB. if kit
+      INTEGER, INTENT(in), OPTIONAL ::   kt_offset    ! time offset in units of timesteps
       !                                               ! is present then units = subcycle timesteps.
       !                                               ! kt_offset = 0 => get data at "now" time level
@@ -106 +105 @@
       ! Initialise data arrays once for all from initial conditions where required
       !---------------------------------------------------------------------------
-      IF( kt == nit000 .AND. .NOT.PRESENT(kit) ) THEN
+      IF( kt == nit000 ) THEN
 
          ! Calculate depth-mean currents
@@ -217 +216 @@
          ! read/update all bdy data
          ! ------------------------
-         CALL fld_read( kt, 1, bf_alias, kit = kit, kt_offset = kt_offset )
+         CALL fld_read( kt, 1, bf_alias, kt_offset = kt_offset )
          ! apply some corrections in some specific cases...
          ! --------------------------------------------------
@@ -275 +274 @@
             END DO
          ENDIF   ! ltotvel
-
-         ! update tidal harmonic forcing
-         IF( PRESENT(kit) .AND. nn_dyn2d_dta(jbdy) .GE. 2 ) THEN
-            CALL bdytide_update( kt = kt, idx = idx_bdy(jbdy), dta = dta_alias, td = tides(jbdy),   & 
-               &                 kit = kit, kt_offset = kt_offset )
-         ENDIF
 
          !  atm surface pressure : add inverted barometer effect to ssh if it was read