Hi,
Yes, there was a mistake in the analysis of the results.
I printed out the results at the boundaries, and found that the values at the boundaries are proper, as given in BORD subroutine. The high currents are seen at nodes which follow the boundary nodes towards right. Another point to be noted is that these boundary nodes are in the inter tidal area near to which the high current follows. i.e the value of BOTTOM variable at these nodes is above 0. Can the error be contributed because of these nodes are in the intertidal region?
I am trying a case now where no forcing is given to the boundary nodes which lie in the intertidal zone.
Thanks,
Gourish