| 查看: 1050 | 回复: 2 | |||
[交流]
【求助】紧急求教读取CHGCAR的matlab脚本，谢谢
|
| 现在需要深入研究vasp的CHGCAR，希望用matlab读取它，但对矩阵的读取一直没查到相关资料，求教哪位有读能提供取CHGCAR的matlab小脚本，谢谢。 |
回复此楼» 猜你喜欢
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有7人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有3人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 继续求助Matlab绘图！
已经有8人回复
- 弱弱问一下：stereochemical activity 是什么意思？ 谢谢
已经有5人回复
- 紧急求助翻译一篇德文文献~翻译其中合成的一部分~急求~谢谢~
已经有4人回复
- matlab 计算过程中工作空间变量保存及读取 求解决内存不足的方法
已经有5人回复
- 帮翻译几段文献，请不要在线翻译，谢谢
已经有2人回复
- 用matlab统计图片中某种颜色成分占总图片中所有颜色成分的比率
已经有4人回复
- 【求助】弱弱的问一下Matlab数据提取的问题
已经有5人回复
- 【求助】用matlab编程，需要将数据自检分类（高手进来指点下）
已经有5人回复
- 【求助】CHGCAR用什么来读取
已经有10人回复
- 【求助】matlab如何只读取图像的一部分？
已经有11人回复
- 【求助】有哪位大侠做过次氯酸叔丁酯啊？帮帮我啊！谢谢
已经有6人回复
- 【求助】matlab读取fortran输出的txt文件【已解决】
已经有3人回复
- 【求助】求教matlab解非线性方程组
已经有9人回复
- 【求助】Weickert的各向异性扩散方程的滤波方法的matlab程序代码
已经有8人回复
» 抢金币啦！回帖就可以得到:
-
项目结题验收，希望一切顺利
+1/971
-
急招科研助理/博士后
+2/124
-
美国肯塔基大学Dr. Sheng Tong博士后招聘（空间可控的基因组编辑与肿瘤免疫治疗）
+1/94
-
【推免】】中科院国家级人才团队招收27级推免硕士生
+1/72
-
福建农林大学王志刚教授团队诚聘青年教师、博士后、博士生 - 材料、化学、力学、高分�
+1/41
-
中山大学中法核工程与技术学院放射生物学课题组招收推免生
+1/37
-
山东征女友，坐标济南
+1/36
-
福建农林大学王志刚教授团队诚聘青年教师、博士后、博士生
+1/36
-
中科院大连化学物理研究所刘生忠研究员团队招收联合培养博士
+1/23
-
澳门大学陈修平课题组招收2027年博士
+1/15
-
深圳大学应用技术学院招聘凝聚态物理博士后
+1/15
-
Biofunctional Materials纳米材料赋能肿瘤新型治疗专题征稿（Scopus、DOAJ收录）
+1/13
-
生物功能应用的仿生超湿材料专题征稿（Scopus、DOAJ 收录）
+1/12
-
【生物材料与生物打印专题征稿】Scopus期刊，免APC
+1/5
-
中国科学院院士成会明招收2027届高分子方向推免直博生（截止2026年9月19日）
+1/4
-
北理工集成电路杰青团队 | 诚招科助理
+1/3
-
口腔特刊客座招募客编Biofunctional Materials牙科生物功能材料专题，Scopus期刊
+1/3
-
信息提供：中国科学院院士成会明招收2027届推免直博生（截止2026年9月19日）
+1/2
-
北京理工大学-集成电路与电子学院杰青团队-招博士后
+1/2
-
北京理工大学-集成电路与电子学院杰青团队-招博士后
+1/1
★ ★
sunyang1988(金币+2): 谢谢帮助。这样直接粘贴可能有格式错误，如果方便请上传到网盘 2011-01-16 14:10:22
雁儿霏霏(金币+20): 谢谢大虾，能给提供个程序文件吗？粘贴有一些格式错误。邮箱：shangyan2009@gmail.com ，谢谢！ 2011-01-16 19:37:24
sunyang1988(金币+2): 谢谢帮助。这样直接粘贴可能有格式错误，如果方便请上传到网盘 2011-01-16 14:10:22
雁儿霏霏(金币+20): 谢谢大虾，能给提供个程序文件吗？粘贴有一些格式错误。邮箱：shangyan2009@gmail.com ，谢谢！ 2011-01-16 19:37:24
|
% Read in CHG file to 3D matrix. fid = fopen('CHG','r'); % Read POSCAR part of CHG file: sysname = fgetl(fid); lconst = fscanf(fid,'%f',1); basisvec = fscanf(fid,'%f',[3,3])'; temp1 = fgetl(fid); temp2 = fgetl(fid); species = sscanf(temp2,'%i'); Natoms = sum(species); temp3 = fgetl(fid); atompos_f3 = fscanf(fid,'%f',[3,Natoms])'; volume = abs(dot(basisvec(1,  ,cross(basisvec(2,,basisvec(3,)))*lconst^3;% Read charge density data. (i,j,k) value in i+(j+k*Ny)*Nx, with coordinate % system defined by basis vectors!: gridsize = fscanf(fid,'%i',3); Nx = gridsize(1); Ny = gridsize(2); Nz = gridsize(3); temp4 = fscanf(fid,'%f',[Nx*Ny*Nz,1]); density_sub_h1_surf=reshape(temp4,Nx,Ny,Nz)/volume; % % Magnitization (if present): % temp5 = fscanf(fid,'%i',3); % temp6 = fscanf(fid,'%f',[Nx*Ny*Nz,1]); % magn_lowacc_t2_N�蕆eshape(temp6,Nx,Ny,Nz)/volume; % minval = min(min(min(density))); % maxval = max(max(max(density))); % For rectangular unit cell, unit axes: % Plot isosurface: % prect = patch(isosurface(density,7.5), 'FaceColor', 'red', 'EdgeColor', 'none'); % isonormals(density,prect); % %axis tight; % camlight; lighting phong; % Non-rectangular (general) unit cell: % Need to created position 3D matrices xtemp = zeros(Nx*Ny*Nz,1); ytemp = zeros(Nx*Ny*Nz,1); ztemp = zeros(Nx*Ny*Nz,1); for k = 0:Nz-1 for j = 0:Ny-1 for i = 0:Nx-1 xtemp(1+i+(j+k*Ny)*Nx) = i/Nx*basisvec(1,1)+j/Ny*basisvec(2,1)+k/Nz*basisvec(3,1); ytemp(1+i+(j+k*Ny)*Nx) = i/Nx*basisvec(1,2)+j/Ny*basisvec(2,2)+k/Nz*basisvec(3,2); ztemp(1+i+(j+k*Ny)*Nx) = i/Nx*basisvec(1,3)+j/Ny*basisvec(2,3)+k/Nz*basisvec(3,3); end end end X = reshape(xtemp,Nx,Ny,Nz)*lconst; Y = reshape(ytemp,Nx,Ny,Nz)*lconst; Z = reshape(ztemp,Nx,Ny,Nz)*lconst; z = squeeze(Z(1,1, );ave_z_h1 = squeeze(sum(sum(density_sub_h1_surf,1),2))/(Nx*Ny); %(sqrt(2)*Nx*Ny)*lconst^2; figure %subplot(211) plot(z,ave_z_h1,'k') % figure % p = patch(isosurface(X,Y,Z,density,7.5), 'FaceColor', 'red', 'EdgeColor', 'none'); % isonormals(density,p); % % maxx = max(basisvec(:,1))*lconst; % % minx = min(basisvec(:,1))*lconst; % % maxy = max(basisvec(:,2))*lconst; % % miny = min(basisvec(:,2))*lconst; % % maxz = max(basisvec(:,3))*lconst; % % minz = min(basisvec(:,3))*lconst; % % axis([minx maxx miny maxy minz maxz]) % camlight; lighting phong; % Extract charge density along certain lines: ------------------------------ |
2楼2011-01-16 12:51:04
3楼2011-12-30 22:59:32
