Scheduled Downtime
On Friday 21 April 2023 @ 5pm MT, this website will be down for maintenance and expected to return online the morning of 24 April 2023 at the latest

Radiance Data Assimilation in 4dvar

This post was from a previous version of the WRF&MPAS-A Support Forum. New replies have been disabled and if you have follow up questions related to this post, then please start a new thread from the forum home page.

Wilson

New member
Hello, in WRFDA (V3.8.1), I use 4DVAR to run Radiance Data Assimilation.Although I have opened the deviation correction, the deviation correction does not work.OMB With BC is the same as OMB No BC.But when I use 3d variational time, the deviation correction works fine.
&wrfvar1
var4d=true,
var4d_lbc=false,
var4d_bin=3600,
var4d_bin_rain=21600,
print_detail_obs=.true.,
print_detail_grad=.true.,
/
&wrfvar2
/
&wrfvar3
num_fgat_time=7,
/
&wrfvar4
use_amsuaobs=.true.
/
&wrfvar5
check_max_iv=true,
/
&wrfvar6
max_ext_its=1,
ntmax=1,
orthonorm_gradient=true,
/
&wrfvar7
cv_options=3
as1(3) = 0.1,0.4,0.6
as2(3) = 0.1,0.4,0.6
as3(3) = 0.1,0.4,0.6
as4(3) = 0.1,0.4,0.6
as5(3) = 0.1,0.4,0.6
/
&wrfvar8
/
&wrfvar9
/
&wrfvar10
test_transforms=false,
test_gradient=false,
/
&wrfvar11
/
&wrfvar12
/
&wrfvar13
/
&wrfvar14
rtminit_nsensor=1,
rtminit_platform=1
rtminit_satid=19
rtminit_sensor=3
thinning=true,
thinning_mesh=60.0,
qc_rad=false,
write_iv_rad_ascii=true,
write_oa_rad_ascii=true,
rtm_option=2,
only_sea_rad=false,
crtm_cloud=true,
use_varbc=true,
use_crtm_kmatrix =true,
freeze_varbc=false,
varbc_factor=1,
/
&wrfvar15
/
&wrfvar16
/
&wrfvar17
/
&wrfvar18
analysis_date=2015-08-03_06:00:00.0000,
/
&wrfvar19
/
&wrfvar20
/
&wrfvar21
time_window_min=2015-08-03_06:00:00.0000,
/
&wrfvar22
time_window_max=2015-08-03_12:00:00.0000,
/
&wrfvar23
/
&time_control
run_hours=6,
start_year=2015,
start_month=08,
start_day=03,
start_hour=06,
end_year=2015,
end_month=08,
end_day=03,
end_hour=12,
interval_seconds=21600,
/
&fdda
/
&domains
time_step=180,
e_we=121,
e_sn=121,
e_vert=42,
num_metgrid_levels=32,
p_top_requested=1000,
interp_type=2,
t_extrap_type=2,
dx=30000,
dy=30000,
/
&dfi_control
/
&tc
/
&physics
mp_physics=2,
ra_lw_physics=1,
ra_sw_physics=1,
radt=30,
sf_sfclay_physics=1,
sf_surface_physics=2,
bl_pbl_physics=1,
cu_physics=1,
cudt=0,
num_soil_layers=4,
mp_zero_out=2,
co2tf=0,
/
&scm
/
&dynamics
w_damping=1,
diff_opt=1,
km_opt=4,
dampcoef=0.2,
time_step_sound=4,
base_temp=290.0,
use_baseparam_fr_nml=true,
/
&bdy_control
specified=true,
real_data_init_type=3,
/
&grib2
/
&fire
/
&namelist_quilt
/
&perturbation
trajectory_io=true,
enable_identity=false,
jcdfi_use=false,
jcdfi_diag=1,
jcdfi_penalty=1000.0,
var4d_detail_out=true,
/
 
hello.
Is this problem solved?
When I run the 3Dvar, I got the ‘VARBC: 0 active observations’ line in the rsl.out.0000 file.
 
Top