| 查看: 1602 | 回复: 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赋值？单精度不是最多七位吗，多了有意义吗？ |
回复此楼» 猜你喜欢
博士读完未来一定会好吗
已经有6人回复
-
小论文投稿
已经有3人回复
-
Bioresource Technology期刊，第一次返修的时候被退回好几次了
已经有9人回复
-
心脉受损
已经有3人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有8人回复
-
申请2026年博士
已经有6人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
求 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人回复
jackplay
金虫 (小有名气)
- 应助: 4 (幼儿园)
- 金币: 1120.8
- 红花: 1
- 帖子: 128
- 在线: 87.6小时
- 虫号: 2945895
- 注册: 2014-01-24
- 性别: GG
- 专业: 分析化学
7楼2016-01-09 08:21:30
2楼2012-05-14 11:26:38
liming9398
金虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 780.9
- 散金: 89
- 红花: 1
- 帖子: 361
- 在线: 1009.2小时
- 虫号: 610070
- 注册: 2008-09-23
- 专业: 代数学
3楼2012-05-14 11:51:31
4楼2012-05-14 13:47:26