尝试解决Multiwfn基于sTDA的光谱图绘制问题

大村驴 · 发表于 Post on 2025-11-30 20:58:48

使用Multiwfn基于sTDA的输出绘制ECD图，当激发态数量超过10000时，列编号会变成“****”，导致Multiwfn无法读取，只能手动列编辑替换。

为此尝试稍微改了一下Multiwfn代码，主要是在spectrum.f90文件里加了两行

复制代码

...
do i=1,numdata
!Note: The last four columns of tda.dat correspond to f_length, f_velocity, R_length, R_velocity
! read(10,*) inouse,datax(i),fl,fv,Rl,Rv
read(10,*) dummy, datax(i), fl, fv, Rl, Rv ! add line
if (ispectrum==3) then
...

复制代码

修改后再次编译，即可正常读取超过10000行的tda.dat。

对比两张图，修改前200nm附近缺了几个激发态，修改后已补齐。
修改后的代码已上传。

spectrum.zip (61.1 KB, 下载次数 Times of downloads: 1)

tda.zip (293.09 KB, 下载次数 Times of downloads: 2)

sobereva · 发表于 Post on 2025-12-1 05:43:38

下次更新的Multiwfn会顺带解决此问题
更简单来说，read(10,*) inouse,datax(i),fl,fv,Rl,Rv改成 read(10,*) c80tmp,datax(i),fl,fv,Rl,Rv即可

大村驴 · 发表于 Post on 2025-12-1 10:55:23

sobereva 发表于 2025-12-1 05:43
下次更新的Multiwfn会顺带解决此问题
更简单来说，read(10,*) inouse,datax(i),fl,fv,Rl,Rv改成 read(10,* ...

谢谢老师！

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[波函数分析交流] 尝试解决Multiwfn基于sTDA的光谱图绘制问题