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

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

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r12072_MERGE_OPTION2_2019/src/TOP/trcbc.F90

    r12246 r12256  
    364364      IF( PRESENT(jit) ) THEN  
    365365         ! 
    366          ! OPEN boundary conditions (use pt_offset=1. as they are applied at the end of the step) 
     366         ! BDY: use pt_offset=0.5 as applied at the end of the step and fldread is referenced at the middle of the step 
    367367         IF( nb_trcobc > 0 ) THEN 
    368368           if (lwp) write(numout,'(a,i5,a,i10)') '   reading OBC data for ', nb_trcobc ,' variable(s) at step ', kt 
    369            CALL fld_read( kt=kt, kn_fsbc=1, sd=sf_trcobc, kit=jit, pt_offset=1.) 
     369           CALL fld_read( kt=kt, kn_fsbc=1, sd=sf_trcobc, kit=jit, pt_offset = 0.5_wp ) 
    370370         ENDIF 
    371371         ! 
     
    384384      ELSE 
    385385         ! 
    386          ! OPEN boundary conditions (use pt_offset=1. as they are applied at the end of the step) 
     386         ! BDY: use pt_offset=0.5 as applied at the end of the step and fldread is referenced at the middle of the step 
    387387         IF( nb_trcobc > 0 ) THEN 
    388388           if (lwp) write(numout,'(a,i5,a,i10)') '   reading OBC data for ', nb_trcobc ,' variable(s) at step ', kt 
    389            CALL fld_read( kt=kt, kn_fsbc=1, sd=sf_trcobc, pt_offset=1.) 
     389           CALL fld_read( kt=kt, kn_fsbc=1, sd=sf_trcobc, pt_offset = 0.5_wp ) 
    390390         ENDIF 
    391391         ! 
Note: See TracChangeset for help on using the changeset viewer.