NCL v6.3.0发布后,发现新增了advect_variable函数,这对于一些复杂诊断方程的计算无疑非常方便,于是决定用ncl来计算Q1.没想到advect_variable函数与垂直积分函数一样,运行后告诉我找不到这个函数,这什么情况,新feature没加好就发布还说有了有了……坑爹啊。
后来用center_finite_diff_n中央差分函数配合完成,计算效率感觉比之前错误得用EC monthly数据算时高非常多(主要是short2flt函数慢),何况NCEP-DOE这套资料还是daily,1度和2.5度分辨率差这么多?没仔细算数据规模,不知道为什么……
写程序时基本没有出什么错误,对这类诊断方程编程时,最需要注意的就是统一国际单位制,量纲要分析正确。还有平流项是严重的非线性项,一定要用daily算再平均到monthly哦~
第一次出图,1979年年均菲律宾海一带Q1与Day200时的Q1,一看就是对的吧:-)
Code放在了github的repo上,供参考
CD的文章默认都是这个。改下设置就可以
打开PDF文件,依次点击“编辑”、“首选项”、“文档”(种类)中选择“从不”选项并确定,默认的为“使用 PDF/A 模式下
参考资料:http://blog.sciencenet.cn/blog-211702-415720.html
又到了批量下数据的时候,linux批量wget真是顺畅无比,脚本直接用C/C++语法即可。
#!/bin/bash for ((i=1979;i<2014;i++)) { wget xxxxx.$i.xxxx }
批量解压文件,依然用循环解决
for tar in *.tar.gz; do tar -xvf $tar; done