Help Needed with Observation Nudging in WRF Using Lidar Data

[Erenw]

Hi everyone,

I'm trying to boost the accuracy of my WRF simulation by using observation nudging with some lidar data I've got. Since the OBSGRID.exe tool isn't suitable for this type of data, I took a shot at creating the OBS_DOMAIN101 file myself using Python as per the observation nudging guide.

I've attached the Python script and the output file OBS_DOMAIN101 here. I've run the obs_nudge with what I think is the right namelist setup, but it seems like the obs data isn't making a difference—the results are the same as without nudging, I think no obs data is read.

Is there anyone who can take a look and tell me if I'm doing something wrong? Any help would be really appreciated!

Thanks a lot!
''
&time
auxinput11_interval_s = 1, 1, 1,
auxinput11_end_h = 48, 48, 24,
auxinput11_end_h = 48,

&fdda
obs_nudge_opt = 1, 1, 1,
max_obs = 15000,
fdda_start = 0., 0., 0.,
fdda_end = 99999., 99999., 99999.,
obs_nudge_wind = 1, 1, 1,
obs_coef_wind = 6.0E-4, 6.0E-4, 6.0E-4,
obs_nudge_temp = 1, 1, 1,
obs_coef_temp = 6.0E-4, 6.0E-4, 6.0E-4,
obs_nudge_mois = 1, 1, 1,
obs_coef_mois = 6.0E-4, 6.0E-4, 6.0E-4,
obs_rinxy = 240., 240., 240.,
obs_rinsig = 0.1,
obs_twindo = 0.6666667,0.6666667,0.6666667,
obs_npfi = 10,
obs_ionf = 180, 180, 2,2,
obs_idynin = 0,
obs_dtramp = 40.,
obs_prt_freq = 10, 10, 10,10,
obs_prt_max = 10
obs_ipf_errob = .true.
obs_ipf_nudob = .true.
obs_ipf_in4dob = .true.
obs_ipf_init = .true.
/
''

 

[kwerner]

Hi,
Can you also share your namelist.oa file, namelist.input file, as well as your met_em* files? I'd like to try to run a test using your file and your specific set-up. The met_em* files are likely going to be too large to attach, so see the home page of this forum for instructions on sharing large files. Thanks!

 

[Erenw]

Thanks, I made it work. i correct some format errors like SPACEs strictly adhering the guide. Now i am trying to integrate more observation stations and looking forward to assessing the potential improvements.

 

[kwerner]

That's great news! I'm so glad you are able to move forward. Thanks for the update!