- Timestamp:
- 2014-04-06T17:28:25+02:00 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_CNRS0_NOC1_LDF/NEMOGCM/NEMO/OPA_SRC/ZDF/zdftke.F90
r4147 r4616 99 99 # include "vectopt_loop_substitute.h90" 100 100 !!---------------------------------------------------------------------- 101 !! NEMO/OPA 4.0 , NEMO Consortium (2011)101 !! NEMO/OPA 3.7 , NEMO Consortium (2014) 102 102 !! $Id$ 103 103 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) … … 256 256 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 257 257 ! en(bot) = (rn_ebb0/rau0)*0.5*sqrt(u_botfr^2+v_botfr^2) (min value rn_emin) 258 ! CDIR NOVERRCHK258 !!bfr - commented area 259 259 !! DO jj = 2, jpjm1 260 !CDIR NOVERRCHK261 260 !! DO ji = fs_2, fs_jpim1 ! vector opt. 262 261 !! ztx2 = bfrua(ji-1,jj) * ub(ji-1,jj,mbku(ji-1,jj)) + & … … 291 290 END DO 292 291 ! ! finite LC depth 293 # if defined key_vectopt_loop294 DO jj = 1, 1295 DO ji = 1, jpij ! vector opt. (forced unrolling)296 # else297 292 DO jj = 1, jpj 298 293 DO ji = 1, jpi 299 # endif300 294 zhlc(ji,jj) = fsdepw(ji,jj,imlc(ji,jj)) 301 295 END DO 302 296 END DO 303 297 zcof = 0.016 / SQRT( zrhoa * zcdrag ) 304 !CDIR NOVERRCHK305 298 DO jk = 2, jpkm1 !* TKE Langmuir circulation source term added to en 306 !CDIR NOVERRCHK 307 DO jj = 2, jpjm1 308 !CDIR NOVERRCHK 299 DO jj = 2, jpjm1 309 300 DO ji = fs_2, fs_jpim1 ! vector opt. 310 301 zus = zcof * SQRT( taum(ji,jj) ) ! Stokes drift … … 425 416 END DO 426 417 ELSEIF( nn_etau == 3 ) THEN !* penetration belox the mixed layer (HF variability) 427 !CDIR NOVERRCHK428 418 DO jk = 2, jpkm1 429 !CDIR NOVERRCHK 430 DO jj = 2, jpjm1 431 !CDIR NOVERRCHK 419 DO jj = 2, jpjm1 432 420 DO ji = fs_2, fs_jpim1 ! vector opt. 433 421 ztx2 = utau(ji-1,jj ) + utau(ji,jj) … … 513 501 zmxlm(:,:,jpk) = rmxl_min ! last level set to the interior minium value 514 502 ! 515 !CDIR NOVERRCHK516 503 DO jk = 2, jpkm1 ! interior value : l=sqrt(2*e/n^2) 517 !CDIR NOVERRCHK518 504 DO jj = 2, jpjm1 519 !CDIR NOVERRCHK520 505 DO ji = fs_2, fs_jpim1 ! vector opt. 521 506 zrn2 = MAX( rn2(ji,jj,jk), rsmall ) … … 588 573 END DO 589 574 END DO 590 !CDIR NOVERRCHK591 575 DO jk = 2, jpkm1 592 !CDIR NOVERRCHK 593 DO jj = 2, jpjm1 594 !CDIR NOVERRCHK 576 DO jj = 2, jpjm1 595 577 DO ji = fs_2, fs_jpim1 ! vector opt. 596 578 zemlm = MIN ( zmxld(ji,jj,jk), zmxlm(ji,jj,jk) ) … … 612 594 ! ! Vertical eddy viscosity and diffusivity (avmu, avmv, avt) 613 595 ! !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 614 !CDIR NOVERRCHK615 596 DO jk = 1, jpkm1 !* vertical eddy viscosity & diffivity at w-points 616 !CDIR NOVERRCHK617 597 DO jj = 2, jpjm1 618 !CDIR NOVERRCHK619 598 DO ji = fs_2, fs_jpim1 ! vector opt. 620 599 zsqen = SQRT( en(ji,jj,jk) )
Note: See TracChangeset
for help on using the changeset viewer.