凝固计算后，提取厚度，总是写入文件不成功

huangxin211

我在凝固计算后，提取厚度，想写入文件，总是不成功，但不知道错误在哪，谁能帮我看一下，谢谢了
这样的一个1/4圆的一个二维非稳态凝固。需要得到如图所示的那条红线随时间的温度变化，那条红线上的温度小于某一值即表示该处为凝固坯壳开始位置，获得该点的位置信息，写入文件。我写的如下：/post1
/uis,msgpop,3
*set,x1,0
*set,x2,0.19
*set,y1,0
*set,y2,0
*set,xdivd,190+1
*dim,x,array,86                        !定义一个临时存储节点坐标2的数组
*dim,cout,array,xdivd
*dim,tem,array,86                          !定义一个临时存储节点温度2的数组
*dim,xy,array,xdivd,2                          !定义一个存储节点坐标1的数组
*dim,NodeTemp,array,xdivd                      !定义一个存储节点温度1的数组
*dim,xx,array,86            !转换为mm单位存储节点X坐标              
*dim,ttem,array,86                  
*cfopen,huaigang,dat                                !输出的文件名为：huaigang.dat
!写TecPlot文件的文件头
*vwrite
('TITLE=counter')
*vwrite
('VARIABLES= "Time","X"')
*do,n,1,86,1  !整个计算时间
set,,,,,1

     path,1,2,,xdivd-1                               !定义路径            
     ppath,1,,x1,y2
     ppath,2,,x2,y2
     pdef,temp,temp
     paget,pathdat,table  
        *do,k,1,xdivd,1   
               xy(k,1)=pathdat(k,4,1)                    !在xy(k,1)中存储点的X坐标
               xy(k,2)=y2                                     !在xy(k,2)中存储点的Y坐标
               NodeTemp(k)=pathdat(k,5,1)                 !在NodeTemp(k)中存储点的温度
         *enddo
      padele,all
      pathdat(1,1,1)=
         *do,k,1,xdivd,1
              *if,NodeTemp(k),lt,1420,and,NodeTemp(k),gt,0,then   !获取有效节点
                tem(n)=NodeTemp(k)
                x(n)=xy(k,1)
               *endif
          *enddo

*vwrite,xdivd
('ZONE I='f6.1', f=point')
      xx(n)=x(n)*1000
      ttem(n)=tem(n)
      *vwrite,time,xx(n)
(3f12.6)
xy(1,1)=
NodeTemp(1)=
tem(1)=
x(1)=
*CFCLOSE
*enddo


      
                     
这之中有做过调整，但总是不能成功写入，谁能帮我看下，谢谢了

lansejixian

请问楼主你解决了没有，最近我也再用数据提取，也是总不成功，希望能从你这学习

huangxin211

lansejixian 发表于 2013-5-13 16:55
请问楼主你解决了没有，最近我也再用数据提取，也是总不成功，希望能从你这学习 ...

还没有，不知道什么原因