请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 8657|回复: 5

WRF跑ERA-Interim

[复制链接]

新浪微博达人勋

发表于 2018-3-19 15:27:24 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 小莹子 于 2019-4-17 11:28 编辑

      之前一直使用的是FNL数据,但在与实测数据做点对点的验证的时候,发现误差有些大。在看了一些资料后了解到ERA-Interim数据的精度更高。现在在尝试ERA-Interim数据。
1.ERA-Interim数据介绍
http://blog.sciencenet.cn/blog-430991-888904.html
2.ERA-Interim数据下载
根据论坛上http://bbs.06climate.com/forum.php?mod=viewthread&tid=13288下载ERA资料http://apps.ecmwf.int/datasets/数据说明http://bbs.06climate.com/forum.php?mod=viewthread&tid=57876
      nc文件不能直接使用,WPS中可以直接处理各种GRIB文件的背景场(ungrib),转化成统一的中间格式数据(FILE:xxxx),再插值成初始场和边界条件,最后驱动模式。如果要用nc格式的背景场,就必须自己将其转换成中间格式。中间格式参考http://www2.mmm.ucar.edu/wrf/OnL ... IM_files/IM_wps.htm
      在下载之后发现数据的命名并不是规则的格式,比如说ERA-Int-2007*之类的,而是一大堆乱码,_mars-atls13-a562cefde8a29a7288fa0b8b7f9413f7-_kFohx.grib。一开始以为是自己下载的问题,毕竟对于长时间序列的数据而言,如果一次性下载许多个文件的话,这样的命名只会增加负担。在交流群STARMAN咨询过后,被解释说这是常态,这是ec的mars系统根据用户需要定制的文件。可以自己命名。
