学习WRFDA安装测试官网观测数据（这里还没有加入辐射数据同化）

任性小王子

Samantha 发表于 2025-2-20 09:01
还有想请教楼主一个WRFDA的问题。就是我在做地面观测数据前处理的时候，下载了区域范围内的数据处理成lit ...

我没使用./obsproc.exe处理过,不过可以对照https://www2.mmm.ucar.edu/wrf/us ... _namelist_variables看一看

任性小王子

任性小王子 发表于 2025-2-20 16:45
抱歉这么晚回复，export WRFPLUS_DIR=${your_source_code_dir}/WRFPLUS，这个执行了吗？再试一下编译，如 ...

https://www2.mmm.ucar.edu/wrf/us ... l/Help/littler.html

Samantha

export已经加到bashrc里了，重新又试了一遍，还是同样的错误

任性小王子

Samantha 发表于 2025-2-20 21:51
export已经加到bashrc里了，重新又试了一遍，还是同样的错误

抱歉这么晚回复，安装的时候把你的辐射传输模块的路径也加进去，重新编一下。

任性小王子

Samantha 发表于 2025-2-20 21:51
export已经加到bashrc里了，重新又试了一遍，还是同样的错误

这是我的编译文件开头，你缺少辐射模块呐。

任性小王子

Samantha 发表于 2025-2-20 09:01
还有想请教楼主一个WRFDA的问题。就是我在做地面观测数据前处理的时候，下载了区域范围内的数据处理成lit ...

我之前跑过一个示例，是看视频（youtube)操作的，如果可以，你在youtube上搜一下wrfda，能搜到的那个示例。这是namelist.obsproc
&record1
obs_gts_filename = 'SURFACE_OBS:2021040500',
obs_err_filename = 'obserr.txt',
gts_from_mmm_archive = .true.,
/

&record2
time_window_min  = '2021-04-04_23:00:00',
time_analysis    = '2021-04-05_00:00:00',
time_window_max  = '2021-04-05_01:00:00',
/

&record3
max_number_of_obs        = 400000,
fatal_if_exceed_max_obs  = .TRUE.,
/

&record4
qc_test_vert_consistency = .TRUE.,
qc_test_convective_adj   = .TRUE.,
qc_test_above_lid        = .TRUE.,
remove_above_lid         = .false.,
domain_check_h           = .true.,
Thining_SATOB            = .false.,
Thining_SSMI             = .false.,
Thining_QSCAT            = .false.,
calc_psfc_from_qnh       = .true.,
/

&record5
print_gts_read           = .TRUE.,
print_gpspw_read         = .TRUE.,
print_recoverp           = .TRUE.,
print_duplicate_loc      = .TRUE.,
print_duplicate_time     = .TRUE.,
print_recoverh           = .TRUE.,
print_qc_vert            = .TRUE.,
print_qc_conv            = .TRUE.,
print_qc_lid             = .TRUE.,
print_uncomplete         = .TRUE.,
/

&record6
ptop =  1000.0,
base_pres       = 100000.0,
base_temp       = 290.0,
base_lapse      = 50.0,
base_strat_temp = 215.0,
base_tropo_pres = 20000.0
/

&record7
IPROJ = 1,
PHIC  = 33.00,
XLONC = -79.00,
TRUELAT1= 30.0,
TRUELAT2= 60.0,
MOAD_CEN_LAT = 33.00,
STANDARD_LON = -79.00,
/

&record8
IDD    =   1,
MAXNES =   1,
NESTIX =  15000,
NESTJX =  15000,
DIS    =  1,
NUMC   =    1,
NESTI  =    1,
NESTJ  =    1,
/

&record9
PREPBUFR_OUTPUT_FILENAME = 'prepbufr_output_filename',
PREPBUFR_TABLE_FILENAME = 'prepbufr_table_filename',
OUTPUT_OB_FORMAT = 2
use_for          = '3DVAR',
num_slots_past   = 3,
num_slots_ahead  = 3,
write_synop = .true.,
write_ship  = .true.,
write_metar = .true.,
write_buoy  = .true.,
write_pilot = .true.,
write_sound = .true.,
write_amdar = .true.,
write_satem = .true.,
write_satob = .true.,
write_airep = .true.,
write_gpspw = .true.,
write_gpsztd= .true.,
write_gpsref= .true.,
write_gpseph= .true.,
write_ssmt1 = .true.,
write_ssmt2 = .true.,
write_ssmi  = .true.,
write_tovs  = .true.,
write_qscat = .true.,
write_profl = .true.,
write_bogus = .true.,
write_airs  = .true.,
/
这是wrfda中namelist.input

&wrfvar1
var4d=false,
print_detail_grad=false,
/
&wrfvar2
/
&wrfvar3
ob_format=1,
/
&wrfvar4
/
&wrfvar5
/
&wrfvar6
max_ext_its=1,
ntmax=50,
orthonorm_gradient=true,
/
&wrfvar7
cv_options=3,
/
&wrfvar8
/
&wrfvar9
/
&wrfvar10
test_transforms=false,
test_gradient=false,
/
&wrfvar11
/
&wrfvar12
/
&wrfvar13
/
&wrfvar14
/
&wrfvar15
/
&wrfvar16
/
&wrfvar17
/
&wrfvar18
analysis_date="2021-04-05_00:00:00.0000",
/
&wrfvar19
/
&wrfvar20
/
&wrfvar21
time_window_min="2021-04-04_23:00:00.0000",
/
&wrfvar22
time_window_max="2021-04-05_01:00:00.0000",
/
&time_control
start_year=2021,
start_month=04,
start_day=05,
start_hour=00,
end_year=2021,
end_month=04,
end_day=05,
end_hour=06,
/
&fdda
/
&domains
e_we=150,
e_sn=130,
e_vert=45,
dx=15000,
dy=15000,
/
&dfi_control
/
&tc
/
&physics
mp_physics=-1,
ra_lw_physics=-1,
ra_sw_physics=-1,
radt=15,
sf_sfclay_physics=-1,
sf_surface_physics=-1,
bl_pbl_physics=-1,
cu_physics=-1,
cudt=0,
num_soil_layers=4,
mp_zero_out=2,
co2tf=0,
/
&scm
/
&dynamics
/
&bdy_control
/
&grib2
/
&fire
/
&namelist_quilt
/
&perturbation
/
这是namelist

Samantha

任性小王子 发表于 2025-2-24 10:18
这是我的编译文件开头，你缺少辐射模块呐。

感谢楼主！我看了下wrfda有自带crtm辐射模块的所以没有安装，想问下必须安装rttov吗？

Samantha

任性小王子 发表于 2025-2-24 10:41
我之前跑过一个示例，是看视频（youtube)操作的，如果可以，你在youtube上搜一下wrfda，能搜到的那个示例 ...

这个问题已解决，是我之前调试其他观测数据的时候增加了一行把数据关掉了，谢谢楼主分享

任性小王子

Samantha 发表于 2025-2-25 10:53
感谢楼主！我看了下wrfda有自带crtm辐射模块的所以没有安装，想问下必须安装rttov吗？

抱歉，这两天一直在外面，身边没电脑。嗯嗯，得有辐射模块，crtm或者rttov 得有，把环境路径加上。

任性小王子

Samantha 发表于 2025-2-25 17:52
这个问题已解决，是我之前调试其他观测数据的时候增加了一行把数据关掉了，谢谢楼主分享

客气啦，客气啦，已解决就好{:eb511:}