在huang子程序中进行晶体取向更新

    前期介绍过使用umat中特定步输出对应的miller指数并结合matlab程序对变形后的模型进行极图绘制的方法，

参考（https://www.jishulink.com/video/c198695），

类似的处理方式也可以参考技术邻的这篇帖子的处理方法

参考（https://www.jishulink.com/post/1926107），

以及购买某鱼上的插件，然而这类似的实现方案通常比较繁琐，且容易出错。同时对于多滑移系统处理方案之间并不是通用的，因此，在程序内对取向进行更新，同时以状态变量的形式进行保存，可以很容易避免这类问题。由于网上更新取向的通常只是公式的简单介绍，这里尝试结合lingzhi的镁合金更新取向的方案对黄永刚的原始程序进行修改，并以状态变量的形式保存。将三个欧拉角分别存储为SDV123，SDV124，SDV125。

模拟FCC沿着ND方向压缩50%之后的取向分布，并与之前介绍的基于matlab脚本生成欧拉角的方法对比。

输入的初始极图

基于Matlab脚本实现的极图

写入fortran状态变量对应的极图

可以看到两者几乎一致，因此更推荐在Fortran程序内部进行取向更新

这里展示部分修改的Fortran程序，大家可以尝试参考lingzhi的程序在黄程序的基础上进行取向更新的加入

优先建议大家自行尝试修改加入，操作出错无法解决可以取技术邻找到我修改后的Fortran程序，链接如下：  

https://www.jishulink.com/post/1936806