- Timestamp:
- 2021-06-18T15:21:42+02:00 (3 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2021/ticket2680_C1D_PAPA/src/OCE/DOM/domzgr.F90
r14433 r15020 340 340 ! ! N.B. top k-index of W-level = mikt 341 341 ! ! bottom k-index of W-level = mbkt+1 342 DO_2D( 1, 0, 1, 0 )342 DO_2D( 0, 0, 0, 0 ) 343 343 miku(ji,jj) = MAX( mikt(ji+1,jj ) , mikt(ji,jj) ) 344 344 mikv(ji,jj) = MAX( mikt(ji ,jj+1) , mikt(ji,jj) ) … … 349 349 END_2D 350 350 ! converte into REAL to use lbc_lnk ; impose a min value of 1 as a zero can be set in lbclnk 351 zk(:,:) = REAL( miku(:,:), wp ) ; CALL lbc_lnk( 'domzgr', zk, 'U', 1.0_wp ) ; miku(:,:) = MAX( NINT( zk(:,:) ), 1 ) 352 zk(:,:) = REAL( mikv(:,:), wp ) ; CALL lbc_lnk( 'domzgr', zk, 'V', 1.0_wp ) ; mikv(:,:) = MAX( NINT( zk(:,:) ), 1 ) 353 zk(:,:) = REAL( mikf(:,:), wp ) ; CALL lbc_lnk( 'domzgr', zk, 'F', 1.0_wp ) ; mikf(:,:) = MAX( NINT( zk(:,:) ), 1 ) 354 ! 355 zk(:,:) = REAL( mbku(:,:), wp ) ; CALL lbc_lnk( 'domzgr', zk, 'U', 1.0_wp ) ; mbku(:,:) = MAX( NINT( zk(:,:) ), 1 ) 356 zk(:,:) = REAL( mbkv(:,:), wp ) ; CALL lbc_lnk( 'domzgr', zk, 'V', 1.0_wp ) ; mbkv(:,:) = MAX( NINT( zk(:,:) ), 1 ) 351 DO_2D( 0, 0, 0, 0 ) 352 zk(ji,jj) = REAL( miku(ji,jj), wp ) 353 END_2D 354 CALL lbc_lnk( 'domzgr', zk, 'U', 1.0_wp ) 355 miku(:,:) = MAX( NINT( zk(:,:) ), 1 ) 356 357 DO_2D( 0, 0, 0, 0 ) 358 zk(ji,jj) = REAL( mikv(ji,jj), wp ) 359 END_2D 360 CALL lbc_lnk( 'domzgr', zk, 'V', 1.0_wp ) 361 mikv(:,:) = MAX( NINT( zk(:,:) ), 1 ) 362 363 DO_2D( 0, 0, 0, 0 ) 364 zk(ji,jj) = REAL( mikf(ji,jj), wp ) 365 END_2D 366 CALL lbc_lnk( 'domzgr', zk, 'F', 1.0_wp ) 367 mikf(:,:) = MAX( NINT( zk(:,:) ), 1 ) 368 ! 369 DO_2D( 0, 0, 0, 0 ) 370 zk(ji,jj) = REAL( mbku(ji,jj), wp ) 371 END_2D 372 CALL lbc_lnk( 'domzgr', zk, 'U', 1.0_wp ) 373 mbku(:,:) = MAX( NINT( zk(:,:) ), 1 ) 374 375 DO_2D( 0, 0, 0, 0 ) 376 zk(ji,jj) = REAL( mbkv(ji,jj), wp ) 377 END_2D 378 CALL lbc_lnk( 'domzgr', zk, 'V', 1.0_wp ) 379 mbkv(:,:) = MAX( NINT( zk(:,:) ), 1 ) 357 380 ! 358 381 END SUBROUTINE zgr_top_bot
Note: See TracChangeset
for help on using the changeset viewer.