| 查看: 159 | 回复: 0 | |||
| 当前主题已经存档。 | |||
benbei626825金虫 (小有名气)
|
[交流]
【讨论】新手上路，这是编写的fcc晶体结构，看看对不对
|
||
因为刚开始学习编写程序，所以想让大家帮忙看看，写对了没有？自己觉得编的有些繁琐，看能不能编写简单点？谢谢了啊 CUA表示的是晶格常数，R1数组就表示坐标这是100 方向的 DO 10 K=-8,10,1 DO 20 J=-10,10,1 DO 30 I=-10,10,1 N(1)=N(1)+1 IF (MOD(K,2)==1) THEN R1(1,N(1),1)=(I-1)*CUA IF (MOD(J,2)==0) THEN R1(1,N(1),1)=(I-0.5)*CUA END IF ELSE R1(1,N(1),1)=(I-0.5)*CUA IF (MOD(J,2)==0) THEN R1(1,N(1),1)=(I-1)*CUA END IF END IF R1(1,N(1),2)=(J-1)*CUA/2.0 R1(1,N(1),3)=(K-1)*CUA/2.0 30 CONTINUE 20 CONTINUE 10 CONTINUE 110方向的 DO 80 K=-10,10,1 DO 90 J=-10,10,1 DO 100 I=-10,10,1 N(2)=N(2)+1 IF (MOD(K,2)==1) THEN R1(2,N(2),1)=(I-1)*CUA/SQRT(2） ELSE R1(2,N(2),1)=I*CUA/SQRT(2.) END IF R1(2,N(2),2)=(J-1)*CUA R1(2,N(2),3)=(K-1)*CUA/SQRT(2.) 100 CONTINUE 90 CONTINUE 80 CONTINUE 111方向的 DO 150 K=-10，10,1 DO 160 J=-10,10,1 DO 170 I=-10,10,1 N(3)=N(3)+1 IF (MOD(K,3)==1) THEN IF (MOD(J,2)==1) THEN R1(3,N(3),1)=(I-1)*CUA/SQRT(2.0) R1(3,N(3),2)=(J-1)*SQRT(6.0)*CUA/4.0 ELSE R1(3,N(3),1)=(2*I-1)*CUA*SQRT(2.0)/4.0 R1(3,N(3),2)=(J-1)*CUA*SQRT(6.0)/4.0 END IF ELSEIF (MOD(K,3)==2) THEN IF (MOD(J,2)==1) THEN R1(3,N(3),1)=(2*I-1)*CUA*SQRT(2.0)/4.0 R1(3,N(3),2)=(J-1)*CUA*SQRT(6.0)/4+SQRT(6.0)*CUA/12.0 ELSE R1(3,N(3),1)=(I-1)*CUA/SQRT(2.0) R1(3,N(3),2)=(J-1)*CUA*SQRT(6.0)/4.0+SQRT(6.0)*CUA/12.0 END IF ELSE IF (MOD(J,2)==1)THEN R1(3,N(3),1)=(I-1)*CUA*SQRT(2.0)/2.0 R1(3,N(3),2)=(J-1)*CUA*SQRT(6.0)/4+CUA/SQRT(6.0) ELSE R1(3,N(3),1)=(2*I-1)*CUA*SQRT(2.0)/4.0 R1(3,N(3),2)=(J-1)*CUA*SQRT(6.0)/4.0+CUA/SQRT(6.0 END IF R1(3,N(3),3)=(K-1)*CUA/SQRT(3.0) END IF 170 CONTINUE 160 CONTINUE 150 CONTINUE [ Last edited by benbei626825 on 2009-3-9 at 22:07 ] |
回复此楼» 猜你喜欢
上海大学实验技术岗位非升即走
已经有9人回复
-
云南大学材料与能源学院解琳课题组钙钛矿博士招生
已经有4人回复
-
博士申请
已经有6人回复
-
评审有感
已经有11人回复
-
面上本子正文33页，违规吗？会被低分嘛？
已经有10人回复
-
国自然上会要求
已经有9人回复
-
今年审到国自然15份，谈谈感受
已经有16人回复
-
考博自荐
已经有6人回复
-
青C资助名额大幅增加！
已经有16人回复
-
重磅！青年科学基金项目（C类）资助增幅预计超过50%
已经有10人回复
30