[ncl-talk] Plotting Lat/Lon Points on a map
Tim Melino
melino33 at gmail.com
Sun May 28 10:13:11 MDT 2017
Thanks Barry and Dennis for your help on this!
Dennis, I added your my other question was why is the point in NCL plotted
south of the North Carolina border with NCL? The actual location is in
Virginia, I noticed this was occurring in the image that you sent me too.
[image: Inline image 1]
On Sun, May 28, 2017 at 12:00 PM, Barry Lynn <barry.h.lynn at gmail.com> wrote:
> Hi Dennis:
>
> Thank you for helping. Sorry I was unavailable to (try to) do more.
>
> Barry
>
> On Sun, May 28, 2017 at 6:25 PM, Dennis Shea <shea at ucar.edu> wrote:
>
>> See attached. Please look carefully at the output produced by
>> printVarSummary and printMinMax
>>
>> Note:
>>
>> [1] GRIB is a moving target. 6.3.0 did not have the look-up tabl builtin
>> while NCL 6.4.0 does.
>>
>> [2] the grib file has 'no report' coded as as -1 ... the script sets
>> those values to _FillValue
>>
>> [3] The grid is *rectilinear* :
>> In netCDF speak, a one dimensional variable where the variable name is
>> the same as the variable dimension name is a coordinate variable.
>> lat_0(lat_0) and lon_0(lon_0)
>>
>> Good Luck
>>
>>
>>
>> On Sun, May 28, 2017 at 10:21 AM, Barry Lynn <barry.h.lynn at gmail.com>
>> wrote:
>>
>>> Tim:
>>>
>>> If you can wait a day or two, I suggest you write to:
>>>
>>> Carrie.Langston at noaa.gov
>>>
>>> and ask Him/Her what is the format of the gridded data and if they have
>>> an NCL program to read it.
>>>
>>> It's looks like very useful data once you "crack the code."
>>>
>>> Barry
>>>
>>> On Sun, May 28, 2017 at 5:09 PM, Tim Melino <melino33 at gmail.com> wrote:
>>>
>>>> Barry,
>>>>
>>>> The source of the data is here:
>>>>
>>>> http://www.nssl.noaa.gov/projects/mrms/operational/tables.php
>>>>
>>>>
>>>>
>>>> Tim
>>>>
>>>>
>>>>
>>>> On Sun, May 28, 2017 at 10:03 AM, Barry Lynn <barry.h.lynn at gmail.com>
>>>> wrote:
>>>>
>>>>> Tim:
>>>>>
>>>>> If you write out the latitude and longitude information is is
>>>>> rectilinear?
>>>>>
>>>>> When trying to sort this out, you should first check if the lat/lon
>>>>> can be plotted by a single coordinate or two. If two, then there is a
>>>>> chance that your grid is not rectilinear, and you will have to account for
>>>>> the distortion in the grid.
>>>>>
>>>>> If there is distortion, then you need to account for it (e.g.,
>>>>> Lambert, curvilinear.
>>>>>
>>>>> I can only offer advice at the moment...
>>>>>
>>>>> What is the data source for the data?
>>>>>
>>>>> Barry
>>>>>
>>>>> On Sun, May 28, 2017 at 4:52 PM, Tim Melino <melino33 at gmail.com>
>>>>> wrote:
>>>>>
>>>>>> Hey Barry,
>>>>>>
>>>>>> The only information about the projection was: *grid_type :
>>>>>> Latitude/Longitude. *I attached the file!
>>>>>>
>>>>>> Tim
>>>>>>
>>>>>> On Sun, May 28, 2017 at 9:00 AM, Barry Lynn <barry.h.lynn at gmail.com>
>>>>>> wrote:
>>>>>>
>>>>>>> Tim:
>>>>>>>
>>>>>>> Please remind me what the name of the file is.
>>>>>>>
>>>>>>> Also, was there any information at the top of the printing about the
>>>>>>> projection?
>>>>>>>
>>>>>>> Yes, the markers would (could) look off.
>>>>>>>
>>>>>>> On Sun, May 28, 2017 at 3:32 PM, Tim Melino <melino33 at gmail.com>
>>>>>>> wrote:
>>>>>>>
>>>>>>>> Hi Barry,
>>>>>>>>
>>>>>>>> Here is the information when printing (a). Would the map projection
>>>>>>>> also cause a marker to be printed in the wrong location?
>>>>>>>>
>>>>>>>> Variable: a
>>>>>>>> Type: file
>>>>>>>> filename: test
>>>>>>>> path: ./test.grib2
>>>>>>>> file global attributes:
>>>>>>>> dimensions:
>>>>>>>> lat_0 = 30
>>>>>>>> lon_0 = 40
>>>>>>>> variables:
>>>>>>>> float VAR_209_3_34_P0_L102_GLL0 ( lat_0, lon_0 )
>>>>>>>> center : US NOAA Office of Oceanic and Atmospheric
>>>>>>>> Research
>>>>>>>> production_status : Research products
>>>>>>>> long_name : unknown variable name
>>>>>>>> units : unknown
>>>>>>>> _FillValue : 1e+20
>>>>>>>> grid_type : Latitude/longitude
>>>>>>>> parameter_discipline_and_category : 209, 3
>>>>>>>> parameter_template_discipline_category_number : (
>>>>>>>> 0, 209, 3, 34 )
>>>>>>>> level_type : Specific altitude above mean sea level (m)
>>>>>>>> level : 500
>>>>>>>> forecast_time : 0
>>>>>>>> forecast_time_units : minutes
>>>>>>>> initial_time : 07/09/2016 (04:46)
>>>>>>>>
>>>>>>>> float lat_0 ( lat_0 )
>>>>>>>> long_name : latitude
>>>>>>>> grid_type : Latitude/Longitude
>>>>>>>> units : degrees_north
>>>>>>>> Dj : 0.01
>>>>>>>> Di : 0.009999676
>>>>>>>> Lo2 : 280.735
>>>>>>>> La2 : 36.695
>>>>>>>> Lo1 : 280.345
>>>>>>>> La1 : 36.405
>>>>>>>>
>>>>>>>> float lon_0 ( lon_0 )
>>>>>>>> long_name : longitude
>>>>>>>> grid_type : Latitude/Longitude
>>>>>>>> units : degrees_east
>>>>>>>> Dj : 0.01
>>>>>>>> Di : 0.009999676
>>>>>>>> Lo2 : 280.735
>>>>>>>> La2 : 36.695
>>>>>>>> Lo1 : 280.345
>>>>>>>> La1 : 36.405
>>>>>>>>
>>>>>>>>
>>>>>>>> On Sun, May 28, 2017 at 12:48 AM, Barry Lynn <
>>>>>>>> barry.h.lynn at gmail.com> wrote:
>>>>>>>>
>>>>>>>>> Hi:
>>>>>>>>>
>>>>>>>>> When the data is skewed, it is often because the projection you're
>>>>>>>>> reading is not the one your plotting.
>>>>>>>>>
>>>>>>>>> Do you know the projection of the data you're using?
>>>>>>>>>
>>>>>>>>> If you print(a) where a is the file you specify from addfile, you
>>>>>>>>> can see the variables (and hopefully projection) indicated within.
>>>>>>>>>
>>>>>>>>> Barry
>>>>>>>>>
>>>>>>>>> On Sat, May 27, 2017 at 4:20 PM, Tim Melino <melino33 at gmail.com>
>>>>>>>>> wrote:
>>>>>>>>>
>>>>>>>>>> Hi,
>>>>>>>>>>
>>>>>>>>>> I am attempting to plot a subset of the MRMS data from NCEP. The
>>>>>>>>>> latitude and longitude I am plotting is :
>>>>>>>>>>
>>>>>>>>>> LAT: 36.5459089
>>>>>>>>>> LON: -79.4585847
>>>>>>>>>>
>>>>>>>>>> I attached an image of the output along with the grib file I am
>>>>>>>>>> using for input. I am noticing a couple things first the data appears to
>>>>>>>>>> be skewed, so I am not sure if it is an issue with the map projection.
>>>>>>>>>> Second the point which should be plotted in Virginia is actually plotted in
>>>>>>>>>> North Carolina, the only thing that I have been able to do if offset the
>>>>>>>>>> marker to compensate for the issue. Does anyone have any ideas on how to
>>>>>>>>>> fix these issues? Any ideas would be greatly appreciated!
>>>>>>>>>>
>>>>>>>>>> Thanks!
>>>>>>>>>> Tim
>>>>>>>>>>
>>>>>>>>>> Code:
>>>>>>>>>>
>>>>>>>>>> load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
>>>>>>>>>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
>>>>>>>>>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
>>>>>>>>>>
>>>>>>>>>> begin
>>>>>>>>>>
>>>>>>>>>> lati=LAT
>>>>>>>>>> loni=LON
>>>>>>>>>> date=DATE
>>>>>>>>>> address=ADDRESS
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>> ; Convert the Lat/Lon to a Float
>>>>>>>>>> latc= stringtofloat(lati)
>>>>>>>>>> lonc = stringtofloat(loni)
>>>>>>>>>>
>>>>>>>>>> ; Set the Map Bounds
>>>>>>>>>> minlat = latc - .12
>>>>>>>>>> maxlat = latc + .12
>>>>>>>>>> minlon = lonc - .15
>>>>>>>>>> maxlon = lonc + .15
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>> ;========================
>>>>>>>>>> ; get list of all files and open as "one big file"
>>>>>>>>>> ;========================
>>>>>>>>>> a = addfile( "./" + lati + "_" + loni + ".grib2","r")
>>>>>>>>>>
>>>>>>>>>> type = "png"
>>>>>>>>>> wks = gsn_open_wks(type,"./" + lati + "" + loni + ".png")
>>>>>>>>>>
>>>>>>>>>> ; Set some Basic Plot options
>>>>>>>>>> res = True
>>>>>>>>>> res at gsnDraw = False ; do not draw the
>>>>>>>>>> plot
>>>>>>>>>> res at gsnFrame = False ; do not advance
>>>>>>>>>> the frame
>>>>>>>>>> pltres = True
>>>>>>>>>> pltres at NoTitles = True
>>>>>>>>>> pltres at CommonTitle = True
>>>>>>>>>> pltres at PlotTitle = ""
>>>>>>>>>> pltres at PanelPlot = True
>>>>>>>>>> pltres at FramePlot = False
>>>>>>>>>>
>>>>>>>>>> ; Set the Colormap
>>>>>>>>>> gsn_define_colormap(wks,"hail")
>>>>>>>>>> lat = a->lat_0
>>>>>>>>>> lon = a->lon_0
>>>>>>>>>> water = a->VAR_209_3_34_P0_L102_GLL0(:,:)
>>>>>>>>>> hail = water *0.0393701;
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>> ; Plotting options for Hail
>>>>>>>>>> res at gsnMaximize = True
>>>>>>>>>> res at mpGeophysicalLineThicknessF = 4.0
>>>>>>>>>> ;res at mpProjection = "LambertConformal"
>>>>>>>>>> res at mpGridLineThicknessF = 0.5
>>>>>>>>>> res at mpLimbLineThicknessF = 4.0
>>>>>>>>>> res at mpNationalLineThicknessF = 4.0
>>>>>>>>>> res at mpUSStateLineThicknessF = 4.0
>>>>>>>>>> res at mpProvincialLineThicknessF = 4.0
>>>>>>>>>> res at mpOutlineBoundarySets = "AllBoundaries"
>>>>>>>>>> res at mpDataBaseVersion = "Ncarg4_1"
>>>>>>>>>> res at mpDataSetName = "Earth..4"
>>>>>>>>>> res at mpFillOn = False ; fill continents
>>>>>>>>>> res at mpPerimOn = False ; no box around map
>>>>>>>>>> res at tiMainOn = True
>>>>>>>>>> res at tiMainFontColor = "Black"
>>>>>>>>>> res at tiMainString =""
>>>>>>>>>> res at tiXAxisString = ""
>>>>>>>>>> res at tiYAxisString = ""
>>>>>>>>>> res at gsnLeftString = "Date: " + date ; add the
>>>>>>>>>> gsn titles
>>>>>>>>>> res at gsnRightString = ""
>>>>>>>>>> res at gsnStringFontHeightF = 0.016
>>>>>>>>>> res at cnInfoLabelOrthogonalPosF = 0.07 ; offset second label
>>>>>>>>>> information
>>>>>>>>>> res at gsnContourLineThicknessesScale = 0.001
>>>>>>>>>> res at cnFillOn = True
>>>>>>>>>> res at lbTitleOn = True ; remove
>>>>>>>>>> field name from label bar
>>>>>>>>>>
>>>>>>>>>> ; Set the Map Bounds
>>>>>>>>>> res at sfYArray = lat
>>>>>>>>>> res at sfXArray = lon
>>>>>>>>>> res at mpLimitMode = "LatLon"
>>>>>>>>>> res at mpMinLatF = minlat
>>>>>>>>>> res at mpMaxLatF = maxlat
>>>>>>>>>> res at mpMinLonF = minlon
>>>>>>>>>> res at mpMaxLonF = maxlon
>>>>>>>>>> res at tmXBOn = False
>>>>>>>>>> res at tmYLOn = False
>>>>>>>>>> res at tmXTOn = False
>>>>>>>>>> res at tmYROn = False
>>>>>>>>>> res at gsnAddCyclic = False
>>>>>>>>>> res at tfDoNDCOverlay = False ; do
>>>>>>>>>> not transform data
>>>>>>>>>> res at cnLevelSelectionMode = "ExplicitLevels" ; set
>>>>>>>>>> explicit contour levels
>>>>>>>>>> res at cnLevels = (/ 0.5,0.75, 1.00, 1.25, 1.50, 1.75,
>>>>>>>>>> 2.00, 2.50, 3.00, 3.50, 4.00, 4.50, 5.00/) ; set levels
>>>>>>>>>> res at tiMainFontHeightF = 0.018 ; font
>>>>>>>>>> height of title
>>>>>>>>>> res at lbAutoManage = False ; we
>>>>>>>>>> control label bar
>>>>>>>>>> res at pmLabelBarDisplayMode = "Always" ; turns on
>>>>>>>>>> label bar
>>>>>>>>>> res at lbOrientation = "Horizontal" ; ncl
>>>>>>>>>> default is vertical
>>>>>>>>>> res at pmLabelBarSide = "Bottom" ; default
>>>>>>>>>> is right
>>>>>>>>>> res at lbLabelStride = 1 ; skip
>>>>>>>>>> every other label
>>>>>>>>>> res at pmLabelBarWidthF = 0.7 ; default
>>>>>>>>>> is shorter
>>>>>>>>>> res at pmLabelBarHeightF = 0.1 ; default
>>>>>>>>>> is taller
>>>>>>>>>> res at lbLabelFontHeightF = .012 ; default
>>>>>>>>>> is HUGE
>>>>>>>>>> res at lbPerimOn = False ; default
>>>>>>>>>> has box
>>>>>>>>>> res at lbTitleString = "Hail Size (Inches)"
>>>>>>>>>> res at lbTitleFontHeightF = 0.012
>>>>>>>>>> plot = gsn_csm_contour_map(wks,hail,res)
>>>>>>>>>>
>>>>>>>>>> ; Setup the point
>>>>>>>>>> mstring1 = "z"
>>>>>>>>>> fontnum1 = 35
>>>>>>>>>> xoffset1 = .4
>>>>>>>>>> yoffset1 = .7
>>>>>>>>>> ratio1 = 0.0
>>>>>>>>>> size1 = 1.0
>>>>>>>>>> angle1 = 0.0
>>>>>>>>>>
>>>>>>>>>> Tm_index = NhlNewMarker(wks, mstring1, fontnum1, xoffset1,
>>>>>>>>>> yoffset1, ratio1, size1, angle1)
>>>>>>>>>>
>>>>>>>>>> ; Plot the sites and label them
>>>>>>>>>> gres = True
>>>>>>>>>> gres at gsMarkerIndex = Tm_index ; 16
>>>>>>>>>> gres at gsMarkerSizeF = 0.015 ; .009
>>>>>>>>>> gres at gsMarkerColor = "black"
>>>>>>>>>> gres at txFont ="helvetica-bold"
>>>>>>>>>> gres at txFontHeightF = "25"
>>>>>>>>>>
>>>>>>>>>> site = (/address/)
>>>>>>>>>> lats = (/latc/)
>>>>>>>>>> lons = (/lonc/)
>>>>>>>>>>
>>>>>>>>>> dot = gsn_add_polymarker(wks,plot,lons,lats,gres)
>>>>>>>>>> tlats =lats + .03
>>>>>>>>>> tlons =lons
>>>>>>>>>> text = gsn_add_text(wks,plot,site,tlons,tlats ,gres)
>>>>>>>>>>
>>>>>>>>>> ; Draw the graphics
>>>>>>>>>> draw(plot)
>>>>>>>>>> frame(wks)
>>>>>>>>>>
>>>>>>>>>> end
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>> _______________________________________________
>>>>>>>>>> ncl-talk mailing list
>>>>>>>>>> ncl-talk at ucar.edu
>>>>>>>>>> List instructions, subscriber options, unsubscribe:
>>>>>>>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> --
>>>>>>>>> Barry H. Lynn, Ph.D
>>>>>>>>> Senior Lecturer,
>>>>>>>>> The Institute of the Earth Science,
>>>>>>>>> The Hebrew University of Jerusalem,
>>>>>>>>> Givat Ram, Jerusalem 91904, Israel
>>>>>>>>> Tel: 972 547 231 170
>>>>>>>>> Fax: (972)-25662581
>>>>>>>>>
>>>>>>>>> C.E.O, Weather It Is, LTD
>>>>>>>>> Weather and Climate Focus
>>>>>>>>> http://weather-it-is.com
>>>>>>>>> Jerusalem, Israel
>>>>>>>>> Local: 02 930 9525
>>>>>>>>> Cell: 054 7 231 170
>>>>>>>>> Int-IS: x972 2 930 9525
>>>>>>>>> US 914 432 3108 <(914)%20432-3108>
>>>>>>>>>
>>>>>>>>
>>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> --
>>>>>>> Barry H. Lynn, Ph.D
>>>>>>> Senior Lecturer,
>>>>>>> The Institute of the Earth Science,
>>>>>>> The Hebrew University of Jerusalem,
>>>>>>> Givat Ram, Jerusalem 91904, Israel
>>>>>>> Tel: 972 547 231 170
>>>>>>> Fax: (972)-25662581
>>>>>>>
>>>>>>> C.E.O, Weather It Is, LTD
>>>>>>> Weather and Climate Focus
>>>>>>> http://weather-it-is.com
>>>>>>> Jerusalem, Israel
>>>>>>> Local: 02 930 9525
>>>>>>> Cell: 054 7 231 170
>>>>>>> Int-IS: x972 2 930 9525
>>>>>>> US 914 432 3108 <(914)%20432-3108>
>>>>>>>
>>>>>>
>>>>>>
>>>>>
>>>>>
>>>>> --
>>>>> Barry H. Lynn, Ph.D
>>>>> Senior Lecturer,
>>>>> The Institute of the Earth Science,
>>>>> The Hebrew University of Jerusalem,
>>>>> Givat Ram, Jerusalem 91904, Israel
>>>>> Tel: 972 547 231 170
>>>>> Fax: (972)-25662581
>>>>>
>>>>> C.E.O, Weather It Is, LTD
>>>>> Weather and Climate Focus
>>>>> http://weather-it-is.com
>>>>> Jerusalem, Israel
>>>>> Local: 02 930 9525
>>>>> Cell: 054 7 231 170
>>>>> Int-IS: x972 2 930 9525
>>>>> US 914 432 3108 <(914)%20432-3108>
>>>>>
>>>>
>>>>
>>>
>>>
>>> --
>>> Barry H. Lynn, Ph.D
>>> Senior Lecturer,
>>> The Institute of the Earth Science,
>>> The Hebrew University of Jerusalem,
>>> Givat Ram, Jerusalem 91904, Israel
>>> Tel: 972 547 231 170
>>> Fax: (972)-25662581
>>>
>>> C.E.O, Weather It Is, LTD
>>> Weather and Climate Focus
>>> http://weather-it-is.com
>>> Jerusalem, Israel
>>> Local: 02 930 9525
>>> Cell: 054 7 231 170
>>> Int-IS: x972 2 930 9525
>>> US 914 432 3108 <(914)%20432-3108>
>>>
>>> _______________________________________________
>>> ncl-talk mailing list
>>> ncl-talk at ucar.edu
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>
>>>
>>
>
>
> --
> Barry H. Lynn, Ph.D
> Senior Lecturer,
> The Institute of the Earth Science,
> The Hebrew University of Jerusalem,
> Givat Ram, Jerusalem 91904, Israel
> Tel: 972 547 231 170
> Fax: (972)-25662581
>
> C.E.O, Weather It Is, LTD
> Weather and Climate Focus
> http://weather-it-is.com
> Jerusalem, Israel
> Local: 02 930 9525
> Cell: 054 7 231 170
> Int-IS: x972 2 930 9525
> US 914 432 3108 <(914)%20432-3108>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170528/8d2ce2d4/attachment.html
-------------- next part --------------
A non-text attachment was scrubbed...
Name: sample_point.png
Type: image/png
Size: 95010 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170528/8d2ce2d4/attachment.png
More information about the ncl-talk
mailing list