- Timestamp:
- 2020-11-12T19:13:23+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/tickets_icb_1900/src/OCE/ICB/icbutl.F90
r13374 r13786 95 95 CALL lbc_lnk_icb( 'icbutl', ssu_e, 'U', -1._wp, 1, 1 ) 96 96 CALL lbc_lnk_icb( 'icbutl', ssv_e, 'V', -1._wp, 1, 1 ) 97 CALL lbc_lnk_icb( 'icbutl', ua_e , 'U', -1._wp, 1, 1 )98 CALL lbc_lnk_icb( 'icbutl', va_e , 'V', -1._wp, 1, 1 )97 CALL lbc_lnk_icb( 'icbutl', ua_e , 'U', -1._wp, 1, 1 ) 98 CALL lbc_lnk_icb( 'icbutl', va_e , 'V', -1._wp, 1, 1 ) 99 99 #if defined key_si3 100 100 hi_e(1:jpi, 1:jpj) = hm_i (:,:) … … 274 274 ! since we're looking for four T points containing quadrant we're in of 275 275 ! current T cell 276 kii = MAX(0, INT( pi ))277 kij = MAX(0, INT( pj )) ! T-point276 kii = MAX(0, INT( pi )) 277 kij = MAX(0, INT( pj )) ! T-point 278 278 zwi = pi - REAL(kii,wp) 279 279 zwj = pj - REAL(kij,wp) 280 280 CASE ( 'U' ) 281 281 kii = MAX(0, INT( pi-0.5_wp )) 282 kij = MAX(0, INT( pj )) ! U-point282 kij = MAX(0, INT( pj )) ! U-point 283 283 zwi = pi - 0.5_wp - REAL(kii,wp) 284 284 zwj = pj - REAL(kij,wp) 285 285 CASE ( 'V' ) 286 kii = MAX(0, INT( pi ))286 kii = MAX(0, INT( pi )) 287 287 kij = MAX(0, INT( pj-0.5_wp )) ! V-point 288 288 zwi = pi - REAL(kii,wp) … … 896 896 897 897 SUBROUTINE test_icb_utl_getkb 898 !!---------------------------------------------------------------------- 899 !! *** FUNCTION test_icb_utl_getkb *** 900 !! 901 !! ** Purpose : Test routine icb_utl_getkb, icb_utl_zavg 902 !! ** Methode : Call each subroutine with specific input data 903 !! What should be the output is easy to determined and check 904 !! if NEMO return the correct answer. 905 !! ** Comments : not called, if needed a CALL test_icb_utl_getkb need to be added in icb_step 906 !!---------------------------------------------------------------------- 898 907 INTEGER :: ikb 899 908 REAL(wp) :: zD, zout
Note: See TracChangeset
for help on using the changeset viewer.