| 查看: 1928 | 回复: 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赋值？单精度不是最多七位吗，多了有意义吗？ |
回复此楼» 猜你喜欢
今天系统多次维护，明天很可能放榜！
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有74人回复
-
岩土工程学报什么时候才能终审完呐
已经有4人回复
-
时间戳又变了8-15
已经有25人回复
-
2027广东省杰青
已经有6人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有9人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有26人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有16人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有12人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求 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
