ASD 降尺度方案

一、       安装MATLAB7.0（高版本会出现.p文件中的函数不兼容问题）

MATLAB7.0版本安装包及其在win10系统下的安装方法见：https://blog.csdn.net/feecooling/article/details/7525140（注意安装完成后将matlab属性中的兼容性改为Windows Vista service pack 2）

安装ASD程序包，安装包见附件。

二、       准备数据

注：本文中GCM数据采用加拿大开发的CanESM2模型输出的数据，下载链接：http://climate-scenarios.canada.ca/?page=pred-canesm2

 研究区气象站点实测数据；NCEP再分析资料；GCM模拟的当前和未来气候条件下的输出数据。

 注：气象站点的日最高、日最低气温和日降水实测资料采用中国气象局气象信息中心实测数据再预处理后的数据；NCEP再分析资料可能会与GCM模式输出的数据分辨率不一致，所以建议采用双线性插值对齐分别率；站点的实测数据时间跨度可能会出现于NCEP、GCMs的时间跨度不一致的情况，因此需要对原始数据进行裁剪对齐数据操作。（附件fixnumber为实现批量裁剪原始数据长度的MATLAB代码）

      三、运行ASD

1.打开MATLAB，把current directory设置为ASD程序所安装的文件路径

2.在command window里输入ASD，打开ASD主页面(见下图)

3.点击setting配置NCEP索引文件：（见下图）模型会根据这个文本去寻找输入的NCEP路径下（NCEP predictors path:）的各气候因子文件。因此，需要将该txt文件名修改的与NCEP各气候因子的文件名一致。

4.点击setting配置GCM索引文件：（见下图）模型会根据这个文本去寻找输入的GCM路径下（GCM predictors path-current/future:）的各气候因子文件。因此，需要将该txt文件名修改的与GCM各气候因子的文件名一致。（附件中indexlist.m是实现读取文件夹内的文件名制作txt文件的程序，rename.m是实现将GCM-future/RCP各情景下的文件名修改成对应的historical中气候因子文件名的程序。其目的是统一current和future路径下个气候因子文件名一致，使得ASD模型可以根据GCM索引文件找到各气候因子的数据。）

5.设置模型参数：（见下图）

基准期，初始时间；未来时期，初始时间；率定期，初始时间；率定期，结束时间；验证期，开始时间；验证期，结束时间；年内日数（注：CanESM2未来时间输出的每年为365天）；未来时期分段（此处将其分为了近中远三个时段分别以2020,2041,2071为初始点）；Bias参数为模型偏差；variance inflation factor为方差膨胀系数，其表示因子间的共线性。一般VIF越大，因子间相关性越大；Percentile，为计算百分位数，一般用以分析极值模拟情况；miss value，ASD模型会将小于miss value这个值的数据自动识别为缺失值；Number of scenarios表示模拟的次数。

         6.设置因子挑选方法、回归方法及阈值：（见下图）

图示中选择的为MLR(多元线性)stepwise(逐步)回归方法；因子挑选显著性水平设置为0.05，最大挑选因子数5，条件模拟（由于降水是条件事件所以选择的是条件模拟；模拟气温时此处改为unconditional）；统计各月参数情况monthly；阈值1，即数据中降水小于1mm的默认为没有发生降水事件。

        7.率定模型：

      8.对历史期进行模拟：

    9.对未来情况进行模拟：

   10.分析模拟情况：

     在graphs选项中可以查看各类图形分析模拟情况。

因子间相关性分析。

NCEP各月均值模拟情况……还有很多图，我们就不一一列举啦！

希望对各位有帮助，也欢迎各位多多交流。