目录

1 简介

2 下载和安装

3 气象数据的获取

4 模式使用

4.1单条后向轨迹

4.1.1建立轨迹

4.1.2运行轨迹

4.1.3轨迹绘图

4.1.4输出轨迹

4.2聚类后向轨迹

4.2.1重新建立轨迹

4.2.2计算日轨迹

4.2.3数据集合

4.2.4运行轨迹聚类

4.2.5运行结果图例

4.3浓度扩散模拟

参考文献

　　HYsplit (HYbrid Single-Particle Lagrangian Integrated Trajectory) ，混合单粒子拉格朗日综合轨迹模型是一个完整的系统，可使用烟团或颗粒方法计算轨迹复杂的弥散和沉积模拟。它由一个模块化的库结构组成，每个主要应用程序均具有主要程序：轨迹和空气浓度。

　　下载网址：https://ready.arl.noaa.gov/HYSPLIT.php，进入网址后，点击2个链接，第一个是HYSPLIT安装包，如Download Public (unregistered) Version，第二个是图形实用程序，如Graphical Utilities，图形实用程序应该在HYSPLIT之前安装。

　　　　　　　　　　　　　  图1 下载页面

    　　安装结束后，界面如图2所示。

　　　　　　　　　  图2  HYSPLIT主界面

　　要计算后退气流轨迹，最重要的是气象数据，一般常用的气象数据为GDAS1。下载地址为：https://nomads.ncep.noaa.gov/pub/data/nccf/com/hysplit/prod/gdas.202012/。

      文件格式为“gdas1.英文月份缩写+年份.第几周”。例如，“gdas1.dec20.w1”指的是2020年12月第1周的气象数据。其他Hysplit支持的气象数据格式详见： https://ready.arl.noaa.gov/archives.php。 　　/data2arl目录中提供了所有以模型可以读取的格式创建气象数据文件的Fortran源代码。

      点击Trajectory,在下拉菜单上点击Setup Run,进行参数设置。

                                          图3  参数设置选择界面

     设置完以上5项后，点击下方save进行保存。

                              图4  参数设置界面

        等到100%后Exit。若不运行，则代表轨迹建立时有误，需返回修改。一般是因为时间和气象文件不一致导致的。

                                          图5  运行轨迹界面和结果

      点击Trajectory,在下拉菜单上点击Display,然后再点击Trajectory进行绘图参数设置。

                                    图6  轨迹绘图界面选择

    设置完以上参数，调整输出图形的大小，调整好后Execute Display。

                              图7  轨迹绘图参数设置界面

        图形使用HYSPLIT自带GS软件查看，并且使用GS转换成WINDOWS可读的格式。

                                            图8  输出轨迹图

        在后向轨迹的基础上进行轨迹分型。分型只能进行1层轨迹的分型，所以要重新进行轨迹的建立。做分型时垂直高度的类型是使用数据“input model data”的分层。

                     图9  轨迹参数设置图

           输入年、月、起始日、持续时间段。选择好后Execute Script。

                                            图10  Special Runs/Daily参数设置图

                                   图11  生成轨迹结果图

        输出的日计算结果tdump*在C:\hysplit4\working\ 或C:\hysplit4\cluster\ working\下，每个文件大约8K，将计算的结果放在一个自建的指定文件夹中，如 20201205，放在C:\hysplit4\cluster\working\下。

                              图12  轨迹输出保存路径

      (1)Time interval 是指轨迹点数据的时间间隔，一般都为1h。输入路径，即上一步输出的文件的自建路径；剩下路径不需要改。Special Runs→Clustering→Standard。

                             图13 轨迹聚类设置

  (2)  点击Make INFILE，每次更改文件内容或输入路径，都需点击此键进行识别，类似“刷新”。若建立轨迹时修改了输出文件名，此处也许需修改成对应的文件名，并点击“Note”进行识别。

                            图14  轨迹聚类设置参数图

（3）点击Run cluster annalysis，如下图所示。

                      图15  轨迹聚类分析结果图

（4）点击“Disply total spatial variance”查看折线图拐点所对的横坐标值；也可通过点击“run”。如下图。Run之后得到了可能的分型条数及变化的可能性（数值越小，变化越少，越好）如：下图所示 18条分型变化最小，但是条数太多不易分析，因此选择7条分型（可变性为92.42）为最佳。

                       图16  轨迹聚类个数求解图

（5）将step4中的分型数改为选定好的聚类数量7，依次点击Note，Run得到分型结果如下：

                           图17  轨迹聚类结果图

（6）分别绘出平均图和分形图，“Display Means”为聚类结果，“Display Clusters”为输入文件包含的所有轨迹图。

        轨迹聚类均值图如图16，所有轨迹聚类图如图17。

                                 图18  轨迹聚类均值图

                                          图19  所有轨迹聚类图

        （1）在浓度菜单下，有8个选项：快速开始（Quick Start），运行设置（Setup Run），运行模型（Run Model），显示（Display），实用程序（Utility）,特殊运行（Special Runs）,多处理器（Multi-Processor）,浓度菜单帮助（Concentration Help）。         （2）运行设置（Setup Run）菜单提出了与轨迹菜单类似的起始信息要求。还有三个子菜单：污染物-可用于设置排放速率，排放的持续时间和开始时间；网格-设置浓度输出网格的位置，分辨率，水平和平均时间；和沉积-设置每种污染物的特性。单击“检索”，输入示例预配置的控制文件的名称：sample_conc，然后单击“确定”。关闭数据输入小部件后，单击“保存”，设置菜单将关闭。

                                     图20  浓度扩散运行设置图

                  图21  浓度扩散污染物、网格和沉积设置图

                        图22  浓度扩散污染物设置图

      （3）从主浓度菜单选项卡中选择运行模型（Run Model），该模型会将设置配置复制到模型的输入CONTROL文件中并开始模型仿真。消息将出现在标准输出上，显示计算进度或计算完成后。请耐心等待，因为浓度计算可能比轨迹计算花费更长的时间。单击退出以关闭消息窗口。此时，二进制浓度输出文件已准备好转换为图形显示。

                              图23 浓度扩散运行结果图

     （4）单击“显示/浓度”（Display/Concentration），然后选择“轮廓”（Contours）以运行一个特殊程序，该程序将二进制浓度文件转换为适合打印的Postscript文件concplot.ps。显示窗口小部件包含用于不同污染物（如果已定义）的多个选项，数据网格，级别和轮廓选项。这些将在图形部分中详细讨论。对于此示例，接受默认值，然后单击“执行显示”。下图显示了释放一小时后的12小时平均空气浓度。

                                                        图24 浓度扩散显示设置图

                               图25 浓度扩散显示结果图

1.https://ready.arl.noaa.gov/HYSPLIT.php （官网） 

2.https://zhuanlan.zhihu.com/p/34885340 （速成指南—Hysplit后推气流轨迹聚类分析 - 知乎）

3.https://wenku.baidu.com/view/fa95028627fff705cc1755270722192e4436585f.html（HYSPLIT功能简介和使用说明解析）