Changeset 888 for trunk/NEMO/LIM_SRC_3/limrhg.F90
- Timestamp:
- 2008-04-11T19:05:03+02:00 (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/LIM_SRC_3/limrhg.F90
r869 r888 16 16 USE dom_oce 17 17 USE dom_ice 18 USE sbc_ice ! Surface boundary condition: ice fields 18 19 USE ice 19 20 USE iceini … … 40 41 !!---------------------------------------------------------------------- 41 42 !! LIM 3.0, UCL-LOCEAN-IPSL (2008) 42 !! $ Header: /home/opalod/NEMOCVSROOT/NEMO/LIM_SRC/limrhg.F90,v 1.5 2005/03/27 18:34:42 opalod Exp $43 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt43 !! $ Id: $ 44 !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) 44 45 !!---------------------------------------------------------------------- 45 46 … … 268 269 / ( e2t(ji,jj+1) + e2t(ji,jj) + epsd ) 269 270 ! 270 u_oce1(ji,jj) = u_ io(ji,jj)271 v_oce2(ji,jj) = v_ io(ji,jj)271 u_oce1(ji,jj) = u_oce(ji,jj) 272 v_oce2(ji,jj) = v_oce(ji,jj) 272 273 273 274 ! Ocean has no slip boundary condition 274 v_oce1(ji,jj) = 0.5*( (v_ io(ji,jj)+v_io(ji,jj-1))*e1t(ji,jj) &275 & +(v_ io(ji+1,jj)+v_io(ji+1,jj-1))*e1t(ji+1,jj)) &275 v_oce1(ji,jj) = 0.5*( (v_oce(ji,jj)+v_oce(ji,jj-1))*e1t(ji,jj) & 276 & +(v_oce(ji+1,jj)+v_oce(ji+1,jj-1))*e1t(ji+1,jj)) & 276 277 & /(e1t(ji+1,jj)+e1t(ji,jj)) * tmu(ji,jj) 277 278 278 u_oce2(ji,jj) = 0.5*((u_ io(ji,jj)+u_io(ji-1,jj))*e2t(ji,jj) &279 & +(u_ io(ji,jj+1)+u_io(ji-1,jj+1))*e2t(ji,jj+1)) &279 u_oce2(ji,jj) = 0.5*((u_oce(ji,jj)+u_oce(ji-1,jj))*e2t(ji,jj) & 280 & +(u_oce(ji,jj+1)+u_oce(ji-1,jj+1))*e2t(ji,jj+1)) & 280 281 & / (e2t(ji,jj+1)+e2t(ji,jj)) * tmv(ji,jj) 281 282 282 283 ! Wind stress. 283 ztagnx = ( 1. - zfrld1(ji,jj) ) * gtaux(ji,jj)284 ztagny = ( 1. - zfrld2(ji,jj) ) * gtauy(ji,jj)284 ztagnx = ( 1. - zfrld1(ji,jj) ) * utaui_ice(ji,jj) 285 ztagny = ( 1. - zfrld2(ji,jj) ) * vtaui_ice(ji,jj) 285 286 286 287 ! Computation of the velocity field taking into account the ice internal interaction. … … 621 622 zdummy = zindb * vt_i(ji,jj) / MAX(at_i(ji,jj) , 1.0e-06 ) 622 623 IF ( zdummy .LE. 5.0e-2 ) THEN 623 u_ice(ji,jj) = u_ io(ji,jj)624 v_ice(ji,jj) = v_ io(ji,jj)624 u_ice(ji,jj) = u_oce(ji,jj) 625 v_ice(ji,jj) = v_oce(ji,jj) 625 626 ENDIF ! zdummy 626 627 END DO
Note: See TracChangeset
for help on using the changeset viewer.