Re: Axis units in wrf_user_intrp3d()

I think (I heard from someone) that '0-100' means 'starting point-ending
point'.
Linear interpolation routine makes the total length 100 (just arbitrary).
For example, if your vertical domain is 0-1000km and you pick up the
domain of interpolated variable a(49:99) out of a(0:99), it covers from
500 km to 1000km. You need to change interpolated vertical coordinate
fit to your case. Although you started from sigma coordinate,
linear-interpolation routine ends up with height coordinate. I played
with idealized WRF case. Sorry, if it is not a good example.

Si-Wan

Gerry Creager wrote:
> Hmmm. The '0-100' sounds like a variation on sigma levels.
>
> IIRC, there's a wrf_Height.ncl example that should allow you to
> generate what you're looking for. I've gotta get "good" with it but
> it's in the example docs.
>
> gerry
>
> Don Morton wrote:
>> Thanks all, I'll try to come up with a very simple
>> case study tonight or tomorrow and submit it.
>>
>> I have 75 levels in my model, so the 0-100 on
>> the left doesn't jive with that. I did do the printVarSummary(z)
>> and saw my units of "m"
>>
>> Again though, rather than confusing everyone, I'll just make
>> up a simple example tonight.
>>
>> On 9/25/07, Mary Haley <haley_at_ucar.edu> wrote:
>>> Don,
>>>
>>> The units that are appearing in the top right corner of your
>>> plot come from the metadata that is attached to
>>> the "z" array that you are plotting.
>>>
>>> If you type "printVarSummary(z)" (or whatever the name of your
>>> variable is), you will see what the units are.
>>>
>>> Perhaps the "units" are not being reset after the calculation is done?
>>> I'm not familiar enough with wrf_user_intrp3d to know if it resets the
>>> units, but I'll include wrfhelp_at_ucar.edu on this one to see if they
>>> can include more information.
>>>
>>> --Mary
>>>
>>>
>>> On Tue, 25 Sep 2007, Don Morton wrote:
>>>
>>>> Howdy,
>>>>
>>>> First, I'd really like to thank you all for the great suite
>>>> of routines for the WRF output files. These are turning
>>>> out to be very useful.
>>>>
>>>> I've created a plot that shows a vertical cross section
>>>> of the u-component of the wind over a defined transect
>>>> (graphic attached). Just in case you're curious, the transect
>>>> is roughly west-to-east, the west end being at Missoula
>>>> International Airport and the east end being in Hellgate
>>>> Canyon on the east side of Missoula, where Interstate
>>>> 90 enters the city. I'm trying to assess WRF's capture
>>>> of the famous easterly "Hellgate Winds" that we often
>>>> get in the morning (sort of a "drainage flow" from the
>>>> east), and you can see a hint of it in this graphic.
>>>>
>>>> But, on to the question - I don't understand what the units
>>>> are on the vertical axis. I've passed in an array of "z" values
>>>> to the function, and I've verified that the values I pass in
>>>> are reasonable, with units of meter. I "suppose" the units
>>>> could be km, and it looks like that would be an appropriate
>>>> unit on the horiz axis, too (the distance is about right).
>>>> But, I'm wondering if someone can verify.
>>>>
>>>> Thanks,
>>>>
>>>> Don
>>>> --
>>>> Arctic Region Supercomputing Center
>>>> http://weather.arsc.edu/
>>>>
>>
>>
>
Received on Tue Sep 25 2007 - 13:12:57 MDT