Changeset 5179 for trunk/NEMOGCM/NEMO/LIM_SRC_3/limvar.F90
- Timestamp:
- 2015-03-26T14:04:33+01:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/LIM_SRC_3/limvar.F90
r5167 r5179 238 238 DO jj = 1, jpj 239 239 DO ji = 1, jpi 240 rswitch = MAX( 0._wp , SIGN( 1._wp , vt_i(ji,jj) - epsi 20 ) )241 tm_i(ji,jj) = tm_i(ji,jj) + ( 1._wp - rswitch ) * rt0 +&242 & rswitch * t_i(ji,jj,jk,jl) * v_i(ji,jj,jl) &243 & / ( REAL(nlay_i,wp) * MAX( vt_i(ji,jj) , epsi20 ) )244 245 246 247 END DO240 rswitch = MAX( 0._wp , SIGN( 1._wp , vt_i(ji,jj) - epsi10 ) ) 241 tm_i(ji,jj) = tm_i(ji,jj) + r1_nlay_i * rswitch * ( t_i(ji,jj,jk,jl) - rt0 ) * v_i(ji,jj,jl) & 242 & / MAX( vt_i(ji,jj) , epsi10 ) 243 END DO 244 END DO 245 END DO 246 END DO 247 tm_i = tm_i + rt0 248 248 ! 249 249 END SUBROUTINE lim_var_glo2eqv … … 396 396 DO jj = 1, jpj 397 397 DO ji = 1, jpi 398 rswitch = MAX( 0._wp , SIGN( 1._wp , vt_i(ji,jj) - epsi06 ) ) 399 tm_i(ji,jj) = tm_i(ji,jj) + rswitch * t_i(ji,jj,jk,jl) * v_i(ji,jj,jl) & 400 & * r1_nlay_i / MAX( vt_i(ji,jj) , epsi06 ) 401 END DO 402 END DO 403 END DO 404 END DO 398 rswitch = MAX( 0._wp , SIGN( 1._wp , vt_i(ji,jj) - epsi10 ) ) 399 tm_i(ji,jj) = tm_i(ji,jj) + r1_nlay_i * rswitch * ( t_i(ji,jj,jk,jl) - rt0 ) * v_i(ji,jj,jl) & 400 & / MAX( vt_i(ji,jj) , epsi10 ) 401 END DO 402 END DO 403 END DO 404 END DO 405 tm_i = tm_i + rt0 405 406 406 407 END SUBROUTINE lim_var_icetm
Note: See TracChangeset
for help on using the changeset viewer.