Welcome to LZN's Blog!
Wind extinguishes a candle but energizes fire.-
CESM RESTFORMAT: NUMBER OF RECORDS ON RESTART FILE NOT SUPPORTED
FS in OUC is trying to setup CESM fully coupled Pre-industrial (B1850 compset) cases while faced with some problem.
In the run log, we found:
ERROR: restformat: number of records on restart file not supported
By searching, the error seems to come from the ice restart file.
We then checked the default restart file, which can be located in
ice_in, in${CESM_INPUT}/ice/cice/iced.0001-01-01.gx1v6_20080212. The file is zero in size. Thus, it must be created by the model, and the model cannot read anything from it. We re-download the corresponding file and the problem solved.Updated 2018-12-26
-
Nuuk in the air
看了下github，已经三个周没有contribution了。上周AGU，应该是毕业前最后一次出国开会，早就准备，YQ也一并同去，两人第一次在国外开会overlap，一起逛了DC的两个博物馆，方尖碑，林肯纪念堂和国会山，看到了梵高自画像的原作，YQ很开心，好久没有见她这么长时间持续的开心了，我也很欣慰。
返程从NY直飞CN，整航班的纨绔美本，着装相貌言谈举止，耳闻目睹间与影视作品中对接角色相差不大。想这代人与他们父辈之间已经判若云泥，略感唏嘘。
航班起飞之后直冲正北，冰洋航线，第一次飞入北极圈。飞西海岸都是走白令海测地线，阿留申群岛掠过，并没有进入北极圈。如今看这长虹贯日般的冰洋航线，想到进入极夜区或许能看到极光，不禁一阵兴奋。
进入北极圈后很快就进入极夜区，与飞越晨昏线的情况类似，似入黄昏一般逐渐暗淡，夕阳在右舷窗外，下面是时而包块时而条状的极地层状云。时值上弦月，斜挂西天，银光下泻，照亮一片云缝之中格陵兰的冰川和海面。
除了引擎的轰鸣，窗外静谧地像油画一般，完全感受不到这片静谧之中是彻骨寒冷。朦胧中忽然发现一片灯光，飞行信息上显示是格陵兰的首府，努克。呵，从这个高度望去，这个小镇的灯火如美国商铺门口的OPEN彩灯一般精巧，不知道这里的人们过着怎样的生活，与我9年前期盼的山顶测站相比如何？有生之年若有机会，愿来这里守候极光。
Updated 2018-12-18
-
Mask region in the grid system according to shapefile
Here is a demand from YQ that the area of Shanghai need to be marked out from the WRF/CMAQ grid system. Previously, we have tried to mask out the whole region, e.g. the PRD. Now, we need to read the shapefile, and mask the region (closed polyline) by the given name.
A NCL Page talking about the shapefile in NCL. However, we do not need to know deep in the GIS architecture. The only thing we need to know is how to mask the region by named polylines.
Still, the
shapefile_mask_dataworks for this. Just an optional attribute is enough.opt = True opt@debug = True opt@shape_var = "NAME_1" opt@shape_names = "Shanghai"ncl_filedump -vcan be used to iterate all possible values forshape_var. Original script is available on the github.Updated 2018-11-13