| 查看: 935 | 回复: 2 | |||
[交流]
【求助】一个fortran程序转换为VC++ 程序 已有2人参与
|
|
这个一个关于建立体心立方结构的小程序，用fortran编的，现在想用VC++ 6.0编写，可是我编写后（我没用动态数组，这个应该问题不大吧？）能够编译，但是运行不出结果，不晓得哪有出错 麻烦各位大侠帮忙看看，能不能编写看看。。。或者看看我编的哪错了。。。十分感谢！ 原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 !定义三个一维动态数组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)) !给动态数组分配内存 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); } |
回复此楼» 猜你喜欢
081700，311，求调剂
已经有15人回复
-
一志愿北京化工085600 310分求调剂
已经有18人回复
-
085600材料与化工301分求调剂院校
已经有4人回复
-
材料与化工371求调剂
已经有14人回复
-
336材料与化工085600求调剂
已经有7人回复
-
材料334求调剂
已经有18人回复
-
331求调剂
已经有8人回复
-
332求调剂
已经有17人回复
-
一志愿南京航空航天大学 材料与化工329分求调剂
已经有4人回复
-
材料专硕322
已经有7人回复
2楼2010-06-04 10:30:24
stonesu
铁虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1079.9
- 散金: 55
- 帖子: 120
- 在线: 81.7小时
- 虫号: 852187
- 注册: 2009-09-20
- 专业: 金属功能材料
3楼2010-06-23 18:29:43
