Re: Verifying a NARR grid value

How about if I wanted to print the grid number that is located at the
polymarker (x,y)?
With the script below I get the error:
fatal:Number of dimensions on right hand side do not match number of
dimension in left hand side
fatal:Execute: Error occurred at or near line 4335 in file
$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl
fatal:Execute: Error occurred at or near line 135 in file 20080522-2100.ncl

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"

;############################

begin

    fili =
"/home/jlahowet/THOR/test2/narr/u-v_comp/narr-a_221_20080522_2100_000.grb"
    f = addfile (fili, "r")

;############################
; Print
;############################

    names = getfilevarnames(f)
    print(names)
    atts = getfilevaratts(f,names(0))
    dims = getfilevardims(f,names(0))
    print(atts)
    print(dims)

;############################
; Variables
;############################

    upt = f->USTM_221_HTGY
    vpt = f->VSTM_221_HTGY
    lat2d = f->gridlat_221
    lon2d = f->gridlon_221

    upt_at_lat2d = lat2d
    upt_at_lon2d = lon2d
    vpt_at_lat2d = lat2d
    vpt_at_lon2d = lon2d

;############################

    ;wks = gsn_open_wks("x11","lcnative")
    wks = gsn_open_wks("ps","lcnative")
    gsn_define_colormap(wks,"gui_default")

;############################
; Resources
;############################

    res = True
    res_at_tiMainString = "Native Lambert Conformal"

    res_at_mpLimitMode = "Corners"
    res_at_mpLeftCornerLatF = 37.0
    res_at_mpLeftCornerLonF = -103.0
    res_at_mpRightCornerLatF = 40.0
    res_at_mpRightCornerLonF = -98.0

    res_at_mpProjection = "LAMBERTCONFORMAL"
    res_at_mpLambertParallel1F = lon2d_at_Latin1
    res_at_mpLambertParallel2F = lon2d_at_Latin2
    res_at_mpLambertMeridianF = lat2d_at_Lov

; print(lon2d_at_Latin1)
; print(lon2d_at_Latin2)
; print(lon2d_at_Lov)

    res_at_pmTickMarkDisplayMode = "Always"
    res_at_mpFillOn = False
    res_at_mpOutlineDrawOrder = "PostDraw"
    res_at_mpOutlineBoundarySets = "GeophysicalAndUSStates"
    res_at_mpGridAndLimbOn = True

    res_at_tfDoNDCOverlay = False

    res_at_cnFillOn = True
    res_at_cnLinesOn = False
    res_at_gsnSpreadColors = True
    res_at_gsnAddCyclic = False
    res_at_gsnFrame = False
    res_at_gsnDraw = False

;############################

    plot = gsn_csm_contour_map(wks,upt,res)
    ;plot = gsn_csm_contour_map(wks,vpt,res)

;############################
; Contour
;############################

    cnres = True
    cnres_at_tfDoNDCOverlay = True ; don't transform
    cnres_at_cnInfoLabelOn = False ; turn off info label
    cnres_at_cnLineLabelsOn = True ; turn off line labels
    cnres_at_tiMainString = "Contour on Contour: Native Grid" ; title
    cnres_at_gsnDraw = False ; don't draw yet
    cnres_at_gsnFrame = False ; don't advance frame
yet

;############################

    ;plotcn = gsn_contour(wks,vpt,cnres)
    ;overlay(plot,plotcn)

;############################
; Plot marker
;############################

    x = -101.01
    y = 39.60

    polyres = True
    polyres_at_gsMarkerIndex = 16
    polyres_at_gsMarkerSizeF = 7.0
    polyres_at_gsMarkerColor = "black"
    draw(plot)
    gsn_polymarker(wks,plot,x,y,polyres)

;############################
; lat/lon Grid
;############################

    gridres = True
    gridres_at_gsMarkerIndex = 2
    gridres_at_gsMarkerSizeF = 10.0
    gridres_at_gsMarkerColor = "white"
    gsn_polymarker(wks,plot,ndtooned(lon2d),ndtooned(lat2d),gridres)

;############################
; Grid Value at point
;############################

    textres = True
    textres_at_txFontHeightF = 0.01 ; May need to change this
    textres_at_txFontColor = "black"
    vpt_text = ndtooned(vpt) + "" ; Convert to string
    ;vpt_text = sprintf("%7.2f",ndtooned(vpt))
    ;gsn_text(wks,plot,vpt_text,ndtooned(lon2d),ndtooned(lat2d),textres)
    gsn_text(wks,plot,vpt_text,ndtooned(x),ndtooned(y),textres)

    frame(wks)

end

On Thu, Sep 11, 2008 at 10:20 AM, Mary Haley <haley_at_ucar.edu> wrote:

