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

Surface data skipped in wrfvar_output

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.


New member
I would like to ask your assistance for my 3D-VAR data assimilation experiment. All elements of my experiment seem ok, the log files tell appropriate work, and wrfvar_output generated correctly.
I have checked the increments (difference between first guess and 3D-VAR output file) and they also seem correct on higher level. BUT on surface level something may be wrong. :cry:

Considering T2m there is slight and singular increment field, moreover U,V fields are identical in both files (no increments), however Jo calculation for synop_u and synop_v appears in log files:

jo_synop 0.211335854820E+03
jo%synop_u 0.107599699190E+03
jo%synop_v 0.186768876058E+01
jo%synop_t 0.101610705821E+03
jo%synop_p 0.257418691903E+00
jo%synop_q 0.342357448447E-03
jo_airep 0.362222449565E+02
jo%airep_u 0.113257960778E+02
jo%airep_v 0.000000000000E+00
jo%airep_t 0.248964488787E+02
jo%airep_q 0.000000000000E+00
jo_profiler 0.346260589824E+02
jo%profiler_u 0.261924596285E+02
jo%profiler_v 0.843359935392E+01

My first idea was that surface observation assimilation results are not written out, but I did not find this option in namelists. What could be wrong?
I have put FG, namelist, and output files together:

Thank you for your help in advance,
I am having the same issue (perhaps). Fields like T2 and U10 are not being updated. But the model levels seem to be reasonably updated. How can I get the derived near-surface fields to update?
Please see the description about update_sfcdiags in var/README.namelist.

It is better to let the model diagnose and output those surface variables from the first couple time steps.
For example, in WRF namelist.input,
auxhist6_outname = "./surface_d<domain>_<date>"
iofields_filename = "./surface_out.txt"
io_form_auxhist6 = 2
auxhist6_begin_m = 3 ! set a proper value according to your time step
auxhist6_end_m = 6 ! set a proper value according to your time step

Edit surface_out.txt file to include the following one line.