<div dir="ltr"><div class="gmail_default" style="font-size:small">Kerwyn,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">If you set both the lat2d/lon2d special attributes *and* you set tfDoNDCOverlay to True, then you are giving conflicting information to the gsn_csm_xxx_map scripts.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">When you set res@tfDoNDCOverlay to True, you are telling NCL that your data is already projected into some map projection, and NCL assumes that you have set all the necessary resources for this map projection. You should NOT set lat2d/lon2d when tfDoNDCOverlay is set to True.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">To fix your script, try commenting out all "tfDoNDCOverlay" settings.</div><div class="gmail_default" style="font-size:small"><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Aug 1, 2016 at 6:24 PM, Kerwyn Texeira <span dir="ltr"><<a href="mailto:ktish86@gmail.com" target="_blank">ktish86@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div><div>Hi ncl-talk,<br><br></div>I made all the changes but I'm still unable to plot a zoomed in an overlay of both winds, temperature, and geopotential height.<br><br></div>I'm not getting any errors though but it just will not plot or plot right. I have attached to figure.<br><br></div><div>Whatever help you can offer me will be greatly appreciated. <br><br></div><div>Thanks!!!<br></div><div><br></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>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>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 = "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@lbTitlePosition = "Right"<br>res@lbTitleFontHeightF = 0.015<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@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>
<br>_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>