Changeset 11424 for branches/UKMO/dev_r5518_GO6_package_FOAMv14_STOPACK/NEMOGCM/NEMO/OPA_SRC/ZDF/zdftke.F90
- Timestamp:
- 2019-08-08T18:08:53+02:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_GO6_package_FOAMv14_STOPACK/NEMOGCM/NEMO/OPA_SRC/ZDF/zdftke.F90
r11417 r11424 358 358 zwlc = zind * rn_lc0(ji,jj) * zus * SIN( rpi * fsdepw(ji,jj,jk) / zhlc(ji,jj) ) 359 359 ! ! TKE Langmuir circulation source term 360 en(ji,jj,jk) = en(ji,jj,jk) + rdt * MAX(0.,1._wp - fr_i(ji,jj) )* ( zwlc * zwlc * zwlc ) / &360 en(ji,jj,jk) = en(ji,jj,jk) + rdt * ( 1._wp - fr_i(ji,jj) )* ( zwlc * zwlc * zwlc ) / & 361 361 & zhlc(ji,jj) * wmask(ji,jj,jk) * tmask(ji,jj,1) 362 362 … … 487 487 DO ji = fs_2, fs_jpim1 ! vector opt. 488 488 en(ji,jj,jk) = en(ji,jj,jk) + rn_efr0(ji,jj) * en(ji,jj,1) * EXP( -fsdepw(ji,jj,jk) / htau(ji,jj) ) & 489 & * MAX(0.,1._wp - fr_i(ji,jj) ) * wmask(ji,jj,jk) * tmask(ji,jj,1)489 & * ( 1._wp - fr_i(ji,jj) ) * wmask(ji,jj,jk) * tmask(ji,jj,1) 490 490 END DO 491 491 END DO … … 496 496 jk = nmln(ji,jj) 497 497 en(ji,jj,jk) = en(ji,jj,jk) + rn_efr0(ji,jj) * en(ji,jj,1) * EXP( -fsdepw(ji,jj,jk) / htau(ji,jj) ) & 498 & * MAX(0.,1._wp - fr_i(ji,jj) ) * wmask(ji,jj,jk) * tmask(ji,jj,1)498 & * ( 1._wp - fr_i(ji,jj) ) * wmask(ji,jj,jk) * tmask(ji,jj,1) 499 499 END DO 500 500 END DO … … 512 512 zdif = rhftau_scl * MAX( 0._wp, zdif + rhftau_add ) ! apply some modifications... 513 513 en(ji,jj,jk) = en(ji,jj,jk) + zbbrau(ji,jj) * zdif * EXP( -fsdepw(ji,jj,jk) / htau(ji,jj) ) & 514 & * MAX(0.,1._wp - fr_i(ji,jj) ) * wmask(ji,jj,jk) * tmask(ji,jj,1)514 & * ( 1._wp - fr_i(ji,jj) ) * wmask(ji,jj,jk) * tmask(ji,jj,1) 515 515 END DO 516 516 END DO … … 885 885 ALLOCATE( rn_efr0 (jpi,jpj) ) ; rn_efr0 = rn_efr 886 886 887 IF( nn_etau == 2 ) THEN 888 ierr = zdf_mxl_alloc() 889 nmln(:,:) = nlb10 ! Initialization of nmln 890 ENDIF 887 891 888 892 IF( nn_etau /= 0 .and. nn_htau == 2 ) THEN
Note: See TracChangeset
for help on using the changeset viewer.