[ncl-talk] Wind speed interpolation about ground level

[Dennis Shea]

Really, somebody on the wrfforum or wrfhelp at ucar.edu should answer
this... not generic ncl-talk. They are the WRF data experts.

---
I *speculate* that the 'height_levels' to which you are interpolating
are relative to sea-level. They are not
'height_levels_above_the_current_terrain_elevation' or
'height_levels_above_the_current_model_height'.

Try:

height_levels = (/ 50., 80., 1000, 2000. /)

I *speculate* that for 1000, you will get the eastern half of the land
filled out and for 2000 most of the area will be filled.

---

Subsequently, maybe the following would work

  thgt = wrf_user_getvar(a,"ter",-1)   ; terrain height
or
  mhgt = wrf_user_getvar(a,"z",-1)    ; full model height

Then in the 'height' loop add the 'thgt' or 'zhgt' to the current grid
point value.

HTH



On Tue, Dec 15, 2015 at 4:07 AM, mehmet ali tor <torm732 at gmail.com> wrote:
> Hello,
>
> Iam using a script which enables interpolate wind speed 50m and 80m. My
> problem is that my output shows only sea surface areas. Land surfaces are
> looking blank. When i make interpolation 50m and 2000m, i can see datas on
> the land surface at 2000m map. 50m map don't show datas on land area.
>
> My code hasn't got masking variables. I think, the problem is about wrong
> vertical coordinat calculation about only sea level. I want to revise my
> code about terrain height. I tried some fixing codes but they are not
> helpfully.
>
> You can look my code this question link =
> http://forum.wrfforum.com/viewtopic.php?f=32&t=9340
>
> My code is at the code panel which is under the question.
>
> Thank you all,
>
> Regards.
>
> --
> Mehmet Ali TOR
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>