Changeset 7170


Ignore:
Timestamp:
2016-11-03T06:08:40+01:00 (4 years ago)
Author:
gm
Message:

#1751 - branch SIMPLIF_6_aerobulk: traadv_fct & zdftke =⇒ back to previous revision

Location:
branches/2016/dev_r6711_SIMPLIF_6_aerobulk/NEMOGCM/NEMO/OPA_SRC
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • branches/2016/dev_r6711_SIMPLIF_6_aerobulk/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_fct.F90

    r7163 r7170  
    728728               ! 
    729729               IF( tmask(ji,jj,jk+1) == 0._wp) THEN   ! Switch to second order centered at bottom 
    730                   zwd (ji,jj,jk) = 2._wp 
    731                   zwi (ji,jj,jk) = 0._wp     ! 1  ici ou à l'autre 
    732                   zws (ji,jj,jk) = 0._wp     !    car 1 dans la diag inferieur au fond et superior en surf 
    733                   zwrm(ji,jj,jk) = 3._wp * pt_in(ji,jj,jk)   
     730                  zwd (ji,jj,jk) = 1._wp 
     731                  zwi (ji,jj,jk) = 0._wp 
     732                  zws (ji,jj,jk) = 0._wp 
     733                  zwrm(ji,jj,jk) = 0.5 * ( pt_in(ji,jj,jk-1) + pt_in(ji,jj,jk) )     
    734734               ENDIF 
    735735            END DO 
  • branches/2016/dev_r6711_SIMPLIF_6_aerobulk/NEMOGCM/NEMO/OPA_SRC/ZDF/zdftke.F90

    r7163 r7170  
    558558         zmxlm(:,:,1) = rn_mxl0 
    559559      ENDIF 
    560        
    561 !!gm  copy from GLS: 
    562 !      ! Set surface roughness length 
    563 !      SELECT CASE ( nn_z0_met ) 
    564 !      ! 
    565 !      CASE ( 0 )             ! Constant roughness           
    566 !         zhsro(:,:) = rn_hsro 
    567 !      CASE ( 1 )             ! Standard Charnock formula for surface roughness 
    568 !         zhsro(:,:) = MAX(rsbc_zs1 * ustars2(:,:), rn_hsro) 
    569 !  with:           rsbc_zs1  = rn_charn/grav 
    570 !      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 SELECT 
    575 !!gm end 
    576  
    577  
    578  
    579560      ! 
    580561      DO jk = 2, jpkm1              ! interior value : l=sqrt(2*e/n^2) 
     
    807788 
    808789      IF( ln_mxl0 ) THEN 
    809          IF(lwp) WRITE(numout,*) 
    810790         IF(lwp) WRITE(numout,*) '   use a surface mixing length = F(stress) :   set rn_mxl0 = rmxl_min' 
    811791         rn_mxl0 = rmxl_min 
Note: See TracChangeset for help on using the changeset viewer.