求助大神！！！小弟初学lammps，想模拟液体粘度，看到“compute viscosity/cos command”命令可以做相关计算，但是不清楚具体使用方法，请各位大神指点迷津。
下面是in文件
units                                           real                        #单位
dimension                                3                                #维度
boundary                        p p p                #边界类型

pair_style                                lj/cut 10.0
atom_style                      full                                                                #原子类型
bond_style                                 harmonic                                #键类型－弹簧势
angle_style                         harmonic                                #角类型
dihedral_style                 harmonic                                #二面角类型

read_data                                216-PEG200-NVE1.data  #读入ｄａｔａ文件
timestep        0.1                                                #时间步长１ｆｓ
neighbor       4.0 bin                                #近邻边界２埃
neigh_modify     every 10 delay 0 check yes                #每１０步更新一次
thermo                100 #
fix                                                        1 all nve

run              20000

compute  cos all viscosity/cos
variable V equal c_cos[7]
variable A equal 0.02E-5
variable density equal density
variable lz equal lz
variable reciprocalViscosity equal v_V/${A}/v_density*39.4784/v_lz/v_lz*100

fix 2 all accelerate/cos 0.02e-5
run  5000
print            "reciprocalViscosity : ${reciprocalViscosity }"
执行以后提示：
compute  cos all viscosity/cos
ERROR: Unrecognized compute style 'viscosity/cos' (../modify.cpp:1206)
Last command: compute  cos all viscosity/cos
已经安装了官方说明书里提示的“USER-MISC package” 返回小木虫查看更多