- 积分
- 1625
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2021-4-7
- 最后登录
- 1970-1-1
![[哲理中的不平凡] 粉丝数:21 微博数:102 新浪微博达人勋](source/plugin/sina_login/img/light.png)
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
经过几周的研究,关于flexpart对于水汽追踪的工作总算是完成了,因为在气象家园也找到了特别有用的帖子,这里也非常感谢那些乐于分享的前辈,所以自己也想写个经验分享贴,希望能帮到到正在学习flexpart的小伙伴。
(1)首先是flexpart的安装,当然网上能搜到比较多的帖子,这里比较推荐FLEXPART安装配置笔记(成功案例)——筑梦之路_筑梦之路的博客-CSDN博客这个帖子比较全面,可以做个参考。但是有一个需要注意的点是,在安装flexpart之前,eccodes和grib_api可以二选一,如果用grib_api安装了flexpart,之后如果模式跑的时候遇到问题(比如识别不了grib文件)可以选择eccodes;
(2)安装好之后就是设置参数,就可以开始模拟,这里推荐看一下flexpart的指导手册,很详细,只要自己仔细看,里面的说明比较全面;
(3)最后是模式数据的后处理,如果是做水汽追踪的话,要画水汽轨迹,要用输出的part*文件(当然这个文件在默认的情况下是二进制文件,但是如果自己会Fortran的话可以去src文件夹下的partoutput.f90[应该是这个脚本]修改一下输出格式,可以设置输出格式为txt)【脚本我放在附件了】;
(4)画完水汽轨迹想继续做水汽轨迹聚类的话,可以把轨迹文件处理成nc文件,因为轨迹文件是递增的,若处理成nc肯定会存在缺测,在聚类前可以把nc数据插值填充,之后可以用kmeans聚类方法,除此之外也可以参考这个库【https://github.com/MasterOnDust/ ... aster/Example.ipynb】。如果把轨迹直接处理成列表【不存在缺测】,可以试试traclus-python这个库;
(5)聚类结束之后就能计算E-P,这个建议还是尽量多看看文献。【有问题的可以留言交流的】
最后希望气象家园能越来越好!
|
|