前期介紹過使用umat中特定步輸出對應的miller指數并結合matlab程序對變形后的模型進行極圖繪制的方法，

參考（http://www.yqgqt.org.cn/video/c198695），

類似的處理方式也可以參考技術鄰的這篇帖子的處理方法

參考（http://www.yqgqt.org.cn/post/1926107），

以及購買某魚上的插件，然而這類似的實現方案通常比較繁瑣，且容易出錯。同時對于多滑移系統處理方案之間并不是通用的，因此，在程序內對取向進行更新，同時以狀態變量的形式進行保存，可以很容易避免這類問題。由于網上更新取向的通常只是公式的簡單介紹，這里嘗試結合lingzhi的鎂合金更新取向的方案對黃永剛的原始程序進行修改，并以狀態變量的形式保存。將三個歐拉角分別存儲為SDV123，SDV124，SDV125。

模擬FCC沿著ND方向壓縮50%之后的取向分布，并與之前介紹的基于matlab腳本生成歐拉角的方法對比。

輸入的初始極圖

基于Matlab腳本實現的極圖

寫入fortran狀態變量對應的極圖

可以看到兩者幾乎一致，因此更推薦在Fortran程序內部進行取向更新

這里展示部分修改的Fortran程序，大家可以嘗試參考lingzhi的程序在黃程序的基礎上進行取向更新的加入

優先建議大家自行嘗試修改加入，操作出錯無法解決可以取技術鄰找到我修改后的Fortran程序，鏈接如下：  

http://www.yqgqt.org.cn/post/1936806