学而时习

[#67] 7Timer! still not available for operating [EN]

Dear all, sorry for the inconvenience. Some unexpected bugs come out and we have to delay on using the 35km datasets.

New ArP is completed, but the main processing software -- GrADS, seems to have encounted a big problem and cannot put out the indexing map file.

傻傻的电脑

  为了第三代7Timer!的自动运行,我花了两小时啃入了BASH的门。我的程序中某一个地方要循环体,需要定义一个循环相加的变量。于是我定义
TEST = 9
  下一个循环中它要是12。根据学C、jsp、php之类的经验,我毫不犹豫地写
TEST = $TEST+3
echo $TEST
  各位猜怎么着?电脑居然显示
9+3
  我狂晕!觉得又好气又好笑。当然,这等问题本来不应该拿到桌面上说,不然肯定被老手笑话你这个菜鸟不见世面,可我的确狂笑了1分钟:电脑啊电脑,你究竟还是只能听人指令的傻东西!哈哈哈哈哈哈哈哈!
  后来查了一下,正确的写法应该是
let "TEST = $TEST+3"
echo $TEST
  电脑就能够显示“12”了。

[#66] Run for Feb.21 delay [EN]

As we are now upgrading the 7Timer! system from 70km resolution to 35km resolution, auto-run shedule for today is delayed.

News will be posted here if there are any.

[#65] Major achievement on 35km-res [中文/EN]

  各位好。我高兴地向大家宣布一个好消息:7Timer!在采用35公里精度GFS数据上取得重要进展。预计在这两天就可以将所用的数据由70公里精度(1度)提升到35公里精度(0.5度),第三代7Timer!触手可及。所取得进展的详细内容:

  • 在本地完成编译cnvgrib模块并在服务器上测试成功。该模块能够将GRIB2格式的35公里精度数据编译为GRIB1格式,供GrADS处理及分析;
  • 服务器上测试grib2ctl模块成功,它可以生成35公里精度数据的ctl文件,供GrADS索引;
  • 服务器上测试get_gfs模块成功,可以有选择地获取远端GFS数据,节省下载时间和运算量。

  要全面采用35公里精度的数据,还有以下几个问题要解决:

  • ArP改为调度get_gfs模块获取远端数据,并需要调度cnvgrib模块来编译数据;
  • 由于数据名称发生改变,需要对gs文件组的相应部分进行修改。

  升级35公里精度以后,预计将有以下计划:

  • 由于精度提高,拟重新开放地点预测功能(7Timer!-Panel系列),并提供视宁度预测等;
  • 开发民用版本(7Timer!-Weatherman系列);
  • 测试WRF-NMM,对精度的进一步提高进行调研。

  由于升级问题,今天7Timer!的运作可能被取消,请各位留意进一步信息。

Dear all, we have a major achievement on the next generation of 7Timer!, it may be applied for use in these two days. By then, the resolution of 7Timer! will be raised from 70km to 35km. We're going to use cnvgrib to convert files from GRIB2 format to GRIB1 format, so GrADS can read and proceed them. We would use "fast download GFS files" technique to require the data files faster. As the resolution raised, 7Timer-Panel might be re-operate to public.

Today's run may be canceled or delayed due to the major operation, we would post news here if there are any.

[#64] New ArP version released [EN]

Hi all, a more powerful version of ArP (Auto-run Program) is just released. It could give more information about ArP, include current status, progress, or whether it has been started on time. Also the program will now auto inform me when the program crashes or fail to operate. You could see the ArP panel on the first page of 7Timer!.

ArP version now is 2.2-s8.

第四代7Timer!—— 5公里精度?

  野心太大,有时真是一件坏事!嘿嘿!

  昨天在7Timer!-INFO N.63里兴致勃勃地发布了GRIB2数据成功解码的消息,又顺利地把首页的图弄成了Mollweide投影,已经觉得是“重要成就”了,一个晚上都乐得屁颠屁颠的。可今天太阳还没落,GRIB2 decode技术好像已经是恐龙时代的事物一样。为什么呀?今天仔细一看才知道,NCEP打算今年8月1日开始全面转用GRIB2格式的数据,GRIB1到了那时候就不再发布。也就是说,7Timer!最迟在今年8月1日就得开始采用0.5deg(今天换算了一下,0.5deg应该是35km。55km是对应赤道的)数据了。今天才知道可以用NCEP自己发布的cnvgrib(http://www.nco.ncep.noaa.gov/pmb/codes/GRIB2/)把GRIB2转为GRIB1,再用GrADS来处理。至此,7Timer!采用35km数据理论上已经扫清了最后一步障碍了。采用GRIB2数据的7Timer!将是第三代7Timer!,其精度可以大大提高。如果第三代7Timer!应用的话,那我的地点预测就可以恢复了。这项业务将命名为7Timer!-Panel2+,而准备投入民用的7Timer!-Weatherman(暂定名)也在开发中。

  最不妙的就是很快又发现,大名鼎鼎的WRF模型原来是可以用GFS/ETA的数据来算出精度更高的预报。其实么,本来是想再看看MM5的可行性。最后发现MM5是可以做到的,因为现在我对Linux/Unix下的操作也算基本入了门,而MM5原来用的是NCEP Reanalysis的数据(还可以用国内的T213?不管,反正不好找到)。用了MM5,精度可以再上一个档次,就可以和北美CSC叫板了,正是我梦寐以求的啊。然而MM5毕竟要把现在的平台全部推倒重来,代价似乎太高了。用WRF的话,可以顺用GFS数据,只是用WRF内核来处理,然后再用WRF2GrADS转到GrADS下处理,现在的平台还可以用上约一半(不知道我理解有没有错误)。精度呢?哇哇哇,WRF可是中小尺度预测的,目标是1-10km。想象一下,5km精度的7Timer!,让人口水都要流出来了!

  但精度做到这么高,代价肯定很高—— 服务器绝对吃不消。我的设想是,先对一些小地区(比如,北京上海广州)进行中小尺度的预测,如果效果不错,说不定可以向崔博士那边申请一台专门的机子……?哈哈,别想那么远了。现在第一步就是要上手WRFSI…… 万里长征第一步还是要迈出去才行啊。

  至于原来的计划么?新版本ArP?唉,后延后延…… 又一个星期没出去了,还是要找天出去溜溜啊。各位,罗马不是一天建成的,但高精度的第四代7Timer!已经在地平线上出现了……

认真研究Kreutz

  最近要做的事情真多,不过也犯不着和别人说“忙”,因为这些事情么,都是召之即来的,要是不喜欢可以拖延,十分自由。

  昨天拿到天爱《增刊》,发现里面介绍有位MM弄Kreutz轨道的。呵呵,牛啊,不看都不知道北方牛人原来还不少。自我批评一下,应该向dreamh学习,d调。

  话说回来到Kreutz,一直很低调地研究克族的性质,现在也可公开一些了。地面研究虽然不是没有可能,但看来难度还是很高的,决定先准确掌握Kreutz性质再说。今天翻到Sekanina的一篇重要论文:http://www.journals.uchicago.edu/ApJ/journal/issues/ApJ/v607n1/59666/59666.html,因为实在太长又太复杂,决定慢慢啃,顺便放在这里看看谁有兴趣也来研究。晚上,朱头把孙同学的论文发了来,可谓雪中送炭哈!原来她把轨道根数已经算了个平均值出来。这些天我还偷懒不想算呢…… 这下有人送上门来了。孙同学还列出了若干重要数据,明天慢慢研究!现在先睡觉。

[#63] GRIB2 decoded succeed, recent progress (中文/EN)

  各位好。非常高兴地向大家宣布一个十分重要的进展:今天经过反复试验,通过degrib的方法顺利对GRIB2格式的数据进行了解码。这些数据的精度达0.5度(约55公里),比现在7Timer!所用的GRIB1格式的数据精确8倍。附测试图片,2月4日北京时间20时的GFS自带云量输出图,可以对比一下新格式数据的效果。

  由于degrib对GRIB2的解码仍然不是非常好用,因此短期内新数据还不会正式投入使用。计划继续对7Timer!-ArP进行升级改造,提供更为智能的ArP运行信息。

[[Image:Other/20070204a.png|500|GRIB1 format (1deg.)]]

[[Image:Other/20070204b.png|500|GRIB2 format (0.5deg.)]]

Hi all, glad to announce a good news: the GRIB2 data has been successfully decoded. GRIB2 format raises the resolution to ~55km, which is eight time better than the currently-used GRIB1 format. Test images are shown here.

As the DEGRIB is not quite smoothly working with GrADS, I've no plan to put the GRIB2 data into use in near future. These days I plan to continue on evaluating 7Timer!-ArP.

7Timer! encounted problem

The server moved to a new IP address and this led to another headache problem. This address has a very slow speed to GFS data server, and the ArP cannot get the data (error: connect time out). So 7Timer! suspended since Feb.1. I'm contacting the administrator for a solution.

7Timer!恢复/GRIB2进展

  恶补了php+shell课程之后,终于依靠自己把7Timer!弄好了,现在又对php+shell大感兴趣,终于觉得linux好处多多!最妙的就是发现php、shell原来和我学的GrADS下的fortran有些相似,学起来就更有底气了。于是决定拿ArP开刀,把lzq的代码优化一下,增添一些有趣的新功能,然后再自己写几个有意思的php小程序。

  本来打算修好7Timer!之后,就完善一下ArP,然后开始弄CMCN。没想到又给GRIB2缠上了,因为NCEP的GFS新产品,精度竟然达到0.5度,比现在我用的1度的好多了。(虽然只是差了2倍,可在3维空间里就是2^3=8!所以要努力争取啊)

  可是么,好东西就是难用上。其实高考结束我就发现它了,但就是因为GrADS不支持GRIB2才放手的。今天忙碌了一个下午,哎…… 这个GrADS也真是慢腾腾,用了2年了,还是1.8-s11和1.9-b4两个版本,估计要等到猴年马月才会支持GRIB2。和我们的班主任王老师讨教了一阵,于是现在有2种办法:

  • grib2ctl+wgrib2尝试。可这个grib2ctl真是很抠门,在linux+perl下才能运行。LASG给的exe转译文件也没法和wgrib2结合起来。wgrib2么?倒也是可以,但是要在linux环境下自己编译和安装。本人太懒,不大想走这条路,于是留下几个网站:
  • degrib。现在就在尝试,因为看来要可行的多。其实它就是先把GRIB2转成GRIB1然后再处理,但如果真的能弄成,估计要耗掉一大堆系统资源…… lamost估计要给弄垮了说。不管,我先试试看。

  于是现在事情看起来简单多了。在服务器上安一个degrib,然后就用新学的php+shell调度下载GRIB2的数据然后degrib成GRIB1,再由你宰割吧。哈哈哈哈!不过,其他的task就又要delay了……