Ik probeer de vergelijkingen er ook nog even netje bij te maken, maar dat is nog best lastig. Willem 22 Dec 2008 17:10 (CET)
∂WX,Y∂t⋅DX⋅DYX=ux,Y⋅Wx,Y*DYx−ux+1,Y⋅Wx+1,Y⋅DYx+1+vX,y*WX,y*DX−vX,y+1*WX,y+1*DX
∂ux,Y∂t⋅DX⋅DYx⋅WxY=impulsnormaal+impulstangentieel+zwaartekracht+kromming+coriolis+getijkracht+frictie
impulsnormaal=Wx−1,Y⋅ux−1,y2⋅DYx−1−Wx+1,Y⋅ux+1,y2⋅DYx+12
impulstangentieel=Wx,y⋅vx,y2⋅DX−Wx+1,y⋅ux+1,y2⋅DX
dwdt[x][y]=u[x][y]*uW[x][y]*dyo[x]−u[x+1][y]*uW[x+1][y]*dyo[x+1]+(v[x][y]*vW[x][y]−v[x][y+1]*vW[x][y+1])*DX;W[x][y]+=dwdt[x][y]*DT/DX/dym[x];dudt[x][y]=(uW[x−1][y]*u[x−1][y]*u[x−1][y]*dyo[x−1]−uW[x+1][y]*u[x+1][y]*u[x+1][y]*dyo[x+1])/2.0+//impulstransport.5*g*(H[x−1][y]−H[x][y])*dyo[x]*uW[x][y]+//potentieleenergie(LOW[x][y]*(v[x][y]+v[x−1][y])*(u[x][y]+u[x][y−1])−ROW[x][y]*(v[x][y+1]+v[x−1][y+1])*(u[x][y+1]+u[x][y]))/4.0*DX+//tangentieelimpulstransport−u[x][y]/R*dyo[x]*DX*vopu[x][y]*sino[x]+//curvature−2*omega*vopu[x][y]*sino[x]*uW[x][y]*DX*dyo[x]+//corioliskrachtuGM[x][y]*uW[x][y]*DX*dyo[x]+//maangetijkracht−u[x][y]*0.01*DX*dyo[x];//frictiedvdt[x][y]=(vW[x][y−1]*v[x][y−1]*v[x][y−1]−vW[x][y+1]*v[x][y+1]*v[x][y+1])/2.0*DX+.5*g*(H[x][y−1]−H[x][y])*vW[x][y]*DX+(LOW[x][y]*(u[x][y]+u[x][y−1])*(v[x][y]+v[x−1][y])*dyo[x]−LBW[x][y]*(u[x+1][y]+u[x+1][y−1])*(v[x+1][y]+v[x][y])*dyo[x+1])/4.0+uopv[x][y]/R*dym[x]*DX*uopv[x][y]*sinm[x]+2*omega*uopv[x][y]*sinm[x]*vW[x][y]*DX*dym[x]+vGM[x][y]*vW[x][y]*DX*dym[x]+−v[x][y]*0.01*DX*dym[x];dudt[x][y]/=uW[x][y]*DX*dyo[x];dvdt[x][y]/=vW[x][y]*DX*dym[x];