Border Condition Algorithm in Python

1

I can not run this in my python, anyone who can help me? The problem is in border conditions.

for it in range(2,nt):
    uu=u*u;
    uv=u*v;
    vv=v*v;
    hbu=(h-b)*u;
    hbv=(h-b)*v;
    gh=g*h;
    # Avançar o sinal em todo o domínio excepto nas fronteiras
    for ix in range(2,Nx-1):
        for iy in range(2,Ny-1):
            uP[ix,iy] = (1/4. * (u[ix-1,iy] + u[ix+1,iy] + u[ix,iy-1] + u[ix,iy+1])) - dt2dx * (uu[ix+1,iy] - uu[ix-1,iy]) - dt2dy * (uv[ix,iy+1] - uv[ix,iy-1]) + dt2dx * (g*h[ix+1,iy] - g*h[ix-1,iy])
            vP[ix,iy] = (1/4. * (v[ix-1,iy] + v[ix+1,iy] + v[ix,iy-1] + v[ix,iy+1])) - dt2dx * (uv[ix+1,iy] - uv[ix-1,iy]) - dt2dy * (vv[ix,iy+1] - vv[ix,iy-1]) + dt2dy * (gh[ix,iy+1] - gh[ix,iy-1])
            hP[ix,iy] = (1/4. * (h[ix-1,iy] + h[ix+1,iy] + h[ix,iy-1] + h[ix,iy+1])) - dt2dx * (hbu[ix+1,iy] - hbu[ix-1,iy]) - dt2dy * (hbv[ix,iy+1] - hbv[ix,iy-1]) 

  #condições fronteira
    hP[1,:]=hP[2,:];
    hP[Nx,:]=hP[Nx-1,:];
    hP[:,1]=hP[:,2];
    hP[:,Ny]=hP[:,Ny-1];

    uP[1,:]=0;
    uP[Nx,:]=0;
    uP[:,1]=uP[:,2];
    uP[:,Ny]=uP[:,Ny-1];

    vP[1,:]=vP[2,:];
    vP[Nx,:]=vP[Nx-1,:];
    vP[:,1]=0;
    vP[:,Ny]=0;

     # actualizar o sinal
    h[:]=hP[:]; u[:]=uP[:]; v[:]=vP[:];

    hStation[it] = h[ixStation, iyStation]
    
asked by anonymous 21.04.2017 / 21:06

0 answers