- Timestamp:
- 2011-09-27T14:33:01+02:00 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2011/UKMO_MERCATOR_obc_bdy_merge/NEMOGCM/NEMO/OPA_SRC/OBC/obcdta.F90
r2831 r2865 55 55 CONTAINS 56 56 57 SUBROUTINE obc_dta( kt, jit )57 SUBROUTINE obc_dta( kt, jit, time_offset ) 58 58 !!---------------------------------------------------------------------- 59 59 !! *** SUBROUTINE obc_dta *** … … 69 69 INTEGER, INTENT( in ) :: kt ! ocean time-step index 70 70 INTEGER, INTENT( in ), OPTIONAL :: jit ! subcycle time-step index (for timesplitting option) 71 INTEGER, INTENT( in ), OPTIONAL :: time_offset ! time offset in units of timesteps. NB. if jit 72 ! is present then units = subcycle timesteps. 73 ! time_offset = 0 => get data at "now" time level 74 ! time_offset = -1 => get data at "before" time level 75 ! time_offset = +1 => get data at "after" time level 76 ! etc. 71 77 !! 72 78 INTEGER :: ib_obc, jfld, jstart, jend, ib, ii, ij, ik, igrd ! local indices … … 202 208 IF( PRESENT(jit) ) THEN 203 209 ! Update barotropic boundary conditions only 204 ! jit is optional argument for fld_read 205 IF( nn_dyn2d(ib_obc) .gt. 0 .and. nn_dyn2d_dta(ib_obc) .eq. 1) THEN206 IF( nn_ tides(ib_obc) .eq. 1 ) THEN210 ! jit is optional argument for fld_read and tide_update 211 IF( nn_dyn2d(ib_obc) .gt. 0 ) THEN 212 IF( nn_dyn2d_dta(ib_obc) .eq. 2 ) THEN ! tidal harmonic forcing ONLY: initialise arrays 207 213 dta_obc(ib_obc)%ssh(:) = 0.0 208 214 dta_obc(ib_obc)%u2d(:) = 0.0 209 215 dta_obc(ib_obc)%v2d(:) = 0.0 210 ELSE 216 ENDIF 217 IF( nn_dyn2d_dta(ib_obc) .eq. 1 .or. nn_dyn2d_dta(ib_obc) .eq. 3 ) THEN ! update external data 211 218 jend = jstart + 2 212 CALL fld_read( kt=kt, kn_fsbc=1, sd=bf(jstart:jend), map=nbmap_ptr(jstart:jend), jit=jit )219 CALL fld_read( kt=kt, kn_fsbc=1, sd=bf(jstart:jend), map=nbmap_ptr(jstart:jend), jit=jit, time_offset=time_offset ) 213 220 ENDIF 214 ENDIF215 IF( nn_tides(ib_obc) .gt. 0 ) THEN216 CALL tide_update( kt=kt, jit=jit, idx=idx_obc(ib_obc), dta=dta_obc(ib_obc), td=tides(ib_obc) )221 IF( nn_dyn2d_dta(ib_obc) .ge. 2 ) THEN ! update tidal harmonic forcing 222 CALL tide_update( kt=kt, idx=idx_obc(ib_obc), dta=dta_obc(ib_obc), td=tides(ib_obc), jit=jit, time_offset=time_offset ) 223 ENDIF 217 224 ENDIF 218 225 ELSE 219 IF( nb_obc_fld(ib_obc) .gt. 0 ) THEN 220 jend = jstart + nb_obc_fld(ib_obc) - 1 221 CALL fld_read( kt=kt, kn_fsbc=1, sd=bf(jstart:jend), map=nbmap_ptr(jstart:jend), timeshift=1 ) 222 ENDIF 223 IF( nn_tides(ib_obc) .eq. 1 ) THEN 226 IF( nn_dyn2d(ib_obc) .gt. 0 .and. nn_dyn2d_dta(ib_obc) .eq. 2 ) THEN ! tidal harmonic forcing ONLY: initialise arrays 224 227 dta_obc(ib_obc)%ssh(:) = 0.0 225 228 dta_obc(ib_obc)%u2d(:) = 0.0 226 229 dta_obc(ib_obc)%v2d(:) = 0.0 227 230 ENDIF 228 IF( nn_tides(ib_obc) .gt. 0 ) THEN 229 !!! THINK ABOUT kt, jit VALUES !!! 230 CALL tide_update( kt=kt, jit=0, idx=idx_obc(ib_obc), dta=dta_obc(ib_obc), td=tides(ib_obc) ) 231 IF( nb_obc_fld(ib_obc) .gt. 0 ) THEN ! update external data 232 jend = jstart + nb_obc_fld(ib_obc) - 1 233 CALL fld_read( kt=kt, kn_fsbc=1, sd=bf(jstart:jend), map=nbmap_ptr(jstart:jend), time_offset=time_offset ) 234 ENDIF 235 IF( nn_dyn2d(ib_obc) .gt. 0 .and. nn_dyn2d_dta(ib_obc) .ge. 2 ) THEN ! update tidal harmonic forcing 236 CALL tide_update( kt=kt, idx=idx_obc(ib_obc), dta=dta_obc(ib_obc), td=tides(ib_obc), time_offset=time_offset ) 231 237 ENDIF 232 238 ENDIF … … 238 244 239 245 IF( ln_full_vel_array(ib_obc) .and. & 240 & ( nn_dyn2d_dta(ib_obc) .eq. 1 .or. nn_dyn 3d_dta(ib_obc) .eq. 1 ) ) THEN246 & ( nn_dyn2d_dta(ib_obc) .eq. 1 .or. nn_dyn2d_dta(ib_obc) .eq. 3 .or. nn_dyn3d_dta(ib_obc) .eq. 1 ) ) THEN 241 247 242 248 igrd = 2 ! zonal velocity … … 326 332 #endif 327 333 ) 334 IF(nn_dta(ib_obc) .gt. 1) nn_dta(ib_obc) = 1 328 335 END DO 329 336 … … 333 340 nb_obc_fld(:) = 0 334 341 DO ib_obc = 1, nb_obc 335 IF( nn_dyn2d(ib_obc) .gt. 0 .and. nn_dyn2d_dta(ib_obc) .eq. 1) THEN342 IF( nn_dyn2d(ib_obc) .gt. 0 .and. ( nn_dyn2d_dta(ib_obc) .eq. 1 .or. nn_dyn2d_dta(ib_obc) .eq. 3 ) ) THEN 336 343 nb_obc_fld(ib_obc) = nb_obc_fld(ib_obc) + 3 337 344 ENDIF … … 408 415 ! Only read in necessary fields for this set. 409 416 ! Important that barotropic variables come first. 410 IF( nn_dyn2d(ib_obc) .gt. 0 .and. nn_dyn2d_dta(ib_obc) .eq. 1) THEN411 412 IF( nn_dyn2d(ib_obc) .ne. jp_frs .and. nn_tides(ib_obc) .ne. 1) THEN417 IF( nn_dyn2d(ib_obc) .gt. 0 .and. ( nn_dyn2d_dta(ib_obc) .eq. 1 .or. nn_dyn2d_dta(ib_obc) .eq. 3 ) ) THEN 418 419 IF( nn_dyn2d(ib_obc) .ne. jp_frs ) THEN 413 420 jfld = jfld + 1 414 421 blf_i(jfld) = bn_ssh … … 419 426 ENDIF 420 427 421 IF( .not. ln_full_vel_array(ib_obc) .and. nn_tides(ib_obc) .ne. 1) THEN428 IF( .not. ln_full_vel_array(ib_obc) ) THEN 422 429 423 430 jfld = jfld + 1 … … 449 456 ! baroclinic velocities 450 457 IF( ( nn_dyn3d(ib_obc) .gt. 0 .and. nn_dyn3d_dta(ib_obc) .eq. 1 ) .or. & 451 ( ln_full_vel_array(ib_obc) .and. nn_dyn2d(ib_obc) .gt. 0 .and. nn_dyn2d_dta(ib_obc) .eq. 1 ) ) THEN 458 & ( ln_full_vel_array(ib_obc) .and. nn_dyn2d(ib_obc) .gt. 0 .and. & 459 & ( nn_dyn2d_dta(ib_obc) .eq. 1 .or. nn_dyn2d_dta(ib_obc) .eq. 3 ) ) ) THEN 452 460 453 461 jfld = jfld + 1 … … 583 591 584 592 IF (nn_dyn2d(ib_obc) .gt. 0) THEN 585 IF( nn_dyn2d_dta(ib_obc) .eq. 0 .or. ln_full_vel_array(ib_obc) .or. nn_tides(ib_obc) .eq. 1) THEN593 IF( nn_dyn2d_dta(ib_obc) .eq. 0 .or. nn_dyn2d_dta(ib_obc) .eq. 2 .or. ln_full_vel_array(ib_obc) ) THEN 586 594 IF( nn_dyn2d(ib_obc) .eq. jp_frs ) THEN 587 595 ilen0(1:3) = nblen(1:3) … … 614 622 ENDIF 615 623 IF ( ( nn_dyn3d(ib_obc) .gt. 0 .and. nn_dyn3d_dta(ib_obc) .eq. 1 ).or. & 616 & ( ln_full_vel_array(ib_obc) .and. nn_dyn2d(ib_obc) .gt. 0 ) ) THEN 624 & ( ln_full_vel_array(ib_obc) .and. nn_dyn2d(ib_obc) .gt. 0 .and. & 625 & ( nn_dyn2d_dta(ib_obc) .eq. 1 .or. nn_dyn2d_dta(ib_obc) .eq. 3 ) ) ) THEN 617 626 jfld = jfld + 1 618 627 dta_obc(ib_obc)%u3d => bf(jfld)%fnow(:,1,:) … … 646 655 ilen0(1:3) = nblenrim(1:3) 647 656 ENDIF 648 ALLOCATE( dta_obc(ib_obc)% ssh(ilen0(1)) )649 ALLOCATE( dta_obc(ib_obc)% u2d(ilen0(1)) )650 ALLOCATE( dta_obc(ib_obc)% v2d(ilen0(1)) )657 ALLOCATE( dta_obc(ib_obc)%frld(ilen0(1)) ) 658 ALLOCATE( dta_obc(ib_obc)%hicif(ilen0(1)) ) 659 ALLOCATE( dta_obc(ib_obc)%hsnif(ilen0(1)) ) 651 660 ELSE 652 661 jfld = jfld + 1
Note: See TracChangeset
for help on using the changeset viewer.
