;+ ;NAME: ; querynednh ;PURPOSE: ; Query the HI column density of the LAB Survey for given position ;CALLING SEQUENCE: ; result=querynednh(ra,dec) ;EXAMPLE: ; IDL ra=2.018361*15 ; IDL dec=0.56667 ; IDL result=querynednh(ra,dec) ; IDL print,result ; 2.1600001e+20 ; ; IDL result=querynednh(ra,dec,/two) ; IDL print,result ; 2.1600001 ; ;PROCESS: ; webget( ),Euler ;INPUT: ; ra ---- right Longitude in degrees (sclar or veltor) ; dec ---- Latitude in degrees (sclar or veltor) ;OPTIONAL KEYWORD INPUT: ; /nearest ----- if set it, will give the HI column density of the ; nearest position of the selected position. ; default it, will give the HI column density of the ; selected position. ; /two ---- if set it, will give the HI column density in ; Unit :10^20 cm^-2 ; default it, will give the HI column density in cm^-2 ;OUTPUT: ; nh ---- the HI column density of line-of-sight ;METHOD: ; see: ; http://www.astro.uni-bonn.de/%7Ewebrai/english/tools_labsurvey.php ; by integrating over the velocity range from -400 to +400 km/s in ; the LSR frame. ; ; Kalberla,P.M.W.,Burton,W.B.,Hartmann,Dap,Arnal,E.M.,Bajaja,E., ; Morras, Poppel,W.G.L.(2005),AA,440,775 ; ; Bajaja,E.,Arnal, E.M., Larrarte,J.J., Morras, R., Poppel, ; W.G.L., Kalberla, P.M.W. 2005, AA, 440, 767 ; ; Arnal, E. M., Bajaja, E., Larrarte, J. J., Morras, R., ; Poppel, W. G. L. 2000, AAS, 142, 35 ; ; Hartmann Burton 1997, Cambridge University Press, ISBN ; 0521471117 ; ;REVISIONHISTORY: ; ORiginal by DL.Wang,26-Jun-2007 ; chage the weblink as follow: ; http://www.astro.uni-bonn.de/~webaiub/english/tools_labsearch.php ; by DL.Wang,Sep-22-2008 ; add !VERSION.OS check by DL.Wang,Oct-27-2011 ;- ;++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ Function querynednh_lab, ra, dec, nearest=nearest, two=two if N_params() lt 2 then begin print,'Syntax - info = querynednh(ra, dec)' print,' RA (degrees), Dec (degrees) -- search coordinates' endif if !VERSION.OS_family eq 'Windows' then stop euler,ra,dec,glon,glat,select=1 glon0=string(glon) glat0=string(glat) QueryURL = 'http://www.astro.uni-bonn.de/~webaiub/english/tools_labsearch.php?alpha='+strtrim(glon0,2)+'beta='+strtrim(glat0,2) result=webget(QueryURL) na=n_elements(result.text) ear=strarr(na) for i=0L,na-1 do begin ear =result.text print,i,' ',ear endfor t=result.text if not keyword_set(nearest) then begin t =strtrim(t ,1) nh0=strmid(t ,18,5) nh1=strmid(t ,38,2) nh=float(nh0)*10.0D0^float(nh1) endif else begin t =strtrim(t ,1) nh0=strmid(t ,18,5) nh1=strmid(t ,38,2) nh=float(nh0)*10.0D0^float(nh1) endelse if not keyword_set(two) then begin nh=nh endif else begin nh=nh*1.0D-20 endelse return,nh End ;++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ Function querynednh,ra, dec, nearest=nearest, two=two n=n_elements(ra) nh=dblarr(n) if not keyword_set(nearest) and not keyword_set(two) then begin for i=0L,n-1 do begin nh =querynednh_lab(ra ,dec ) endfor endif if not keyword_set(nearest) and keyword_set(two) then begin for i=0L,n-1 do begin nh =querynednh_lab(ra ,dec ,/two) endfor endif if keyword_set(nearest) and not keyword_set(two) then begin for i=0L,n-1 do begin nh =querynednh_lab(ra ,dec ,/nearest) endfor endif if keyword_set(nearest) and keyword_set(two) then begin for i=0L,n-1 do begin nh =querynednh_lab(ra ,dec ,/nearest,/two) endfor endif return,nh End
Pro nvss_read,filename, ra, dec, dist, flux, major, minor, pa, res, $ p_flux, p_ang, field, xpix, ypix, $ nvss, not_nvss, type_trans=type_trans ;+ ;NAME: ; nvss_read ;PURPOSE: ; read data from outfile by NVSSlist ;CLING SEQUENCE: ; nvss_read,filename, ra, dec, dist, flux, major, minor, pa, res, $ ; p_flux, p_ang, field, xpix, ypix, $ ; nvss, not_nvss ;CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC ;CC ;CC Waring: ;CC ;CC This code is used under the condition that one source have ;CC only one matched source in the NVSS catalog ;CC ;CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC ;INPUT: ; filename -- name of NVSSlist outfile ;OUTPUT: ; ra --- RA in h:m:s ; dec --- DEC in d:m:s ; dist --- distance between matched source NVSS source ; in Unit: arcsec ; flux --- flux in Unit: mJy ; major --- ; minor --- ; pa --- ; res --- Residual (Res) code ; P* = high peak, ; R* = high RMS, ; S* = high integral ; p_flux --- ; p_ang --- ; field --- ; xpix --- ; ypix --- ; nvss --- index of detected source ; not_nvss --- index of not detected source ; ;OPTIONAL KEYWORD OUTPUT: ; type_trans --- if set it,dist,flux,p_flux will be float, ; default it,they will be string ; ;EXAMPLE: ; IDL filename='rgb_nvss.dat' ; IDL nvss_read,filename, ra, dec, dist, flux, major, minor, pa, res, ; p_flux, p_ang, field, xpix, ypix, nvss, not_nvss ; ;PROCESS: ; extract_multarr ;REVISION HISTORY: ; Original by DL.Wang,Aug-29-2007 ; ;- nlines = NUMLINES( filename ) if nlines LT 0 then return ;########################## ;## definite arrary keep=lonarr(nlines) index0=lonarr(nlines) radec=strarr(nlines) ra=strarr(nlines) dec=strarr(nlines) dist=strarr(nlines) flux=strarr(nlines) major=strarr(nlines) minor=strarr(nlines) pa=strarr(nlines) res=strarr(nlines) p_flux=strarr(nlines) p_ang=strarr(nlines) field=strarr(nlines) xpix=strarr(nlines) ypix=strarr(nlines) temp = '' openr,lun,filename,/get_lun for j = 0L, nlines-1 do begin readf, lun, temp ;########################################### ;########### blank row pos=strpos(temp,'') if pos eq -1L then begin ; length=strlen(temp) ; if length eq 0L then begin index0 =0L keep =0L ra ='' dec ='' dist ='' flux ='' major ='' minor ='' pa ='' res ='' p_flux ='' p_ang ='' field ='' xpix ='' ypix ='' endif ;########################################### ;######### Source row unpos=strpos(temp,'SOURCE NOT FOUND') if unpos ne -1L then begin index0 =1L keep =0L ra ='' dec ='' dist ='' flux ='' major ='' minor ='' pa ='' res ='' p_flux ='' p_ang ='' field ='' xpix ='' ypix ='' endif rapos=strpos(temp,'|') if rapos ne -1 then begin dat=strsplit(temp,'|',/extract) if n_elements(dat) eq 13 then begin index0 =1L keep =1L ra =dat dec =dat dist =dat flux =dat major =dat minor =dat pa =dat res =dat p_flux =dat p_ang =dat field =dat xpix =dat ypix =dat ;print,ra ;pause endif if n_elements(dat) ne 13 then begin index0 =0L keep =0L ra ='' dec ='' dist ='' flux ='' major ='' minor ='' pa ='' res ='' p_flux ='' p_ang ='' field ='' xpix ='' ypix ='' endif endif endfor free_lun,lun ;###################################################### ;######### pick out needed data nvssindex=where(index0 eq 1L and ra ne ' RA_2000 ' and ra ne ' h m s ',cp) print,'match',cp index0=index0 keep=keep ra=ra dec=dec dist=dist flux=flux major=major minor=minor pa=pa res=res p_flux=p_flux p_ang=p_ang field=field xpix=xpix ypix=ypix nvss=where(index0 eq 1L and keep eq 1L and ra ne '',count) print,'NVSS',count not_nvss=where(index0 eq 1L and keep eq 0L,count1) print,'not NVSS',count1 ra=ra dec=dec dist=dist flux=flux major=major minor=minor pa=pa res=res p_flux=p_flux p_ang=p_ang field=field xpix=xpix ypix=ypix ;##################################################### ;######## data type transformation if keyword_set(type_trans) then begin ;removes leading and trailing blank from the input String for i=0L,n_elements(ra)-1 do begin dist =strtrim(dist ,2) flux =strtrim(flux ,2) p_flux =strtrim(p_flux ,2) endfor dist=float(dist) flux=float(flux) p_flux=float(p_flux) endif ;stop End
Pro nvss_upfile,ra,dec,radius=radius,filename=filename,format=format ;+ ;NAME: ; nvss_upfile ;PURPOSE: ; write a upfile for NVSS matched source ;CALLING SEQUENCE: ; nvss_upfile,ra,dec,filename=filename ; ;INPUT: ; ra ----- RA for source in Unit:degree ; dec ----- DEC for source in Unit:degree ;OPTIONAL KEYWORD INPUT: ; radius ---- matched radius in Unit:arcsec ; default it,it will set it 15 arcsec ; format ---- if set it,it will output according to format ;OPTIONAL KEYWORD OUT: ; filename ---- upfile name ;EXAMPLE: ; IDL nvss_upfile,ra,dec,filename='NVSS_update.dat' ; ;REVISION HISTORY: ; Original by DL.Wang,Aug-30-2007 ;- if ( N_PARAMS() lt 2 ) then begin message,'Syntax:nvss_upfile,ra,dec, ' return endif if not keyword_set(radius) then radius=15 openw,lun,filename,/get_lun n=n_elements(ra) for i=0L,n-1 do begin if not keyword_set(format) then begin printf,lun,adstring(ra ,dec ,1),' ',radius,' ','0' endif else begin printf,lun,adstring(ra ,dec ,1),radius,'0',format=format endelse endfor free_lun,lun End
Pro write_2mass_upfile,ra,dec,radius=radius,name=name, $ filename=filename,format=format ;+ ;NAME: ; write_2mass_upfile ;PURPOSE: ; write a upfile for matched source list,matched with 2MASS ;CALLING SEQUENCE: ; write_2mass_upfile,ra,dec,radius=radius,filename=filename,format=format, ; /name ;INPUT: ; RA ----- RA for matched source list in Unit:degree ; DEC ----- DEC for matched source list in Unit:degree ;OPTIONAL KEYWORD INPUT: ; radius ----- match radius in Unit:arcsec ; default it,radius is 7.0 arcsec ; name ----- if set it,it will print source name. ; default it,it will print source sequence ;OPTIONAL KEYWORD OUTPUT: ; filename ----- name for output file ; default it,filename is '2mass_update.dat' ; format ----- format for output data ; ;EXAMPLE: ; IDL write_2mass_upfile,ra,dec,filename='rgb_2mass_update.dat' ; ;REVISION HISTORY: ; Original by DL,Wang,Aug-17-2007 ;- Npar = N_params() if ( Npar lt 2 ) then message, $ 'ERROR - RA and Declination do not have equal number of elements' if not keyword_set(radius) then radius=7.0 if not keyword_set(format) then format='(1x,A14,1x,F10.6,1x,F10.6,1x,F3.1)' if not keyword_set(filename) then filename='2mass_update.dat' openw,lun,filename,/get_lun printf,lun,'\ Example of cone search' printf,lun,"\EQUINOX = 'J2000.0'" printf,lun,'| objname | ra | dec | radius |' printf,lun,'| string | double | double | double |' printf,lun,'| unit | unit | unit | arcsec |' if not keyword_set(name) then begin for i=0L,n_elements(ra)-1 do begin printf,lun,i+1,ra ,dec ,radius,F=format endfor endif else begin for i=0L,n_elements(ra)-1 do begin printf,lun,name ,ra ,dec ,radius,F=format endfor endelse free_lun,lun End
SAGE 的研发至少通过两个公司对后来的电子商务发展产生了深远的影响，一个是我们所熟悉的 IBM ，另一个是我们不太熟悉的 DEC 。 IBM 走的是大而全的发展方向，而 DEC 则是小而精。前者为美国的众多企业尤其是服务业全面信息化提供了有效但是昂贵的解决方案。而后者则以交互式运算为主要技术特色，通过降低价位和用户共同开发等策略，开创了个人电脑的发展方向，也为早期互联网的发展提供了硬件支持。 从两个公司的历史来看， DEC 的创立是因为晶体管新技术的出现使得用其组装计算机成为可能，所以 DEC 想做的，就是用新技术做新计算机，没有历史包裹需要背负。而 IBM 则是从一个行业（读卡机）涉足另一个行业（计算机），并且已经在所涉足的行业里有了以上一代技术（电子管）开发的产品和客户，所以对 IBM 来说，晶体管的出现是技术的一个升级，而是否完全利用这个新技术取决于很多现有客户市场和产品更新换代的因素。这个区别使得两个公司在这里历史新起点上走了不同的发展道路。 我们先来看看 IBM 。 创建于 1911 年的 IBM 是个独特的公司。美国历史上能够延续超过百年的公司屈指可数，它算一家。尽管它的业务从读卡机到计算机一直是当时时代的高科技，它的最初两代领导人却没有什么高科技背景。老沃森是个做推销员起家的商人，他的儿子小沃森是个富二代，喜欢开飞机。但历史给了 IBM 很多得天独厚的机会。首先是美国政府在十九世纪末由于人口增长导致的管理方面的技术革新需要，使得读卡机（从提高效率来看，相当于那个时代的计算机）的生意蒸蒸日上。老沃森靠此发家。然后，由于两次世界大战和美国的战后技术领先地位，老沃森又通过跟政府打交道建立起的人脉把读卡机把生意扩展到全球（尤其是战后的欧洲）。当小沃森 1956 年接过 CEO 的位置时， IBM 已经在读卡机和计算机领域里成为无可非议的老大。而小沃森接任之后，正赶上前面所提到的美国的企业信息革命，大型机系统在商业领域的应用成为当时的市场焦点。这使得小沃森有了跟老沃森在 50 多年前几乎同样的历史机遇，于是 IBM 再上一个台阶。继续引领计算机产业的龙头老大地位，继续走开发政府和军工项目（降低自我研发成本），然后将成熟技术扩散到民企的路子。整个 50 到 60 年代， IBM 几乎介入了美国所有大型军工项目，从中赚足了钱。这一点从 SAGE 上就可以看出来，仅仅这一个项目就占 IBM 整个公司五十年代利润的 10% 。 IBM 的另一个特色是不只卖产品，而是卖服务。这是从老沃森时代开始留下来的传统。读卡机的制造成本相对昂贵，但对企业来说，从长远看，买还是比租要划算。可以推销员出身的老沃森发现，如果光卖机器本身，对企业来说需要在短时间内支付一大笔钱，对 IBM 来说失去了未来的源源不断的现金流，是个一锤子买卖。另外，读卡机的更新换代也比较快，新的型号出来了，企业可能不愿意再付一大笔钱更新换代。所以还是租赁比较符合当时的双方需求。于是 IBM 就坚持租赁的形式。既然是租赁，就有客服，也就有了提供全套服务，也就是完整解决方案的需求。这既给 IBM 不断的新挑战，也给 IBM 发现新产品和市场，和锁定用户的能力。所以这个策略可谓一石二鸟。 IBM 的租赁策略一直到 80 年代初都屡试不爽，不过其领导层渐渐志满意得，偏离了老沃森策略的深意，认为大型机是其核心竞争力，而忽略了用户真正需要的是服务。这一点直到 90 年代初郭士纳接管 IBM 时才发现，于是 IBM 后来重回正轨，这是后话。 回头再看 DEC 。 创立于 1957 年的 DEC 从创立到产品开发策略再到最后的归宿，跟现在的很多创业公司非常类似。它的两个创始人奥尔森和安德森都在麻省理工的林肯实验室工作，曾经参与开发了 SAGE 的交互式雷达监控系统以及在这之前空军需要的模拟飞行系统。我们前面讲过，雷达信号被转换成数字信号，传送到指挥中心时， AN/FSQ-7 计算机将其显示在雷达监控屏幕上，而指挥人员一边监控屏幕，一边用光笔点触代表敌机的信号，然后再通过一系列屏幕互动指挥拦截。这是一个典型的交互式系统。而当时计算机的主流并不是交互式系统，因为交互式系统将大量的机时耗费在与人互动的等待过程中。这对每月花几十万美元租赁大型机的公司来说，得不偿失。所以当时往往是人一排排站在计算机旁边，等待递交计算任务。可是交互式系统对非商业领域却有很大的需求。比如当时的麻省理工很多学生宁愿使用比较慢的交互式系统来处理他们的数据而不用更快的 IBM 的计算机（其实这是所有面向个人的计算技术的共性，也是苹果公司成功的原因）。这些经验和观察都给了 DEC 的两个创立人足够的信心来开发交互式计算机。 可是两人很快发现，当时的投资人不太愿意投资给小计算机公司。因为这些小公司往往不久就垮掉。于是他们不得不修改他们的创业计划，把成立交互式计算机公司改成成立交互式计算机“部件”公司。这样投资规模小，并且既可以日后生产交互式计算机也可以首先为其他计算机公司提供晶体管部件。于是他们以 70% 的股份换取了七万美元的投资，成立了 DEC 。 投资规模和发展方向都决定了 DEC 跟 IBM 走了两条不同的路。因为投资规模小，他们只生产以晶体管为基础计算机核心部件，把外围的东西留给用户自己解决（就好比只卖主板，声卡，显卡，显示器等等都让用户自己去解决）。而且他们只卖不租。通过提供详细的技术手册的方式，让用户有最大限度的，按照自己的需求，来利用 DEC 提供的计算机。这种模式适合的是对计算机技术有相当了解的用户，而大学和研发机构恰恰是这种用户。这一产品策略衍生出了后来 70 年代的个人计算机市场发展模式以及黑客文化。 DEC 专注的是以晶体管为基础部件的交互式计算机市场。其运算速度虽然比 IBM 的大型机慢，但对用户来说不是问题，所以市场是有的，关键是如何扩大这个市场。而要扩大市场，在价格上要有竞争力。要在价格上有竞争力，就要在设计上省钱。于是在设计上（或者说计算机架构上）， DEC 选择了用软件替代部分硬件功能的办法，这样牺牲了一部分效率但降低了成本，这个选择导致大型机与小型机以及后来个人计算机设计上的一个核心区别。 这里的替代主要是在硬件设计的字节长度方面。 IBM 的大型机字节一般至少是 32 位，而 DEC 的最初小型机字节只有 12 位（后来的英特尔系列处理器字节稍长一些，固定为 16 位）。字节贯穿计算机信息处理的几个关键步骤，包括数据的存储，计算和寻址。长的字节跟长的带宽类似，可以在相同时间内更快速地处理数据，而短的字节在这些方面则需要通过软件方法进行转换。简单的讲，一个 9 位数的整数，如果用 32 位字节的计算机来存储，就可以直接从一个地址上读取；而如果用 16 位字节的计算机来读取，因为字节长度不够，就需要至少两次从不同地址读取。这些就导致了软件复杂度增加，而且运算速度也相应减慢。不过好处是硬件费用降低，而速度的减慢比例对交互式运算的用户来说不是问题。在 50 年代末和 60 年代初，硬件部分，尤其是中央处理器和磁芯内存占了计算机生产成本的大头。所以 DEC 这么做，可以将一台计算机的成本控制在一万美元以内。另外，他们把使用晶体管组装计算机的过程标准化，然后通过机器自动装配，这样在提高效率的同时，减少了人工操作的误差。还有很重要的一点就是， DEC 将 IBM 大型机必备的数据处理通道（相当于第二个 CPU ，主要是管理数据的输入输出）拿掉，让内存可以与外围设备直接交换数据（我们现在称其为 DMA ，直接内存访问）。而一个数据处理通道的成本就相当于 DEC 后来一台计算机的价格。就这样， DEC 设计的交互式计算机成功的降低了成本。该公司在推出的第一版计算机 PDP-1 售价是 12 万美元。等到 1965 年推出后来最流行的 PDP-8 计算机时，销售价格不到 1 万美元 。值得一提的是， DEC 的交互式小型机的首批用户除了大学和研究机构之外，还有地处麻省剑桥的 BBN 咨询公司，该公司在六十年代研发了互联网前身阿帕网。 PDP-8小型机
http://www.nature.com/news/index.html US Congress passes strategic science bill Nature News 22 Dec 2010 African elephants are two distinct species Nature News 21 Dec 2010 Love thy lab neighbour Nature News 21 Dec 2010 Integrity policy unveiled at last Nature News 20 Dec 2010 UK science faces facilities freeze Nature News 20 Dec 2010 comments No evidence of time before Big Bang 54 comments 10 December 2010 Arsenic-eating microbe may redefine chemistry of life 33 comments 02 December 2010 Women scientists must speak out 30 comments 08 December 2010 Telomerase reverses ageing process 20 comments 28 November 2010 Optical wing generates lift from light 15 comments 05 December 2010
标签: Dec
