Changeset 12256
- Timestamp:
- 2019-12-16T13:52:49+01:00 (4 years ago)
- Location:
- NEMO/branches/2019/dev_r12072_MERGE_OPTION2_2019/src
- Files:
-
- 5 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r12072_MERGE_OPTION2_2019/src/OCE/BDY/bdydta.F90
r12246 r12256 75 75 CONTAINS 76 76 77 SUBROUTINE bdy_dta( kt , pt_offset)77 SUBROUTINE bdy_dta( kt ) 78 78 !!---------------------------------------------------------------------- 79 79 !! *** SUBROUTINE bdy_dta *** … … 85 85 !!---------------------------------------------------------------------- 86 86 INTEGER, INTENT(in) :: kt ! ocean time-step index 87 REAL(wp),INTENT(in), OPTIONAL :: pt_offset ! time offset in units of timesteps88 87 ! 89 88 INTEGER :: jbdy, jfld, jstart, jend, ib, jl ! dummy loop indices … … 210 209 ! read/update all bdy data 211 210 ! ------------------------ 212 CALL fld_read( kt, 1, bf_alias, pt_offset = pt_offset ) 211 ! BDY: use pt_offset=0.5 as applied at the end of the step and fldread is referenced at the middle of the step 212 CALL fld_read( kt, 1, bf_alias, pt_offset = 0.5_wp ) 213 213 214 214 ! 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 342 342 ! 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 ) 344 345 ENDIF 345 346 ENDIF -
NEMO/branches/2019/dev_r12072_MERGE_OPTION2_2019/src/OCE/SBC/fldread.F90
r12251 r12256 163 163 164 164 IF( nn_components == jp_iam_sas ) THEN ; zt_offset = REAL( nn_fsbc, wp ) 165 ELSE 165 ELSE ; zt_offset = 0. 166 166 ENDIF 167 167 IF( PRESENT(pt_offset) ) zt_offset = pt_offset -
NEMO/branches/2019/dev_r12072_MERGE_OPTION2_2019/src/OCE/step.F90
r12246 r12256 114 114 IF( ln_tide ) CALL tide_update( kstp ) ! update tide potential 115 115 IF( ln_apr_dyn ) CALL sbc_apr ( kstp ) ! atmospheric pressure (NB: call before bdy_dta which needs ssh_ib) 116 IF( ln_bdy ) CALL bdy_dta ( kstp , pt_offset = 1. )! update dynamic & tracer data at open boundaries116 IF( ln_bdy ) CALL bdy_dta ( kstp ) ! update dynamic & tracer data at open boundaries 117 117 IF( ln_isf ) CALL isf_stp ( kstp ) ! ice shelf/ocean boundary condition 118 118 CALL sbc ( kstp ) ! Sea Boundary Condition (including sea-ice) -
NEMO/branches/2019/dev_r12072_MERGE_OPTION2_2019/src/SAS/step.F90
r12246 r12256 96 96 ! From SAS: ocean bdy data are wrong (but we do not care) and ice bdy data are OK. 97 97 ! This is not clean and should be changed in the future. 98 IF( ln_bdy ) CALL bdy_dta ( kstp, pt_offset=1. )! update dynamic & tracer data at open boundaries98 IF( ln_bdy ) CALL bdy_dta( kstp ) ! update dynamic & tracer data at open boundaries 99 99 ! ==> 100 100 CALL sbc ( kstp ) ! Sea Boundary Condition (including sea-ice) -
NEMO/branches/2019/dev_r12072_MERGE_OPTION2_2019/src/TOP/trcbc.F90
r12246 r12256 364 364 IF( PRESENT(jit) ) THEN 365 365 ! 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 367 367 IF( nb_trcobc > 0 ) THEN 368 368 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 ) 370 370 ENDIF 371 371 ! … … 384 384 ELSE 385 385 ! 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 387 387 IF( nb_trcobc > 0 ) THEN 388 388 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 ) 390 390 ENDIF 391 391 !
Note: See TracChangeset
for help on using the changeset viewer.