[ncl-talk] get sensible heat(HFX) values from wrfout at specific height.

Bill Ladwig ladwig at ucar.edu
Mon Jul 25 10:02:39 MDT 2016


Hi Ananya,

While this may be different for LES runs, I'm pretty sure that the HFX
variable is still going to be the upward heat flux at the surface, which is
a 2D variable.  In order to do any vertical interpolation, you would need a
three dimensional heat flux variable.  You're trying to use a
one-dimensional interpolation routine, which is really meant for one
vertical column in a three-dimensional variable.  If you had a 3D heat flux
variable, you would want to use wrf_interp_3d_z or wrf_user_intrp3d.

However, since all you have is the 2D heat flux at the surface, I think
you're going to have to resort to the radiation equations for calculating
the 20 m, 30 m, and 100 m values.  Off the top of my head, I don't think
NCL has anything built in to do this for you, so you'll probably have to
implement this yourself.  It's been a long time since I took radiation, so
if all this is new to you, start with the concept of 'optical depth' and go
from there.  Maybe others have better suggestions here.

Bill

On Mon, Jul 25, 2016 at 1:51 AM, ananya sritawan <ananbeam76 at gmail.com>
wrote:

> Dear all
>
> I have been working on WRF-LES and now I am trying to get a data at a
> specific height says 20m 30m 100m.
> The data that I want to obtain is sensible heat flux from my wrfout at 20m
> 30m and 100m. Can we do that? Therefore I tried to look for the ncl script
> to do so, I have tried many set of interpolate script but I cannot
> accomplished my goal. Then I found one script which have been successfully
> run by other user, then I modified the script according to my need with my
> very little skill of coding. When I ran the script, I've got an error about
> the dimension. something like this : (when I use  wrf_interp_1d)
>
> fatal:wrf_interp_1d: The v_in, z_in, and z_out arrays must be the same
> number of dimensions
>
> fatal:["Execute.c":8575]:Execute: Error occurred at or near line 51 in
> file fillvalue100
>
> which I have tried to check how to make it into the same number of
> dimension, but I didn't found one although I still quite not sure its
> dimension process. I also attached the script and my nameless input
> below(in the script "fillvalue100" for getting data at 100m).
>
> I am working with ncl-6.3.0 on Darwin Kernel mac os x
>
> I'm a newbie for both wrf and ncl, so any suggestion would be greatly
> appreciate. I really need help.
>
> Best regard,
>
>
> --
> Ananya SRITAWAN
> Hydro-Environment Laboratory
> Graduate School of Environmental Studies, Tohoku University.
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20160725/89dcddd6/attachment.html 


More information about the ncl-talk mailing list