和友人吃饭聊天,谈起一个米籍华人教授。友人说,他被判刑了,你不知道吗?这是 2 年前的事了!被判刑一年多,现在应该出来了。 为什么呀? 因为他用自己公司的名义申请了 NASA 的经费,说主要的研究在这个公司做。而实际上,他的公司根本没有做研究的设备,研究都是在大学做的,却以公司做研究的名义向 NASA 要钱,还不跟 NASA 说公司是他们两口子的。 另外,申请书说的是她老婆要指导研究生做实验。事实上他老婆一天都没去指导过。 他还说,这个教授“太抠”,给学生钱太少,发论文写报告只有他们两口子,没学生的事 …… NASA 以没有用 NASA 的经费做该做的研究为由,起诉他们诈骗,要让他们还钱。他老婆也被判了刑。 我知道这个教授非常有成就,在业界非常有名,太可惜了…… 【友人让我参见如下网络新闻】 http://www.mcall.com/news/breaking/mc-lehigh-professor-nasa-fraud-verdict-20151120-story.html http://www.mcall.com/news/breaking/mc-lehigh-professor-yujie-ding-nasa-fraud-sentencing-20160928-story.html 联想到, 我们自己在使用科研经费的时候, 是都用在申请书中指定的科研上了吗?昨天审计的还问:你派学生参加的这个会和项目有关吗?你参加的这个会和项目有关吗?绩效发给的这两个老师,是课题组成员吗?问得我有些心虚。虽然其中一部分经费是拨给别人的,里面的花费有些不合理,我也得担当,把账调出去。
JPL 行星精密历表的使用 上海航天技术研究院 李云飞 (blitheli@gmail.com) 当计算太阳系某行星相对太阳或者相对其他行星的位置和速度的时候,在精度不太高的情况下,可以用行星的平均根数解析的得出,但是如果需要计算其高精度的位置和速度的时候,通常采用 NASA( 美国航空航天局 ) 的 JPL( 喷气推进实验室 ) 给出的行星精密历表 (DE200,DE403,DE405) 。 在 JPL 的太阳系动力学网页 (http://ssd.jpl.nasa.gov/) 上有太阳系行星历表计算的一般精度和高精度计算方法和子程序。 本文以行星历表的高精度计算 (DE405) 为例,阐述如何从 JPL 的 FTP 上下载相关数据,如何验证子程序的正确性,以及最后如何使用子程序。 本示例中以 DE405 为例,下载的子程序为 FORTRAN 语言,对于其他的历表 ( 如 DE200) 和语言 (C++) ,过程和以下相同,只要将相应的代码对应修改即可。 详细的过程可参考 http://ssd.jpl.nasa.gov/?planet_eph_export ,作者的操作系统为 Windows XP ,对于 UNIX 操作系统,请阅读上行中的链接中的说明。由于下载数据后要对其验证子程序做一些修改,而读者在初次读其验证子程序时难免有些疑惑,因此下面直接给出步骤: 1. 登录 JPL 的 ftp(ftp://ssd.jpl.nasa.gov/), 下载以下数据: 一共 9 个文件: /pub/eph/planets/usrguide(需要仔细阅读) /pub/eph/planets/fortran/testeph.f /pub/eph/planets/fortran/selcon.f /pub/eph/planets/fortran/asc2eph.f /pub/eph/planets/fortran/binmerge.f /pub/eph/planets/fortran/binshort.f /pub/eph/planets/ascii/DExxx/ascpYYYY.XXX (多个文件) /pub/eph/planets/ascii/DExxx/header.XXX /pub/eph/planets/test-data/testpo.XXX 将上述文件放在同一文件夹下,为了方便下面描述,设文件夹名字为 D:\DE405 。 请仔细阅读 usrguide 文件,看不明白也行,只要按照以后步骤做即可。 上述文件中, ascpYYYY.XXX 系列文件就是我们所需要的精密历表的数值文件。其中 XXX 表示星历表编号,本示例中为 405 ; YYYY 表示年份,如果你只需要公元 2000 到公元 2060 年之间的行星位置数据,那么只下载 ascp2000.405 到 ascp2040.405 的文件即可。 2. 打开“开始”菜单,点击“运行”,进入窗口 DOS 命令行界面 , 并进入 DE405 文件夹下。 首先你的系统内要有 fortran 的编译器 ( 如 Compaq visual fortran 6.5 版 ) ;那么在命令行界面输入 : copy header.405+ascp2000.405+ascp2020.405+ascp2040.405 infile.405 此命令会产生一个 infile.405 的文件。继续在命令行界面输入 : DF asc2eph.f ( df 表示对 .f 文件使用 fortran 进行编译,生成 exe 文件。若在 linux 系统下,命令为 gfortran ;在 windows 下,为 ifort 。因为我是第一次用编译器来做,为 win7 下安装 vs2010+studio XE2011 。其对应的编译器的命令串口可以在开始菜单中找到,根据自己计算机的位数, 32 或 64 进行使用,如图 asc2eph infile.405 ( 此命令会产生一个名为 JPLEPH 的文件 ) 3. 打开 ”testeph.f” 文件,修改程序中以下部分: 在子程序 ”FSIZER3” 中,将 ”NRECL” 设为 4, 将 ”NAMFIL” 设为 'JPLEPH' 。将 ”KSIZE” 设为 2036 ,对于其它的历表 ( 如 DE200), ”KSIZE” 的值有所不同,参考程序中旁边的提示。 在子程序 ”STATE” 中,将语句 ”CALL FSIZER3(NRECL,KSIZE,NRFILE,NAMFIL)” 最前端的 ’C’ 去掉。关闭 ”testeph.f” 。 4. 继续在命令行界面中输入: DF testeph.f testeph testpo.405 此时界面中会输出一系列的常用参数,以及和标准数据的比较结果,其中对应 ”difference” 的位置的数据全是 ”0.****E-13”, 则表示 ”testeph.f” 程序修改的正确,否则则会出现 '***** WARNING : next difference = 1.D-13 *****' 的字样,那么表示程序修改错误,需要继续修改。如果出现错误,则尝试将步骤 3 中的 ”NRECL” 设为 1 重新试试。如果仍不行,那么可能由于系统的原因,请读者自己按照其英文的说明文档自己做,或者来信与笔者交流 ( blitheli@gmail.com ) 。 5. 打开 ”testeph.f” 文件,将子程序 ’FSIZER3’( 包括它 ) 以下所有的子程序全部拷贝到到 ’selcon.f’ 文件中。 6. 至此,精密历表的前期工作已经完成,读者所需要的计算行星位置速度的高精度星历表子程序全部在文件 ’selcon.f’ 中,读者在实际应用中可以直接调用 ’selcon.f’ 文件中的相关子程序 ( 相应的 ’JPLEPH’ 数据文件要存在 ) 。如果需要的精密历表的时间与上面的时间范围不同,读者可以自己转换所需要时间段的 ( 如公元 2020-2100 年间的精密历表 ) 星历表数值文件,具体方法同本文的步骤类似。 7. 关于 ’selcon.f’ 文件中的行星历表的子程序名称及其接口,请参考附录 A 。 附录 A 精密历表子程序 在最后的 ’selcon.f’ 文件中,包含着读者所需要的几个有关精密历表的子程序,下面给出各个子程序的简介及其接口参数说明。 : 给定某一儒略日时刻,得到目标行星相对某行星中心的位置和速度 : 获取所有的有关星历表的常数 : 获得读者所需的星历表某一常数 : 同子程序 PLEPH ,但输入的儒略日时刻更加精确 l PLEPH ( ET, NTARG, NCENT, RRD ) Input: ET : 儒略日 NTARG : 目标行星的编号 NCENT : 中心行星的编号 Output: RRD(6) : 目标行星相对中心行星的位置和速度 ; 单位 : for nutations, d(psi), d(eps), d(psi)-dot, d(eps)-dot ; for librations, (Euler angles and rates, w.r.t. the ephemeris reference frame) ---------------------------------------------------------------------------- 关于目标行星和中心行星的编号,它们的含义如下 : 1 = MERCURY 2 = VENUS 3 = EARTH 4 = MARS 5 = JUPITER 6 = SATURN 7 = URANUS 8 = NEPTUNE 9 = PLUTO 10 = MOON 11 = SUN 12 = SOLAR-SYSTEM BARYCENTER 13 = EARTH-MOON BARYCENTER 14 = NUTATIONS (LONGITUDE AND OBLIQ) 15 = LIBRATIONS, IF ON EPH FILE (IF NUTATIONS ARE WANTED, SET NTARG = 14. FOR LIBRATIONS, SET NTARG = 15. SET NCENT=0.) l DPLEPH ( ET2, NTARG, NCENT, RRD ) Input: ET2(2) : 通常此向量有以下三种用法: 1 最简单的用法,将儒略日赋值给 ET2(1), 将 ET2(1) 置 0 ; 2 为了最大的提高计算精度,将 ET2(1) 赋值为离所要计算时刻最近的午夜对应的儒略日,将剩下的儒略日赋值给 ET2(2) ; 3 为了方便,也可将 ET2(1) 赋值为某一固定时刻的儒略日 ( 如开始积分的时刻 ) ,将剩下的儒略日赋值给 ET2(2) ; 其他参数参考 PLEPH ( ET, NTARG, NCENT, RRD ) 。 l CONST(NAM,VAL,SSS,N) Output: NAM(N) : 常数名称向量 VAL(N) : 常数值向量 SSS(3) : sss(1) : 星历表数据对应的起始时刻儒略日 sss(2) : 星历表数据对应的终点时刻儒略日 sss(3) : Chebychev 系数每块数据所包含的天数 N : 常数向量的长度 l selcon(nams,nns,vals) Input: nams : 所要计算常数的名称 nns : 所要计算常数的编号 Output: vals : 所要计算常数的数值 ' -------------------------------------------------------------------------------- 部分常数名称 ( 按编号顺序 ) : DENUM Planetary ephemeris number. LENUM Lunar ephemeris number. TDATEF, TDATEB Dates of the Forward and Backward Integrations CLIGHT Speed of light (km/s). AU Number of kilometers per astronomical unit. EMRAT Earth-Moon mass ratio. GMi GM for ith planet . GMB GM for the Earth-Moon Barycenter . GMS Sun (= k**2) . X1, ..., ZD9 Initial conditions for the numerical integration, given at "JDEPOC", with respect to "CENTER". JDEPOC Epoch (JED) of initial conditions, normally JED 2440400.5. CENTER Reference center for the initial conditions. (Sun: 11, Solar System Barycenter: 12) RADi Radius of ith planet . *** 上面, 表示双精度类型数据 (double precision); 表示整型类型数据 ; 表示字符型数据。 *** 某一格林尼治时刻所对应的儒略日 ( 例如 2008-Jan-20 15:00:00 时刻对应的儒略日为 2454486.125) 的计算方法的子程序见我个人主页中的 ’ UTC 到 Julday 转 换 ’ 文件 ( http://blitheli.googlepages.com/programs ) ;也可用 JPL 网站上现成的程序来进行转换 ( http://ssd.jpl.nasa.gov/tc.cgi#top ) 。 附录 B JPL 精密历表分类 l DE200 : ( 包含章动,但不含岁差 ) 对应的历元起始时刻为: JED 2305424.5 (1599 DEC 09) to JED 2513360.5 (2169 MAR 31) 。 This ephemeris has been the basis of the Astronomical Almanac since 1984. It is based upon the dynamical equator and equinox of J2000 (see Standish, 1982 and Standish, 1990). l DE403 : ( 包含岁差和章动 ) 对应的历元起始时刻为: JED 2433264.5 (1949 DEC 14) to JED 2469808.5 (2050 JAN 02) Tied to the International Celestial Reference Frame through comparison of UTPM estimates from Lunar Laser Ranging and VLBI l DE405 : ( 包含岁差和章动 ) 对应的历元起始时刻为: JED 2305424.50 (1599 DEC 09) to JED 2525008.50 (2201 FEB 20) Tied to the International Celestial Reference Frame through VLBI observations of the Magellan spacecraft in orbit around Venus. DE405 was created in May-June 1997. l DE406 : 长期星历表 ( 不含岁差章动 ) 对应的历元起始时刻为: JED 0624976.50 (-3001 FEB 04) to 2816912.50 (+3000 MAY 06) This is the same ephemeris as DE405, though the accuracy of the interpolating polynomials has been lessened (interpolation on the 64-day mesh points remains exact, however). For DE406/LE406, the interpolating accuracy is no worse than 25 meters for any planet and no worse than 1 meter for the moon. DE406 requires about 10 megabytes for each 300-year block. l DE410 : Specialzed ephemeris used for Mars Exploration Rover navigation (DE409 is identical except for the reference vale uf GM for Mars) Tied to ICRF by VLBI observations of Mars Global Surveyor and Mars Odyssey. Created 24 April 2003 l DE413 : A special ephemeris to update the orbit of Pluto to aid in planning for an accolutation of a relatively bright start by Pluto's satellite Charon on 11 July 2005. Created 4 November 2004 l DE414 : An ephemeris fit to ranging data from MGS and Odyssey through 2003 in addition to many other data types for all planets. Integration covers 1599 to 2201. Some aspects are documented in Alex S Konopliv et al., Icarys vol 182, pp 23-50 (2006). Created May 2005
权利人真需要保留所有权利吗 ? 著作权法保护作者的两类权利:人身权和财产权。其中人身权包括发表权、署名权、修改权、保护作品完整权等精神权利,财产权则主要是以复制权为核心的从作品中获得经济利益的权利。在开放存取出现以前,唯一的版权声明方式是保留所有权利( All Rights Reserved ),即未经著作权人同意,使用者不能随意对作品进行传播、复制、下载、打印等处理。这种对任何作品都声明保留所有权利,已经被滥用到任何人都可能触犯版权法的地步,而且使得很多优秀作品得不到广泛传播,无法实现效益最大化。作者发表论文,尤其是科技学术论文,很少出于考虑为获得直接的经济报酬,目的主要还是在于学术交流,让更多同行了解和认可,共同发展本学科领域的知识。因此,现在越来越多的学者开始意识到自己并非需要保留所有权利,相反他们更愿意选择保留部分权利( Some Rights Reserved )或不保留权利( No Rights Reserved )。在 OA 出版模式下,作者选择的是保留对作品的部分权利,同时通过授权许可协议(最常用的是创作共享许可协议)对使用者尽可能的开放其他权利,一般而言作者保留的是精神权利,放弃的是财产权利。 笔者主张和呼吁,为扩大作品在网络环境下的传播效益最大化,版权人应舍得放弃,采用保留部分权利的版权声明,适当让渡和放弃部分权利,于人于己均好。