<div dir="ltr"><div>Hi Alexander,<br><br><br></div><div>It worked. I catch what you are saying too. <br><br></div><div>Thanks,<br></div><div>Kerwyn<br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Aug 2, 2016 at 3:20 PM, Alexander Schaefer <span dir="ltr"><<a href="mailto:aschae11@uncc.edu" target="_blank">aschae11@uncc.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word">Kerwyn,<div><br></div><div>NARR data is 32 km grid spacing. You are plotting 1 degrees by 0.65 degrees. 1 degree is about 111 km so you would only expect to have 4x3 or less wind barbs plotted. </div><div><br></div><div>If you change the contour interval on the temps(0.25 plotted here) then you don’t have a single value covering the plot. This is the same issue with the heights. I set it to 5 for the spacing a height line fits within your plot space.</div><div><br></div><div>You can see the wind barbs you had, with 1 geopotential line and some colored contours for your temps.</div><div><br></div><div>Plotting such a small area of interest with coarse data what is making the plot difficult.</div><div><br></div><div>I hope that offers some insight,</div><div>-Alex</div><div><br></div><div><img src="cid:9F4202AE-4467-4705-9E48-8737CA5DDD9E" height="539" width="539"></div><div><br></div><div><br><div><blockquote type="cite"><div><div class="h5"><div>On Aug 2, 2016, at 5:43 PM, Kerwyn Texeira <<a href="mailto:ktish86@gmail.com" target="_blank">ktish86@gmail.com</a>> wrote:</div><br></div></div><div><div><div class="h5"><div dir="ltr"><div><div><div>Hi ncl-talk,<br><br></div>I have made changes to my NARR script but it is still not plotting a zoomed in overlay plot. I tried several times playing around with the script by commenting and adding stuff but I'm not seeing much of any progress. <br><br></div>When I commented out all the tfDoNDCOverlay and left lat2d and lon2d uncommented, I got the entire USA plotted with both temp, geopotential height and winds overlayed.<br><br></div>When I commented out lat2d and lon2d and left tfDoNDCOverlay uncommented, I got mainly a plain figure with one or two wind barbs. (temp_narr) figure attached.<br><div><div><br></div><div>Kindly help, please. I do appreciate the help that I'm getting but I'm just not quite reach my ultimate goal which is a zoomed in overlay plot of temps, winds, and geopotential height. <br><br></div><div>Thanks<br></div><div><br></div><div>Some outputs:<br><br>Variable: temp<br>Type: float<br>Total Size: 386692 bytes<br> 96673 values<br>Number of Dimensions: 2<br>Dimensions and sizes: [gridx_221 | 277] x [gridy_221 | 349]<br>Coordinates:<br>Number Of Attributes: 15<br> lv_ISBL3 : 600<br> sub_center : The North American Regional Reanalysis (NARR) Project<br> center : US National Weather Service - NCEP (WMC)<br> long_name : Temperature<br> units : K<br> _FillValue : 1e+20<br> coordinates : gridlat_221 gridlon_221<br> level_indicator : 100<br> grid_number : 221<br> parameter_table_version : 131<br> parameter_number : 11<br> model : North American Regional Reanalysis (NARR)<br> forecast_time : 0<br> forecast_time_units : hours<br> initial_time : 01/28/2014 (00:00)<br><br>Variable: u_wind<br>Type: float<br>Total Size: 386692 bytes<br> 96673 values<br>Number of Dimensions: 2<br>Dimensions and sizes: [gridx_221 | 277] x [gridy_221 | 349]<br>Coordinates:<br>Number Of Attributes: 15<br> lv_ISBL3 : 600<br> sub_center : The North American Regional Reanalysis (NARR) Project<br> center : US National Weather Service - NCEP (WMC)<br> long_name : u-component of wind<br> units : m/s<br> _FillValue : 1e+20<br> coordinates : gridlat_221 gridlon_221<br> level_indicator : 100<br> grid_number : 221<br> parameter_table_version : 131<br> parameter_number : 33<br> model : North American Regional Reanalysis (NARR)<br> forecast_time : 0<br> forecast_time_units : hours<br> initial_time : 01/28/2014 (00:00)<br><br>Variable: v_wind<br>Type: float<br>Total Size: 386692 bytes<br> 96673 values<br>Number of Dimensions: 2<br>Dimensions and sizes: [gridx_221 | 277] x [gridy_221 | 349]<br>Coordinates:<br>Number Of Attributes: 15<br> lv_ISBL3 : 600<br> sub_center : The North American Regional Reanalysis (NARR) Project<br> center : US National Weather Service - NCEP (WMC)<br> long_name : v-component of wind<br> units : m/s<br> _FillValue : 1e+20<br> coordinates : gridlat_221 gridlon_221<br> level_indicator : 100<br> grid_number : 221<br> parameter_table_version : 131<br> parameter_number : 34<br> model : North American Regional Reanalysis (NARR)<br> forecast_time : 0<br> forecast_time_units : hours<br> initial_time : 01/28/2014 (00:00)<br><br>Variable: hgt<br>Type: float<br>Total Size: 386692 bytes<br> 96673 values<br>Number of Dimensions: 2<br>Dimensions and sizes: [gridx_221 | 277] x [gridy_221 | 349]<br>Coordinates:<br>Number Of Attributes: 15<br> lv_ISBL3 : 600<br> sub_center : The North American Regional Reanalysis (NARR) Project<br> center : US National Weather Service - NCEP (WMC)<br> long_name : Geopotential height<br> units : gpm<br> _FillValue : 1e+20<br> coordinates : gridlat_221 gridlon_221<br> level_indicator : 100<br> grid_number : 221<br> parameter_table_version : 131<br> parameter_number : 7<br> model : North American Regional Reanalysis (NARR)<br> forecast_time : 0<br> forecast_time_units : hours<br> initial_time : 01/28/2014 (00:00)<br><br></div><div>My Script:<br><br>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"<br>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"<br>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"<br><br>begin<br>a = addfile("./narr-a_221_20140128_0000_000.grb", "r")<br><br>wks_type = "png"<br>wks_type@wkWidth = 2500<br>wks_type@wkHeight = 2500<br><br>wks= gsn_open_wks(wks_type, "temp_narr")<br>gsn_define_colormap(wks,"matlab_jet")<br><br>lat2d = a->gridlat_221<br>lon2d = a->gridlon_221<br>hgt = a->HGT_221_ISBL({600},:,:)<br>temp = a->TMP_221_ISBL({600},:,:)<br>pres = a->PRES_221_HYBL<br>u_wind = a->U_GRD_221_ISBL({600},:,:)<br>v_wind = a->V_GRD_221_ISBL({600},:,:)<br>lev = a->lv_ISBL3 ; 29 levels<br><br>printVarSummary(temp)<br>printVarSummary(u_wind)<br>printVarSummary(v_wind)<br>printVarSummary(hgt)<br><br>spd = (u_wind*u_wind + v_wind*v_wind)^(0.5) ; m/s<br><br>pres = pres*0.01 ; converts to hpa<br><br>u_wind = u_wind*1.94384449<br>v_wind = v_wind*1.94384449<br><br>temp1 = temp - 273<br><br>temp1@lat2d = lat2d<br>temp1@lon2d = lon2d<br><br>hgt@lat2d = lat2d<br>hgt@lon2d = lon2d<br><br>;pres@lon2d = lon2d<br>;pres@lat2d = lat2d<br><br>u_wind@lon2d = lon2d<br>u_wind@lat2d = lat2d<br><br>spd = spd*1.94384449<br><br>spd@lon2d = lon2d<br>spd@lat2d = lat2d<br><br><br>res = True<br>res@gsnDraw = False<br>res@gsnFrame = False<br>res@cnLineLabelsOn = False<br>res@cnFillOn = True<br>res@tiMainString = "Temperature (C) at 600hpa on Jan 28 at 00:00UTC"<br>res@tiMainFont = "helvetica"<br>res@tiMainFontHeightF = 0.02<br>res@pmTickMarkDisplayMode = "Always"<br>res@mpProjection = "CylindricalEquidistant"<br>res@mpDataBaseVersion = "MediumRes"<br>res@mpOutlineOn = "True"<br>res@lbOrientation = "Vertical"<br>res@tiMainOffsetYF = -0.03<br>res@mpFillOn = False<br>res@mpOutlineOn = True<br>;res@tfDoNDCOverlay = True<br>res@mpMinLatF = 37.85<br>res@mpMaxLatF = 38.50<br>res@mpMinLonF = -120.0<br>res@mpMaxLonF = -119.0<br>res@gsnLeftString = "Geopotential Height (m) and Winds (kts)"<br>res@gsnStringFontHeightF = 0.020<br>res@gsnRightString = ""<br>res@gsnMaximize = True<br>res@mpShapeMode = "FreeAspect"<br>res@lbTitleString = "Temperature (C)"<br>res@lbTitleDirection = "Across"<br>res@lbTitleAngleF = 90<br>res@lbTitlePosition = "Right"<br>res@lbTitleFontHeightF = 0.020<br>res@gsnAddCyclic = False<br>;res@cnLevelSelectionMode = "ManualLevels"<br>;res@cnMinLevelValF = <br>;res@cnMaxLevelValF =<br>;res@cnLevelSpacingF =<br><br>res1 = True<br>res1@gsnDraw = False<br>res1@gsnFrame = False<br>res1@vcWindBarbLineThicknessF = 7.0<br>res1@vcRefLengthF =0.018<br>res1@vcRefMagnitudeF = 10<br>;res1@tfDoNDCOverlay = True<br>res1@vcMinDistanceF = 0.017<br>res1@vcGlyphStyle = "WindBarb"<br>res1@gsnRightString =""<br>res1@gsnLeftString = ""<br>res1@vcRefAnnoOn = False<br>res1@gsnAddCyclic = False<br><br>res2 = True<br>res2@gsnDraw = False<br>res2@gsnFrame = False<br>res2@cnLineColor = "Black"<br>res2@gsnContourLineThicknessesScale = 5.0<br>res2@gsnAddCyclic = False<br>;res2@tfDoNDCOverlay = True<br>res2@gsnLeftString = ""<br>res2@gsnRightString = ""<br>res2@cnInfoLabelOn = False<br>;res2@cnLevelSpacingF = 2<br>res2@cnLineLabelBackgroundColor = "white"<br>res2@cnLineLabelDensityF = 1<br><br>geo = gsn_csm_contour(wks, hgt, res2)<br><br>t600 = gsn_csm_contour_map(wks, temp1, res)<br><br>vector = gsn_csm_vector(wks, u_wind, v_wind, res1)<br><br>overlay(t600, geo)<br>overlay(t600, vector)<br><br>draw(t600)<br>frame(wks)<br><br>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<br>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<br> <br>end<br><br><br></div><div><br></div></div></div>
</div></div><span><temp_narr.png></span>_______________________________________________<br>ncl-talk mailing list<br><a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>List instructions, subscriber options, unsubscribe:<br><a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br></div></blockquote></div><br></div></div></blockquote></div><br></div>