Changeset 4861 for branches/2014/dev_r4765_CNRS_agrif
- Timestamp:
- 2014-11-18T12:26:48+01:00 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_r4765_CNRS_agrif/NEMOGCM/NEMO/OPA_SRC/ZDF/zdftke.F90
r4789 r4861 852 852 CALL tke_rst( nit000, 'READ' ) !* read or initialize all required files 853 853 ! 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 ! 854 877 END SUBROUTINE zdf_tke_init 855 878
Note: See TracChangeset
for help on using the changeset viewer.