Changeset 12255 for NEMO/branches/2019/dev_r11943_MERGE_2019/src/TOP/trcbc.F90

Timestamp:

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

Author:

smasson

Message:

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

File:

• NEMO/branches/2019/dev_r11943_MERGE_2019/src/TOP/trcbc.F90 (modified) (2 diffs)

NEMO/branches/2019/dev_r11943_MERGE_2019/src/TOP/trcbc.F90

@@ -367 +367 @@
       IF( PRESENT(jit) ) THEN 
          !
-         ! OPEN boundary conditions (use pt_offset=1. as they are applied at the end of the step)
+         ! BDY: use pt_offset=0.5 as applied at the end of the step and fldread is referenced at the middle of the step
          IF( nb_trcobc > 0 ) THEN
            if (lwp) write(numout,'(a,i5,a,i10)') '   reading OBC data for ', nb_trcobc ,' variable(s) at step ', kt
-           CALL fld_read( kt=kt, kn_fsbc=1, sd=sf_trcobc, kit=jit, pt_offset=1.)
+           CALL fld_read( kt=kt, kn_fsbc=1, sd=sf_trcobc, kit=jit, pt_offset = 0.5_wp )
          ENDIF
          !
@@ -387 +387 @@
       ELSE
          !
-         ! OPEN boundary conditions (use pt_offset=1. as they are applied at the end of the step)
+         ! BDY: use pt_offset=0.5 as applied at the end of the step and fldread is referenced at the middle of the step
          IF( nb_trcobc > 0 ) THEN
            if (lwp) write(numout,'(a,i5,a,i10)') '   reading OBC data for ', nb_trcobc ,' variable(s) at step ', kt
-           CALL fld_read( kt=kt, kn_fsbc=1, sd=sf_trcobc, pt_offset=1.)
+           CALL fld_read( kt=kt, kn_fsbc=1, sd=sf_trcobc, pt_offset = 0.5_wp )
          ENDIF
          !