位移差？应力差？利用fish语言实现结果文件间运算。

;建立模型并储存结果文件model newzone create brick size 20 20 20 group 'a'zone group 'b' range position-x 5 15 position-z 5 15zone gridpoint initialize pore-pressure 1e3 model save 'a'zone gridpoint initialize pore-pressure 0 range group 'a'model save 'b';设置文件读写参数&获取数组所需sizemodel restore 'a.f3sav'fish def setup    io_read = 0     io_write = 1    io_ascii = 1    arr_size = gp.numend@setup;提取结果a中的孔隙压力值并写出到文件fish def pp_out    array id1(arr_size)    array pp1(arr_size)    n = 0    loop foreach gp1 gp.list        n = n+1        id1(n) = gp.id(gp1)        pp1(n) = gp.pp(gp1)    endloop    status = file.open('id_list',io_write,io_acsii)    status = file.write(id1,string(arr_size))    status = file.close    status = file.open('pp_list',io_write,io_acsii)    status = file.write(pp1,string(arr_size))    status = file.closeend@pp_out;设置文件读写参数&获取数组所需sizemodel restore 'b.f3sav'fish def setup    io_read = 0    io_write = 1    io_ascii = 1    arr_size = gp.numend@setup;读取结果a的孔隙压力值并与结果b的孔隙压力值相减fish def pp_minus    array id(arr_size)    array pp(arr_size)    status = file.open('id_list',io_read,io_acsii)    status = file.read(id,arr_size)    status = file.close    status = file.open('pp_list',io_read,io_acsii)    status = file.read(pp,arr_size)    status = file.close    loop i(1,arr_size)        gp = gp.find(id(i))        gp.pp(gp) = pp(i) - gp.pp(gp)    endloopend@pp_minusmodel save 'c'