| 查看: 4260 | 回复: 11 | ||
[求助]
差分电荷密度的计算 已有2人参与
|
||
| 生成的CHGCAR文件网格下面的数据太大，我把三个系统的数据分别导入到excel中相减，但由于数据太多，excel中只能导入一部分，无法全部倒入或不能导入，这样肯定结果不全，我想问前辈有什么方法吗？谢谢！ |
回复此楼» 猜你喜欢
今天系统多次维护，明天很可能放榜！
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有74人回复
-
岩土工程学报什么时候才能终审完呐
已经有4人回复
-
时间戳又变了8-15
已经有25人回复
-
2027广东省杰青
已经有6人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有9人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有26人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有16人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有12人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于差分密度研究电荷转移中的坐标问题
已经有3人回复
- 【求教】电子密度，差分密度，原子电荷这三个之间有什么关系？
已经有7人回复
- vasp如何计算差分电荷密度
已经有7人回复
- 如何做电子密度差分图
已经有16人回复
- 【求助】求助差分电荷密度的脚本！
已经有45人回复
- 【求助】做差分电荷密度
已经有5人回复
- 【求助】差分电荷密度
已经有19人回复
hn19870519
木虫 (正式写手)
- 1ST强帖: 1
- 应助: 45 (小学生)
- 金币: 1225.2
- 散金: 1861
- 红花: 11
- 帖子: 379
- 在线: 203.1小时
- 虫号: 867945
- 注册: 2009-10-10
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2011-10-22 23:22:16
hn19870519
木虫 (正式写手)
- 1ST强帖: 1
- 应助: 45 (小学生)
- 金币: 1225.2
- 散金: 1861
- 红花: 11
- 帖子: 379
- 在线: 203.1小时
- 虫号: 867945
- 注册: 2009-10-10
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★
youzhizhe(金币+1): 建议直接上传附件。谢谢交流。 2011-10-22 23:43:38
youzhizhe(金币+1): 建议直接上传附件。谢谢交流。 2011-10-22 23:43:38
|
subroutine chgdiff !################################################ !DMol3 total electron density !(1p,e12.5) ! 9.000 8.200 9.400 90.000 90.000 90.000 ! 45 41 47 ! 1 -23 22 -20 21 -22 25 !################################################# !initializing all paremeters implicit none integer :: NGX,NGY,NGZ real*8 :: amod,bmod,cmod,alpha,beta,gamma,VOL real*8 :: scale real*8,allocatable :: a(  ,b(,c(real*8,allocatable :: atomnum( real*8,allocatable :: chgdata1(:,:, ,chgdata2(:,:,integer :: i,j,k,n,m,iserror,p open(11,file='CHGCAR1',form='formatted',status='old',iostat=iserror) open(13,file='CHGCAR2',form='formatted',status='old',iostat=iserror) if(iserror>0) then write(*,*)"The CHGCAR file does not exist, please try again!" stop endif !读入原子种类数量 write(*,*)"---------------------------------" write(*,*)"Please input atomic species amount: " read(*,*)n !read initialized paremeters and discard head sections read(11,*) read(11,*)scale allocate(a(3)) allocate(b(3)) allocate(c(3)) allocate(atomnum(n)) read(11,*)a(1:3) read(11,*)b(1:3) read(11,*)c(1:3) read(11,*)atomnum(1:n) read(11,*) do i=1,sum(atomnum(1:n))+1 read(11,*) enddo read(11,*)NGX,NGY,NGZ !compute angle call getmod(a,scale,amod) call getmod(b,scale,bmod) call getmod(c,scale,cmod) alpha=acos(dot_product(a,c)*(scale**2)/(amod*cmod))*180/3.1415926 beta=acos(dot_product(b,c)*(scale**2)/(bmod*cmod))*180/3.1415926 gamma=acos(dot_product(a,b)*(scale**2)/(amod*bmod))*180/3.1415926 !charge part open(12,file='ms.grd',status='replace') write(12,'(A)')"DMol3 spin density" write(12,'(A)')'(1p,e12.5)' write(12,'(6f8.3)')amod,bmod,cmod,alpha,beta,gamma write(12,'(3I5)')NGX-1,NGY-1,NGZ-1 write(12,'(7I5)')1,-(NGX/2-1),NGX/2,-(NGY/2-1),NGY/2,-(NGZ/2-1),NGZ/2 allocate(chgdata1(NGX,NGY,NGZ)) allocate(chgdata2(NGX,NGY,NGZ)) !--将电荷密度数据存入数组chgdata中---- read(11,*)chgdata1 do i=1,57 read(13,*) enddo read(13,*)chgdata2 !------------------------------------- VOL=scale*scale*scale*(a(1)*b(2)*c(3)+a(2)*b(3)*c(1)+a(3)*b(1)*c(2)& &-a(1)*b(3)*c(2)-a(2)*b(1)*c(3)-a(3)*b(2)*c(1)) write(*,*)VOL do i=1,NGZ do j=1,NGY do k=1,NGX write(12,'(6E12.5)')(chgdata1(k,j,i)-chgdata2(k,j,i))/VOL enddo enddo enddo write(*,*)"Write charge density successfully!!!" end subroutine |
3楼2011-10-22 23:23:51
hn19870519
木虫 (正式写手)
- 1ST强帖: 1
- 应助: 45 (小学生)
- 金币: 1225.2
- 散金: 1861
- 红花: 11
- 帖子: 379
- 在线: 203.1小时
- 虫号: 867945
- 注册: 2009-10-10
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2011-10-22 23:25:26
我对程序很无奈啊！我确实不会，希望大侠拯救一下下啊！！！ 5楼2011-10-23 11:39:06
6楼2011-10-23 11:42:42
|
本帖内容被屏蔽 |
7楼2011-10-24 09:00:29
8楼2011-10-24 10:44:21
sgquan
铁杆木虫 (著名写手)
- 应助: 46 (小学生)
- 金币: 6381.4
- 散金: 1692
- 红花: 22
- 帖子: 2964
- 在线: 1377.4小时
- 虫号: 1328308
- 注册: 2011-06-22
- 性别: GG
- 专业: 冶金物理化学与冶金原理
9楼2011-10-24 11:18:08
10楼2014-01-13 11:03:39
25
