if( abs(yc(i,j,bi,bj)*deg2rad) .lt. pio3 ) then
sinarg = 3.*yc(i,j,bi,bj)*deg2rad/2.
ypart = yc(j)
yc(1)=2500.
yc(j)=yc(j-1)+degY
print*,'yc(1), yc(Ny): ',yc(1), yc(Ny)
double precision depth(Nx,Ny),xc(Nx),yc(Ny),degX,degY