| 查看: 1656 | 回复: 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赋值？单精度不是最多七位吗，多了有意义吗？ |
回复此楼» 猜你喜欢
求助一下有机合成大神
已经有3人回复
-
拟解决的关键科学问题还要不要写
已经有8人回复
-
最失望的一年
已经有13人回复
-
存款400万可以在学校里躺平吗
已经有30人回复
-
求推荐英文EI期刊
已经有5人回复
-
请教限项目规定
已经有4人回复
-
国自然申请面上模板最新2026版出了吗？
已经有20人回复
-
26申博
已经有3人回复
-
基金委咋了？2026年的指南还没有出来？
已经有10人回复
-
基金申报
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求 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人回复
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
5楼2012-05-15 15:02:48
liming9398
金虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 780.9
- 散金: 89
- 红花: 1
- 帖子: 361
- 在线: 1009.2小时
- 虫号: 610070
- 注册: 2008-09-23
- 专业: 代数学
6楼2012-05-15 15:16:59
jackplay
金虫 (小有名气)
- 应助: 4 (幼儿园)
- 金币: 1120.8
- 红花: 1
- 帖子: 128
- 在线: 87.6小时
- 虫号: 2945895
- 注册: 2014-01-24
- 性别: GG
- 专业: 分析化学
7楼2016-01-09 08:21:30
43fd6ys
木虫 (正式写手)
- 应助: 37 (小学生)
- 金币: 2050.4
- 红花: 7
- 帖子: 379
- 在线: 147小时
- 虫号: 3935320
- 注册: 2015-06-22
- 专业: 零件成形制造
8楼2016-01-09 21:20:09