- Timestamp:
- 2017-12-21T12:03:00+01:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/ROMS_WAD_7832/NEMOGCM/NEMO/OPA_SRC/BDY/bdydyn2d.F90
r8981 r9150 292 292 ! & zssh(ii ,ij-1) * tmask(ii ,ij-1,1) + & 293 293 ! & zssh(ii ,ij+1) * tmask(ii ,ij+1,1) 294 zcoef = bdytmask(ii-1,ij) + bdytmask(ii+1,ij) + bdytmask(ii,ij-1) + bdytmask(ii,ij+1) 295 zssh(ii,ij) = zssh(ii-1,ij ) * bdytmask(ii-1,ij ) + & 296 & zssh(ii+1,ij ) * bdytmask(ii+1,ij ) + & 297 & zssh(ii ,ij-1) * bdytmask(ii ,ij-1) + & 298 & zssh(ii ,ij+1) * bdytmask(ii ,ij+1) 299 zssh(ii,ij) = ( zssh(ii,ij) / MAX( 1, zcoef) ) * tmask(ii,ij,1) 294 !CEODORIG zcoef = bdytmask(ii-1,ij) + bdytmask(ii+1,ij) + bdytmask(ii,ij-1) + bdytmask(ii,ij+1) 295 !CEODORIG zssh(ii,ij) = zssh(ii-1,ij ) * bdytmask(ii-1,ij ) + & 296 !CEODORIG & zssh(ii+1,ij ) * bdytmask(ii+1,ij ) + & 297 !CEODORIG & zssh(ii ,ij-1) * bdytmask(ii ,ij-1) + & 298 !CEODORIG & zssh(ii ,ij+1) * bdytmask(ii ,ij+1) 299 !CEODORIG zssh(ii,ij) = ( zssh(ii,ij) / MAX( 1, zcoef) ) * tmask(ii,ij,1) 300 !CEOD Just look at diagonals bdytmaks will only be non zero for point required. 301 zssh(ii,ij) = zssh( ii-1, ij-1 ) * bdytmask( ii-1, ij-1) + & 302 & zssh( ii+1, ij+1 ) * bdytmask( ii+1, ij+1) + & 303 & zssh( ii+1, ij-1 ) * bdytmask( ii+1, ij-1) + & 304 & zssh( ii-1, ij+1 ) * bdytmask( ii-1, ij+1) 300 305 ELSE 301 306 ip = bdytmask(ii+1,ij ) - bdytmask(ii-1,ij )
Note: See TracChangeset
for help on using the changeset viewer.