Changeset 13295 for NEMO/trunk/src/ICE/icethd.F90
- Timestamp:
- 2020-07-10T20:24:21+02:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/ICE/icethd.F90
r13226 r13295 109 109 zu_io(:,:) = u_ice(:,:) - ssu_m(:,:) 110 110 zv_io(:,:) = v_ice(:,:) - ssv_m(:,:) 111 DO_2D _00_00111 DO_2D( 0, 0, 0, 0 ) 112 112 zfric(ji,jj) = rn_cio * ( 0.5_wp * & 113 113 & ( zu_io(ji,jj) * zu_io(ji,jj) + zu_io(ji-1,jj) * zu_io(ji-1,jj) & … … 115 115 END_2D 116 116 ELSE ! if no ice dynamics => transmit directly the atmospheric stress to the ocean 117 DO_2D _00_00117 DO_2D( 0, 0, 0, 0 ) 118 118 zfric(ji,jj) = r1_rho0 * SQRT( 0.5_wp * & 119 119 & ( utau(ji,jj) * utau(ji,jj) + utau(ji-1,jj) * utau(ji-1,jj) & … … 126 126 ! Partial computation of forcing for the thermodynamic sea ice model 127 127 !--------------------------------------------------------------------! 128 DO_2D _11_11128 DO_2D( 1, 1, 1, 1 ) 129 129 rswitch = tmask(ji,jj,1) * MAX( 0._wp , SIGN( 1._wp , at_i(ji,jj) - epsi10 ) ) ! 0 if no ice 130 130 ! … … 196 196 ! select ice covered grid points 197 197 npti = 0 ; nptidx(:) = 0 198 DO_2D _11_11198 DO_2D( 1, 1, 1, 1 ) 199 199 IF ( a_i(ji,jj,jl) > epsi10 ) THEN 200 200 npti = npti + 1
Note: See TracChangeset
for help on using the changeset viewer.