- Timestamp:
- 2016-07-01T18:02:45+02:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/OPA_SRC/ZDF/zdftke_crs.F90
r6101 r6772 260 260 DO ji = 2, jpi_crs-1 ! vector opt. 261 261 DO jk = 2, jpkm1 ! from the surface to the bottom : lup 262 zmxld_crs(ji,jj,jk) = MIN( zmxld_crs(ji,jj,jk-1) + e3t_crs(ji,jj,jk-1), zmxlm_crs(ji,jj,jk) )262 zmxld_crs(ji,jj,jk) = MIN( zmxld_crs(ji,jj,jk-1) + fse3t_crs(ji,jj,jk-1), zmxlm_crs(ji,jj,jk) ) 263 263 END DO 264 264 DO jk = jpkm1, 2 , -1 ! from the bottom to the surface : ldown 265 zmxlm_crs(ji,jj,jk) = MIN( zmxlm_crs(ji,jj,jk+1) + e3t_crs(ji,jj,jk+1), zmxlm_crs(ji,jj,jk) )265 zmxlm_crs(ji,jj,jk) = MIN( zmxlm_crs(ji,jj,jk+1) + fse3t_crs(ji,jj,jk+1), zmxlm_crs(ji,jj,jk) ) 266 266 END DO 267 267 END DO … … 312 312 DO ji = 2, jpi_crs-1 ! vector opt. 313 313 DO jk = 2, jpkm1 314 zcoef = avm_crs(ji,jj,jk) * 2._wp * e3w_crs(ji,jj,jk) *e3w_crs(ji,jj,jk)314 zcoef = avm_crs(ji,jj,jk) * 2._wp * fse3w_crs(ji,jj,jk) * fse3w_crs(ji,jj,jk) 315 315 ! ! shear 316 316 zdku = avmu_crs(ji-1,jj,jk) * ( un_crs(ji-1,jj,jk-1) - un_crs(ji-1,jj,jk) ) * ( ub_crs(ji-1,jj,jk-1) - ub_crs(ji-1,jj,jk) ) &
Note: See TracChangeset
for help on using the changeset viewer.