爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 19428|回复: 27

wrf同化卫星辐射率资料出现问题

[复制链接]
发表于 2012-5-9 10:16:42 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 小虹 于 2012-5-9 15:28 编辑

大家好,我是一个wrf初学者,最近在学习卫星辐射率同化,按照说明书上转wrf提供的例子,按照下列步骤进行

1) ln -sf /home/name/wrf3.2.1/WRFDA/run/LANDUSE.TBL ./LANDUSE.TBL
2) ln -sf /home/name /wrf3.2.1/WRFDA/var/run/be.dat.cv3 ./be.dat
3) ln -sf /home/name /wrf3.2.1/WRFDA/var/run/radiance_info ./radiance_info
4) ln -sf /home/name/wrf3.2.1/WRFDA/var/run/VARBC.in ./VARBC.in
5)ln -sf /home/name/wrf3.2.1/WRFDA/var/run/crtm_coeffs ./crtm_coeffs
6)cp /home/name/wrf3.2.1/WRFDA/var/test/radiance/namelist.input ./namelist.input
7) ln -sf /home/name/wrf3.2.1/WRFDA/var/da/da_wrfvar.exe ./da_wrfvar.exe




然后将它提供的wrfinput_d01复制到此目录下,改名为fg,将amsua.bufr,amsub.bufr也复制到此目录下,这些文件都用的是testdata,绝对没问题的
namelist.input设置如下

&wrfvar1
print_detail_grad=.false.,
/
&wrfvar2
/
&wrfvar3
ob_format=1
&wrfvar4
use_amsuaobs=.true.
use_amsubobs=.true.
/
&wrfvar5
/
&wrfvar6
/
&wrfvar7
/
&wrfvar8
/
&wrfvar9
/
&wrfvar10
/
&wrfvar11
/
&wrfvar12
/
&wrfvar13
/
&wrfvar14
rtminit_nsensor=6,
rtminit_platform=1,1,1,1,1,1,
rtminit_satid=15,16,18,15,16,17,
rtminit_sensor=3,3,3,4,4,4,
thinning_mesh=120.0,120.0,120.0,120.0,120.0,120.0,
thinning=true,
qc_rad=true,
write_iv_rad_ascii=false,
write_oa_rad_ascii=true,
rtm_option=2,
only_sea_rad=false,
use_varbc=true,
crtm_atmosphere=6,
/
&wrfvar15
/
&wrfvar16
/
&wrfvar17
/
&wrfvar18
analysis_date="2008-02-05_12:00:00.0000",
/
&wrfvar19
/
&wrfvar20
/
&wrfvar21
time_window_min="2008-02-05_11:00:00.0000",
/
&wrfvar22
time_window_max="2008-02-05_13:00:00.0000",
/
&wrfvar23
/
&time_control
start_year=2008,
start_month=02,
start_day=05,
start_hour=12,
end_year=2008,
end_month=02,
end_day=05,
end_hour=12,
/
&fdda
/
&domains
e_we=90,
e_sn=60,
e_vert=41,
dx=60000,
dy=60000,
/
&dfi_control
/
&tc
/
&physics
mp_physics=3,
ra_lw_physics=1,
ra_sw_physics=1,
radt=60,
sf_sfclay_physics=1,
sf_surface_physics=1,
bl_pbl_physics=1,
cu_physics=1,
cudt=5,
num_soil_layers=5,
mp_zero_out=2,
co2tf=0,
/
&scm
/
&dynamics
/
&bdy_control
/
&grib2
/
&fire
/
&namelist_quilt
/


运行 ./da_wrfvar.exe   rsl.out.0000 file  最后出现如下错误

**************BUFR ARCHIVE LIBRARY ABORT*****************
BUFRLIB: OPENBF - ERROR READING INPUT FILE CONNECTED TO UNIT  99 WHEN CHECKING
FOR 'BUFR' IN FIRST 4 BYTES OF RECORD            
**************BUFR ARCHIVE LIBRARY ABORT*****************



请教大家是怎么回事啊?有谁知道的告诉我一声吧,非常感谢啦~


密码修改失败请联系微信:mofangbao
 楼主| 发表于 2012-5-13 11:18:48 | 显示全部楼层
其实还是testdata格式和自己机器格式不一致造成的,从ssrc.c重新将卫星原始数据处理了一下,再run就可以了,给遇到同样错误的同学提个醒~
密码修改失败请联系微信:mofangbao
发表于 2012-5-31 21:46:48 | 显示全部楼层
怎么处理阿?介绍下。
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2012-6-10 16:18:47 | 显示全部楼层

就是一个可以将卫星格式进行转化的工具,和cword一样的功能 ssrc.c (3.74 KB, 下载次数: 68)
密码修改失败请联系微信:mofangbao
发表于 2012-6-6 11:52:23 | 显示全部楼层
ssrc.c是什么文件,在哪里?
密码修改失败请联系微信:mofangbao
发表于 2012-6-12 13:34:27 | 显示全部楼层
正在做卫星资料同化工作,一起学习
密码修改失败请联系微信:mofangbao
发表于 2012-11-9 10:29:28 | 显示全部楼层
我也是在做卫星资料同化工作,一起学习
密码修改失败请联系微信:mofangbao
发表于 2013-5-9 12:47:41 | 显示全部楼层
通过ssrc.c进行编辑了 还是不对
密码修改失败请联系微信:mofangbao
发表于 2013-5-14 22:40:29 | 显示全部楼层
是不是ARW WRF手册上说的,对BUFR数据先要用一个C程序装换,将本来在IBM平台上编码的数据改变成我们自己平台上的数据格式啊。
俺看这个手册好几遍了,还是不会跑,郁闷中。。。
密码修改失败请联系微信:mofangbao
发表于 2014-12-14 19:28:01 | 显示全部楼层
本帖最后由 lee309 于 2014-12-14 19:29 编辑
小虹 发表于 2012-5-13 11:18
其实还是testdata格式和自己机器格式不一致造成的,从ssrc.c重新将卫星原始数据处理了一下,再run就可以了 ...

问下怎么生成自己的be.dat文件?也就是 background error statistics file

密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

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

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

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