- Timestamp:
- 2016-05-22T20:11:06+02:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/dev_r6409_SIMPLIF_2_usrdef/NEMOGCM/NEMO/OPA_SRC/SBC/sbcisf.F90
r6140 r6596 536 536 END SUBROUTINE sbc_isf_cav 537 537 538 538 539 SUBROUTINE sbc_isf_gammats(pgt, pgs, pqhisf, pqwisf ) 539 540 !!---------------------------------------------------------------------- … … 635 636 636 637 !! compute eta* (stability parameter) 637 zetastar = 1._wp / ( SQRT(1._wp + MAX(zxsiN * zustar(ji,jj) / ( ABS(ff (ji,jj)) * zmols * zRc ), 0.0_wp)))638 zetastar = 1._wp / ( SQRT(1._wp + MAX(zxsiN * zustar(ji,jj) / ( ABS(ff_f(ji,jj)) * zmols * zRc ), 0._wp))) 638 639 639 640 !! compute the sublayer thickness … … 641 642 642 643 !! compute gamma turb 643 zgturb = 1._wp / vkarmn * LOG(zustar(ji,jj) * zxsiN * zetastar * zetastar / ( ABS(ff (ji,jj)) * zhnu )) &644 zgturb = 1._wp / vkarmn * LOG(zustar(ji,jj) * zxsiN * zetastar * zetastar / ( ABS(ff_f(ji,jj)) * zhnu )) & 644 645 & + 1._wp / ( 2 * zxsiN * zetastar ) - 1._wp / vkarmn 645 646
Note: See TracChangeset
for help on using the changeset viewer.