| 查看: 2078 | 回复: 7 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
liming9398金虫 (正式写手)
|
[求助]
numerical recipes in fortan 代码的疑惑 已有1人参与
|
|
|
SUBROUTINE flmoon(n,nph,jd,frac) INTEGER jd,n,nph REAL frac,RAD PARAMETER (RAD=3.14159265/180.) INTEGER i REAL am,as,c,t,t2,xtra c=n+nph/4. t=c/1236.85 t2=t**2 as=359.2242+29.105356*c am=306.0253+385.816918*c+0.010730*t2 jd=2415020+28*n+7*nph xtra=0.75933+1.53058868*c+(1.178e-4-1.55e-7*t)*t2 if(nph.eq.0.or.nph.eq.2)then xtra=xtra+(0.1734-3.93e-4*t)*sin(RAD*as)-0.4068*sin(RAD*am) else if(nph.eq.1.or.nph.eq.3)then xtra=xtra+(0.1721-4.e-4*t)*sin(RAD*as)-0.6280*sin(RAD*am) else pause 'nph is unknown in flmoon' endif if(xtra.ge.0.)then i=int(xtra) else i=int(xtra-1.) endif jd=jd+i frac=xtra-i return END fotran77b标准谢的，为什么要用八位或者九位的常数给real数据RAD，as,am赋值？单精度不是最多七位吗，多了有意义吗？ |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
延安大学回应“副教授留绝笔信”：确已离世，正核实举报内容
已经有6人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有5人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求 Risk, Hazards & Crisis in Public Policy期刊缩写
已经有6人回复
- 应虫友要求，共享《Catalysis in C1 Chemistry》
已经有292人回复
- 【转载】物理著作 The Theory of Transport Phenomena in Solids
已经有86人回复
- Nuclear-Power-in-China-Market-Outlook-to-2020
已经有19人回复
- Postdoctoral Position in the MSE in NC state university
已经有15人回复
- 求分子动力学入门fortran程序代码
已经有1人回复
- 急问：连续资助项目的申请代码1必须与原青年科学基金项目保持一致吗？
已经有8人回复
- International Journal for Numerical Methods in Engineering二审时间
已经有13人回复
- International Journal for Numerical Methods in Engineering 审稿周期
已经有10人回复
- Numerical Recipes in Fortran
已经有181人回复
- 【原创+分享】书写金Au(111)面块体电极的代码
已经有4人回复
43fd6ys
木虫 (正式写手)
- 应助: 37 (小学生)
- 金币: 2050.4
- 红花: 7
- 帖子: 379
- 在线: 147小时
- 虫号: 3935320
- 注册: 2015-06-22
- 专业: 零件成形制造
8楼2016-01-09 21:20:09
2楼2012-05-14 11:26:38
liming9398
金虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 780.9
- 散金: 89
- 红花: 1
- 帖子: 361
- 在线: 1010.1小时
- 虫号: 610070
- 注册: 2008-09-23
- 专业: 代数学
3楼2012-05-14 11:51:31
4楼2012-05-14 13:47:26