其中变量下载参照http://bbs.06climate.com/forum.php?mod=viewthread&tid=30962
三维数据(例如,气压层的数据)
· 温度(130.128)
· 风场u,v分量(131.128, 132.128)
· 位势高度(129.128对应位势)
· 相对湿度(157.128/比湿 (133.128)
二维数据
· 地面气压 (134.128)
· 平均海平面气压(151.128)
· 地表温度(235.128)
· 2米温度 (167.128, 2m 露点温度168.128)
· 2米相对湿度或比湿
· 10风场u,v分量(165.128, 166.128)  
       如果需要用到masked,建议提供LANDSEA场。  
      如果在进行物理过程参数化时需要用到Noah LSM陆面模式,需要提供两层Soil Temperature (139,170,183,236.128)和Soil       Moisture(39,40,41,42.128)。
      如果模式运行时需要更新SST,需要提供SST(34.128)。
      建议提供Soil Height(可选),这样就可以调整土壤温度达到与WRF模式地形高度向适应。
      建议提供SNOW(141.128场。
      建议对高纬度冬季个例,提供SEAICE (31.128)
http://bbs.06climate.com/forum.php?mod=viewthread&tid=57876数据介绍

3.数据与处理
根据帖子http://bbs.06climate.com/forum.php?mod=viewthread&tid=30962上介绍的方法,

编辑namelist.wps:             vi namelist.wps
运行geogrid:                    ./geogrid.exe

链接Vtable:                     ln -sf ungrib/Variable_Tables/Vtable.ECMWF Vtable
链接高空数据:               ./link_grib.csh ./usr/local/data/ERA/ERA-pl201703.grib
编辑namelist.wps:            prefix = '3D',
运行ungrib:                    ./ungrib.exe
链接地面数据:              ./link_grib.csh ./usr/local/data/ERA/ERA-sfc201703.grib
编辑namelist.wps:           prefix = 'SFC',
运行ungrib:                   ./ungrib.exe
metgrib部分
编辑namelist.wps里metgrid部分: fg_name = '3D','SFC',  
运行                              ./metgrid.exe

但是在./metgrid.exe时,出现错误:
[root@localhost WPS]# ./metgrid.exe
Processing domain 1 of 1
Processing 2017-01-01_00
    3D
    SFC
ERROR: The mandatory field UU was not found in any input data.
Warning: ieee_inexact is signaling
FORTRAN STOP

在/../../WPS文件夹内查看数据,发现3D:2017-01**和SFC:2017-01**数据是已经生成的
在将陆面数据换成FNL数据之后,发现./metgrid.exe时可以跑下去的。纠其原因还要回到最初的数据下载,变量没有下载完全。重新下载多变量数据,然后./metgrid.exe

链接Geopotential.grib 和 Land-sea mask.grib后,
[root@localhost WPS]# ./metgrid.exe
Processing domain 1 of 1
Processing 1989-01-01_12
    PFILE
WARNING: Couldn't open file PFILE:1989-01-01_12 for input.

ERROR: The mandatory field TT was not found in any input data.
Warning: ieee_inexact is signaling
FORTRAN STOP


vi namelist.input
./real.exe
flag_soil_layers read from met_em file is  1
Max map factor in domain 1 =  1.13. Scale the dt in the model accordingly.
p_top_requested =     5000.000   
allowable grid%p_top in data   =     80000.00
   
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE:  <stdin>  LINE:    1137
p_top_requested < grid%p_top possible from data
-------------------------------------------
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE:  <stdin>  LINE:    1137
p_top_requested < grid%p_top possible from data
解决方法:将namelist.input中p_top_requested=5000改为80000.因为我下载的800hpa的探空数据。
之后出现问题:Could not find level of surface
思考了一下,明白了探空数据不能只单独下载某一层数据,而是要从地面下载到你需要的大气层高度
./real.exe
./wrf.exe成功!!!

上面用的是monthly means of daily means数据,换成daily试跑。
./real.exe
d01 2016-01-01_00:00:00 Total post number of sea ice location changes (water to land) =     11
the grid%landmask and soil/veg cats do not match
i,j=           20            4

grid%landmask=    0.000000   
grid%ivgtyp=           17
grid%isltyp=            6
iswater=           17
grid%tslb=    0.000000        0.000000        0.000000        0.000000   
grid%sst=    0.000000   
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE:  <stdin>  LINE:    3273
mismatch_landmask_ivgtyp
-------------------------------------------
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE:  <stdin>  LINE:    3273
mismatch_landmask_ivgtyp
第一直觉是与时间无关的两个物理变量中的land-sea mask.grib没有更新成功有关。查找相关帖子,看到别人出现的类似问题是将

namelist.input中surface_input_source 是3调成1,打开自己的namelist.input发现没有这一设置,手动添加,./real.exe成功!!!


在链接多月份数据的时候总会出现

ERROR: Data not found: 2017-03-01_00:00:00.0000
以为要分别链接,但数据跑起来的时候,只会承认最新链接的数据。
解决方法,直接用*将这些数据全部链接了,肯定还有简单便捷的方法


下载的数据只下载到了500hpa,然后我将p_top_requested=5000(50hpa),改为50000(500hpa)。./real.exe出现错误
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE:  <stdin>  LINE:    5696
troubles, could not find trapping x locations

-------------------------------------------
查看帖子说是"网格设置方面的问题"或者"p_top_requested这一项,检查一下,改一下试试;或者可能是设置的格点超过了数据的范围。有必要提醒一下的是,一般下载的选择的“经纬度”实际上是一个梯形,而wrf模式的区域是一个“正方形”

ERA数据下载可以设置范围,大大减少数据量。但在./metgrid.exe时出现如下错误
WARNING: Field PSFC has missing values at level 200100 at (i,j)=(1,63)
ERROR: Missing values encountered in interpolated fields. Stopping

查看帖子,是驱动范围小于实际范围导致的。
变量名可以到这个网页查找 http://apps.ecmwf.int/codes/grib/param-db



但因为数据下载变量多,
而且勾选起来的话特别麻烦,所以选择批量下载。
python调用ECMWF欧洲中心api下载数据https://blog.csdn.net/ginynu/article/details/54990330







2018-05-18 11-13-22 的屏幕截图.png
2018-05-18 11-13-38 的屏幕截图.png
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-3-27 09:42:18 | 显示全部楼层
本帖最后由 小莹子 于 2019-5-21 22:03 编辑

http://blog.sina.com.cn/s/articlelist_3852831855_0_1.html
https://www.guokr.com/blog/441223/
http://bbs.06climate.com/forum.php?mod=viewthread&tid=24226
下载SRTM数据放到设定的geog_data_path下,
设定geog_data_res = “srtm”,
再在GEOGRID.TBL.ARW中加入interp_option = srtm:four_pt和rel_path= srtm:srtm/
这样就OK了

https://apps.ecmwf.int/datasets/ ... ll-moda/levtype=pl/
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2020-1-2 16:21:28 | 显示全部楼层
请问楼主在遇到这个问题时如何解决?
Processing domain 1 of 1
    LSM:1989-01-01_12
    Z:1989-01-01_12
Processing 1989-01-01_12
    3D
WARNING: Couldn't open file 3D:1989-01-01_12 for input.
    SFC
WARNING: Couldn't open file SFC:1989-01-01_12 for input.
ERROR: The mandatory field TT was not found in any input data.
application called MPI_Abort(MPI_COMM_WORLD, 0) - process 0
[unset]: write_line error; fd=-1 buf=:cmd=abort exitcode=0
:
system msg for write_line failure : Bad file descriptor
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2021-11-18 13:13:55 | 显示全部楼层
ssn_h 发表于 2020-1-2 16:21
请问楼主在遇到这个问题时如何解决?
Processing domain 1 of 1
    LSM:1989-01-01_12

请问这个问题是如何解决的?
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2022-6-17 00:29:20 | 显示全部楼层
您好,我在执行real.exe时出现了如下错误:troubles, could not find trapping x locations。更改ptop和网格依旧报错,求问大神还有什么解决办法嘛?
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2023-5-4 14:58:08 | 显示全部楼层
我也是troubles, could not find trapping x locations。更改ptop和网格依旧报错,请教一下,你解决了吗?谢谢
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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