悼念

国内是地震悼念日的第二天了,今天一个朋友约了些在这里的中国人,在他家后院举行了一个很小的悼念会。大家叠了一些纸鹤画了一些家电烧了,希望不幸离开的人在另一个世界可以一路走好。

回到家,见LP在线。说到我退税拿回来一些钱,她说那再捐点吧。我对上会的中国红十字基金会是彻底失望了(见:搞不懂的红十字会基金会),给他们写信反映可能存在的问题,已经过去好几天,没有任何回复。而且最近中国红十字会闹出挺不让人开心的事情。前面朋友给介绍过乐施会,连接到乐施会美国,又捐了一点钱。希望这点点钱真能到达需要的人手上。

Difference between Radar precipitation data stage 2 and stage 4

Zhuotong Nan ([email protected])

Stage II

Be produced directly from the real-time gauge and WSR-88D radar data received at NCEP.

raingauges:

  • from ~1450 ASOS sites from near-real time METAR file
  • and from ~5500 HADS automated gauge reports via the GEOS data Collection Platform.

Radar:

  • Hourly real-time Digital Precipitation Arrays (DPAs) from 140 ConUS WSR-88D radar sites.

The primary input to the ETA/EDAS precipitation assimilation before 2003.

Stage IV

A mosaic of regional multi-sensor analysis produced by NWS River Forecast Centers (RFCs) and benefits from the RFCs’ manual quality control step

mosaicked from the regional multi-sensor 1h and 6h analyses produced by the 12 ConUS RFCs.

The primary input to the ETA/EDAS precipitation assimilation after 2003. Stage II is only used when Stage IV is unavailable.

In summary,

Since Stage IV benefits from manual QC at the local RFCs, it is generally better than Stage II, which has no manual QC. There are a number of other issues for your consideration:

1. Coverage period: Stage II began to be archived at NCAR on 1 May 1996. Stage IV did not began to be archived until 1 Jan 2002. (NOTE: we have data from 1997 on for OHIO river basin. I guess this statement is for entire continuous continental, rather than for an individual RFC. -Zhuotong)

2. Do you need hourly or 6-hourly analysis? For 6-hourlies, Stage IV is the one to use (if you don’t need the pre-2002 data). Stage IV hourlies have some coverage problem since we do not receive the local analyses from all the RFCs for all hours.

3. How soon do you need the analysis to be available? Stage II hourly analysis is made at 35 minutes after the top of the hour, then re-made twice at 6 hours and 18 hours later. Stage VI runs each hour at approximately 35 minutes past the top of the hour.

4. If I choose Stage II, which type should I use? We recommend the multi-sensor analysis with the exception of early spring – early Aug 2000, when the radar coverage was bad. You might want to use the gauge-only analysis for that time period.

Validation examples:

clip_image004 clip_image006

figure 1 24 hour precipitation accumulative amount on OCT 10, 2004, Stage II (left) and Stage IV (right).

clip_image008clip_image010

figure 2 Left: Stage II; Right: Stage VI. On Jun 28 2005. Here we can see the variance between them.

References:

1. https://nanzt.info/wp-content/uploads/2008/05/stage2-4.19hydro.pdf
2. http://www.emc.ncep.noaa.gov/mmb/ylin/pcpanl/QandA/
3. http://www.emc.ncep.noaa.gov/mmb/ylin/pcpverif/daily/

患难见真情

媳妇看了我的BLOG说,怎么才捐了这么一点,比你妈还少。她自然得意了一下,说她在单位捐了300。不过我奇怪的是,我妈居然也要捐款?记忆中父母是十分省的人,比如从温州到兰州,他们都是坐汽车30几个小时过来的,舍不得吃舍不得穿,对他们自己已经到了一种吝啬的地步。不过对儿子很大方,每次回家回程的机票都是他们给买的,从机场到家里还有一些路,也总是让打车回家或者找人开车去接我们,我们呆在老家的日子总是每天做大鱼大肉,我念大学的时候,每年寒假回去都能胖上10斤左右。我记忆中,他们对别人也不是那种省得花钱的人。