>
> Jamie,
>
> I believe the reason your marker is not showing up is because you are
> using a longitude value that is outside the range of the map area you
> are displaying.
>
> Your latitude/longitude area is defined by:
>
> res_at_mpLeftCornerLonF = -103.0
> res_at_mpRightCornerLonF = -96.0
> res_at_mpLeftCornerLatF = 36.0
> res_at_mpRightCornerLatF = 43.0
>
> and you are trying to plot the point:
>
> x = 98
> y = 56
>
> Longitude 98 is well outside this area, and latitude 56 is as well.
>
> By "grid" do you mean the lat/lon grid? If so, you could plot it using
> the code you've set up below:
>
> polyres = True
> polyres_at_gsMarkerIndex = 16
> polyres_at_gsMarkerSizeF = 10.0
> polyres_at_gsMarkerColor = "white" ; The (/ /) is not necessary
> draw(plot)
> gsn_polymarker(wks,plot,ndtooned(lon2d),ndtooned(lat2d),polyres)
>
> If you want to print the grid value at the lat/lon locations, then
> you would use something like:
>
> textres = True
> textres_at_txFontHeightF = 0.01 ; May need to change this
> textres_at_txFontColor = "white"
> draw(plot)
> vpt_text = ndtooned(vpt) + "" ; Convert to string
> gsn_text(wks,plot,vpt_text,ndtooned(lon2d),ndtooned(lat2d),textres)
>
> You can format the vpt strings using sprintf:
>
> vpt_text = sprintf("%5.2f",ndtooned(vpt))
>
> Of course, you may need to change the "5.2" format to accommodate
> your data range.
>
> --Mary
>
>
> On Wed, 10 Sep 2008, Jamie Lahowetz wrote:
>
> I'm looking for a way to verify that my algorithm is getting the right
>> grid
>> points.
>> Is there a way to display the grid and the number that that grid
>> represents?
>> Right now, The point I try to mark is not showing up and I used the same
>> scheme as before, any ideas?
>> Also, can I show the exact value that is represented by a grid position?
>> Thanks for any help.
>>
>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
>>
>> ;############################
>>
>> begin
>>
>> fili =
>> "/home/jlahowet/THOR/test2/narr/u-v_comp/narr-a_221_20080522_2100_000.grb"
>> f = addfile (fili, "r")
>>
>> ;############################
>> ; Print
>> ;############################
>>
>> names = getfilevarnames(f)
>> print(names)
>>
>> atts = getfilevaratts(f,names(0))
>> dims = getfilevardims(f,names(0))
>> print(atts)
>> print(dims)
>>
>> ;############################
>> ; Variables
>> ;############################
>>
>> upt = f->USTM_221_HTGY
>> vpt = f->VSTM_221_HTGY
>> lat2d = f->gridlat_221
>> lon2d = f->gridlon_221
>>
>> upt_at_lat2d = lat2d
>> upt_at_lon2d = lon2d
>> vpt_at_lat2d = lat2d
>> vpt_at_lon2d = lon2d
>>
>> ;############################
>>
>> wks = gsn_open_wks("x11","lcnative")
>> gsn_define_colormap(wks,"gui_default")
>>
>> ;############################
>> ; Resources
>> ;############################
>>
>> res = True
>> res_at_tiMainString = "Native Lambert Conformal"
>>
>> res_at_mpLimitMode = "Corners"
>> res_at_mpLeftCornerLatF = 36.0
>> res_at_mpLeftCornerLonF = -103.0
>> res_at_mpRightCornerLatF = 43.0
>> res_at_mpRightCornerLonF = -96.0
>>
>> res_at_mpProjection = "LAMBERTCONFORMAL"
>> res_at_mpLambertParallel1F = lon2d_at_Latin1
>> res_at_mpLambertParallel2F = lon2d_at_Latin2
>> res_at_mpLambertMeridianF = lat2d_at_Lov
>>
>> ; print(lon2d_at_Latin1)
>> ; print(lon2d_at_Latin2)
>> ; print(lon2d_at_Lov)
>>
>> res_at_pmTickMarkDisplayMode = "Always"
>> res_at_mpFillOn = False
>> res_at_mpOutlineDrawOrder = "PostDraw"
>> res_at_mpOutlineBoundarySets = "GeophysicalAndUSStates"
>> res_at_mpGridAndLimbOn = True
>>
>> res_at_tfDoNDCOverlay = False
>>
>> res_at_cnFillOn = True
>> res_at_cnLinesOn = False
>> res_at_gsnSpreadColors = True
>> res_at_gsnAddCyclic = False
>> res_at_gsnFrame = False
>> res_at_gsnDraw = False
>>
>> ;############################
>>
>> ;plot = gsn_csm_contour_map(wks,upt,res)
>> plot = gsn_csm_contour_map(wks,vpt,res)
>>
>> ;############################
>> ; Plot marker
>> ;############################
>>
>> ;x = 39.60
>> ;y = -101.01
>> x = 98
>> y = 56
>> polyres = True
>> polyres_at_gsMarkerIndex = 16
>> polyres_at_gsMarkerSizeF = 10.0
>> polyres_at_gsMarkerColor = (/"white"/)
>> draw(plot)
>> gsn_polymarker(wks,plot,x,y,polyres)
>> frame(wks)
>>
>> end
>>
>>
>>
>> --
>> Jamie Ryan Lahowetz
>> University of Nebraska - Lincoln
>> Graduate Student - Geosciences
>> 402.304.0766
>> jlahowe2_at_bigred.unl.edu
>>
>>

-- 
Jamie Ryan Lahowetz
University of Nebraska - Lincoln
Graduate Student - Geosciences
402.304.0766
jlahowe2_at_bigred.unl.edu