Changeset 12256 for NEMO/branches/2019/dev_r12072_MERGE_OPTION2_2019/src/OCE/BDY/bdydta.F90

Timestamp:

2019-12-16T13:52:49+01:00 (4 years ago)

Author:

smasson

Message:

rev12232_dev_r12072_MERGE_OPTION2_2019: bugfix in fldread offset in bdy, see #2263. all sette tests ok

File:

• NEMO/branches/2019/dev_r12072_MERGE_OPTION2_2019/src/OCE/BDY/bdydta.F90 (modified) (4 diffs)

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

@@ -75 +75 @@
 CONTAINS
 
-   SUBROUTINE bdy_dta( kt, pt_offset )
+   SUBROUTINE bdy_dta( kt )
       !!----------------------------------------------------------------------
       !!                   ***  SUBROUTINE bdy_dta  ***
@@ -85 +85 @@
       !!----------------------------------------------------------------------
       INTEGER, INTENT(in)           ::   kt           ! ocean time-step index 
-      REAL(wp),INTENT(in), OPTIONAL ::   pt_offset    ! time offset in units of timesteps
       !
       INTEGER ::  jbdy, jfld, jstart, jend, ib, jl    ! dummy loop indices
@@ -210 +209 @@
          ! read/update all bdy data
          ! ------------------------
-         CALL fld_read( kt, 1, bf_alias, pt_offset = pt_offset )
+         ! BDY: use pt_offset=0.5 as applied at the end of the step and fldread is referenced at the middle of the step
+         CALL fld_read( kt, 1, bf_alias, pt_offset = 0.5_wp )
 
          ! apply some corrections in some specific cases...
@@ -341 +341 @@
          ELSE ! Add tides if not split-explicit free surface else this is done in ts loop
             !
-            CALL bdy_dta_tides( kt=kt, pt_offset=pt_offset )
+            ! 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
+            CALL bdy_dta_tides( kt=kt, pt_offset = 1._wp )
          ENDIF
       ENDIF