- Timestamp:
- 2011-07-11T12:53:56+02:00 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2011/UKMO_MERCATOR_obc_bdy_merge/NEMOGCM/NEMO/OPA_SRC/DYN/dynspg_ts.F90
r2724 r2797 25 25 USE domvvl ! variable volume 26 26 USE zdfbfr ! bottom friction 27 USE obcdta ! open boundary condition data28 USE obcfla ! Flather open boundary condition29 27 USE dynvor ! vorticity term 30 28 USE obc_oce ! Lateral open boundary condition 31 USE obc_par ! open boundary condition parameters 32 USE bdy_oce ! unstructured open boundaries 33 USE bdy_par ! unstructured open boundaries 34 USE bdydta ! unstructured open boundaries 35 USE bdydyn ! unstructured open boundaries 36 USE bdytides ! tidal forcing at unstructured open boundaries. 29 USE obcdta ! open boundary condition data 30 USE obcdyn2d ! open boundary conditions on barotropic variables 37 31 USE lib_mpp ! distributed memory computing library 38 32 USE lbclnk ! ocean lateral boundary conditions (or mpp link) … … 352 346 zssh_sum(:,:) = sshn (:,:) 353 347 354 #if defined key_obc355 ! set ssh corrections to 0356 ! ssh corrections are applied to normal velocities (Flather's algorithm) and averaged over the barotropic loop357 IF( lp_obc_east ) sshfoe_b(:,:) = 0.e0358 IF( lp_obc_west ) sshfow_b(:,:) = 0.e0359 IF( lp_obc_south ) sshfos_b(:,:) = 0.e0360 IF( lp_obc_north ) sshfon_b(:,:) = 0.e0361 #endif362 363 348 ! ! ==================== ! 364 349 DO jn = 1, icycle ! sub-time-step loop ! (from NOW to AFTER+1) … … 367 352 IF( jn == 1 ) z2dt_e = rdt / nn_baro 368 353 369 ! !* Update the forcing (OBC , BDYand tides)354 ! !* Update the forcing (OBC and tides) 370 355 ! ! ------------------ 371 IF( lk_obc ) CALL obc_dta_bt ( kt, jn ) 372 IF( lk_bdy ) CALL bdy_dta_fla( kt, jn+1, icycle ) 356 IF( lk_obc ) CALL obc_dta ( kt, jit=jn ) 373 357 374 358 ! !* after ssh_e … … 384 368 ! 385 369 #if defined key_obc 386 ! ! OBC : zhdiv must be zero behind the open boundary 387 !! mpp remark: The zeroing of hdiv can probably be extended to 1->jpi/jpj for the correct row/column 388 IF( lp_obc_east ) zhdiv(nie0p1:nie1p1,nje0 :nje1 ) = 0.e0 ! east 389 IF( lp_obc_west ) zhdiv(niw0 :niw1 ,njw0 :njw1 ) = 0.e0 ! west 390 IF( lp_obc_north ) zhdiv(nin0 :nin1 ,njn0p1:njn1p1) = 0.e0 ! north 391 IF( lp_obc_south ) zhdiv(nis0 :nis1 ,njs0 :njs1 ) = 0.e0 ! south 392 #endif 393 #if defined key_bdy 394 zhdiv(:,:) = zhdiv(:,:) * bdytmask(:,:) ! BDY mask 370 zhdiv(:,:) = zhdiv(:,:) * obctmask(:,:) ! OBC mask 395 371 #endif 396 372 ! … … 489 465 ! !* domain lateral boundary 490 466 ! ! ----------------------- 491 ! ! Flather's boundary condition for the barotropic loop : 492 ! ! - Update sea surface height on each open boundary 493 ! ! - Correct the velocity 494 495 IF( lk_obc ) CALL obc_fla_ts ( ua_e, va_e, sshn_e, ssha_e ) 496 IF( lk_bdy .OR. ln_tides ) CALL bdy_dyn_fla( sshn_e ) 467 ! OBC open boundaries 468 IF( lk_obc .OR. ln_tides ) CALL obc_dyn2d( kt, sshn_e ) 497 469 ! 498 470 CALL lbc_lnk( ua_e , 'U', -1. ) ! local domain boundaries
Note: See TracChangeset
for help on using the changeset viewer.