爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 65331|回复: 117

WRF模式的后处理程序-可生成GrADS和Micaps格式

  [复制链接]

新浪微博达人勋

0
早起挑战累计收入
发表于 2011-11-6 17:31:38 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
转载自国家气象中心的网站,作者说明非商业用途可以传播和使用。
下面是一个使用说明:

  1.                       wrfcnv程序说明(2008年11月4日)
  2. 开发: 邓莲堂  
  3. 单位: 国家气象中心
  4. 发布: 2008年11月
  5. EMAIL: lotusjune@163.com
  6. ========================================================================
  7.          1. 程序特点
  8.          2. 程序编译
  9.          3. 程序使用
  10.          4. 提示信息
  11.          5. 程序修改
  12. ========================================================================
  13. wrfcnv是WRF模式的后处理程序. 本程序是在参考WRF模式官方网站后处理程序和
  14. 其他模式的后处理程序的基础上开发的. 程序完全免费,任何人和单位可以自由使
  15. 用于非商业目的活动中. 程序是在平时工作中逐渐开发,逐步改进的, 因此本程序
  16. 可能存在许多小问题. 非常欢迎任何人和单位把在使用过程中发现的问题和对程
  17. 序的改进的信息反馈给我,但对程序不完善的问题造成的过错和损失不承担任何
  18. 责任.
  19. ========================================================================
  20. 1. 程序特点:
  21. (1) 处理WRF模式输出的GRIB1和NETCDF两种输入数据格式
  22. (2) 输出MICAPS和GRADS两种数据格式
  23. (3) 模式输入数据可以不进行水平插值处理,也可以进行水平插值处理,并能输出指
  24.     定的经纬度范围的数据
  25. (4) 模式输入数据在垂直方向可以不插值处理,也可以进行等压面和等高面插值处理.
  26. (5) 可以为每个输出变量指定输出层数和具体的层次信息.
  27. (6) 对于量级较小的变量(如涡度,散度等),程序会自动进行量级缩放,并在变量的单
  28.     位中给出相应的缩放信息.
  29. (7) 可以输出MICAPS的第2类格式,第3类格式,第4类格式和第11类格式.
  30. ========================================================================
  31. 2.  程序编译:
  32.     本程序是WRF模式的后处理程序, 输入数据的程序接口是采用了WRF模式的函数库,
  33.     因此在编译程序前, 必须先编译好WRF模式.
  34. (1)下载文件包 wrfcnv.tar.gz 文件,并解压
  35.     > gzip -d wrfcnv.tar.gz
  36.     > tar -xf wrfcnv.tar.gz
  37. (2)进入到程序所在目录wrfcnv,修改编译配置文件makefile.
  38.    修改makefile文件中的关键词 WRF_DIR, 将其指定为WRF模式的安装路径.   
  39.    > cd wrfcnv
  40.    > vi makefile
  41.    makefile 的内容示例如下:
  42.    WRF_DIR         =       /u/denglt/chlafs/WRFV2
  43.    WRF_INCLUDE     =       -I$(WRF_DIR)/external/io_netcdf \
  44.                            -I$(WRF_DIR)/external/io_grib_share \
  45.                            -I$(WRF_DIR)/external/io_grib1 \
  46.                            -I$(WRF_DIR)/external/io_int \
  47.                            -I$(WRF_DIR)/inc
  48.    include $(WRF_DIR)/configure.wrf
  49.    
  50.    ...
  51. (3)输入命令进行编译.
  52.    > make all
  53. (4)编译成功后,会在当前目录下生成可知性文件wrfcnv.exe.
  54.    > ls *.exe
  55.    -rwxr-xr-x 1 denglt wrf_xp 1569601 Nov 04 02:24 wrfcnv.exe
  56. ========================================================================   
  57. 3. 程序的使用
  58. (1) 进入到程序所在目录
  59.    
  60.     > cd wrfcnv

  61. (2) 编辑参数配置文件namelist.post, 具体的参数配置信息可以参考
  62.     namelist.readme文件
  63.     > vi namelist.post
  64.    
  65. (3) 运行程序.  
  66.   
  67.     > wrfcnv.exe
  68.   
  69.     当输入的数据格式是GRIB1时,在运行的当前路径下必须存在gribinfo.txt
  70.     和gribmap.txt两个文件.
  71.    
  72. ========================================================================
  73. 4. 提示信息
  74. (1) 关于垂直层次:
  75.    a. 插值到等高面时,要求以千米为单位。
  76.    b. 插值到等压面时,要求以百帕为单位。
  77. (2) 参数配置文件中有一个控制逻辑型参数overwrite,当它取为.true.时,则输
  78.     出结果文件名如果有重复,则进行覆盖; 否则输出新的文件名。
  79.    
  80. (3)  参数配置文件中针对于MICAPS数据格式的变量变量的等值线间隔值和加粗值
  81.      等信息, 在输出GRADS数据格式时不起作用.
  82. (4)  输出的变量名不区分大小写, 但只有数据文件中出现的变量名和一些特定的
  83.      诊断变量名才有效. 这些诊断变量名包括:
  84.      PRES    = " - ; 气压",
  85.      Z       = " - ; 高度",
  86.      TC      = " - ; 温度",
  87.      T2C     = " - ; 2米温度",
  88.      UA      = " - ; 东西风"      ! 水平风东西分量
  89.      UA10M   = " - ; 10米东西风", ! 10米高度水平风东西分量
  90.      VA      = " - ; 南北风",     ! 水平风南北分量
  91.      VA10M   = " - ; 10米南北风", ! 10米高度水平风南北分量
  92.      WA      = " - ; 垂直速",
  93.      THETA   = " - ; 位温",      
  94.      WIND    = " - ; 水平风",   
  95.      WIND10M = " - ; 10米风",  
  96.      RH      = " - ; 相对湿度",
  97.      RH2     = " - ; 2米相对湿度",
  98.      TD      = " - ; 露点温度",
  99.      TD2     = " - ; 2米露点温度",
  100.      T_TD    = " - ; 温度露点差",
  101.      T2_TD2  = " - ; 2米温度露点差",
  102.      THETASE = " - ; 假相当位温",      
  103.      DBZ     = " - ; 反射率",
  104.      CAPE    = " - ; 对流有效未能",
  105.      VIS     = " - ; 水平能见度",
  106.      PW      = " - ; 总可降水量",      
  107.      SSP     = " - ; 饱和比湿",
  108.      MPV     = " - ; 湿位涡, 未经检验, 不建议使用",
  109.      PMSL    = " - ; 海平面气压",      
  110.      DIV     = " - ; 散度",
  111.      VOR     = " - ; 涡度",  
  112. ========================================================================
  113. 5. 程序修改:
  114. 2007.03.15
  115. 1. 修改了xy2latlon.F中find_point_map程序中的错误。
  116.    地图投影的标准经度给错了,给成了中心经度。
  117. 2. 修改了wrf2micaps.F中求水平风场(包括三维风场和10米风场)的程序错误。
  118.    风场数据应该是用wrf_get_interp_3dvar函数获取的结果,而不是用get_wrf_u和get_wrf_v
  119.    函数获取的u和v数值。因此改用wrf_get_interp_3dvar函数分别获取u和v数值。
  120. 2007.05.24
  121. 1. 修改了read_namelist_util.F中第1112行格式语句因为不同编译平台而造成的错误。?
  122.   
  123. 20080828
  124. 1. 风场插值到等经纬度坐标时进行了旋转变换处理
  125. 2. 输出部分进行了模块话处理,将GRADS和MICAPS分别进行了模块封装。
  126. 20081022
  127. 1. 修改了垂直插植程序,对温度,高度,风场,湿度,气压等量进行特殊外插处理
  128.    在类型变量var中添加了垂直外插类型interp,对应值为:
  129.    1:温度 2:风场 3:湿度 4:气压 5:高度 6:垂直速度 9:默认变量 0:不插值
  130. 2. 改正了MICAPS第二类格式中经纬度以及地形信息的错误
  131. 20081025
  132. 1. 改正了micaps读取小于1的变量,如涡度、散度时,不能正确显示的错误
  133. 20081129
  134. 1. 修改了编译文件makefile,适用于pgf90编译器
  135. 20081129
  136. 1. 采用了固定温度直减率,修改了温度的垂直外插时过大或者过小的问题。
  137. 20081225
  138. 1. 添加了风场(WIND)的GrADS输出, 输出形式为全风速和风向  
  139. 2. 添加了2米温度的摄氏温度输出
  140. 3. 修改了MICAPS格式第3蕾格式中降水类型数据头信息输出错误
  141.   编译时,可能会由于编译器不同,造成编译不通过的问题。可能要将库函数的顺序重排一下。
  142. 200908
  143.   将二维变量的levidx改成1, 但会影响MICAPS的输出,要在module_output.F中进行修改


复制代码

wrfcnv3.2.1.tar_.gz__0.zip (257.08 KB, 下载次数: 1039)

评分

参与人数 1金钱 +2 收起 理由
meehooqq + 2

查看全部评分

密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2011-11-6 18:03:24 | 显示全部楼层
哇,果然好东西,果断顶
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2011-11-6 18:58:47 | 显示全部楼层
清风继续给力了,哈哈
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2011-11-6 19:46:39 | 显示全部楼层
绝对的好东西啊,哈哈
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2011-11-8 18:14:35 | 显示全部楼层
楼主是好人
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2011-11-17 22:26:17 | 显示全部楼层
楼主,你太好了!!!!!!!!!!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2011-11-18 16:55:45 | 显示全部楼层
我按上面的方法发现用不了啊。能具体点吗?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

0
早起挑战累计收入
 楼主| 发表于 2011-11-18 16:59:46 | 显示全部楼层
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2011-11-18 17:01:27 | 显示全部楼层
mofangbao 发表于 2011-11-18 16:59
里面有个说明的,我不做wrf,没法给你建议了

依然谢谢了啊。说明看了,还是找不到合理的解决办法。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2011-12-9 23:36:53 | 显示全部楼层
棒!! 感謝分享!!!
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表