<div dir="ltr"><div><div><div><div><div><div><div><div>Hi ncl-talk,<br><br></div>I&#39;m trying to plot and overlay of (NARR) North American Regional Reanalysis temperature, geopotential height contours, and wind barbs at 600mb but I&#39;m running into some problems.<br><br></div>1.  Where on the script should I mention the pressure level eg 600hpa so that both the temperature, wind barbs and geopotential heights can be plotted at this level that level?<br><br></div>2. How do I zoom into an area?  I specified the min lat, lon and max lat, lon but it seems not to zoom into the area of interest?<br><br></div>3. I do not understand why I&#39;m getting the following error:<br><br> (0)     gsn_csm_contour_map_ce: Fatal: the input data array must be 1D or 2D<br>fatal:Illegal right-hand side type for assignment<br>fatal:[&quot;Execute.c&quot;:8578]:Execute: Error occurred at or near line 113 in file narr_temp.ncl<br><br></div>My script:<br><br>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl&quot;<br>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl&quot;<br>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl&quot;<br><br>begin<br>a = addfile(&quot;./narr-a_221_20140128_0000_000.grb&quot;, &quot;r&quot;)<br><br>wks_type = &quot;png&quot;<br>wks_type@wkWidth = 2500<br>wks_type@wkHeight = 2500<br><br>wks= gsn_open_wks(wks_type, &quot;temp_narr&quot;)<br>gsn_define_colormap(wks,&quot;matlab_jet&quot;)<br><br>lat2d = a-&gt;gridlat_221<br>lon2d = a-&gt;gridlon_221<br>hgt = a-&gt;HGT_221_ISBL<br>temp = a-&gt;TMP_221_ISBL<br>pres = a-&gt;PRES_221_HYBL<br>u_wind = a-&gt;U_GRD_221_ISBL<br>v_wind = a-&gt;V_GRD_221_ISBL<br>lev    = a-&gt;lv_ISBL3            ; 29 levels<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>temp = temp - 273<br><br>temp@lat2d = lat2d<br>temp@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>v_wind@lat2d = lat2d<br>v_wind@lon2d = lon2d<br><br><br>res = True<br>res@gsnDraw = False<br>res@gsnFrame = False<br>res@cnLineLabelsOn = False<br>res@cnFillOn = True<br>res@tiMainString = &quot;Temperature (C) at 700hpa on Jan 28 at 00:00UTC&quot;<br>res@tiMainFont = &quot;helvetica&quot;<br>res@tiMainFontHeightF = 0.02<br>res@pmTickMarkDisplayMode = &quot;Always&quot;<br>res@mpProjection = &quot;CylindricalEquidistant&quot;<br>res@mpOutlineOn = &quot;True&quot;<br>res@lbOrientation = &quot;Vertical&quot;<br>res@tiMainOffsetYF = -0.03<br>res@mpFillOn  = False<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 = &quot;Geopotential Height (m) and Winds (kts)&quot;<br>res@gsnStringFontHeightF = 0.020<br>res@gsnRightString = &quot;&quot;<br>res@gsnMaximize = True<br>res@mpShapeMode = &quot;FreeAspect&quot;<br>res@lbTitleString = &quot;Temperature (C)&quot;<br>res@lbTitlePosition = &quot;Right&quot;<br>res@lbTitleFontHeightF = 0.020<br>res@gsnAddCyclic = False<br>;res@cnLevelSelectionMode = &quot;ManualLevels&quot;<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@vcMinDistanceF = 0.05<br>res1@vcGLyphStyle = &quot;WindBarb&quot;<br>res1@gsnRightString =&quot;&quot;<br>res1@vcRefAnnoOn = False<br>res1@gsnAddCyclic = False<br><br>res2 = True<br>res2@gsnDraw = False<br>res2@gsnFrame = False<br>res2@cnLineColor = &quot;Black&quot;<br>res2@gsnContourLineThicknessesScale = 5.0<br>res2@gsnAddCyclic = False<br>res2@tfDoNDCOverlay = True<br>res2@gsnLeftString = &quot;&quot;<br>res2@gsnRightString = &quot;&quot;<br>res2@cnInfoLabelOn = False<br>;res2@cnLevelSpacingF = 2<br>res2@cnLineLabelBackgroundColor = &quot;white&quot;<br>res2@cnLineLabelDensityF = 1<br><br>geo = gsn_csm_contour(wks, hgt, res2)<br><br>temp = gsn_csm_contour_map(wks, temp, res)<br><br>vector = gsn_csm_vector(wks, u_wind, v_wind, res1)<br><br>overlay(temp, geo)<br>overlay(temp, vector)<br><br>draw(temp)<br>frame(wks)<br><br>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<br>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<br> <br>end<br><br></div>Whatever help anyone can give to me will be greatly appreciated. I thank you very much in advance!!!.<br></div></div></div>