 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 930 | 回復: 2 | |||
_yixiao鐵蟲 (小有名氣)
|
[交流]
【求助】一個fortran程序轉換為VC++ 程序 已有2人參與
|
|
這個一個關于建立體心立方結構的小程序，用fortran編的，現(xiàn)在想用VC++ 6.0編寫，可是我編寫后（我沒用動態(tài)數組，這個應該問題不大吧？）能夠編譯，但是運行不出結果，不曉得哪有出錯 麻煩各位大俠幫忙看看，能不能編寫看看。。。或者看看我編的哪錯了。。。十分感謝！ 原fortran程序： program structure implicit none integer :: i,ix,iy,iz,ncellx,ncelly,ncellz,iatoms,natoms,ibas !整數 real*8 :: xunit(2),yunit(2),zunit(2),a !三個一維的實數數組 real*8,allocatable,dimension(  :: xbas,ybas,zbas !定義三個一維動態(tài)數組a=2.8664d0 ! Fe的晶格常數 ！ 單胞原子坐標 xunit(1)=0.0d0 yunit(1)=0.0d0 zunit(1)=0.0d0 xunit(2)=0.50d0 yunit(2)=0.50d0 zunit(2)=0.50d0 ！ 在xyz方向重復次數 write(*,*) 'Numbers of unit cells in x,y,z?' read (*,*) ncellx, ncelly, ncellz iatoms=0 natoms=2 * ncellx * ncelly * ncellz allocate(xbas(natoms)) !給動態(tài)數組分配內存 allocate(ybas(natoms)) allocate(zbas(natoms)) do ix=1,ncellx do iy=1,ncelly do iz=1,ncellz do ibas=1,2 iatoms=iatoms+1 xbas(iatoms)=(float(ix-1)+xunit(ibas))*a ybas(iatoms)=(float(iy-1)+yunit(ibas))*a zbas(iatoms)=(float(iz-1)+zunit(ibas))*a enddo enddo enddo enddo ！ 輸出一個xyz文件 open(11,file='Fe.xyz',status='new') do i=1,iatoms write(11,10)i,xbas(i),ybas(i),zbas(i) enddo 10 format(i6,1x,7(f10.5,1x),4(i3,1x)) end program 我用VC++ 6.0改的程序： void CBCC_CrystalDlg::OnCalculate() { // TODO: Add your control notification handler code here FILE * fileFe_xyz = fopen("Fe_xyz.txt", "w+"  ;fprintf(fileFe_xyz, "i xbas ybas zbas\n" ;a=2.8664; xunit[0]=0.0; yunit[0]=0.0; zunit[0]=0.0; xunit[1]=0.50; yunit[1]=0.50; zunit[1]=0.50; ncellx=M; ncelly=M; ncellz=M; iatoms=0; natoms=2*ncellx*ncelly*ncellz; for(int i=0;i xbas=0; ybas=0; zbas=0; } for(ix=1;i<=ncellx;ix++) { for(iy=1;iy<=ncelly;iy++) { for(iz=1;iz<=ncellz;iz++) { for(ibas=0;ibas<2;ibas++) { xbas[iatoms]=(float(ix-1)+xunit[ibas])*a; ybas[iatoms]=(float(iy-1)+yunit[ibas])*a; zbas[iatoms]=(float(iz-1)+zunit[ibas])*a; iatoms=iatoms+1; } } } } for(i=0;i fprintf(fileFe_xyz, "%6d %10.5f %10.5f %10.5f\n", i,xbas,ybas,zbas); } fclose(fileFe_xyz); } |
回復此樓» 猜你喜歡
化學求調劑
已經有5人回復
-
一志愿中國石油大學（華東） 本科齊魯工業(yè)大學
已經有3人回復
-
332求調劑
已經有4人回復
-
265求調劑
已經有9人回復
-
328求調劑，英語六級551，有科研經歷
已經有8人回復
-
求調劑
已經有3人回復
-
278求調劑
已經有6人回復
-
307求調劑
已經有10人回復
-
279分求調劑 一志愿211
已經有12人回復
-
297求調劑
已經有9人回復
2樓2010-06-04 10:30:24
stonesu
鐵蟲 (小有名氣)
- 應助: 0 (幼兒園)
- 金幣: 1079.9
- 散金: 55
- 帖子: 120
- 在線: 81.7小時
- 蟲號: 852187
- 注冊: 2009-09-20
- 專業(yè): 金屬功能材料
3樓2010-06-23 18:29:43
相關版塊跳轉
我要訂閱樓主 _yixiao 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿中國石油大學（華東） 本科齊魯工業(yè)大學 +3 | 石能偉 2026-03-17 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 278求調劑 +6 | 煙火先于春 2026-03-17 | 6/300 |
|
|
|
[考研] 280求調劑 +7 | 咕嚕曉曉 2026-03-18 | 8/400 |
|
|
|
[考研] 一志愿西南交大，求調劑 +5 | 材化逐夢人 2026-03-18 | 5/250 |
|
|
|
[考研] 288求調劑 +16 | 于海海海海 2026-03-19 | 16/800 |
|
|
|
[考研] 一志愿 南京航空航天大學大學 ，080500材料科學與工程學碩 +5 | @taotao 2026-03-20 | 5/250 |
|
|
|
[考研] 一志愿中國海洋大學，生物學，301分，求調劑 +5 | 1孫悟空 2026-03-17 | 6/300 |
|
|
|
[考研] 0703化學調劑 +10 | 妮妮ninicgb 2026-03-15 | 14/700 |
|
|
|
[考研] 復試調劑 +4 | z1z2z3879 2026-03-14 | 6/300 |
|
|
|
[考研] 0703化學調劑 +4 | 18889395102 2026-03-18 | 4/200 |
|
|
|
[考研] 085600材料與化工求調劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|
|
|
[考研] 344求調劑 +6 | knight344 2026-03-16 | 7/350 |
|
|
|
[考研] 311求調劑 +6 | 26研0 2026-03-15 | 6/300 |
|
|
|
[考研] 304求調劑 +12 | 小熊joy 2026-03-14 | 13/650 |
|
|
|
[考研] 0703化學調劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
|
[考研] 283求調劑 +3 | 聽風就是雨； 2026-03-16 | 3/150 |
|
|
|
[考研] 304求調劑 +4 | ahbd 2026-03-14 | 4/200 |
|
|
|
[考研] 318求調劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
|
[考研] 0703化學調劑 290分有科研經歷，論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
|
[考研] 復試調劑 +3 | 呼呼？~+123456 2026-03-14 | 3/150 |
|