求助 haitang大神的 CAD to surf to flca中

B530-2

table.data文件都是零 并且只能算到table132 而且报错  下面是我的命令流  错误代码
大神帮看看 帮我修改一下



surf to table命令
;---------------------------------------------------------------
def setup
a_size = 100000
IO_READ = 0
IO_WRITE = 1
IO_FISH = 0
IO_ASCII = 1
IN_FILE_NAME = suferfile ;the data file from sufer!
end
setup
;----------------SOMETHING SHOULD BE PREDIFINED----------------
;----------------DEFINE THE FUNCTIN TO LOAD DATAFILE---------------
def read_data
        array x_cor(a_size)   ;DIFINE ARRAY TO STORE THE X_COR
        array y_cor(a_size)   ;DIFINE ARRAY TO STORE THE Y_COR
        array zg_cor(a_size)  ;DIFINE ARRAY TO STORE THE Z_COR of Ground
  array zr_cor(a_size)  ;DIFINE ARRAY TO STORE THE Z_COR of Rock
  array sData(a_size)   ;DEFINE ARRAY TO STORE FILE INFO
        array sOut(a_size)    ;DEFINE ARRAY TO STORE FILE TO BE SAVED(Optional,for debug)
;------------------------------------------------------------------------------------
        status = open(IN_FILE_NAME, IO_READ, IO_ASCII)
        status = read(sData, a_size)       
        ReciveNum = status
        if status >= a_size then
                oo = out('Array is not big enough! Please Adjust the dimension of the array!')
                command
                   pause
                endcommand
        endif
        status = close
;------------------------------------------------------------------------------------
        loop n (1,ReciveNum)
                caseof pre_parse(sData(n), 1)
                case 1  
                        x_cor(n)=parse(sData(n), 1)
                case 2
                        x_cor(n)=parse(sData(n), 1)
                case 3
                        oo = out('Error,input x_cor')
                        exit
                case 0
                        oo = out('Error,input x_cor')
                        exit
                endcase
   ;--------------------
                caseof pre_parse(sData(n), 2)
                case 1  
                        y_cor(n)=parse(sData(n), 2)
                case 2
                        y_cor(n)=parse(sData(n), 2)
                case 3
                        oo = out('Error,input y_cor')
                        exit
                case 0
                        oo = out('Error,input y_cor')
                        exit
                endcase
   ;--------------------
                caseof pre_parse(sData(n), 3)
                case 1  
                        zg_cor(n)=parse(sData(n), 3)
                case 2
                        zg_cor(n)=parse(sData(n), 3)
                case 3
                        oo = out('Error,input z_cor')
                        exit
                case 0
                        oo = out('Error,input z_cor')
                        exit
                endcase
         ;--------------------
                caseof pre_parse(sData(n), 4)
                case 1  
                        zr_cor(n)=parse(sData(n), 4)
                case 2
                        zr_cor(n)=parse(sData(n), 4)
                case 3
                        oo = out('Error,input z_cor')
                        exit
                case 0
                        oo = out('Error,input z_cor')
                        exit
                endcase
        endloop       
end

















command   


;*************************************************************
;*************************************************************
;*********   Origin version from Starsmoon *******************
;*************************************************************
;*********        Adapted by Haitang       *******************
;*************************************************************
;*************************************************************
new
def para_set
suferfile = 'Sufer_Data.dat'   ;File from Sufer which contains the layers' coords.
Tablefile = 'Table_data.dat'   ;File which contains the table.
AnsysFile = 'Ansys_Model.dat'  ;Ansys  Model
D3DECFile = 'D3dec_Model.dat'  ;3dec   Model
FlacFile  = 'Flac3D_Model.dat' ;Flac3D Model
Col_Num   = 100                ;number of elements in x-direction
Row_Num   = 66                ;number of elements in y-direction
dx_size   = 22
dy_size   = 22
z_base    = 0
z_size    = 10
xor       = 0.0
yor       = 0.0
n_zon_col = Col_Num - 1
n_zon_row = Row_Num - 1
xdiv      = 1
ydiv      = 1
zdiv      = 2
end
para_set

