- Timestamp:
- 2016-05-22T20:11:06+02:00 (8 years ago)
- Location:
- branches/2016/dev_r6409_SIMPLIF_2_usrdef/NEMOGCM/NEMO/OPA_SRC/SBC
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/dev_r6409_SIMPLIF_2_usrdef/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_if.F90
r6140 r6596 134 134 ! # qns = zqrp -2(-4) watt/m2 if climatological ice and opa ice (zfr_obs=1, fr_i=1) 135 135 ! (-2=arctic, -4=antarctic) 136 zqi = -3. + SIGN( 1. e0, ff(ji,jj) )136 zqi = -3. + SIGN( 1._wp, ff_f(ji,jj) ) 137 137 qns(ji,jj) = ( ( 1.- zfr_obs ) * qns(ji,jj) & 138 138 & + zfr_obs * fr_i(ji,jj) * zqi ) * tmask(ji,jj,1) & -
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.