- Timestamp:
- 2010-10-13T10:48:48+02:00 (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/DEV_r2106_LOCEAN2010/NEMO/OPA_SRC/DYN/dynldf_bilap.F90
r1152 r2240 162 162 zuf(ji,jj,jk) = fmask(ji,jj,jk) * ( zcv(ji+1,jj ) - zcv(ji,jj) & 163 163 & - zcu(ji ,jj+1) + zcu(ji,jj) ) & 164 #if defined key_zco165 & / ( e1f(ji,jj)*e2f(ji,jj) )166 #else167 164 & * fse3f(ji,jj,jk) / ( e1f(ji,jj)*e2f(ji,jj) ) 168 #endif169 165 END DO 170 166 END DO … … 173 169 DO jj = 1, jpjm1 174 170 DO ji = 1, fs_jpim1 ! vector opt. 175 #if defined key_zco176 zlu(ji,jj,jk) = e2u(ji,jj) * zlu(ji,jj,jk)177 zlv(ji,jj,jk) = e1v(ji,jj) * zlv(ji,jj,jk)178 #else179 171 zlu(ji,jj,jk) = e2u(ji,jj) * fse3u(ji,jj,jk) * zlu(ji,jj,jk) 180 172 zlv(ji,jj,jk) = e1v(ji,jj) * fse3v(ji,jj,jk) * zlv(ji,jj,jk) 181 #endif182 173 END DO 183 174 END DO … … 186 177 DO jj = 2, jpj 187 178 DO ji = fs_2, jpi ! vector opt. 188 #if defined key_zco189 zbt = e1t(ji,jj) * e2t(ji,jj)190 #else191 179 zbt = e1t(ji,jj) * e2t(ji,jj) * fse3t(ji,jj,jk) 192 #endif193 180 zut(ji,jj,jk) = ( zlu(ji,jj,jk) - zlu(ji-1,jj ,jk) & 194 181 & + zlv(ji,jj,jk) - zlv(ji ,jj-1,jk) ) / zbt … … 210 197 DO jj = 2, jpjm1 211 198 DO ji = fs_2, fs_jpim1 ! vector opt. 212 #if defined key_zco213 ze2u = e2u(ji,jj)214 ze2v = e1v(ji,jj)215 #else216 199 ze2u = e2u(ji,jj) * fse3u(ji,jj,jk) 217 200 ze2v = e1v(ji,jj) * fse3v(ji,jj,jk) 218 #endif219 201 ! horizontal biharmonic diffusive trends 220 202 zua = - ( zuf(ji ,jj,jk) - zuf(ji,jj-1,jk) ) / ze2u &
Note: See TracChangeset
for help on using the changeset viewer.