Ignore:
Timestamp:
2019-12-16T13:52:25+01:00 (11 months ago)
Author:
smasson
Message:

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

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r11943_MERGE_2019/src/OCE/BDY/bdydta.F90

    r12250 r12255  
    7575CONTAINS 
    7676 
    77    SUBROUTINE bdy_dta( kt, Kmm, pt_offset ) 
     77   SUBROUTINE bdy_dta( kt, Kmm ) 
    7878      !!---------------------------------------------------------------------- 
    7979      !!                   ***  SUBROUTINE bdy_dta  *** 
     
    8686      INTEGER, INTENT(in)           ::   kt           ! ocean time-step index  
    8787      INTEGER, INTENT(in)           ::   Kmm          ! ocean time level index 
    88       REAL(wp),INTENT(in), OPTIONAL ::   pt_offset    ! time offset in units of timesteps 
    8988      ! 
    9089      INTEGER ::  jbdy, jfld, jstart, jend, ib, jl    ! dummy loop indices 
     
    211210         ! read/update all bdy data 
    212211         ! ------------------------ 
    213          CALL fld_read( kt, 1, bf_alias, pt_offset = pt_offset, Kmm = Kmm ) 
     212         ! BDY: use pt_offset=0.5 as applied at the end of the step and fldread is referenced at the middle of the step 
     213         CALL fld_read( kt, 1, bf_alias, pt_offset = 0.5_wp, Kmm = Kmm ) 
    214214         ! apply some corrections in some specific cases... 
    215215         ! -------------------------------------------------- 
     
    341341         ELSE ! Add tides if not split-explicit free surface else this is done in ts loop 
    342342            ! 
    343             CALL bdy_dta_tides( kt=kt, pt_offset=pt_offset ) 
     343            ! 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 
     344            CALL bdy_dta_tides( kt=kt, pt_offset = 1._wp ) 
    344345         ENDIF 
    345346      ENDIF 
Note: See TracChangeset for help on using the changeset viewer.