New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 4861 – NEMO

Changeset 4861


Ignore:
Timestamp:
2014-11-18T12:26:48+01:00 (9 years ago)
Author:
rblod
Message:

dev_r4765_CNRS_agrif : missing coeff for agrif_tke

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2014/dev_r4765_CNRS_agrif/NEMOGCM/NEMO/OPA_SRC/ZDF/zdftke.F90

    r4789 r4861  
    852852      CALL tke_rst( nit000, 'READ' )  !* read or initialize all required files 
    853853      ! 
     854      IF(.NOT. Agrif_Root()) THEN 
     855         wei3d(:,:,:) = 1. 
     856         DO jk = 1, jpkm1 
     857            DO jj = 2, jpjm1 
     858               DO ji = 2, jpim1 
     859                  wei3d(ji,jj,jk) =   & 
     860                     &   1.*tmask(ji-1,jj-1,jk) + 2.*tmask(ji,jj-1,jk) + 1.*tmask(ji+1,jj-1,jk)& 
     861                     & + 2.*tmask(ji-1,jj  ,jk) + 4.*tmask(ji,jj  ,jk) + 2.*tmask(ji+1,jj  ,jk)& 
     862                     & + 1.*tmask(ji-1,jj+1,jk) + 2.*tmask(ji,jj+1,jk) + 1.*tmask(ji+1,jj+1,jk) 
     863                  wei3d(ji,jj,jk) = tmask(ji,jj,jk) / MAX( 1., wei3d(ji,jj,jk) ) 
     864               END DO 
     865            END DO 
     866         END DO 
     867         CALL lbc_lnk( wei3d, 'T', 1. ) 
     868 
     869         wmix(:,:) = 0. 
     870         wmix(mi0(2):mi1(jpiglo-1),mj0(2):mj1(jpjglo-1)) = 1. 
     871         wmix(mi0(6):mi1(jpiglo-5),mj0(6):mj1(jpjglo-5)) = 0.75 
     872         wmix(mi0(7):mi1(jpiglo-6),mj0(7):mj1(jpjglo-6)) = 0.5 
     873         wmix(mi0(8):mi1(jpiglo-7),mj0(8):mj1(jpjglo-7)) = 0.25 
     874         wmix(mi0(9):mi1(jpiglo-8),mj0(9):mj1(jpjglo-8)) = 0. 
     875      END IF 
     876      !   
    854877   END SUBROUTINE zdf_tke_init 
    855878 
Note: See TracChangeset for help on using the changeset viewer.