<div dir="ltr">Hi, <div><br></div><div>I am attempting to plot a subset of the MRMS data from NCEP.  The latitude and longitude I am plotting is : </div><div><br></div><div>LAT: <span style="color:rgb(0,0,0);font-family:&quot;Times New Roman&quot;;font-weight:bold;text-align:-webkit-center">36.5459089</span></div><div>LON: <span style="color:rgb(0,0,0);font-family:&quot;Times New Roman&quot;;font-weight:bold;text-align:-webkit-center">-79.4585847</span></div><div><br></div><div>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! </div><div><br></div><div>Thanks! </div><div>Tim </div><div><br></div><div>Code: </div><div><br></div><div><div>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl&quot;</div><div>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl&quot;</div><div>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl&quot;</div><div><br></div><div>begin</div><div><br></div><div>lati=LAT</div><div>loni=LON</div><div>date=DATE</div><div>address=ADDRESS</div><div><br></div><div><br></div><div>; Convert the Lat/Lon to a Float</div><div>latc= stringtofloat(lati)</div><div>lonc = stringtofloat(loni)</div><div><br></div><div>; Set the Map Bounds</div><div>minlat = latc - .12</div><div>maxlat = latc + .12</div><div>minlon = lonc - .15</div><div>maxlon = lonc + .15</div><div><br></div><div><br></div><div>  ;========================</div><div>  ; get list of all files and open as &quot;one big file&quot;</div><div>  ;========================</div><div>  a = addfile( &quot;./&quot; + lati + &quot;_&quot; + loni + &quot;.grib2&quot;,&quot;r&quot;)</div><div><br></div><div>  type = &quot;png&quot;</div><div>  wks = gsn_open_wks(type,&quot;./&quot; + lati + &quot;&quot; + loni + &quot;.png&quot;)</div><div><br></div><div>; Set some Basic Plot options</div><div>  res = True</div><div>  res@gsnDraw      =  False                   ; do not draw the plot</div><div>  res@gsnFrame     =  False                   ; do not advance the frame</div><div>  pltres = True</div><div>  pltres@NoTitles = True</div><div>  pltres@CommonTitle = True</div><div>  pltres@PlotTitle = &quot;&quot;</div><div>  pltres@PanelPlot = True</div><div>  pltres@FramePlot = False</div><div><br></div><div>; Set the Colormap</div><div>  gsn_define_colormap(wks,&quot;hail&quot;)</div><div>    lat  = a-&gt;lat_0</div><div>    lon  = a-&gt;lon_0</div><div>    water = a-&gt;VAR_209_3_34_P0_L102_GLL0(:,:)</div><div>    hail = water  *0.0393701;</div><div><br></div><div><br></div><div>  ; Plotting options for Hail</div><div>  res@gsnMaximize = True</div><div>  res@mpGeophysicalLineThicknessF = 4.0</div><div>  ;res@mpProjection      = &quot;LambertConformal&quot;</div><div>  res@mpGridLineThicknessF        = 0.5</div><div>  res@mpLimbLineThicknessF        = 4.0</div><div>  res@mpNationalLineThicknessF    = 4.0</div><div>  res@mpUSStateLineThicknessF     = 4.0</div><div>  res@mpProvincialLineThicknessF  = 4.0</div><div>  res@mpOutlineBoundarySets = &quot;AllBoundaries&quot;</div><div>  res@mpDataBaseVersion     = &quot;Ncarg4_1&quot;</div><div>  res@mpDataSetName         = &quot;Earth..4&quot;</div><div>  res@mpFillOn            = False         ; fill continents</div><div>  res@mpPerimOn           = False        ; no box around map</div><div>  res@tiMainOn = True</div><div>  res@tiMainFontColor = &quot;Black&quot;</div><div>  res@tiMainString =&quot;&quot;</div><div>  res@tiXAxisString   = &quot;&quot;</div><div>  res@tiYAxisString   = &quot;&quot;</div><div>  res@gsnLeftString   = &quot;Date: &quot; + date               ; add the gsn titles</div><div>  res@gsnRightString  = &quot;&quot;</div><div>  res@gsnStringFontHeightF = 0.016</div><div>  res@cnInfoLabelOrthogonalPosF = 0.07  ; offset second label information</div><div>  res@gsnContourLineThicknessesScale = 0.001</div><div>  res@cnFillOn = True</div><div>  res@lbTitleOn = True                                 ; remove field name from label bar</div><div><br></div><div> ; Set the Map Bounds</div><div>        res@sfYArray = lat</div><div>        res@sfXArray = lon</div><div>        res@mpLimitMode                 = &quot;LatLon&quot;</div><div>        res@mpMinLatF                   = minlat</div><div>        res@mpMaxLatF                   = maxlat</div><div>        res@mpMinLonF                   = minlon</div><div>        res@mpMaxLonF                   = maxlon</div><div>        res@tmXBOn = False</div><div>        res@tmYLOn = False</div><div>        res@tmXTOn = False</div><div>        res@tmYROn = False</div><div>        res@gsnAddCyclic                = False</div><div>        res@tfDoNDCOverlay              = False            ; do not transform data</div><div>        res@cnLevelSelectionMode = &quot;ExplicitLevels&quot;   ; set explicit contour levels</div><div>        res@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</div><div>        res@tiMainFontHeightF     = 0.018               ; font height of title</div><div>        res@lbAutoManage          = False             ; we control label bar</div><div>        res@pmLabelBarDisplayMode = &quot;Always&quot;          ; turns on label bar</div><div>        res@lbOrientation         = &quot;Horizontal&quot;      ; ncl default is vertical</div><div>        res@pmLabelBarSide        = &quot;Bottom&quot;          ; default is right</div><div>        res@lbLabelStride         = 1                 ; skip every other label</div><div>        res@pmLabelBarWidthF      = 0.7               ; default is shorter</div><div>        res@pmLabelBarHeightF     = 0.1               ; default is taller</div><div>        res@lbLabelFontHeightF    = .012              ; default is HUGE</div><div>        res@lbPerimOn             = False             ; default has box</div><div>        res@lbTitleString = &quot;Hail Size (Inches)&quot;</div><div>        res@lbTitleFontHeightF = 0.012</div><div>        plot = gsn_csm_contour_map(wks,hail,res)</div><div><br></div><div>; Setup the point</div><div> mstring1 = &quot;z&quot;</div><div> fontnum1 = 35</div><div> xoffset1 = .4</div><div> yoffset1 = .7</div><div> ratio1   = 0.0</div><div> size1    = 1.0</div><div> angle1   = 0.0</div><div><br></div><div> Tm_index = NhlNewMarker(wks, mstring1, fontnum1, xoffset1, yoffset1, ratio1, size1, angle1)</div><div><br></div><div>            ; Plot the sites and label them</div><div>                gres = True</div><div>                gres@gsMarkerIndex = Tm_index ; 16</div><div>                gres@gsMarkerSizeF = 0.015 ; .009</div><div>                gres@gsMarkerColor = &quot;black&quot;</div><div>                gres@txFont =&quot;helvetica-bold&quot;</div><div>                gres@txFontHeightF = &quot;25&quot;</div><div><br></div><div>site = (/address/)</div><div>lats = (/latc/)</div><div>lons = (/lonc/)</div><div><br></div><div>dot = gsn_add_polymarker(wks,plot,lons,lats,gres)</div><div>tlats =lats + .03</div><div>tlons =lons</div><div>text = gsn_add_text(wks,plot,site,tlons,tlats ,gres)</div><div><br></div><div>; Draw the graphics</div><div>  draw(plot)</div><div>  frame(wks)</div><div><br></div><div> end</div></div><div><br></div></div>