24小时热门版块排行榜    

查看: 387  |  回复: 1

zml90h

银虫 (小有名气)

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

下面一段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
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

做一个拥有哲学头脑的技术人
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

zml90h

银虫 (小有名气)

如何放弃求助啊 问题已经解决了
做一个拥有哲学头脑的技术人
2楼2013-12-19 11:19:41
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zml90h 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料工程专硕调剂 +5 204818@lcx 2026-03-17 6/300 2026-03-18 22:55 by 204818@lcx
[考研] 321求调剂 +3 何润采123 2026-03-18 3/150 2026-03-18 21:27 by li123456789.
[考研] 0854可跨调剂,一作一项核心论文五项专利,省、国级证书40+数一英一287 +8 小李0854 2026-03-16 8/400 2026-03-18 14:35 by 搏击518
[考研] 288求调剂,一志愿华南理工大学071005 +4 ioodiiij 2026-03-17 4/200 2026-03-18 12:36 by Linda Hu
[考研] 0703化学求调剂 总分331 +3 ZY-05 2026-03-13 3/150 2026-03-18 10:58 by macy2011
[考研] 0703化学调剂 +3 妮妮ninicgb 2026-03-17 3/150 2026-03-18 10:29 by macy2011
[考研] 211本,11408一志愿中科院277分,曾在中科院自动化所实习 +6 Losir 2026-03-12 7/350 2026-03-17 12:09 by danranxie
[考研] 285化工学硕求调剂(081700) +9 柴郡猫_ 2026-03-12 9/450 2026-03-17 10:18 by Sammy2
[考研] 考研调剂 +3 淇ya_~ 2026-03-17 5/250 2026-03-17 09:25 by Winj1e
[考研] 11408 一志愿西电,277分求调剂 +3 zhouzhen654 2026-03-16 3/150 2026-03-17 07:03 by laoshidan
[考研] 304求调剂 +4 ahbd 2026-03-14 4/200 2026-03-16 16:48 by 我的船我的海
[考研] 283求调剂 +10 小楼。 2026-03-12 14/700 2026-03-16 16:08 by 13811244083
[考研] 085600调剂 +5 漾漾123sun 2026-03-12 6/300 2026-03-16 15:58 by 漾漾123sun
[考研] 070303 总分349求调剂 +3 LJY9966 2026-03-15 5/250 2026-03-16 14:24 by xwxstudy
[考研] 070305求调剂 +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
[考研] 265求调剂 +4 威化饼07 2026-03-12 4/200 2026-03-14 17:23 by userper
[考研] 一志愿哈工大材料324分求调剂 +5 闫旭东 2026-03-14 5/250 2026-03-14 14:53 by 木瓜膏
[考研] 招收0805(材料)调剂 +3 18595523086 2026-03-13 3/150 2026-03-14 00:33 by 123%、
[考研] 295求调剂 +3 小匕仔汁 2026-03-12 3/150 2026-03-13 15:17 by vgtyfty
[考研] 0817化学工程与技术考研312分调剂 +3 T123 tt 2026-03-12 3/150 2026-03-13 10:49 by houyaoxu
信息提示
请填处理意见