24小时热门版块排行榜    

Znn3bq.jpeg
查看: 397  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化工求调剂! +27 RichLi_ 2026-04-06 27/1350 2026-04-08 22:57 by jean5056
[考研] 293调剂 +20 yj1221 2026-04-08 21/1050 2026-04-08 22:37 by 猪会飞
[考研] 085801 总分275 本科新能源 求调剂 +4 bradoner 2026-04-08 5/250 2026-04-08 21:27 by cyh—315
[考研] 273求调剂 +41 麦小叮当 2026-04-06 48/2400 2026-04-08 15:16 by screening
[考研] 一志愿生物与医药,296分,求调剂 +14 66鹿 2026-04-03 16/800 2026-04-08 10:38 by tjzhao
[考研] 生物学363调剂求助 +7 fanzhang6666 2026-04-06 9/450 2026-04-07 17:37 by lijunpoly
[考研] 080500求调剂 +12 黄宇博 2026-04-06 12/600 2026-04-07 12:41 by upczlm1989
[考研] 372分材料与化工(085600)英二数二求调剂 +4 蓝笺片 2026-04-06 4/200 2026-04-07 12:30 by dongzh2009
[考研] 材料调剂 +17 小刘同学吖吖 2026-04-06 18/900 2026-04-07 11:41 by 诗与自由
[考研] 软工学硕299求调剂 +6 useryy 2026-04-07 6/300 2026-04-07 09:50 by vgtyfty
[考研] 求调剂 +4 电气小神童 2026-04-04 6/300 2026-04-07 00:14 by guanxin1001
[考研] 304求调剂 +4 luoye0105 2026-04-05 4/200 2026-04-06 21:05 by 木子君1218
[考研] 22408 331分求调剂 +4 y__1 2026-04-06 4/200 2026-04-06 17:26 by 土木硕士招生
[考研] 269电子信息求调剂,可转专业 +5 独酌wl 2026-04-06 5/250 2026-04-06 17:23 by 土木硕士招生
[考研] 一志愿南航,数一英一学硕317求调剂!! +6 Acaciad 2026-04-04 6/300 2026-04-06 12:13 by 考研学校招点人
[考研] 313求调剂 +5 海日海日 2026-04-04 5/250 2026-04-05 15:52 by jndximd
[考研] 368求调剂 +5 今华习 2026-04-03 7/350 2026-04-04 18:47 by imissbao
[考研] 325求调剂 +4 春风不借意 2026-04-04 4/200 2026-04-04 14:46 by 湘农储能材料
[考研] 一志愿双非085502,267分,过四级求调剂 +3 再忙也要吃饭啊 2026-04-03 3/150 2026-04-04 05:03 by gswylq
[考研] 数二英二348求调剂 +4 hxdzj1 2026-04-03 5/250 2026-04-03 21:25 by zhq0425
信息提示
请填处理意见