Re: Wrfout wind profiler scripting

From: Tardif, Elliot M <elliot.tardif_at_nyahnyahspammersnyahnyah>
Date: Wed Sep 14 2011 - 15:16:55 MDT

Hi Adam, thanks very much for your response! After some further analysis, I think using pressure as my y-axis variable would be more suitable. The big challenge is how to fit the plotting of the pressure variable into the script. I have the plot of the winds and wind speed ("windlayer = gsn_vector_scalar(wks,ugrid(bottom_top|:,Time|:),vgrid(bottom_top|:,Time|:),wspd(bottom_top|:,Time|:),uv_res)"), but there's no room in this particular program to plot the wind barbs on pressure levels across the time period, instead of on sigma levels. I want to plot them like this to show the variation in pressure for a given sigma level across the time period. Will an overlay of some type work in this case?

Thanks,

Elliot Tardif.

Elliot Tardif, Meteorologist I
NC DENR, Division of Air Quality
Planning Section, Attainment Planning Branch
1641 Mail Service Center
Raleigh, NC 27699-1641
Phone:  919-715-9773
Fax  :  919-715-7476
Email:  Elliot.Tardif@ncdenr.gov
Web  :  http://www.ncair.org

*****************************************************************************************
Email correspondence to and from this address is subject to the North Carolina Public Records Law and may be disclosed to third parties unless the content is exempt by statute or other regulation.
*****************************************************************************************

-----Original Message-----
From: ncl-talk-bounces@ucar.edu [mailto:ncl-talk-bounces@ucar.edu] On Behalf Of Adam Phillips
Sent: Tuesday, September 13, 2011 6:23 PM
To: ncl-talk@ucar.edu
Subject: Re: Wrfout wind profiler scripting

Hi Elliot.
Looking at your script I don't see an array that would house the height
in meters of each sigma level. Sigma levels vary in height, so I'm not
sure what you want to use as your heights. Regardless, assuming you have
an array that houses the heights for each sigma level, I think the
easiest thing for you to do would be to set uv_res@vfYArray = to your
height array. Note that dimensions of your height array would have to
match the bottom_top dimension of your ugrid/vgrid/wspd array.

Hope that helps in some way. If not, let the group know..
Adam

On 09/13/2011 01:42 PM, Tardif, Elliot M wrote:
> Good afternoon,
>
> I am attempting to produce wind profiler-type output from wrf output
> files, colorizing the wind barbs by speed. I have successfully done
> something like this using a combination of the wrf_meteo_5.ncl (link:
> http://www.mmm.ucar.edu/wrf/OnLineTutorial/Graphics/NCL/Examples/SPECIAL/wrf_meteo_5.ncl)
> and the barb_4.ncl scripts (link:
> http://www.ncl.ucar.edu/Applications/Scripts/barb_4.ncl), but the y-axis
> is in sigma levels (see attached meteo4a.ps). Instead of sigma levels, I
> would like to change the y-axis to height, in meters, using the height
> of each sigma level for each hour. What do I need to change in my NCL
> script (attached) to make this happen?
>
> Thanks,
>
> Elliot Tardif.
>
> Elliot Tardif, Meteorologist I
>
> NC DENR, Division of Air Quality
>
> Planning Section, Attainment Planning Branch
>
> 1641 Mail Service Center
>
> Raleigh, NC 27699-1641
>
> Phone: 919-715-9773
>
> Fax : 919-715-7476
>
> Email: Elliot.Tardif@ncdenr.gov <mailto:Nick.Witcraft@ncdenr.gov>
>
> Web : http://www.ncair.org <http://www.ncair.org/>
>
> *****************************************************************************************
>
> *Email correspondence to and from this address is subject to the North
> Carolina Public Records Law and may be disclosed to third parties unless
> the content is exempt by statute or other regulation*.
>
> *****************************************************************************************
>
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

-- 
______________________________________________________________
Adam Phillips                                asphilli@ucar.edu
NCAR/Climate and Global Dynamics Division       (303) 497-1726
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 Wed Sep 14 17:17:09 2011

This archive was generated by hypermail 2.1.8 : Thu Sep 22 2011 - 17:12:44 MDT