<html><head><meta http-equiv="Content-Type" content="text/html charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">Kerwyn,<div class=""><br class=""></div><div class="">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 class=""><br class=""></div><div class="">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 class=""><br class=""></div><div class="">You can see the wind barbs you had, with 1 geopotential line and some colored contours for your temps.</div><div class=""><br class=""></div><div class="">Plotting such a small area of interest with coarse data what is making the plot difficult.</div><div class=""><br class=""></div><div class="">I hope that offers some insight,</div><div class="">-Alex</div><div class=""><br class=""></div><div class=""><img height="539" width="539" apple-width="yes" apple-height="yes" apple-inline="yes" id="8C59B74C-50DF-495C-ABBD-804838772546" src="cid:9F4202AE-4467-4705-9E48-8737CA5DDD9E" class=""></div><div class=""><br class=""></div><div class=""><br class=""><div><blockquote type="cite" class=""><div class="">On Aug 2, 2016, at 5:43 PM, Kerwyn Texeira <<a href="mailto:ktish86@gmail.com" class="">ktish86@gmail.com</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div dir="ltr" class=""><div class=""><div class=""><div class="">Hi ncl-talk,<br class=""><br class=""></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 class=""><br class=""></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 class=""><br class=""></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 class=""><div class=""><div class=""><br class=""></div><div class="">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 class=""><br class=""></div><div class="">Thanks<br class=""></div><div class=""><br class=""></div><div class="">Some outputs:<br class=""><br class="">Variable: temp<br class="">Type: float<br class="">Total Size: 386692 bytes<br class=""> 96673 values<br class="">Number of Dimensions: 2<br class="">Dimensions and sizes: [gridx_221 | 277] x [gridy_221 | 349]<br class="">Coordinates:<br class="">Number Of Attributes: 15<br class=""> lv_ISBL3 : 600<br class=""> sub_center : The North American Regional Reanalysis (NARR) Project<br class=""> center : US National Weather Service - NCEP (WMC)<br class=""> long_name : Temperature<br class=""> units : K<br class=""> _FillValue : 1e+20<br class=""> coordinates : gridlat_221 gridlon_221<br class=""> level_indicator : 100<br class=""> grid_number : 221<br class=""> parameter_table_version : 131<br class=""> parameter_number : 11<br class=""> model : North American Regional Reanalysis (NARR)<br class=""> forecast_time : 0<br class=""> forecast_time_units : hours<br class=""> initial_time : 01/28/2014 (00:00)<br class=""><br class="">Variable: u_wind<br class="">Type: float<br class="">Total Size: 386692 bytes<br class=""> 96673 values<br class="">Number of Dimensions: 2<br class="">Dimensions and sizes: [gridx_221 | 277] x [gridy_221 | 349]<br class="">Coordinates:<br class="">Number Of Attributes: 15<br class=""> lv_ISBL3 : 600<br class=""> sub_center : The North American Regional Reanalysis (NARR) Project<br class=""> center : US National Weather Service - NCEP (WMC)<br class=""> long_name : u-component of wind<br class=""> units : m/s<br class=""> _FillValue : 1e+20<br class=""> coordinates : gridlat_221 gridlon_221<br class=""> level_indicator : 100<br class=""> grid_number : 221<br class=""> parameter_table_version : 131<br class=""> parameter_number : 33<br class=""> model : North American Regional Reanalysis (NARR)<br class=""> forecast_time : 0<br class=""> forecast_time_units : hours<br class=""> initial_time : 01/28/2014 (00:00)<br class=""><br class="">Variable: v_wind<br class="">Type: float<br class="">Total Size: 386692 bytes<br class=""> 96673 values<br class="">Number of Dimensions: 2<br class="">Dimensions and sizes: [gridx_221 | 277] x [gridy_221 | 349]<br class="">Coordinates:<br class="">Number Of Attributes: 15<br class=""> lv_ISBL3 : 600<br class=""> sub_center : The North American Regional Reanalysis (NARR) Project<br class=""> center : US National Weather Service - NCEP (WMC)<br class=""> long_name : v-component of wind<br class=""> units : m/s<br class=""> _FillValue : 1e+20<br class=""> coordinates : gridlat_221 gridlon_221<br class=""> level_indicator : 100<br class=""> grid_number : 221<br class=""> parameter_table_version : 131<br class=""> parameter_number : 34<br class=""> model : North American Regional Reanalysis (NARR)<br class=""> forecast_time : 0<br class=""> forecast_time_units : hours<br class=""> initial_time : 01/28/2014 (00:00)<br class=""><br class="">Variable: hgt<br class="">Type: float<br class="">Total Size: 386692 bytes<br class=""> 96673 values<br class="">Number of Dimensions: 2<br class="">Dimensions and sizes: [gridx_221 | 277] x [gridy_221 | 349]<br class="">Coordinates:<br class="">Number Of Attributes: 15<br class=""> lv_ISBL3 : 600<br class=""> sub_center : The North American Regional Reanalysis (NARR) Project<br class=""> center : US National Weather Service - NCEP (WMC)<br class=""> long_name : Geopotential height<br class=""> units : gpm<br class=""> _FillValue : 1e+20<br class=""> coordinates : gridlat_221 gridlon_221<br class=""> level_indicator : 100<br class=""> grid_number : 221<br class=""> parameter_table_version : 131<br class=""> parameter_number : 7<br class=""> model : North American Regional Reanalysis (NARR)<br class=""> forecast_time : 0<br class=""> forecast_time_units : hours<br class=""> initial_time : 01/28/2014 (00:00)<br class=""><br class=""></div><div class="">My Script:<br class=""><br class="">load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"<br class="">load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"<br class="">load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"<br class=""><br class="">begin<br class="">a = addfile("./narr-a_221_20140128_0000_000.grb", "r")<br class=""><br class="">wks_type = "png"<br class="">wks_type@wkWidth = 2500<br class="">wks_type@wkHeight = 2500<br class=""><br class="">wks= gsn_open_wks(wks_type, "temp_narr")<br class="">gsn_define_colormap(wks,"matlab_jet")<br class=""><br class="">lat2d = a->gridlat_221<br class="">lon2d = a->gridlon_221<br class="">hgt = a->HGT_221_ISBL({600},:,:)<br class="">temp = a->TMP_221_ISBL({600},:,:)<br class="">pres = a->PRES_221_HYBL<br class="">u_wind = a->U_GRD_221_ISBL({600},:,:)<br class="">v_wind = a->V_GRD_221_ISBL({600},:,:)<br class="">lev = a->lv_ISBL3 ; 29 levels<br class=""><br class="">printVarSummary(temp)<br class="">printVarSummary(u_wind)<br class="">printVarSummary(v_wind)<br class="">printVarSummary(hgt)<br class=""><br class="">spd = (u_wind*u_wind + v_wind*v_wind)^(0.5) ; m/s<br class=""><br class="">pres = pres*0.01 ; converts to hpa<br class=""><br class="">u_wind = u_wind*1.94384449<br class="">v_wind = v_wind*1.94384449<br class=""><br class="">temp1 = temp - 273<br class=""><br class="">temp1@lat2d = lat2d<br class="">temp1@lon2d = lon2d<br class=""><br class="">hgt@lat2d = lat2d<br class="">hgt@lon2d = lon2d<br class=""><br class="">;pres@lon2d = lon2d<br class="">;pres@lat2d = lat2d<br class=""><br class="">u_wind@lon2d = lon2d<br class="">u_wind@lat2d = lat2d<br class=""><br class="">spd = spd*1.94384449<br class=""><br class="">spd@lon2d = lon2d<br class="">spd@lat2d = lat2d<br class=""><br class=""><br class="">res = True<br class="">res@gsnDraw = False<br class="">res@gsnFrame = False<br class="">res@cnLineLabelsOn = False<br class="">res@cnFillOn = True<br class="">res@tiMainString = "Temperature (C) at 600hpa on Jan 28 at 00:00UTC"<br class="">res@tiMainFont = "helvetica"<br class="">res@tiMainFontHeightF = 0.02<br class="">res@pmTickMarkDisplayMode = "Always"<br class="">res@mpProjection = "CylindricalEquidistant"<br class="">res@mpDataBaseVersion = "MediumRes"<br class="">res@mpOutlineOn = "True"<br class="">res@lbOrientation = "Vertical"<br class="">res@tiMainOffsetYF = -0.03<br class="">res@mpFillOn = False<br class="">res@mpOutlineOn = True<br class="">;res@tfDoNDCOverlay = True<br class="">res@mpMinLatF = 37.85<br class="">res@mpMaxLatF = 38.50<br class="">res@mpMinLonF = -120.0<br class="">res@mpMaxLonF = -119.0<br class="">res@gsnLeftString = "Geopotential Height (m) and Winds (kts)"<br class="">res@gsnStringFontHeightF = 0.020<br class="">res@gsnRightString = ""<br class="">res@gsnMaximize = True<br class="">res@mpShapeMode = "FreeAspect"<br class="">res@lbTitleString = "Temperature (C)"<br class="">res@lbTitleDirection = "Across"<br class="">res@lbTitleAngleF = 90<br class="">res@lbTitlePosition = "Right"<br class="">res@lbTitleFontHeightF = 0.020<br class="">res@gsnAddCyclic = False<br class="">;res@cnLevelSelectionMode = "ManualLevels"<br class="">;res@cnMinLevelValF = <br class="">;res@cnMaxLevelValF =<br class="">;res@cnLevelSpacingF =<br class=""><br class="">res1 = True<br class="">res1@gsnDraw = False<br class="">res1@gsnFrame = False<br class="">res1@vcWindBarbLineThicknessF = 7.0<br class="">res1@vcRefLengthF =0.018<br class="">res1@vcRefMagnitudeF = 10<br class="">;res1@tfDoNDCOverlay = True<br class="">res1@vcMinDistanceF = 0.017<br class="">res1@vcGlyphStyle = "WindBarb"<br class="">res1@gsnRightString =""<br class="">res1@gsnLeftString = ""<br class="">res1@vcRefAnnoOn = False<br class="">res1@gsnAddCyclic = False<br class=""><br class="">res2 = True<br class="">res2@gsnDraw = False<br class="">res2@gsnFrame = False<br class="">res2@cnLineColor = "Black"<br class="">res2@gsnContourLineThicknessesScale = 5.0<br class="">res2@gsnAddCyclic = False<br class="">;res2@tfDoNDCOverlay = True<br class="">res2@gsnLeftString = ""<br class="">res2@gsnRightString = ""<br class="">res2@cnInfoLabelOn = False<br class="">;res2@cnLevelSpacingF = 2<br class="">res2@cnLineLabelBackgroundColor = "white"<br class="">res2@cnLineLabelDensityF = 1<br class=""><br class="">geo = gsn_csm_contour(wks, hgt, res2)<br class=""><br class="">t600 = gsn_csm_contour_map(wks, temp1, res)<br class=""><br class="">vector = gsn_csm_vector(wks, u_wind, v_wind, res1)<br class=""><br class="">overlay(t600, geo)<br class="">overlay(t600, vector)<br class=""><br class="">draw(t600)<br class="">frame(wks)<br class=""><br class="">;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<br class="">;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<br class=""> <br class="">end<br class=""><br class=""><br class=""></div><div class=""><br class=""></div></div></div>
<span id="cid:4CA81EA9-2B3F-4DC1-BBD8-56DFC9B7BA8A"><temp_narr.png></span>_______________________________________________<br class="">ncl-talk mailing list<br class=""><a href="mailto:ncl-talk@ucar.edu" class="">ncl-talk@ucar.edu</a><br class="">List instructions, subscriber options, unsubscribe:<br class="">http://mailman.ucar.edu/mailman/listinfo/ncl-talk<br class=""></div></blockquote></div><br class=""></div></body></html>