(i8,2x,i8,2x,i8,2x,i8)") np,ne,nm,nxРwrite(10,"(f8.0,2x,f8.3,2x,f8.3)")((node(i,j),j=1,3),i=1,np)Рwrite(10,"(i8,2x,i8,2x,i8,2x,i8,2x,i8,2x,i8)")((element(i,j),j=1,6),i=1,ne)Рwrite(10,"(i8,2x,i8,2x,i8)")((fix(i,j),j=1,3),i=1,nx)Рwrite(10,"(f8.0,2x,e16.3e3,2x,f8.3,2x,f10.3,2x,f8.3)")(mat(i),i=1,5)Рwrite(10,"(i8,2x,i8,2x,i8)")node_fn,Elem_fn,Face_FnРif (node_fn>0)thenР write(10,"(i8,2x,i8,2x,i8)")((node(i,j),j=1,3),i=1,node_fn) Рend ifРif (face_fn>0)thenР do i=1,face_fnР write(10,"(i8,2x,i8,2x,i8,2x,i8,2x,i8,2x,i8)")(facef(i,j),j=1,6)Р write(10,"(10(i8))")(apply(i,j),j=1,facef(i,2))Р end doРend ifРcall calc()РpauseРendРFEMQ4.FORР SUBROUTINE CALC()Р CHARACTER*hР DIMENSION SK(80000),COOR(2,300),AE(4,11),MEL(5,200),Р & WG(4),JR(2,300),MA(600),R(600),iew(30),STRE(3,200)Р COMMON /CMN1/ NP,NE,NM,NR
全文预览