媳妇看到我的惊讶,说,你居然不知道?你妈说她最乐意干的事就是捐款了,她在家里也有帮助别人,只是没有太多的钱可捐。——我还真不知道,十来年长期在外,对自己父母居然这样不了解。

我妈昨晚看了电视,看到血淋淋的电视镜头,拿出她自己的钱,一定让媳妇要将这钱汇到四川给受灾的人。看她认真的样子,媳妇说,今天要从银行汇,还一定要写上你妈的名字。

灾难是惨重的,我们在悲痛的同时,也看到很多美好的画面:温总理忙碌在第一线,可爱的子弟兵不辞劳苦抢救灾民,单位、个人、海国外捐款捐物,志愿者一波波奔赴灾区,……

平时大家也许对生活现状都不满之外,但所谓爱之深、责之切,患难才见真情!

搞不懂的红十字会基金会

给四川地震救灾捐点钱。前面看到公布出来的中国红十字会的帐号,发现外币开户名都是中文的,想想是无法在国外给这种帐户转钱的,人家又不认中文,我又不知道他的准确英文名是什么,而且也没有更快捷的渠道可以直接从信用卡上划钱。据说可以从它的网站在线捐款,但问题是我居然没法子打开它的网站,是不是人一下子太多给挂了。

从文字城上看到红十字会基金会可以海外在线划款,进入 http://www.crcf.org.cn/donationol/donation.asp,很快填了50USD,看起来好像是委托第三方,要收3.2%的费用,OK,可以接受。一切很顺利。通过得到的捐赠序号,看到,

捐助序列号

捐助者名称

捐赠类别

币种

金额(元)

折合人民币(¥)

是否到帐

DCS0805140xxxxx

Zhuotong Nan

抗震救灾

USD美元

50

413.87

在途

貌似都顺利。到银行帐户上一看,发现

Effective Date

Description

Amount

05/13/2008

VISA – www.crcf.org.cn

$7.42

怎么回事呀?一算这个数字大概相当于人民币50块,怎么前面的美金变成人民币了,郁闷。想着可能是那个基金会的网站与第三方划款平台间的连接有问题,那头的美金设定没有传递到划款平台。

朋友在旁边看了说,你小心点,不要被骗了。我这才发现原来红十字基金会跟中国红十字会不是一回事,好在后来在红十字会官方网站上看到红十字基金会的链接,想来不是骗人,再看到李亚鹏的基金帮过不少人,想来不至于被恶意网站欺骗。不过朋友又说了,你看这些人,连个网站都做的要不连不上,要不连划个款也不对,这些人其实根本不在乎你们捐的这些钱,多点少点根本不管。

想想好像也有点道理,很让人寒心,不禁郁闷。联想到我大学时某天冲动,非拉着一个好朋友去南关什字的献血站去献血,发现门口苍蝇到处飞(堆了医疗垃圾),热情一下子被泼了冷水,吓的跟朋友赶紧就跑回来,从此没再去献过血。

朋友说他要晚些时候捐,也许是教会的,否则他不放心,错了不说,也许捐的钱还到不了需要这笔钱的人手上,被中间截住花掉了,这种新闻不是没有。

很无奈不是,我只有自我安慰,总是对得起自己了。

汶川地震

看到报道,温总理跑四川灾区最前线亲自指挥,真是很感动,看到为国为民的好总理的形象。08年好事不少,奥运会要办了,但好像不好的事情也不少,天灾人祸,又是雪灾又是地震还有列车事故,每每都看到60多岁的温总理忙碌的身影。感动之余,不免感叹除了总理外的大小官员们都到哪去了?我宁愿看到总理坐阵北京,指挥一批得力干将,让他们去最前线冲锋陷阵,以最高的效率去解救受灾的群众。然而现在的情况看来是,总理每事事亲躬,一方面反映总理关心受灾人民心切,另一方面仍然缺少一个高效运作的政府体系,后者对于一个总理来讲更是一个艰巨的任务。看着60来岁还这么操劳的身影,免不了担心他的身体情况,经不起这样的折腾呀。

