24小时热门版块排行榜    

Znn3bq.jpeg
查看: 398  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化工求调剂! +28 RichLi_ 2026-04-06 28/1400 2026-04-09 10:30 by QYQX_123
[考研] 一志愿鲁东大学071000生物学学硕初试分数276求调剂 +3 慕绝cc 2026-04-09 3/150 2026-04-09 09:57 by liuhuiying09
[考研] 286求调剂 +19 Faune 2026-04-08 20/1000 2026-04-09 08:36 by 哦哦123
[考研] 296求调剂 +3 汪!?! 2026-04-08 3/150 2026-04-08 22:00 by zhouyuwinner
[考研] 一志愿华东理工085601材料工程303分求调剂 +15 a1708 2026-04-06 15/750 2026-04-08 16:23 by luoyongfeng
[考研] 085600材料与化工301分求调剂院校 +27 刺痛jk 2026-04-06 28/1400 2026-04-08 16:16 by luoyongfeng
[考研] 263分B区求调剂 +6 李nihao 2026-04-08 6/300 2026-04-08 09:38 by 南开小綦
[考研] 277求调剂 数一104分 +9 瓶子PZ 2026-04-05 14/700 2026-04-07 17:52 by 蓝云思雨
[考研] 316求调剂 +7 yyx想调剂 2026-04-05 7/350 2026-04-07 14:31 by shdgaomin
[考研] 319分085702安全工程求调剂 +6 rious 2026-04-05 6/300 2026-04-07 09:42 by jp9609
[考研] 285求调剂 +15 哦呦呼o 2026-04-04 17/850 2026-04-06 23:02 by chenzhimin
[考研] 一志愿C9的化学工程(085602) 340分,感觉校内调剂无望,求调剂 +12 万事宜臻 2026-04-04 12/600 2026-04-06 07:46 by 无际的草原
[考研] 一志愿北京交通大学材料工程总分358求调剂 +4 cs0106 2026-04-04 4/200 2026-04-05 18:46 by imissbao
[考研] 求调剂 +4 晟功? 2026-04-03 4/200 2026-04-04 21:58 by hemengdong
[考研] 338求调剂 +7 晟功? 2026-04-03 7/350 2026-04-04 20:37 by 蓝云思雨
[考研] 材料295 +13 小英11 2026-04-03 14/700 2026-04-04 09:02 by 来看流星雨10
[考研] 本科985,专业0812分336求调剂 +4 莫莫很行 2026-04-03 4/200 2026-04-03 21:31 by zhq0425
[考研] 315求调剂 +6 顺理成张 2026-04-03 8/400 2026-04-03 14:04 by 百灵童888
[考研] 330求调剂 +3 白神呜呼呼 2026-04-02 3/150 2026-04-03 10:15 by 蓝云思雨
[考研] 一志愿武汉理工0856,初试334 +3 26考研材料 2026-04-02 3/150 2026-04-02 21:22 by dongzh2009
信息提示
请填处理意见