2013年第1期(总241期)大地春秋栏目,“南卓铜:引航西部地学信息化研究”
这名字起的真得让人汗颜,西部研究人才济济,愧不敢当呀。
原题 中科院知识创新工程重要方向项目群“地表过程集成系统研究”各项目顺利通过验收
9月9日,中科院知识创新工程重要方向项目群“地表过程集成系统研究”各项目结题验收会议在地理资源所召开。验收组由8位专家组成:国家自然科学基金委宋长青研究员、冷疏影研究员,中国科学院遥感应用研究所施建成研究员,中国科学院寒区旱区环境与工程研究所李新研究员,中国科学院地理科学与资源研究所周成虎研究员,中国科学院水土保持研究所邵明安研究员,中国科学院新疆生态与地理研究所陈亚宁研究员,中国科学院南京土壤研究所张甘霖研究员。项目组织部门中科院资环局副局长冯仁国、国土与遥感处处长翟金良、副处长段晓男出席了会议。
首先,冯仁国在会上致辞,他介绍了项目的立项背景,强调了中国科学院对纯科学问题研究的重视及投入,简述了地表过程的集成研究的研究背景,同时感谢各项目和课题的承担单位及研究人员的辛勤工作,并对做好课题验收及后续项目验收工作提出了要求。
中科院寒区旱区环境与工程研究所南卓铜研究员、中科院遥感应用研究所陈良富研究员、中科院南京地理与湖泊研究所张奇研究员、中科院地理资源所刘慧研究员先后代表四个项目做了结题报告,向专家组汇报了课题的目标任务及完成的主要成果、应用情况、经费支出等情况,并对专家组的疑问与质询进行了逐一解答。
专家组在听取了项目的汇报,查阅相关技术报告等项目验收材料,并进行质询和讨论后,对四个项目分别给出验收意见。专家组一致认为,各项目都完成了项目任务书规定的研究内容,达到了考核指标要求,验收文档齐全,同意通过验收。
出处:北京地理所
受青海省环保厅委托,我们执行的青海湖流域生态环境数据共享平台的建设,现在已经完成主要的功能,界面也出来了。把数据/元数据正式进去后,就可以较完满的完成项目目标了。
DEEP-Q: Data Center for Eco-Environment Protection in the Qinghai Lake Basin
上午银川某部队水文勘测的M总和W工过来找我们。他们要给叶城—阿里区域边防驻地解决用水问题,通俗的讲就是打个水井。但因为在冻土区可能面临的问题他们没有经验,要求我们给提供一些有关冻土分布、厚度等资料以及冻土对地下水资源影响的有关专家经验。谈话里讲到冻土区打井可能遇到的一些问题,感觉很有意思并且深受启发,有机会可以在这方面深入做一些工作。
比如在多年冻土区打水井,一般需要取冻土层下承压层的地下水,那么如果打穿的话,水在压力作用下会上涌,区别于非冻土区,经过一段时间井里的水可能会结冰,那么井就废了。怎么解决这个问题? 在井下增加加温装置可能不大实用,边防驻地没有大功率发电机,而且费用不菲。一个可能的方案是水冻结前,使用水泵再抽一次水,利用地下水本身的正温防止水井冻结。然而这里还有问题,多长时间抽一次呢? 可能需要做一些实验建一些模型进行计算,得到合理的一个时间间隔。
这些实际的工程需求平时都没有注意到。埋首实验室辛苦工作的同时,如果能结合一些实际需求做一些工作也是很有意思的。
swat 2005源代码是可以从官方网站上下载到的。尝试了gfortran 编译,首先删除了*.mod,使用 gfortran –c modparm.f 重新生成 parm.mod。在main.f里注释了第一行,即include ‘modparm.f’。否则会出现重复定义的问题。最后编译,可以成功,但在运行时出现 Fortran runtime error,大致是从文件里读取end of file之类的错误。
安装了intel fortran 编译器,用的是官方的30天试用版。在visual studio 2010里建立了fortran console application项目,将swat 2005里的全部*.f 加到里面。同样重新生成 parm.mod,并修改了main.f (同上),编译,试运行,可以用。
学生下载了gldas v1 0.25x 0.25度的 3-hourly 数据,据介绍里面是有rainfall数据的。但用grads或者matlab read_grib 读取时,发现没有此项,如图:
这很奇怪,因为相关说明里分明是有rainfall一项的,如在http://hydro1.sci.gsfc.nasa.gov/dods/GLDAS_NOAH025SUBP_3H.info介绍的,
而且对比两者,里面的参数名也不一致。
仔细研究gldas数据对应的元数据,发现每个文件每项里都有对应的最大值最小值,试着读取并对比,发现第7项4LFTX即是rainf 一项,最大值一致。看样子是read_grib等代码是专门for WMO的,跟gldas的描述不一致。所以在读取gldas等数据时,还需要注意了。
在ftp://hydro1.sci.gsfc.nasa.gov/data/s4pa/GLDAS/GLDAS_NOAH025_3H_E1.002/doc/README.GLDAS2.pdf找到parameter table的相关信息,如下图
注意其中的pds ids,rainfall rate是132,对应于图1的第7个记录,#=132。
大概是因为read_grib的参数表是基于ncep的某种格式,读取gldas需要配置自定义的参数表。
本来是要写另一个的,结果三天前换成决策支持系统的。累的呀,不过终于在截止期前也赶出来了。感谢书功、彦博、小余、小曹的大力帮忙,光我一个人那累死也忙不出来的。
这期间还搞了一个三江源生态环境监测的本子,今天315也是截止日期,因为要查新报告和合作协议,今天也是赶出来的,让小余租了个车给西宁那边给送过去了。感谢小余和老杨。
希望过了这阵子可以稍轻松一些了。