Changeset 589 for trunk/NEMO/OPA_SRC
- Timestamp:
- 2007-01-22T15:01:49+01:00 (17 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/TRA/traldf_bilap.F90
r474 r589 90 90 REAL(wp) :: zta, zsa ! temporary scalars 91 91 REAL(wp), DIMENSION(jpi,jpj) :: & 92 zeeu, zeev, zbtr ! 2D workspace arrays 92 zeeu, zeev, zbtr, & ! 2D workspace 93 zlt, zls 93 94 REAL(wp), DIMENSION(jpi,jpj,jpk) :: & 94 zsu, zsv , zlt, zls ! 3D workspace arrays95 zsu, zsv ! 3D workspace 95 96 !!---------------------------------------------------------------------- 96 97 … … 100 101 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~' 101 102 ENDIF 103 102 104 103 105 ! ! =============== … … 160 162 DO jj = 2, jpjm1 161 163 DO ji = fs_2, fs_jpim1 ! vector opt. 162 zlt(ji,jj ,jk) = zbtr(ji,jj) * ( ztu(ji,jj,jk) - ztu(ji-1,jj,jk) + ztv(ji,jj,jk) - ztv(ji,jj-1,jk) )163 zls(ji,jj ,jk) = zbtr(ji,jj) * ( zsu(ji,jj,jk) - zsu(ji-1,jj,jk) + zsv(ji,jj,jk) - zsv(ji,jj-1,jk) )164 zlt(ji,jj) = zbtr(ji,jj) * ( ztu(ji,jj,jk) - ztu(ji-1,jj,jk) + ztv(ji,jj,jk) - ztv(ji,jj-1,jk) ) 165 zls(ji,jj) = zbtr(ji,jj) * ( zsu(ji,jj,jk) - zsu(ji-1,jj,jk) + zsv(ji,jj,jk) - zsv(ji,jj-1,jk) ) 164 166 END DO 165 167 END DO … … 168 170 DO jj = 2, jpjm1 169 171 DO ji = fs_2, fs_jpim1 ! vector opt. 170 zlt(ji,jj,jk) = fsahtt(ji,jj,jk) * zlt(ji,jj,jk) 171 zls(ji,jj,jk) = fsahtt(ji,jj,jk) * zls(ji,jj,jk) 172 END DO 173 END DO 174 ENDDO 175 176 ! Lateral boundary conditions on the laplacian (zlt,zls) (unchanged sgn) 177 CALL lbc_lnk( zlt, 'T', 1. ) ; CALL lbc_lnk( zls, 'T', 1. ) 178 179 DO jk = 1, jpkm1 180 172 zlt(ji,jj) = fsahtt(ji,jj,jk) * zlt(ji,jj) 173 zls(ji,jj) = fsahtt(ji,jj,jk) * zls(ji,jj) 174 END DO 175 END DO 176 177 ! Lateral boundary conditions on the laplacian (zlt,zls) (unchanged sgn) 178 CALL lbc_lnk( zlt, 'T', 1. ) ; CALL lbc_lnk( zls, 'T', 1. ) 179 181 180 ! 2. Bilaplacian 182 181 ! -------------- … … 185 184 DO jj = 1, jpjm1 186 185 DO ji = 1, fs_jpim1 ! vector opt. 187 ztu(ji,jj,jk) = zeeu(ji,jj) * ( zlt(ji+1,jj ,jk) - zlt(ji,jj,jk) )188 zsu(ji,jj,jk) = zeeu(ji,jj) * ( zls(ji+1,jj ,jk) - zls(ji,jj,jk) )189 ztv(ji,jj,jk) = zeev(ji,jj) * ( zlt(ji ,jj+1 ,jk) - zlt(ji,jj,jk) )190 zsv(ji,jj,jk) = zeev(ji,jj) * ( zls(ji ,jj+1 ,jk) - zls(ji,jj,jk) )186 ztu(ji,jj,jk) = zeeu(ji,jj) * ( zlt(ji+1,jj ) - zlt(ji,jj) ) 187 zsu(ji,jj,jk) = zeeu(ji,jj) * ( zls(ji+1,jj ) - zls(ji,jj) ) 188 ztv(ji,jj,jk) = zeev(ji,jj) * ( zlt(ji ,jj+1) - zlt(ji,jj) ) 189 zsv(ji,jj,jk) = zeev(ji,jj) * ( zls(ji ,jj+1) - zls(ji,jj) ) 191 190 END DO 192 191 END DO
Note: See TracChangeset
for help on using the changeset viewer.