版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

liming9398

金虫 (正式写手)

应助: 1 (幼儿园)
金币: 780.9
散金: 89
红花: 1
帖子: 361
在线: 1009.2小时
虫号: 610070
注册: 2008-09-23
专业: 代数学

[求助] 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人回复
2025冷门绝学什么时候出结果已经有6人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有7人回复
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗已经有8人回复
寻求一种能扛住强氧化性腐蚀性的容器密封件已经有5人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有4人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复

» 本主题相关价值贴推荐，对您同样有帮助:

求 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人回复

1楼2012-05-14 11:05:53

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

hklulu

木虫 (正式写手)

000

应助: 18 (小学生)
金币: 2645
散金: 61
红花: 3
帖子: 341
在线: 154.9小时
虫号: 108850
注册: 2005-11-18
专业: 原子和分子物理

【答案】应助回帖

感谢参与，应助指数 +1

哈哈，新手吧。

八位不是指的数字的位数，而是二进制的8B，就是通常用64位二进制数字来表示双精度数字。

此位非彼位。

赞一下

回复此楼

2楼2012-05-14 11:26:38

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

liming9398

金虫 (正式写手)

应助: 1 (幼儿园)
金币: 780.9
散金: 89
红花: 1
帖子: 361
在线: 1009.2小时
虫号: 610070
注册: 2008-09-23
专业: 代数学

可能我没说清楚3.14159265不就是九位的吗，可RAD是real类型的也就是单精度类型的，最多七位有效数字啊，赋值时写这么多位还有意义吗，同样的问题出现在给real类型数据as赋值时表达式中出现29.105356，也超过7位有效数字啊

赞一下

回复此楼

3楼2012-05-14 11:51:31

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

hklulu

木虫 (正式写手)

000

应助: 18 (小学生)
金币: 2645
散金: 61
红花: 3
帖子: 341
在线: 154.9小时
虫号: 108850
注册: 2005-11-18
专业: 原子和分子物理

【答案】应助回帖

奥，当然有意义：
32位机的单双精度的数据：单精度的有效数字：6位；双精度的有效数字：15位
64位机器上有最大18位有效数字和4932为指数的类型。

有些计算需要有效位数1000多位，比如Pi的计算等，这时候就需要自己编写程序来控制有效位数的输出。我记着10年前，我为了在PC上输出20位有效位数，亲自编写了一个多有效位（名为doubledoubledboule类型）的程序，不过现在网上多的很，楼主可以自己找一找，如果感兴趣的话。

赞一下

回复此楼

4楼2012-05-14 13:47:26

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

Carpathia

银虫 (小有名气)

应助: 35 (小学生)
金币: 5.6
散金: 12
红花: 2
帖子: 140
在线: 24.9小时
虫号: 1070633
注册: 2010-08-06
性别: GG
专业: 动力学与控制

【答案】应助回帖

感谢参与，应助指数 +1

建议ls看清楚lz的问题
lz说得没错，给单精度赋值超过位数是没有意义的，同样给双精度数赋值也要加上D0才有意义。Fortran默认的是单精度数。

赞一下

回复此楼

追求理想的道路上东张西望

5楼2012-05-15 15:02:48

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

liming9398

金虫 (正式写手)

应助: 1 (幼儿园)
金币: 780.9
散金: 89
红花: 1
帖子: 361
在线: 1009.2小时
虫号: 610070
注册: 2008-09-23
专业: 代数学

引用回帖:

4楼: Originally posted by hklulu at 2012-05-14 13:47:26:
奥，当然有意义：
32位机的单双精度的数据：单精度的有效数字：6位；双精度的有效数字：15位
64位机器上有最大18位有效数字和4932为指数的类型。

有些计算需要有效位数1000多位，比如Pi的计算等，这时候就需 ...

不太明白你的回复，你解释一下吗？我只是觉得这是本有名的书，似乎不应该犯这种低级错误，或许是作者另有初衷，我没理解到吗？

赞一下

回复此楼

6楼2012-05-15 15:16:59

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