- Timestamp:
- 2012-11-22T18:32:14+01:00 (11 years ago)
- Location:
- branches/2012/dev_NOC_2012_rev3555/NEMOGCM/NEMO/OPA_SRC/DYN
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2012/dev_NOC_2012_rev3555/NEMOGCM/NEMO/OPA_SRC/DYN/dynldf_bilap.F90
r3294 r3634 144 144 145 145 ! Multiply by the eddy viscosity coef. (at u- and v-points) 146 zlu(:,:,jk) = zlu(:,:,jk) * fsahmu(:,:,jk) 147 zlv(:,:,jk) = zlv(:,:,jk) * fsahmv(:,:,jk) 146 zlu(:,:,jk) = zlu(:,:,jk) * ( fsahmu(:,:,jk) * (1-nkahm_smag) + nkahm_smag) 147 148 zlv(:,:,jk) = zlv(:,:,jk) * ( fsahmv(:,:,jk) * (1-nkahm_smag) + nkahm_smag) 148 149 149 150 ! Contravariant "laplacian" … … 200 201 & + ( zut(ji,jj+1,jk) - zut(ji ,jj,jk) ) / e2v(ji,jj) 201 202 ! add it to the general momentum trends 202 ua(ji,jj,jk) = ua(ji,jj,jk) + zua 203 va(ji,jj,jk) = va(ji,jj,jk) + zva 203 ua(ji,jj,jk) = ua(ji,jj,jk) + zua * ( fsahmu(ji,jj,jk)*nkahm_smag +(1 -nkahm_smag )) 204 va(ji,jj,jk) = va(ji,jj,jk) + zva * ( fsahmv(ji,jj,jk)*nkahm_smag +(1 -nkahm_smag )) 204 205 END DO 205 206 END DO -
branches/2012/dev_NOC_2012_rev3555/NEMOGCM/NEMO/OPA_SRC/DYN/dynldf_bilapg.F90
r3294 r3634 414 414 ! II.3 Divergence of vertical fluxes added to the horizontal divergence 415 415 ! --------------------------------------------------------------------- 416 417 IF( kahm == 1 ) THEN 416 IF( (kahm -nkahm_smag) ==1 ) THEN 418 417 ! multiply the laplacian by the eddy viscosity coefficient 419 418 DO jk = 1, jpkm1 … … 430 429 END DO 431 430 END DO 432 ELSEIF( kahm== 2 ) THEN431 ELSEIF( (kahm +nkahm_smag ) == 2 ) THEN 433 432 ! second call, no multiplication 434 433 DO jk = 1, jpkm1
Note: See TracChangeset
for help on using the changeset viewer.