- Timestamp:
- 2020-05-20T16:08:51+02:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/r4.0-HEAD_r12713_clem_dan_fixcpl/src/OCE/SBC/sbccpl.F90
r12948 r12955 1503 1503 INTEGER :: ji, jj ! dummy loop indices 1504 1504 INTEGER :: itx ! index of taux over ice 1505 REAL(wp) :: zztmp1, zztmp2 1505 1506 REAL(wp), DIMENSION(jpi,jpj) :: ztx, zty 1506 1507 !!---------------------------------------------------------------------- … … 1566 1567 p_taui(:,:) = frcv(jpr_itx1)%z3(:,:,1) ! (U,V) ==> (U,V) 1567 1568 p_tauj(:,:) = frcv(jpr_ity1)%z3(:,:,1) 1568 CASE( 'F' )1569 DO jj = 2, jpjm1 ! F ==> (U,V)1570 DO ji = fs_2, fs_jpim1 ! vector opt.1571 p_taui(ji,jj) = 0.5 * ( frcv(jpr_itx1)%z3(ji,jj,1) + frcv(jpr_itx1)%z3(ji ,jj-1,1) )1572 p_tauj(ji,jj) = 0.5 * ( frcv(jpr_ity1)%z3(ji,jj,1) + frcv(jpr_ity1)%z3(ji-1,jj ,1) )1573 END DO1574 END DO1575 1569 CASE( 'T' ) 1576 1570 DO jj = 2, jpjm1 ! T ==> (U,V) 1577 1571 DO ji = fs_2, fs_jpim1 ! vector opt. 1578 p_taui(ji,jj) = 0.5 * ( frcv(jpr_itx1)%z3(ji+1,jj ,1) + frcv(jpr_itx1)%z3(ji,jj,1) ) 1579 p_tauj(ji,jj) = 0.5 * ( frcv(jpr_ity1)%z3(ji ,jj+1,1) + frcv(jpr_ity1)%z3(ji,jj,1) ) 1572 ! take care of the land-sea mask to avoid "pollution" of coastal stress. p[uv]taui used in frazil and rheology 1573 zztmp1 = 0.5_wp * ( 2. - umask(ji,jj,1) ) * MAX( tmask(ji,jj,1),tmask(ji+1,jj ,1) ) 1574 zztmp2 = 0.5_wp * ( 2. - vmask(ji,jj,1) ) * MAX( tmask(ji,jj,1),tmask(ji ,jj+1,1) ) 1575 p_taui(ji,jj) = zztmp1 * ( frcv(jpr_itx1)%z3(ji+1,jj ,1) + frcv(jpr_itx1)%z3(ji,jj,1) ) 1576 p_tauj(ji,jj) = zztmp2 * ( frcv(jpr_ity1)%z3(ji ,jj+1,1) + frcv(jpr_ity1)%z3(ji,jj,1) ) 1580 1577 END DO 1581 1578 END DO 1582 CASE( 'I' ) 1583 DO jj = 2, jpjm1 ! I ==> (U,V) 1584 DO ji = 2, jpim1 ! NO vector opt. 1585 p_taui(ji,jj) = 0.5 * ( frcv(jpr_itx1)%z3(ji+1,jj+1,1) + frcv(jpr_itx1)%z3(ji+1,jj ,1) ) 1586 p_tauj(ji,jj) = 0.5 * ( frcv(jpr_ity1)%z3(ji+1,jj+1,1) + frcv(jpr_ity1)%z3(ji ,jj+1,1) ) 1587 END DO 1588 END DO 1579 CALL lbc_lnk_multi( 'sbccpl', p_taui, 'U', -1., p_tauj, 'V', -1. ) 1589 1580 END SELECT 1590 IF( srcv(jpr_itx1)%clgrid /= 'U' ) THEN1591 CALL lbc_lnk_multi( 'sbccpl', p_taui, 'U', -1., p_tauj, 'V', -1. )1592 ENDIF1593 1581 1594 1582 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.