灾害也反应了太多的事情,地震局的职能多大程度上得到了体现,据说当地的地震监测站都被开发商给拆了,地震消息居然还是美国的地质调查局先于发布。而据朋友观察,当CNN、BBC等外媒网站都第一时间发布这个重大事故,CCTV网站还兴高采烈的庆祝奥运火炬。中国红十字会动员大家捐款,我试图登陆它的网站,却发现总也连接不上。公布的帐号是工商行的,却被网友说,银行还要收他们的手续费。他们肯定会说系统还没有更新过来,然后推上一天两天再公告说免捐款手续费。只是那一两天已经足够将银行本身捐的一千来万给挣回来了还有余利呢,真是两面讨好呀,不用出一分钱,又捐了款,又说是为大家最快时间推出了免费的绿色捐款渠道。

无论如何,看到不足的,我们更是看到了太多感人的内容,全国人民空前的团结起来,论坛上基本都是祝福和关心的帖子,谁说中国人的素质不高了?

Expensive sun?

看到网易上的一个新闻,太搞笑了。成都铁路局贵阳客运段将贵阳翻译成“expensive sun”,警(路)风翻译成“Police road breeze”等低级错误(想乐请见新闻原出处)。这也算了,承认个错误,说不懂英语,重新改过就好了。当记者采访他们的领导,这位可爱的主任居然说,“贵阳客运段的列车也有外国友人乘坐,但此前没有接到来自外国乘客的反映”,“当初这些标牌制作时曾经请教过外语方面的专家,但可能由于中英文理解上有出入,所以出现了以上的问题”。言外之意难道是想说,你们这些人懂不懂,老外和外语专家都看过了,(没问题,你们瞎起哄什么呀)。

老外没报告,我可以理解,大多数老外来中国是找乐子的,看到这么好玩的,这趟也没白来,回去跟朋友吹更有谈资了。另外老外的性格也是不愿意说三道四,更没有义务向你什么政府的报告。只是这些“外语方面的专家”我就奇怪了。任何一个学过英文的,即使理解再有出入,也绝不会翻译成这样,除非根本没看过,又或者那专家根本就不认识中文。说“曾经请教过外语方面的专家”是百分之二百的睁着眼睛说瞎话。这位主任想着这回出漏子了,想给自己扯点借口而已,却不想这借口也太无知了。

读取网站的Alexa排名/Get Alexa ranking data for your site

南卓铜(Zhuotong Nan, [email protected])

由于网站自己设置的网站访问数有时不真实,为了比较网站的访问量,我们一般使用权威的第三方网站来比较访问量。Alexa网站提供被大家认可的排名数据。比如,访问http://www.alexa.com/data/details/traffic_details/westdc.westgis.ac.cn,可以看到“西部数据中心”目前排名访问。

Alexa提供了收费的Web service允许大家使用其数据,大概是每1000次请求0.15美金(见这里)。收费并不高,而且包括众多的功能。

然而作为程序员,有时候宁愿挑战一下自己的能力。比如有没有一种免费而且合法的手段来获取它的排名数据,比如Westdc.westgis.ac.cn目前排名1,080,823里的这个名次(May 06 2008)。

Alexa为了挣钱,使用了一些方法来防止简单的页面数据获取。比如我们看排名的HTML片断:

<span class=”descBold”> &nbsp;<!–Did you know? Alexa offers this data programmatically.  Visit http://aws.amazon.com/awis for more information about the Alexa Web Information Service.–><span class=”c669″>1,</span><span class=”cbf1″>34</span>0<span class=”cd05″>80</span><span class=”c9d1″>,8</span><span class=”c2e8″>23</span></span>

