自从2002年左右用过一段时间Ansys后，一直没再碰过这个软件。最近帮别人计算几个流体力学问题，输出数据再处理时遇到了一些麻烦。数值结果列表中的那些非数字字符可以用命令“/page,34454,132,99999,243”去除。

19    -38.6603000000       21.3985000000       0.00800000000    
20    -38.6603000000       21.0985000000       24.00000000000    
 
  NODE        X                   Y                   Z
21    -38.6603000000       20.7985000000       80.00000000000    
22    -38.6603000000       20.4985000000       56.90320000000    

节点坐标每隔20行就出现“NODE”和“X”等字符。很少使用Ansys，系统学习它命令流太不划算，用熟悉的fortran77处理吧！在后面的代码中，n为节点总数，xy和xy.dat文件分别为处理前后的文件名。用gfortran a.for编译,用./a.out执行。

   implicit real*8(a-h,o-z)
   parameter(n=8095)
   character*4 cn
   character*1 cx,cy,cz
 
   open(1,file='xy')
   open(2,file='xy.dat')
   kk=0
   ii=0
       do 1 i=1,n
   kk=kk+1
   do 2 while(kk.le.20)
   ii=ii+1
   if(ii.le.n) then  
   read(1,*)nnn,x,y,z
   write(2,2012)nnn,x,y,z
   endif
   kk=kk+1
2    continue
   kk=0
   if(ii.le.n) then  
   read(1,*)cn,cx,cy,cz
   endif
1    continue      
   close(1)
   close(2)
2012    format(I7,f9.4,f9.4,f9.4)    
         
   end