************************************2022.3.16*******************************************
该代码在Julia程序中运行,但由于2015年以后就没有更新,不能用目前的Julia版本(1.4.2)运行
有两种方案:用当时的版本(估计为0.6或更早的版本)进行计算,但需要把完全匹配才行,包括一系列的相互依赖的运行包
对代码进行手动更新,已经完成,适配1.4.2
****************************************************************************************
出图需要三个文件:.xyz .vtk .cube
.xyz 定义了原子位置,用openbabel即可,但建议用软件内置的程序
include("src/framework.jl")
replicate_cssr_to_xyz(${yourstructurename})
.vtk 定义了边界
include("src/framework.jl")
write_unitcell_boundary_vtk("${yourstructurename}")
.cube包含能量信息,建议gridspacing取小一些
include("src/energygrid.jl")
writegrid("CH4", "IRMOF-1", "UFF", gridspacing=1.0, cutoff=12.5)
***************************************************************************************
导入到visit后如果发现.xyz文件所展示的原子范围太大,可以将.xyz文件转换为excel文件,再对照.vtk文件对边界的定义进行筛选,再导出为.xyz文件即可