;==============================================================
cal 01_Surfer_To_Table.fis
read_data    ;------------ [ read sufer network information! ]

cal 02_creat_table.fis
creat_table  ;------------ [ creat table " MyTopoTab.dat " ! ]
cal Tablefile

cal 03_Gen_FlacModel.fis

cal 04_Gen_AnsysModel.fis

cal 05_Gen_3DECModel.fis
;===============================================================

new
pause
cal Flac3D_Model.dat


02_creat_table.fis


;----------------CREAT THE TABLE DATA---------------
;In order to employ the topo.fis,several tables are
;required. The style of the table is the same as the
;tables in the file "Topotab.dat"
;---------------------------------------------------
def creat_table
        ;paremeters list
        ;TAB_NUM---the number of the tabs created
        ;        the tab_num here equals to the number
        ;        of colums(y).Ref to the topotab.dat
        ;ROW_NUM,COL_NUM
       
        status = open(Tablefile, IO_WRITE, IO_ASCII)
        Index=1
        ;for debud
        n_zone_col=Col_Num
        n_zone_row=Row_Num

        ;prepare for the data to be exported
        loop ii (1,ROW_NUM)       
                loop kk (1,COL_NUM)
                n=ii*ROW_NUM+kk
                sOut(Index)='Table'+'  '+string(ii)+'  '+string(x_cor(Index)) + '  ' + string(zg_cor(Index))
                Index=Index+1
                endloop
        endloop
       
        startnum = ROW_NUM + 1
        endnum   = ROW_NUM*2
        totalnum = ROW_NUM*COL_NUM
        loop ii (startnum,endnum)       
                loop kk (1,COL_NUM)
                sOut(Index)='Table'+'  '+string(ii)+'  '+string(x_cor(Index - totalnum)) + '  ' + string(zr_cor(Index - totalnum))
                Index=Index+1
                endloop
        endloop               
       
        TotalWriteNum=n_zone_col*n_zone_row*2
        status = write(sOut, TotalWriteNum)
        if status >= TotalWriteNum then
                oo = out(' Bad number of lines')
        endif
        status = close
       
end


