24小时热门版块排行榜

>论坛更新日志 (2975)
>虫友互识 (718)
>文献求助 (156)
>休闲灌水 (133)
>导师招生 (132)
>考博 (123)
>教师之家 (79)
>博后之家 (71)
>硕博家园 (71)
>找工作 (71)
>论文投稿 (68)
>论文道贺祈福 (64)
>基金申请 (61)
>考研 (54)
>公派出国 (49)
>招聘信息布告栏 (31)

返回列表

wangwen552

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 24.9
红花: 1
帖子: 140
在线: 85.2小时
虫号: 2723215
注册: 2013-10-14
性别: MM
专业: 固体力学

[求助] 求助下面一段fortran程序的意思

open(unit=Fout,file='out.dat')
      write(fout,'(A)') 'zone '
      do i=1,counter
            if(z(i)>z0min+5.0.and. z(i)<z0max-5.0) then
            write(fout,'(5(3x e15.6))') x(i),y(i),z(i),vone(i)
   :  ,rot(i)
                     if(vone(i)<5.0) then
                              AllStra=vone(i)+AllStra
                              numall=numall+1
                     endif
            endif
      enddo

      write(fout,'(A)') 'zone '
      do i=1,gbnum
            j=gbid(i)
            if(z(j)>z0min+5.0.and. z(j)<z0max-5.0) then
            write(fout,'(5(3x e15.6))') x(j),y(j),z(j),vone(j)
   :  ,rot(j)
                     if(vone(j)<5.0) then
                              GbStra=vone(j)+GbStra
                              numgb=numgb+1
                     endif
            endif
      enddo
      GbStra=GbStra
      AllStra = AllStra
      print*, '%Gb atoms:', 1.0*numgb/numall
      print*,'--------------------strain %------------------'
      print*,'       Ave          %GB       %ResT'
      write(*,'(3(3x e15.6))') AllStra/numall,GbStra/AllStra,
   :  1.0-GbStra/AllStra

c       print*,'e=',AllStra/numall,'GB%=',GbStra/AllStra,
c    : 'Rest%=', 1.0-GbStra/AllStra

      print*,'--------------------strain ave------------------'
      print*,'       Ave          GBe       ResTe'
      write(*,'(3(3x e15.6))') AllStra/numall,
   : GbStra/numgb,
   : (AllStra-GbStra)/(numall-numgb)
      print*,'-------------normal-strain ave------------------'
      write(*,'(3(3x e15.6))') AllStra/numall,
   : GbStra/numgb/AllStra*numall,
   : (AllStra-GbStra)/(numall-numgb)/AllStra*numall
c       print*,'Av e=',AllStra/numall,'GB=',GbStra/numgb,
c    : 'Rest=',(AllStra-GbStra)/(numall-numgb)

      close(fout)

我只知道这段程序能生成一个out.dat文件，不懂是怎么生成的、如何计算的。求程序大神帮忙

回复此楼