Re: meteogram

Zhan,

As Adam correctly stated, you generally get this message if you are setting the sf/vfX/YArray resources to arrays that are not equal in size to your original data array.

Another possibility is that internally a cyclic point may be getting added to your longitude dimension, although I didn't think the WRF scripts did this.

You can try this:

   uv_res@gsnAddCyclic = False

and ditto for the contour resource list.

--Mary

On Sep 7, 2010, at 1:46 PM, Adam Phillips wrote:

> I am not terribly familiar with WRF output. Does zz (z_point) actually
> vary with time? If not, then zz could be set up as a 1D array as in the
> following:
> zz = z_point(bottom_top|:,Time|0)
>
> If zz does vary with time then both vfYArray and vfXArray must have the
> same dimension size. (bottom_top x Time) You will then want to set
> vfYArray = zz as you are doing now. However, vfXArray is set to be equal
> to taus which is a 1D array of times. You will need to convert taus into
> a 2D array of the same sizes as z_point/zz. (The times will be the same
> for each level.) You can use conform to do this, or set up a quick do loop.
> Adam
>
>
> On 09/07/2010 01:23 PM, Zhan Zhao wrote:
>> Hi Adam,
>>
>> Thanks for your reply. This did not work, and I got error message:
>>
>> warning:ScalarFieldSetValues: 2d coordinate array sfYArray has an
>> incorrect dimension size: defaulting sfYArray
>> warning:VectorFieldSetValues: 2d coordinate array vfYArray has an
>> incorrect dimension size: defaulting vfYArray
>>
>> So maybe there are some problems in the way zz is derived?
>>
>> z = wrf_user_getvar(a,"z",-1) ;
>> z_point = z(:,:,locY,locX)
>>
>> ; Swap the dimensions as we want to plot time on the X axis later
>> tt = tc_point(bottom_top|:,Time|:)
>> ugrid = u_point(bottom_top|:,Time|:)
>> vgrid = v_point(bottom_top|:,Time|:)
>> zz = z_point(bottom_top|:,Time|:)
>>
>> Zhan
>>
>>
>>
>> On Tue, Sep 7, 2010 at 12:12 PM, Adam Phillips <asphilli@ucar.edu
>> <mailto:asphilli@ucar.edu>> wrote:
>>
>> Hi Zhan,
>> You are setting the vfXArray resource, as well as the sfYArray resource.
>> I think you should change this line:
>> uv_res@sfYArray = zz
>> to this:
>> uv_res@vfYArray = zz
>>
>> Let the group know if that doesn't work...
>> Adam
>>
>> On 09/07/2010 12:54 PM, Zhan Zhao wrote:
>>> Hi NCL group,
>>>
>>> I am trying to make meteogram plot following the example
>> wrf_meteo_5.ncl
>>>
>> <http://www.mmm.ucar.edu/wrf/OnLineTutorial/Graphics/NCL/Examples/SPECIAL/wrf_meteo_5.ncl>
>> on
>>>
>>>
>> http://www.mmm.ucar.edu/wrf/OnLineTutorial/Graphics/NCL/Examples/SPECIAL/meteograms.htm.
>>>
>>> I'd like to use Height as the y axis for the vertical profile
>> plot, so I
>>> made some changes to this example (the attachment). Seems the changes
>>> did not work, the y axis is still the vertical level indexes
>> instead of
>>> the height. Could you please give me some suggestions on this?
>>>
>>> Thanks a lot!
>>>
>>> Zhan
>>>
>>>
>>>
>>>
>>> _______________________________________________
>>> ncl-talk mailing list
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>> --
>> __________________________________________________
>> Adam Phillips
>> asphilli@ucar.edu <mailto:asphilli@ucar.edu>
>> National Center for Atmospheric Research tel: (303) 497-1726
>> Climate and Global Dynamics Division fax: (303) 497-1333
>> P.O. Box 3000
>> Boulder, CO 80307-3000 http://www.cgd.ucar.edu/cas/asphilli
>> _______________________________________________
>> ncl-talk mailing list
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>>
>
> --
> __________________________________________________
> Adam Phillips
> asphilli@ucar.edu
> National Center for Atmospheric Research tel: (303) 497-1726
> Climate and Global Dynamics Division fax: (303) 497-1333
> P.O. Box 3000
> Boulder, CO 80307-3000 http://www.cgd.ucar.edu/cas/asphilli
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri Sep 10 12:58:02 2010