连续无雨日数(持续干燥指数)与连续有雨日数（连续湿润指数）的计算软件rclimdex

hillside

        连续无雨日数、持续干燥指数与抗旱天数三个词的基本意思是相同的，只不过使用在不同的场合。
        连续无雨日数是一个直白且常用的术语。持续干燥指数（CDD）近年常用于作为衡量极端气候的一项指标。抗旱天数则是水利部门习用的表达。
       在阈值上，有各种不同的门槛。气象部门常用0.1mm、1mm，也有根据农业干旱情景使用3mm、5mm不等。
       依干旱类型而言，有气象干旱、水文干旱、农业干旱、生态干旱等。

      Cllimdex系列软件(含Rlimdex、Fclimdex、Climedex的Excel宏程序)是受国际认可的权威性软件，可以用来计算连续无雨日数(持续干燥指数)与连续有雨日数及其他相关指标（Climedex的Excel宏程序中无CWD）。具体含义如下：
RX5day 5 日最大降水量                       每月内连续五日的最大降水量 mm
CDD 持续干燥指数                              日降水量<1mm 的最长连续日数 天
CWD 持续湿润指数                             日降水量≥mm 的最大持续日数 天
SDII 普通日降水强度                            降水量≥1mm 的总量与日数之比 mm
R10 强降水日数                                   每年日降水量＞=10mm 的总日数 天
R95pTOT                                             强降水量 95%分位值强降水之和 mm


        Climdex软件本人已经在极端气候标题的帖子中作了介绍，但仅需要计算连续无雨日数与连续有雨日数的网友们不一定能够直接发现，因此，本处特意专门介绍。另外需指出，该软件计算连续无雨日数的程序自带最低阈值是1mm。        在不修改程序的情况下，Rlimdex软件的使用并不复杂，比较简单。基本步骤大体如下：
        1.下载并安装R语言（近于一键式，不需设置，安装方便）；
        2.下载本站帖子已提供的rclimdex软件；
        2.打开R语言窗口；
        3.点击R语言窗口左上角的FILE菜单项下的“Source R code”加载rclimdex;
        4.LOAD  DATA AND RUN QC(数据合理性检验，包括日降水值不能为负，日最低气温不能高于最高气温)。
        5...

       以下二段摘自已在气象家园RLICIMDEX帖子中上传的《内蒙古地区极端气候事件时空变化及其与NDVI的相关性（RClimdex）》硕士论文 ：
        “RClimDex有着比较友好的操作界面，使用者可以在窗口界面的提示和引导下将研究区域的气象记录值文件进行处理和计算，结果将以excel 文件的形式输出。
       在应用RClimDex 处理数据之前，必须确保每个站点的数据以文本格式储存，并且储存的气象数据必须按照年、月、日、24 小时日降水量、日最高气温、日最低气温等顺序排列，各记录项之间通过空格将其隔开。由于研究的气象记录年限跨度很大，难免有些记录值缺测，因此程序要求所有的缺测值都统一设为-99.9。启动R 编辑器，读入RClimDex 程序，在桌面窗口的引导下读入某站点气象数据的文本文件，程序将对文件数据进行检验处理，用于提高所用数据的质量。检验处理主要涉及几个方面的错误记录检查：（1）日最低温>日最高温（2）日降水量<0mm(3)记录值严重偏离本地区气象实际情况，即超出3 倍标准差的值定义为出界值，通过人工检查同相邻站点的记录进行比对，合理的保留，不合理的按缺测进行处理。进行完质量控制处理后就可以对数据进行计算处理，得出一系列的excel 格式的文件结果。”
        关于输入数据格式，Climedex的Excel宏程序可以使用.csv文件，rclimdex必须使用文本格式。        

                                                                       在Rclimdex中修改连续无雨日数降水阈值的方法
       因为RCLIMDEX自带的CDD阈值仅是1mm，而农业干旱等实际研究中存在多种降水阈值的设定。如3mm，5mm等。一篇美国农业干旱文献设为6.35mm（见附件）。
        经过对rclimdex源程序的分析，本人已经了解修改阈值的语句。
        仅需修改3处，其他操作不变：
         1.第300行
          将prcp>=1  改为  prcp>=6.35（或其他因地制宜之干旱指标阈值）
         2.第303行
         将PRCP>=1mm  改为  PRCP>=6.35mm（此句仅用作绘图标题，并不参与计算，不改亦无大碍）     
         3.第1905行
          这是最关键的一句。应将mid[j]<1  改为   mid[j]<6.35

附：气象家园相关帖子：计算27项极端气候指数的权威级软件RCLIMDEX（R语言）
农业干旱的计算方法_D_M_赫什菲尔德.caj (129.42 KB, 下载次数: 103)

2013-8-8 11:37 上传

点击文件名下载附件
rclimdex.r (105.61 KB, 下载次数: 179)

2013-8-8 11:37 上传

点击文件名下载附件

才仔

长知识了，楼主好多气象气候相关知识

又是那隻貓

长知识了  

kongfeng0824

期待更好地更多的资源

霄霄

希望能像楼主一样，能够搞懂R语言的程序，能把这个程序最大化的应用。

qxtlyf

长知识了      

等待为了一切

谢谢楼主

meredithxf

长知识了，得好好学习一下

飞天88

谢谢楼主分享

寻路人

{:eb502:}{:eb502:}