关于fix shake以及compute命令的应用问题

我最近写了一个关于超临界水的in文件，里面用到了fix shake命令，但运行的时候一直报错ERROR: Invalid fix style，我对照manuel检查了命令格式，没有发现问题，禁了fix shake后运行又出现ERROR: Invalid compute style，都没查出问题在哪里，请各位大侠不吝赐教，小女子先谢过了
以下是我的in文件
units real
dimension 3
boundary p p p
atom_style full
pair_style lj/charmm/coul/charmm 7.5 8.0 8.0 10.0
bond_style harmonic
angle_style harmonic 
read_data data.Mol 
angle_coeff 1  55.0000   104.5200 
bond_coeff  1  450.0000   0.9572 
pair_coeff 1 1 0.1521 3.1507
pair_coeff 1 2 0.0836 1.7753
pair_coeff 2 2 0.0460 0.4000
region 1 block INF INF INF INF INF INF 
velocity all create 640 100 
velocity all scale 640 
#group water type 1 2
fix 1 all shake 0.0001 20 0 b 1 a 1
fix 2 all npt temp 640 760 1000 iso 236.92 236.92 1000 drag 1
timestep 0.01
compute 1 all ke