【报告一个Multiwfn/GaussView兼容性问题】关于multiwfn输出Gaussian型cube文件

yjcmwgk

multiwfn输出Gaussian型cube文件的时候，当某个格点函数值大于等于1.00 x 10^100的时候，输出为.cube文件的时候会写成1.00000+100而不是1.00000E+100
比如下面这一行里：

  1.12207E+52  6.41195E+69  4.37285E+87  1.02032+101  1.28909+105  1.67292E+98

红色的两个数字，这是GaussView无法读取的数据。
谢谢sob！！！！！！！
Multiwfn Version 3.6(dev), release date: 2019-May-14

=====我的临时措施=======
直接吧cube文件读进Python里，扫描每个格点的数值，凡是倒数第四位出现+的，一律改成E。比如，1.28909+105改成1.28909E105后，即可被GaussView正常读取。

sobereva

其实这不是bug，而是Fortran本身的输出格式特征
https://software.intel.com/en-us ... mpiler/topic/271144

为了对于数值很大的情况能另其它程序兼容，我已经改了导出cube文件时用的格式了，现在应该没问题了，我已更新了官网上的程序。

让你变成回忆

我也遇到过类似的问题。
我利用Multiwfn计算Overlap矩阵的时候，当计算的结果非常小的时候，程序输出如下：
  592    588 -0.130156D-53 -0.279543D-26  0.212582D-25  0.978887D-29 -0.652195D-26
  593    589 -0.594195D-57 -0.127619D-29  0.978887D-29 -0.183792D-27 -0.297745D-29
  594    590  0.257790-133 0.273744D-39 -0.272875D-38 -0.212888D-41  0.150792D-38
  595    591  0.273744D-39  0.194355D-22 -0.113531D-21 -0.885733D-25  0.627382D-22
  596    592  0.272875D-38  0.113531D-21 -0.654953D-21 -0.517397D-24  0.366482D-21
因为这里是调用了showmatgau子程序来输出的，所以我是通过修改输出的格式（由默认的改为了f14.6的形式），这样就不会输出成科学计数的形式。

sobereva

让你变成回忆 发表于 2019-12-5 08:35
我也遇到过类似的问题。
我利用Multiwfn计算Overlap矩阵的时候，当计算的结果非常小的时候，程序输出如下 ...

输出格式后头加个E3即可解决，比如1PE16.8改成1PE16.8E3，这样指数总是输出成3位且不会把前头的E给顶了

让你变成回忆

sobereva 发表于 2019-12-5 08:41
输出格式后头加个E3即可解决，比如1PE16.8改成1PE16.8E3，这样指数总是输出成3位且不会把前头的E给顶了

好的，我去修改一下。
谢谢sob老师！

yjcmwgk

sobereva 发表于 2019-12-5 07:06
其实这不是bug，而是Fortran本身的输出格式特征
https://software.intel.com/en-us ... mpiler/topic/2711 ...

谢谢sob！是我表述有问题。已经把【报告一个小bug】改成【报告一个Multiwfn/GaussView兼容性问题】