bug: wrf model_interpolate .or. QTY_SURFACE_ELEVATION

[hkershaw-brown]

🐛

Describe the bug

1. Run model_mod_check interpolate QTY_SURFACE_ELEVATION with a location in m
2. What was the expected outcome?
   successful interpolation
3. What actually happened?
   Error code 2

Error Message

 fld 2 is missing F
member     1, ERROR with error code ::            2

Which model(s) are you working with?

wrf

Screenshots

If applicable, add screenshots to help explain your problem.

Start of the if statement:

DART/models/wrf/model_mod.f90

Line 1169 in a32ab03

if(is_vertical(location,"LEVEL")) then

if(is_vertical(location,"LEVEL")) then 
elseif(is_vertical(location,"PRESSURE")) then
elseif(is_vertical(location,"HEIGHT"))
elseif(is_vertical(location,"SURFACE") .or. obs_kind == QTY_SURFACE_ELEVATION) then
elseif(is_vertical(location,"UNDEFINED")) then
endif

This if .or. QTY_SURFACE_ELEVATION will not be entered for level, pressure, height

DART/models/wrf/model_mod.f90

Line 1242 in a32ab03

elseif(is_vertical(location,"SURFACE") .or. obs_kind == QTY_SURFACE_ELEVATION) then

which gives you a failing model_interpolate even though the interpolation has been done for the single level.

Version of DART

Which version of DART are you using?
You can find the version using git describe --tags
v11.19.0-34-g13dec7c94

Have you modified the DART code?

Yes

Build information

Please describe:

1. intel m
2. GNU Fortran (MacPorts gcc13 13.3.0_2+stdlib_flag) 13.3.0

I believe this came in 49c0dd6#diff-e12d14ed6b47c2ee908ccf4ad30e45a718ebd9e26297e77254a978e8804d7403 same as issue #1013 not setting surf_var