关于地面观测资料转little_r格式的疑问

schLltt

小虹 发表于 2014-12-9 21:48
1层的不影响的，建议你先写1个观测，运行成功了再写多个观测的，good luck

好的，thks!

kzpyezi

正好需要这个东西

kzpyezi

正好需要，金钱不够好忧伤

追风的阳光

甜果果 发表于 2014-12-2 16:04
是micaps2的数据吧？那如果转换1的 是不是把相应的变量更换就可以了啊？

请问用micaps，dimond1数据转换为little_R的问题你解决了吗？

甜果果

没有 后来没用wrfda了

追风的阳光

小虹 发表于 2014-12-9 21:48
1层的不影响的，建议你先写1个观测，运行成功了再写多个观测的，good luck

lz能把obsproc.namelist给我参考下吗？另想问一下你的背景协方差文件是怎么求的？

小虹

追风的阳光 发表于 2015-5-12 16:40
lz能把obsproc.namelist给我参考下吗？另想问一下你的背景协方差文件是怎么求的？

&record1
obs_gts_filename = 'ob_tem_afqc20%-kv15',
obs_err_filename = 'obserr.txt',
fg_format='WRF',
/

&record2
time_window_min  = '2008-01-30_17:00:00.0000',
time_analysis    = '2008-01-30_18:00:00.0000',
time_window_max  = '2008-01-30_19:00:00.0000',
/

&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         = .TRUE.,
domain_check_h           = .true.,
Thining_SATOB            = false,
Thining_SSMI             = false,
Thining_QSCAT            = false,
/

&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  = 34.183,
XLONC = 106.016,
TRUELAT1= 30.0,
TRUELAT2= 60.0,
MOAD_CEN_LAT =34.183,
STANDARD_LON = 106.016,
/

&record8
IDD    =   1,
MAXNES =   1,
NESTIX =  200,  200,  136,  181,  211,
NESTJX =  160,  200,  181,  196,  211,
DIS    =  30000,  10.,  3.3,  1.1,  1.1,
NUMC   =    1,    1,   2,     3,    4,
NESTI  =    1,   40,  28,    35,   45,
NESTJ  =    1,   60,  25,    65,   55,
/

&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.,
/                                                                                                                                                

背景误差协方差是利用NMC方法做了一个月的预报统计出来的，具体可以详见WRF的说明书

追风的阳光

小虹 发表于 2015-5-12 21:46
&record1
obs_gts_filename = 'ob_tem_afqc20%-kv15',
obs_err_filename = 'obserr.txt',

同化小白好多不懂，我看了下，1如何要同化多种资料的话，要在Litter_R中多些几个report，后面我发现Litter_R转换程序的输出文件名如果一样的话，每次的输出就会加到原先的输出就会加到原先的输出资料上去，不知道这样可不可行呢？2.我看了一下你的是用了一个时次的观测资料，如果我需要写入多个观测资料的话，namelist中的资料路径是否只写成资料路径就可以了？我看了教程里4dvar的测试数据，里面只输入了一个时次的观测资料，最后却生成了整个事件窗的6个obs_gets文件，这是为什么呢？

小虹

追风的阳光 发表于 2015-5-13 09:31
同化小白好多不懂，我看了下，1如何要同化多种资料的话，要在Litter_R中多些几个report，后面我发现Litte ...

1同化多种资料，将多种资料写到一个little-r文件中即可。
2.需要写入多个观测资料，可以将多个时次的观测资料写入一个little-r文件，然后改动namelist里面的时间，可以得到不同时次的obs_gts_*.3DVAR文件。
3 4DVAR的同化不懂，相关知识不知道，抱歉。

追风的阳光

小虹 发表于 2015-5-13 16:12
1同化多种资料，将多种资料写到一个little-r文件中即可。
2.需要写入多个观测资料，可以将多个时次的观 ...

我后来改好Little_R后，obsproc.exe,提示空间不够，请问你碰到过这个情况吗？