- Timestamp:
- 2016-11-03T06:08:40+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/dev_r6711_SIMPLIF_6_aerobulk/NEMOGCM/NEMO/OPA_SRC/ZDF/zdftke.F90
r7163 r7170 558 558 zmxlm(:,:,1) = rn_mxl0 559 559 ENDIF 560 561 !!gm copy from GLS:562 ! ! Set surface roughness length563 ! SELECT CASE ( nn_z0_met )564 ! !565 ! CASE ( 0 ) ! Constant roughness566 ! zhsro(:,:) = rn_hsro567 ! CASE ( 1 ) ! Standard Charnock formula for surface roughness568 ! zhsro(:,:) = MAX(rsbc_zs1 * ustars2(:,:), rn_hsro)569 ! with: rsbc_zs1 = rn_charn/grav570 ! CASE ( 2 ) ! Roughness formulae according to Rascle et al., Ocean Modelling (2008)571 ! zdep(:,:) = 30.*TANH(2.*0.3/(28.*SQRT(MAX(ustars2(:,:),rsmall)))) ! Wave age (eq. 10)572 ! zhsro(:,:) = MAX(rsbc_zs2 * ustars2(:,:) * zdep(:,:)**1.5, rn_hsro) ! zhsro = rn_frac_hs * Hsw (eq. 11)573 ! !574 ! END SELECT575 !!gm end576 577 578 579 560 ! 580 561 DO jk = 2, jpkm1 ! interior value : l=sqrt(2*e/n^2) … … 807 788 808 789 IF( ln_mxl0 ) THEN 809 IF(lwp) WRITE(numout,*)810 790 IF(lwp) WRITE(numout,*) ' use a surface mixing length = F(stress) : set rn_mxl0 = rmxl_min' 811 791 rn_mxl0 = rmxl_min
Note: See TracChangeset
for help on using the changeset viewer.