求助：关于将输出结果转化为molden正则坐标问题

wang_

       最近想做振动模式分配，但是苦于CRYSTAL的输出文件无法转化为输入，sobereva老师有一篇文章是使CP2K计算的振动模式可以被GaussView观看的程序：MfakeG（http://sobereva.com/656），所想着能否通过这种方法制作molden文件实现
       但是由于我是用CRYSTAL计算的，无法直接得到molden格式的文件，所以想仿照着老师给出的格式进行改写，但是在简振振动坐标[FR-NORM-COORD]这一块碰上了问题，sobereva老师有一个帖子基于fch中的Hessian矩阵计算振动频率的简单程序Hess2freq（http://bbs.keinsci.com/forum.php ... rdertype=1&tid=3368）解释了Hessian矩阵转换为振动频率计算流程，说明了[FR-NORM-COORD]这一块是归一化后的振动坐标，另外，molden的红外强度似乎与GaussianView显示的强度不一样，是有什么关系吗
       我的问题是，我的输出结果附带说明是经过归一化的，但是我计算其平方和结果却不是等于1，我想请教下我如何才能将这个结果转化为molden文件格式的要求呢