- Timestamp:
- 2016-11-04T15:23:32+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/r5936_INGV1_WAVE-coupling/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfqiao.F90
r7167 r7193 31 31 REAL(wp),PUBLIC,ALLOCATABLE,DIMENSION (:,:,:) :: QBv, QBvu, QBvv 32 32 33 !! * Substitutions 34 # include "domzgr_substitute.h90" 33 35 !!---------------------------------------------------------------------- 34 36 !! NEMO/OPA 4.0 , NEMO Consortium (2011) … … 82 84 QBv(ji,jj,jk) = 1.0 * 0.353553 * swh(ji,jj) * tsd2d(ji,jj) * & 83 85 & exp(3.0 * wnum(ji,jj) * & 84 & (-MIN( gdept_0(ji ,jj ,jk) , gdept_0(ji+1,jj ,jk), &85 & gdept_0(ji ,jj+1,jk),gdept_0(ji+1,jj+1,jk))))86 & (-MIN( fsdept(ji ,jj ,jk) , fsdept(ji+1,jj ,jk), & 87 & fsdept(ji ,jj+1,jk) , fsdept(ji+1,jj+1,jk)))) 86 88 END DO 87 89 END DO … … 98 100 DO jj = 1, jpjm1 99 101 DO ji = 1, jpim1 100 QBvu(ji,jj,jk) =0.5 * ( 2. - umask(ji,jj,jk) ) *&101 & ( QBv(ji ,jj,jk) * tmask(ji,jj,jk) &102 & + QBv(ji+1,jj,jk) * tmask(ji+1,jj,jk) ) 103 QBvv(ji,jj,jk) =0.5 * ( 2. - vmask(ji,jj,jk) ) *&104 & ( QBv(ji,jj ,jk) * tmask(ji,jj,jk)&102 QBvu(ji,jj,jk) = 0.5 * umask(ji,jj,jk) * & 103 & ( QBv(ji ,jj,jk) * tmask(ji ,jj,jk) & 104 & + QBv(ji+1,jj,jk) * tmask(ji+1,jj,jk) ) 105 QBvv(ji,jj,jk) = 0.5 * vmask(ji,jj,jk) * & 106 & ( QBv(ji,jj ,jk) * tmask(ji,jj ,jk) & 105 107 & + QBv(ji,jj+1,jk) * tmask(ji,jj+1,jk) ) 106 107 108 END DO 108 109 END DO
Note: See TracChangeset
for help on using the changeset viewer.