GULP采用经典势能模型输出shengBTE计算用力常数文件

仅需在option参数中使用supercell控制2阶力常数的扩胞倍数，以及output shengBTE:

但是注意，gulp这里存在一个bug，输出的2阶力常数文件中超胞中原子次序与phonopy以及shengBTE用不一致：shengBTE中FORCE_CONSTANT中原子次序为：ucatom*uz*uy*ux；但是gulp中为ucatom*ux*uy*uz，对于各向同性系统，则计算没有区别，但是对于各向异性系统，会产生计算错误，需要转换2阶力常数的中的原子次序。