flexpart学习总结

在野靡生

  这几周学习flexpart和flex_extract在气象家园和github上找到了很多有用的经验和代码，遇到了很多错误，得到了一些浅薄的经验，在此做一小结希望能为后来的读者一点启发。
1、如果要安装flexpart和flex_extract  b站有个up主的视频非常详细和简介，比简书的教程好不少（安装这个教程安装容易出现问题），或者直接在github下载上传的Docker，都会事半功倍，虽然安装过程中也能学到一些东西，特别是对初学者而言付出的时间和收益相比性价比非常低。
2、flex_extract目前已经可以通过cdsapi下载最新的era5数据，但最好是使用网线和虚拟机的物理网络连接，否则不保证下载成功，我下载了全球近20几年6小时的era5，有缘的话可以分享给大家。此外，也可以用fnl之类的下载GFS数据作为输入，我没有用这个所以不太了解。
3、flexpart的输入参数不是很复杂，网站里介绍的也很清晰，要改的也就是COMMADN、RELEASES、available之类的几个文件，种类设置为24就是水汽。默认最大粒子数量是十万，如果要修改，在src里面找到includepar_mod.f90文件（我记得应该是这个），修改后重新make ncf=yes即可，很简单。
4、运行完成后output会有一堆二进制文件，也可以修改输出为txt、csv之类的格式，但是没必要，最后都要处理为list的。可以批量读取然后处理为list of arrays，每行都处理为最大长度，不足的长度补齐为空即可，把这个不管存储为nc还是csv还是任何格式，就可以绘制全部的轨迹图了，下一步就可以聚类（这个气象家园已有记录，不赘述），之后的安装公式计算即可，gpt可。
5、有一个我之前的迷惑也记录在这里，我一直认为粒子应该有一个独立id，这样才能绘制轨迹，但二进制文件里面的npoint又只是RELEASES的id，这个问题耽误了我两三天，发现直接按照行来拼接起来就行，这么简单ai完完全全几分钟可以写出代码，居然浪费我这么久。
6、记得多多阅读github，上面有用的代码太多太多了，相比之下气象家园的帮助小了很多。

大部分人不愿意分享自己的经验  导致一个不是很复杂的软件有了较高的门槛，只有师门传承才最简单的学到方法，望大家都进入大组，名师出高徒