Re: Question About gsn_csm_pres_hgt

Hi Elizabeth,
This is unfortunately not a trivial thing to do. We have one possible
suggestion for you:
Try setting pltres@trGridType = "TriangularMesh"

If that does not work please let ncl-talk know.
Adam

On 04/03/2014 09:04 AM, Elizabeth Wilson wrote:
>
> Hello,
>
> I'm using the gsn_csm_pres_hgt function to plot pressure vs. longitude
> for a specific latitude. Right now, I'm just choosing an arbitrary
> latitude to see if the plot works. I do get a plot; however, I want
> the longitude values to appear on the x-axis. Right now, only the
> index values are showing up. I've tried messing with the tickmark
> options but that usually turns them off all together and nothing shows
> up. I'll post my code below and the image that is created. I'm not
> technically getting any errors but I can't get the longitude values to
> appear. Also, as an extension of this question, I've tried messing
> around with coordinate subscripting when using this function instead
> of picking an index; however, I get this error: (gridy_221) of
> (uwind_kts) does not have an associated coordinate variable. I've
> copied over all the coordinates from a previous variable so I wouldn't
> think it should be giving me this error. Any suggestions?
>
> Thanks!
>
> -Elizabeth
>
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
>
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
>
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
>
> load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
>
> begin
>
> ; ######### Open file ###############
>
> f=addfile("/lustre-sas1/ewilson/UPP/postprd/narr-a_221_20121030_0000_000.grb","r")
>
> lat2d=f->gridlat_221(:,:)
>
> lat2d!1="lat"
>
> lat2d!0="lon"
>
> lat=f->gridlat_221(0,:)
>
> lon2d=f->gridlon_221(:,:)
>
> lon2d!1="lat"
>
> lon2d!0="lon"
>
> lon=f->gridlon_221(:,0)
>
> dimy=dimsizes(lat)
>
> dimx=dimsizes(lon)
>
> uwind=f->U_GRD_221_ISBL(:,:,:)
>
> ;print(ind_nearest_coord(35.8189,lat,0))
>
> ;ind= ind(lat.eq.35.8189)
>
> ;print(ind)
>
> ;uwind!0="level"
>
> ;uwind!2="lat2d"
>
> ;uwind!1="lon2d"
>
> ;uwind@coordinates="lon2d lat2d"
>
> ;uwind@lat2d=lat2d
>
> ;uwind@lon2d=lon2d
>
> uwind_kts=uwind*1.944
>
> copy_VarCoords(uwind,uwind_kts)
>
> printVarSummary(uwind_kts)
>
> ;copy_VarMeta(uwind,uwind_kts)
>
> ;uwind_kts!0="level"
>
> ;uwind_kts!2="lat2d"
>
> ;uwind_kts!1="lon2d"
>
> ;uwind_kts@coordinates="lon2d lat2d"
>
> ;uwind_kts@lat2d=lat2d
>
> ;uwind_kts@lon2d=lon2d
>
> wks=gsn_open_wks("pdf","test")
>
> gsn_define_colormap(wks,"WhBlGrYeRe")
>
> pltres=True
>
> pltres@gsnDraw=True
>
> pltres@gsnFrame=True
>
> pltres@gsnAddCyclic=False
>
> plot=gsn_csm_pres_hgt(wks,uwind_kts(:,:,50),pltres)
>
> delete(pltres)
>
> delete(uwind_kts)
>
> end
>
>
>
> _______________________________________________
> 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/staff/asphilli

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Apr 07 11:27:34 2014