• sidama
    了解作者
  • Python
    开发工具
  • 1KB
    文件大小
  • rar
    文件格式
  • 0
    收藏次数
  • 1 积分
    下载积分
  • 1
    下载次数
  • 2020-06-22 18:04
    上传日期
将VASP生成的CHGCAR格式文件合并后转变成Tecplot格式文件
CHGCAR2Tecplot.rar
  • CHGCAR2tecplot.py
    1.3KB
  • chgsum.pl
    2.1KB
内容介绍
import numpy as np import glob # to convert a 'CHGCAR' style file to a tecplot style file filename='CHGCAR' inp=[ line for line in open(filename) if line.strip()] scale=float(inp[1]) cell = np.array([line.split() for line in inp[2:5]], dtype=float) cell *= scale ta = inp[5].split() tb = inp[6].split() Nelem = np.array(tb,dtype=int) Nions = Nelem.sum() NX=int(inp[8+Nions].split()[0]) NY=int(inp[8+Nions].split()[1]) NZ=int(inp[8+Nions].split()[2]) ELF=np.zeros(NX*NY*NZ) inp=[line for line in open(filename) if line.strip()] N=0 for i in range(9+Nions,len(inp)): if inp[i].split()[0].isalpha(): break else: for j in range(len(inp[i].split())): ELF[N]=inp[i].split()[j] N=N+1 ELF=ELF.reshape((NX,NY,NZ,-1),order="F") f=open(str(filename)+'.dat','wt') print >>f, 'TITLE = "Example: Simple 3D-Volume Data"' print >>f, 'VARIABLES = "X", "Y", "Z", "Fz"' print >>f, 'ZONE I=%d, J=%d, K=%d, F=POINT' %(NX,NY,NZ) position=np.zeros(3) for k in range(NZ): for j in range(NY): for i in range(NX): position[0]=(i+0.5)/float(NX) position[1]=(j+0.5)/float(NY) position[2]=(k+0.5)/float(NZ) positionC=np.dot(position,cell) print >>f, '%-12.3f%12.3f%12.3f%12.3f' %(positionC[0],positionC[1],positionC[2],ELF[i,j,k,0])
评论
    相关推荐
    • matlabcnhelp.rar
      matlab中文帮助很难找的,快速下载
    • MobilePolice.rar
      移动警察,车牌识别,车牌定位系统源代码,已经运用在移动车载稽查系统中。
    • SVM(matlab).rar
      支持向量机(SVM)实现的分类算法源码[matlab]
    • svm.zip
      用MATLAB编写的svm源程序,可以实现支持向量机,用于特征分类或提取
    • Classification-MatLab-Toolbox.rar
      模式识别matlab工具箱,包括SVM,ICA,PCA,NN等等模式识别算法,很有参考价值
    • VC++人脸定位实例.rar
      一个经典的人脸识别算法实例,提供人脸五官定位具体算法及两种实现流程.
    • QPSK_Simulink.rar
      QPSK的Matlab/Simulink的调制解调仿真系统,给出接收信号眼图及系统仿真误码率,包含载波恢复,匹配滤波,定时恢复等重要模块,帮助理解QPSK的系统
    • LPRBPDemo2009KV.rar
      车牌识别,神经网络算法,识别率高达95%,识别时间低于80ms。
    • MODULATION.RAR
      这个源程序代码包提供了通信系统中BPSK,QPSK,OQPSK,MSK,MSK2,GMSK,QAM,QAM16等调制解调方式 用matlab的实现,以及它们在AWGN和Rayleigh信道下的通信系统实现及误码率性能
    • algorithms.rar
      十大算法论文,包括遗传算法,模拟退火,蒙特卡罗法等等,对于初学者很有帮助!!