24小时热门版块排行榜

>论坛更新日志 (3561)
>虫友互识 (437)
>休闲灌水 (219)
>导师招生 (199)
>文献求助 (157)
>考博 (60)
>论文道贺祈福 (56)
>博后之家 (53)
>论文投稿 (53)
>硕博家园 (48)
>考研 (44)
>基金申请 (34)
>公派出国 (29)
>招聘信息布告栏 (28)
>找工作 (19)
>教师之家 (15)

返回列表

zml90h

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 108.6
红花: 5
帖子: 209
在线: 51小时
虫号: 1824115
注册: 2012-05-18
性别: GG
专业: 金属材料的制备科学与跨学

[求助] 请问一下，应该从命令窗口输入什么，格式是什么样子的。

下面一段Fortran代码会在命令窗口提示用户输入转轴和转角，但是我不明白输入的格式是什么，因为我刚刚接触到Fortran写的程序，所以看不懂。因此在这里请求各位虫友们帮小弟看看。谢谢了。主要是前几行。

c  define rotation of lattice axes
c
2 write(*,'(a)',advance='NO') ' axis (x/y/z/q) and degree : '
   read(*,*,end=3) axis,degree
   degr = 3.1415926535898d0 * degree / 180.0d0
   if(axis(1:1).eq.'q' .or. axis(1:1).eq.'Q') goto 3
   if(axis(1:1).eq.'x' .or. axis(1:1).eq.'X') then
      rotation(1,1) = one
      rotation(2,1) = zero
      rotation(3,1) = zero
      rotation(1,2) = zero
      rotation(2,2) = dcos(degr)
      rotation(3,2) = - dsin(degr)
      rotation(1,3) = zero
      rotation(2,3) = dsin(degr)
      rotation(3,3) = dcos(degr)
      call matrix_product(tmp_rotation,rotation,all_rotation)
      tmp_rotation = all_rotation
   elseif(axis(1:1).eq.'y' .or. axis(1:1).eq.'Y') then
      rotation(1,1) = dcos(degr)
      rotation(2,1) = zero
      rotation(3,1) = dsin(degr)
      rotation(1,2) = zero
      rotation(2,2) = one
      rotation(3,2) = zero
      rotation(1,3) = - dsin(degr)
      rotation(2,3) = zero
      rotation(3,3) = dcos(degr)
      call matrix_product(tmp_rotation,rotation,all_rotation)
      tmp_rotation = all_rotation
   elseif(axis(1:1).eq.'z' .or. axis(1:1).eq.'Z') then
      rotation(1,1) = dcos(degr)
      rotation(2,1) = - dsin(degr)
      rotation(3,1) = zero
      rotation(1,2) = dsin(degr)
      rotation(2,2) = dcos(degr)
      rotation(3,2) = zero
      rotation(1,3) = zero
      rotation(2,3) = zero
      rotation(3,3) = one
      call matrix_product(tmp_rotation,rotation,all_rotation)
      tmp_rotation = all_rotation
   endif

回复此楼