- Timestamp:
- 2012-11-21T19:01:50+01:00 (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2012/dev_r3604_LEGI8_TAM/NEMOGCM/NEMO/OPATAM_SRC/DYN/dynadv_ubs_tam.F90
r3611 r3627 253 253 ! 254 254 zfv_ftl(ji ,jj ,jk) = ( zfvitl - gamma2 * ( zlv_vutl(ji,jj,jk,2) + zlv_vutl(ji+1,jj ,jk,2) ) ) & 255 & * ( un(ji,jj,jk) + un(ji ,jj+1,jk) - gamma1 * zl_u ) + ( zfvi - gamma2 * ( zlv_vu(ji,jj,jk,2) + zlv_vu(ji+1,jj ,jk,2) ) ) & 255 & * ( un(ji,jj,jk) + un(ji ,jj+1,jk) - gamma1 * zl_u ) & 256 & + ( zfvi - gamma2 * ( zlv_vu(ji,jj,jk,2) + zlv_vu(ji+1,jj ,jk,2) ) ) & 256 257 & * ( un_tl(ji,jj,jk) + un_tl(ji ,jj+1,jk) - gamma1 * zl_utl ) 257 258 zfu_ftl(ji ,jj ,jk) = ( zfujtl - gamma2 * ( zlu_uvtl(ji,jj,jk,2) + zlu_uvtl(ji ,jj+1,jk,2) ) ) & 258 & * ( vn(ji,jj,jk) + vn(ji+1,jj ,jk) - gamma1 * zl_v ) + ( zfuj - gamma2 * ( zlu_uv(ji,jj,jk,2) + zlu_uv(ji ,jj+1,jk,2) ) ) & 259 & * ( vn(ji,jj,jk) + vn(ji+1,jj ,jk) - gamma1 * zl_v ) & 260 & + ( zfuj - gamma2 * ( zlu_uv(ji,jj,jk,2) + zlu_uv(ji ,jj+1,jk,2) ) ) & 259 261 & * ( vn_tl(ji,jj,jk) + vn_tl(ji+1,jj ,jk) - gamma1 * zl_vtl ) 260 262 END DO … … 289 291 DO jj = 2, jpjm1 290 292 DO ji = fs_2, fs_jpim1 291 zfu_uwtl(ji,jj, 1 ) = 2.e0 * ( zfwtl(ji,jj,1) + zfwtl(ji+1,jj ,1) ) * un(ji,jj,1) + 2.e0 * ( zfw(ji,jj,1) + zfw(ji+1,jj ,1) ) * un_tl(ji,jj,1) 292 zfv_vwtl(ji,jj, 1 ) = 2.e0 * ( zfwtl(ji,jj,1) + zfwtl(ji ,jj+1,1) ) * vn(ji,jj,1) + 2.e0 * ( zfw(ji,jj,1) + zfw(ji ,jj+1,1) ) * vn_tl(ji,jj,1) 293 zfu_uwtl(ji,jj, 1 ) = 2.e0 * ( zfwtl(ji,jj,1) + zfwtl(ji+1,jj ,1) ) * un(ji,jj,1) & 294 & + 2.e0 * ( zfw(ji,jj,1) + zfw(ji+1,jj ,1) ) * un_tl(ji,jj,1) 295 zfv_vwtl(ji,jj, 1 ) = 2.e0 * ( zfwtl(ji,jj,1) + zfwtl(ji ,jj+1,1) ) * vn(ji,jj,1) & 296 & + 2.e0 * ( zfw(ji,jj,1) + zfw(ji ,jj+1,1) ) * vn_tl(ji,jj,1) 293 297 END DO 294 298 END DO … … 297 301 DO jj = 2, jpjm1 298 302 DO ji = fs_2, fs_jpim1 ! vector opt. 299 zfu_uwtl(ji,jj,jk) = ( zfwtl(ji,jj,jk)+ zfwtl(ji+1,jj ,jk) ) * ( un(ji,jj,jk) + un(ji,jj,jk-1) ) + ( zfw(ji,jj,jk)+ zfw(ji+1,jj ,jk) ) * ( un_tl(ji,jj,jk) + un_tl(ji,jj,jk-1) ) 300 zfv_vwtl(ji,jj,jk) = ( zfwtl(ji,jj,jk)+ zfwtl(ji ,jj+1,jk) ) * ( vn(ji,jj,jk) + vn(ji,jj,jk-1) ) + ( zfw(ji,jj,jk)+ zfw(ji ,jj+1,jk) ) * ( vn_tl(ji,jj,jk) + vn_tl(ji,jj,jk-1) ) 303 zfu_uwtl(ji,jj,jk) = ( zfwtl(ji,jj,jk)+ zfwtl(ji+1,jj ,jk) ) * ( un(ji,jj,jk) + un(ji,jj,jk-1) ) & 304 & + ( zfw(ji,jj,jk)+ zfw(ji+1,jj ,jk) ) * ( un_tl(ji,jj,jk) + un_tl(ji,jj,jk-1) ) 305 zfv_vwtl(ji,jj,jk) = ( zfwtl(ji,jj,jk)+ zfwtl(ji ,jj+1,jk) ) * ( vn(ji,jj,jk) + vn(ji,jj,jk-1) ) & 306 & + ( zfw(ji,jj,jk)+ zfw(ji ,jj+1,jk) ) * ( vn_tl(ji,jj,jk) + vn_tl(ji,jj,jk-1) ) 301 307 END DO 302 308 END DO
Note: See TracChangeset
for help on using the changeset viewer.