Fortran90 编写理想流体的平面圆柱绕流程序

)*f(zt(i,2))+cc(i,3)*f(zt(i,3))Р?vy(i)=-(bb(i,1)*f(zt(i,1))+bb(i,2)*f(zt(i,2))+bb(i,3)*f(zt(i,3)))Р?v(i)=sqrt(vx(i)**2+vy(i)**2) Р?write(10,'(a,i3,a,f10.6,2x,a,f10.6)'),'第',i,'单元x向流速',vx(i),'y向流速',vy(i) Р?write(11,'(a,i3,a,f10.6)'),'第',i,'个单元流速',v(i) Рend do Рclose(10)Рclose(11)Рallocate(zx(dyz),zy(dyz),zxx(dyz),zyy(dyz),jx(dyz))Рallocate (jtx1(dyz),jty1(dyz),jtx2(dyz),jty2(dyz))Рopen(12,file='夹角.txt') !绘制流场矢量图Рdo i=1,dyz !找到各个单元的中心Р?zx(i)=(x(zt(i,1))+x(zt(i,2))+x(zt(i,3)))/3Р?zy(i)=(y(zt(i,1))+y(zt(i,2))+y(zt(i,3)))/3Р?if(vx(i)==0)then !计算流速与x轴的夹角Р jx(i)=pi/2.0Р?elseР jx(i)=atan(vy(i)/vx(i))Р?end ifР?write(12,'(a,i3,a,f10.6)'),'第',i,'单元流速与x轴夹角',jx(i)Р?zxx(i)=zx(i)+cos(jx(i))*v(i)*0.1 !流线终点x坐标Р?zyy(i)=zy(i)+sin(jx(i))*v(i)*0.1 !流线终点y坐标Р?call line(zx(i),zy(i),zxx(i),zyy(i),ly(:length)) !画流线