[交流] 【求助】一个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);          }

回复此楼