基于mfix-19.1.2版本
方法一：直接输出差值网格固相速度
注：这种方式只适用于garg 2012颗粒差值格式
在DEM中，默认是无法输出固相速度的网格值的：
但是通过搜索des文件夹下V_s关键字可以看到，实际上网格固相速度是会通过差值计算得到的：
并且在计算DEM颗粒脉动速度的地方也用到了网格固相速度：
说明网格固相速度只是默认没有写出，因此只需要人为写出一下即可。
首先找到网格输出的代码位置：
由于MMAX代表TFM的固相数，因此在循环范围部分需要修改，默认只写出TFM的网格固相速度。修改如下：
最后在dat文件里手动添加一下VTK_U_S(VTK, PHASE)这个关键字（且确保颗粒的差值格式为garg 2012）例如：
结果如下：
方法二：保存网格内颗粒的平均速度到网格
网格颗粒速度还可以仿照DES_GRANULAR_TEMPERATURE函数内计算网格拟颗粒温度的方式，对网格内颗粒速度进行统计平均。
按照J. Musser的建议，可以把需要后处理的变量存到ReactionRates这个数组中，虽然这个数组本来是用来后处理反应速率的，不过由于可以直接在gui里设置该数组的个数以及输出频率，而且该变量仅仅用来保存输出，计算中未被使用，因此可以灵活使用，来保存其他需要保存的变量。
例如这里我要保存网格内颗粒Y方向的平均速度V_s，可以先设置一个ReactionRates数组，以备后面使用：
在des_granular_temperature中USE这个数组：
定义一个局部变量用于后面保存网格内颗粒速度的和：
将颗粒速度累加到对应网格内的M相内：
除以网格内M相颗粒个数，得到网格内颗粒的平均速度，并保存到ReactionRates(ijk,1)中。
注意拟颗粒温度由于其定义跟维数有关，因此要多除以一个维数，但是速度不需要：
添加代码如下：
现在，只要调用这个函数，就会往ReactionRates中写入固相速度，因此只要在输出的部分，每次输出之前调用一次DES_GRANULAR_TEMPERATURE就可以了：
重新编译计算可以看到输出Y方向速度：