等参元画三角单元的电位云图.docxVIP

!C-----!等参元画三角单元的电位云图---------------------------------------- !C program EMF2D--finite element method!C!C definition of the variables!C nodei(e),nodej(e),nodem(e):the number of i j m nodes of triangle element e!C l(n):address of diagonal elements of coefficient matrix!C sl:limit of 1-D variable band storage range!C k(i):1-D varible band storage array of coefficient matrix!C u(i):potential matrix!C bnode(i):the number of first kind boundary nodes!C bptal(i):potential of first kind boundary nodes!C cpm:characteristic parameter of medium!C sd:source density!C number1:the last number of elements in medium 1(source area) parameter(lnd=700,lel=1200) common/xy/lastnd,nm(lnd),x(lnd),y(lnd),bnode(lnd),bptal(lnd)/elem/lastel,p(lel),nodei(lel),nodej(lel),nodem(lel) open(2,file=jdf.txt,status=old)!节点数、单元数、 open(3,file=jdxy.txt,status=old) !节点x，y坐标 open(4,file=dyijm.txt,status=old)!单元编号和节点编号 open(7,file=jdu.txt)!节点电位值 open(1,file=ccc.txt)! call readst call chazhi call huitu READ (*,*) end subroutine readst parameter(lnd=700,lel=1200) common/xy/lastnd,nm(lnd),x(lnd),y(lnd),bnode(lnd),bptal(lnd)/elem/lastel,p(lel),nodei(lel),nodej(lel),nodem(lel) read(2,*)lastnd read(2,*)lastel do i=1,lastnd read(3,*)nm(i) read(3,*)x(i) read(3,*)y(i) enddo read(4,*)(p(i),nodei(i),nodej(i),nodem(i),i=1,lastel) read(7,*)(bnode(i),bptal(i),i=1,lastnd) close(2) close(3) close(4) close(7) return end subroutine huitu USE MSFLIB parameter(lnd=700,lel=1200) common/xy/lastnd,nm(lnd),x(lnd),y(lnd),bnode(lnd),bptal(lnd)/elem/lastel,p(lel),nodei(lel),nodej(lel),nodem(lel) TYPE(XYCOORD)::xy TYPE (wxycoord)::wxy CHARACTER tfont*8 CHARACTER string(11)*1 INTEGER n(7) integer kyun INTEGER(2) status REAL(8) wx,wy real(8) p1,p2 real cc,up real umin,umax real u(0:9) write(*,*)请输入或,画云图输，画等位线输 read(*,*)kyun string(10)=0 string(1)=1 string(2)=2 string(3)=3 string(4)=4 string(5)=5 string(6)=6 string(7)=7 string(8)=8 string(9)=9 string(11)=. umin=bptal(1) umax