直接从Web页面拷贝的结果是1,34080,823,而不是正确的1,080,823。这是因为Alexa增加了一些<span>标签来混淆HTML代码,这些<span>的CSS被设置成display:none,所以在浏览器里显示却是正确的。而且这些混淆的<span>标签是随机任何组合的。

解决方案可以从模拟浏览器显示出发,逐步剥离没用的信息,最终获取排名数字。

a. 获取整个HTML源代码;分析获取源代码中有关排名的HTML片断;
b. 下载干扰的CSS表,取得display属性为none的全部css类名;
c. 利用css类名列表,从HTML片断中移去对应的<span>标签和标签内的数字;
d. 移去剩余的HTML标签;
e. 转成数值输出。

以下代码演示了此方法,使用了c# 2.0,在Visual Studio 2005编译运行通过。代码里使用了正则表达式。

/* Purpose: to get Alexa ranking data by using c#
* Author: Zhuotong Nan ([email protected])
* Date: May 06 2008
*/
using System;
using System.Collections.Generic;
using System.Text;
using System.Text.RegularExpressions;

namespace wml.stat
{
class AlexaRanking
{
public static int Rank(string url)
{
int ret = -1;

Uri uri = new Uri(url);
string newUrl = “http://www.alexa.com/data/details/traffic_details/” + uri.Host;
System.Net.WebClient wc = new System.Net.WebClient();
string html=wc.DownloadString(newUrl);

//pattern for obtaining html codes in relation to ranking data
string htmlpattern = @” about the Alexa Web Information Service.–>(.+?)</span><!–“;
string snipet = Regex.Match(html, htmlpattern).Groups[1].Value;

//get css file which store css classes preventing from scrambling
string cssUrl = “http://client.alexa.com/common/css/scramble.css”;
string cssfile = wc.DownloadString(cssUrl);

//css class pattern for getting CSS class listing with no display to the browse
string cssclassPattern=@”.(.*?) {“;
MatchCollection cssmc = Regex.Matches(cssfile, cssclassPattern);
//css classes without display, forming reg patterns
List<string> css_nodisp_patterns = new List<string>();
foreach (Match m in cssmc)
{
css_nodisp_patterns.Add( “<span class=”” + m.Groups[1].Value
+””>.*?</span>”);
}
//remove those classes from html snippet
foreach (string p in css_nodisp_patterns)
{
snipet=Regex.Replace(snipet, p, “”);
}

//see html snippet left
//remove span tags
string tagPattern = “<[^>]*>”;
snipet=Regex.Replace(snipet, tagPattern, “”);

ret = Int32.Parse(snipet, System.Globalization.NumberStyles.AllowThousands);
return ret;
}

static void Main(string[] args)
{
AlexaRanking.Rank(“http://westdc.westgis.ac.cn”);
}
}
}

本文独立实现,但后来google发现有人利用了差不多的方法,只不过在实现上用了PHP,最终产生的结果稍有不同,见 http://plice.net/?p=10

爱,是要表达的

上午跟一个朋友讲家里情况。才想起又是很长时间没给家里打电话了。

每次打过去,我爸接的一般都很快转到我妈手里。我一般只会说,妈,都好着吧?你跟爸的身体都好着吧?她就说,都好都好。然后再问问朵朵的情况,接着说电话费贵,不说了,接着挂掉。前后一般不到3分钟。

我爸妈从不将家里发生的事讲给我听。比如我妈跟我弟吵架,都是从我堂妹或表妹那知晓的。当打回去的时候,问及这些事,总是讲没事,小事,已经过去了。也许真是没事小事,但我知道其实他们有时候是很渴望有人听他们讲。这人不是我。他们总是怕我在外面担心他们的事情——所以从来不跟我讲。

而我过的总是很平淡的生活。科研的事情是没法跟他们讲的,他们搞不清楚。讲到最后,一般他们很茫然的问,那今年挣了多少钱?温州人的思维都是将什么事情都要换算成钱来衡量做事是否值得。自从有了朵朵,情况稍有改善,因为多了朵朵的话题。

所以听到孩子妈妈给她家里打电话,总是有说不尽的话,感觉特是羡慕。也许女孩这点天生比男孩强,女儿跟父母会比较亲昵一些。

有时候想起来给家里挂个电话,但想想没啥话讲,摇摇头作罢。所以发展到现在,除非有什么事情,比如买房子要向他们借款,再比如出国了,才能给家里打个电话,顺便问候一下老人的身体情况——尽管毫无例外他们会告诉我,都好。

以前有时候看电视小说,看到亲人团聚,双方都不吝表达他们的亲情。大学时尝试学学他们,跟妈妈说几句爱她的话,鼓起勇气。妈听了呵呵的说,问我,你是不是读书读傻了?从此,我没敢尝试。父母那辈不会说爱儿女,听着儿女说爱他们也会怪怪的感觉肉麻。虽然他们无时无刻不牵挂着在外的儿子。

也许等朵朵长大,可以大方的跟我和她妈妈说,爸妈我爱你们。

我跟朋友说,我是爱他们的,他们明白。朋友说,那也是需要表达的。

所以决定今天还是打个电话回去报告一下我在这边的生活吧——因为朵朵和她妈妈回兰州了,也算个大事了。

爱国

韩国逮捕了一名21岁的中国留学生,起因是这名留学生在奥运圣火传递期间的集会上暴力攻击韩国民众。我没看到攻击的图片,但我同事告诉说打的头破血流,很是血腥。中国外交部发言人此前就韩国声明要严惩此次集会中的暴力行为发言说,中国留学生只是为了保卫圣火,请韩国考虑实际情况。

这些被逮捕或者将被逮捕的留学生面临的韩国法律的惩罚,可能导致驱逐出境。这后果估计是这些学生事先没有想到的。我不知道这位被逮捕的学生现在是怎么想,会怨恨国家没出面救他吗?

由于奥运会,大家爱国热情大涨,这本是好事,我感觉也是政府在努力倡导的。有了凝聚力是好事前提是爱国得用到合法的事情上去。如果借爱国的名义干违反法律的事情,或者扰乱正常的社会和经济秩序,不管在国内还是国外,估计都是被禁止的。事实上对中国的形象也不见得是好事情。有这样一句话,我可以不赞同你的观点,但坚决悍卫你说话的权利。民主社会里你有集会的权力,同时也赋予别人有不同声音的权力。暴力是不对的。

我并不为留学生感到惋惜。然而,外交部发言人不痛不痒的发言也许是种外交辞令,但总让我不是很开心。但凡说话示弱的只有三种情况,有求于人,有把柄示于人,或者人格或国格的缺失。外国人在国内处在大老爷的位置就很不正常,政府是要代表自己民众的利益。胶济铁路事故后,几个法国人受了轻伤但被迅速转移到北京协和医院,铁道路副部长还亲自去看望,而且全免治疗费,弄的几个法国人很惊讶,说想不到,他们在淄博医院已经受到很好的治疗。然而不和谐的是,有报道指出由于床位紧张,一些中国伤员得不到很好的休养,要转院却面临医院不接收的问题。类似例子比比皆是,外国人可以在中国的土地上恣意妄为,处在高级公民的位置,有时候让外国人自己也惊讶。韩国人在国内也是乱来惯了,作为表达中国的一种态度,我更希望政府就逮捕几个违法的韩国人也赶了出去,未必不是种强势的外交手段。

离别

时间好快。才记得过去机场接她们过来的时候,又到了离别的时间。上午将两位lady送到了机场。有点不大敢回家了,怕面对空荡荡的房间。有时候烦朵朵的不听话或者缠人,但其实如果没有她的声音会少掉好多的乐趣。

然而为了生活总得放弃一些东西,不能总是在一起,好在剩余的时间也不很多了,盼望着5个月后回家再见她们了。

祝旅途一切都顺利了,当妈妈的只有辛苦着多担待一些了。