424549.08894647,3894126.38194189,1350.00000000
424548.61102192,3894124.81079630,1350.00000000
424547.44181536,3894121.25025252,1350.00000000
424546.16305583,3894117.84616687,1350.00000000
424544.66929531,3894113.50779939,1350.00000000
424543.76785871,3894110.80617541,1350.00000000
424543.05361416,3894108.71416784,1350.00000000
424541.39904606,3894103.72850294,1350.00000000
424539.94102897,3894099.21007322,1350.00000000
424539.16387657,3894096.45996792,1350.00000000
424538.22352672,3894092.57258228,1350.00000000
424537.78439481,3894090.37534147,1350.00000000
424537.17790959,3894086.28835521,1350.00000000
424536.88485607,3894084.02371183,1350.00000000
424536.19757849,3894078.27727820,1350.00000000
424535.65498717,3894074.23503084,1350.00000000
424535.10911924,3894070.19048919,1350.00000000
424534.73979305,3894066.16579737,1350.00000000
424534.58997102,3894062.50233530,1350.00000000
424534.59375445,3894058.84266296,1350.00000000
424534.62620822,3894055.20306576,1350.00000000
424534.47746479,3894050.84940094,1350.00000000
424534.05709169,3894047.47992376,1350.00000000
424533.35558478,3894043.32640272,1350.00000000
424532.48942821,3894039.05759295,1350.00000000
424531.81822994,3894034.92529421,1350.00000000
424531.73026828,3894031.20138138,1350.00000000
424532.12508169,3894028.05722389,1350.00000000
424533.11237311,3894024.84449675,1350.00000000
424533.91237282,3894021.74234001,1350.00000000
424536.15817677,3894019.65254441,1350.00000000
424539.30880327,3894023.02813409,1350.00000000
424541.08600010,3894025.72037617,1350.00000000
424542.64044319,3894027.85378932,1350.00000000
424544.82049653,3894030.84764643,1350.00000000
424546.83171226,3894033.20567428,1350.00000000
424548.85210519,3894035.55059568,1350.00000000
424550.78474104,3894038.02084742,1350.00000000
424552.67206441,3894040.55581222,1350.00000000
424554.44983481,3894043.24723515,1350.00000000
424556.08707916,3894046.13935045,1350.00000000
424557.70539552,3894049.05849779,1350.00000000
424559.27897296,3894052.04153902,1350.00000000
424560.86000690,3894055.01393127,1350.00000000
424562.38138895,3894058.07151539,1350.00000000
424563.98880737,3894061.00622663,1350.00000000
424565.70635233,3894063.78366058,1350.00000000
424567.60744151,3894066.29896571,1350.00000000
424569.62209868,3894068.65207864,1350.00000000
424571.50999563,3894071.18622428,1350.00000000
424573.27227949,3894073.89976433,1350.00000000
424574.91067099,3894076.79024133,1350.00000000
424576.41599293,3894079.87076171,1350.00000000
424577.86281014,3894083.03483567,1350.00000000
424579.34518905,3894086.14812216,1350.00000000
424580.89410274,3894089.16638695,1350.00000000
424582.52618491,3894092.06587463,1350.00000000
424584.01544673,3894095.16933129,1350.00000000
424585.56206612,3894098.19087270,1350.00000000
424586.77625042,3894100.47066409,1350.00000000
424003.59730442,3894180.73528976,1350.00000000
424009.39509584,3894175.11354444,1350.00000000
424017.26436529,3894168.83128220,1350.00000000
424025.87167832,3894162.92201994,1350.00000000
424033.74331413,3894156.63786708,1350.00000000
424042.45805271,3894151.08258612,1350.00000000
424052.40098958,3894148.58238781,1350.00000000
424062.22141338,3894153.03601576,1350.00000000
424066.06359666,3894162.73507132,1350.00000000
424068.56379521,3894172.67802180,1350.00000000
424072.40597848,3894182.37707736,1350.00000000
424076.83678028,3894191.45066126,1350.00000000
424079.21834484,3894192.18882243,1350.00000000
424079.74724566,3894194.48416985,1350.00000000
424080.51428478,3894198.68357694,1350.00000000
424081.29852610,3894202.89502913,1350.00000000
424081.79803442,3894205.69121876,1350.00000000
424082.74413585,3894211.23189716,1350.00000000
424083.29328039,3894215.27873311,1350.00000000
424083.73757339,3894219.25215138,1350.00000000
424084.28507963,3894223.29784019,1350.00000000
424085.29704536,3894227.66874643,1350.00000000
424086.17711979,3894230.81930654,1350.00000000
424087.13754479,3894234.67802184,1350.00000000
424088.09452834,3894238.54165207,1350.00000000






Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  1  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  2  0  0
Table  3  0  0
Table  3  0  0
Table  3  0  0
Table  3  0  0
Table  3  0  0
Table  3  0  0
Table  3  0  0
Table  3  0  0
Table  3  0  0
Table  3  0  0
Table  3  0  0
Table  3  0  0
Table  3  0  0
Table  3  0  0
Table  3  0  0
Table  3  0  0
Table  3  0  0
Table  3  0  0
Table  3  0  0
Table  3  0  0
Table  3  0  0
Table  3  0  0
Table  3  0  0
Table  3  0  0
Table  3  0  0
Table  3  0  0
Table  3  0  0
Table  3  0  0
Table  3  0  0
Table  3  0  0
Table  3  0  0
Table  3  0  0
Table  3  0  0
Table  3  0  0
Table  3  0  0
Table  3  0  0
Table  3  0  0
Table  3  0  0
Table  3  0  0
Table  3  0  0
Table  3  0  0
Table  3  0  0
Table  3  0  0
   


报错：mesh primitive does not conform to node numbering convertion