<html><head><meta http-equiv="Content-Type" content="text/html charset=windows-1252"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">Hi Mira,<div class=""><br class=""></div><div class="">You aren’t seeing vectors or map because they a buried under the filled contours. You may want to flip the plots around and do the contours as a map and overlay the vectors on top of the contours.</div><div class=""><br class=""></div><div class="">Upon further analysis of your script, you were doing the overlay on contour plots but then using different contour plots for the panel plot. Below I have altered your script, try it this way and see how it goes.</div><div class=""><br class=""></div><div class="">I am not sure if you have made any more recent changes to this script but this is using the last one you sent.</div><div class=""><br class=""></div><div class="">Hope that helps,</div><div class="">-Alex</div><div class=""><br class=""></div><div class=""><br class=""></div><div class="">;*****************************************<br class="">; ; plot average winds overlayed on geopotental filled contours for the<br class="">years that are extreme<br class="">; the years are found with my matlab script findExtremeYrs.m<br class="">; this particlar script does the years that correspond to high correlation<br class="">with DJF, SE precip and NAO<br class="">;*****************************************<br class="">; the original data goes from 1948 January to April 2015.<br class="">; I will cut it out to 2014 December, so we have something divisible by 12<br class="">so we can do seasonal averages...<br class=""><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="">load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"<br class="">;************************************************<br class="">begin<br class="">;************************************************<br class="">; read in netCDF file s<br class="">;************************************************<br class="">a = addfile("uwnd.mon.mean.alllevels.nc","r") ; u winds<br class="">b = addfile("vwnd.mon.mean.alllevels.nc","r") ; v winds<br class=""><br class="">c = addfile("../Geopotential/hgt.mon.mean.nc","r") ; geopotential heights.<br class=""><br class="">;************************************************<br class="">; read in zonal [u] and meridional [v] winds (July)<br class="">;************************************************<br class=""><br class="">u = a->uwnd(0:803,{500},{45:90},{270:357.5})<br class="">v = b->vwnd(0:803,{500},{45:90},{270:357.5}) ; Get u, v, time (1),level<br class="">(1000hpa),latitude(-90:90) and longitude(0:360) data.<br class="">z = c->hgt(0:803,{500},{45:90},{270:357.5}) ; get geopotenial heights...<br class=""><br class="">printVarSummary(u)<br class="">printVarSummary(v)<br class="">printVarSummary(z)<br class=""><br class="">; Calculate the seasonal averages.<br class="">uDJF = month_to_season(u, "DJF")<br class="">vDJF = month_to_season(v, "DJF")<br class="">zDJF = month_to_season(z, "DJF")<br class=""><br class="">printVarSummary(uDJF)<br class="">printVarSummary(vDJF)<br class="">printVarSummary(zDJF)<br class=""><br class="">; from the matlab script i wrote: findExtremeYrs, i pulled out the extreme<br class="">years (> or < 1std) that i want to average and plot here.<br class=""><br class="">; for ans = 4 (NAO)<br class="">; yearList_hi = 1973 1975 1983 1989 1995 <br class="">2000 2007 2012<br class="">; yearList_lo = 1963 1964 1965 1969 1977 <br class="">1979 1996 1997 2010 2011<br class=""><br class=""><br class="">; this data starts at 1948 (this is index 0), so 1953=5, 1963=15 etc.<br class=""><br class="">uDJF_NAO_hi = uDJF((/25,27,35,41,47,52,59,64/),:,:)<br class="">uDJF_NAO_lo = uDJF((/15,16,17,21,29,31,48,49,62,63/),:,:)<br class=""><br class="">vDJF_NAO_hi = vDJF((/25,27,35,41,47,52,59,64/),:,:)<br class="">vDJF_NAO_lo = vDJF((/15,16,17,21,29,31,48,49,62,63/),:,:)<br class=""><br class="">zDJF_NAO_hi = zDJF((/25,27,35,41,47,52,59,64/),:,:)<br class="">zDJF_NAO_lo = zDJF((/15,16,17,21,29,31,48,49,62,63/),:,:)<br class=""><br class=""><br class="">uAvgTime_hi = dim_avg_n_Wrap(uDJF_NAO_hi,0)<br class="">uAvgTime_lo = dim_avg_n_Wrap(uDJF_NAO_lo,0)<br class=""><br class="">printVarSummary(uAvgTime_hi)<br class="">printVarSummary(uAvgTime_lo)<br class=""><br class="">vAvgTime_hi = dim_avg_n_Wrap(vDJF_NAO_hi,0)<br class="">vAvgTime_lo = dim_avg_n_Wrap(vDJF_NAO_lo,0)<br class=""><br class="">printVarSummary(vAvgTime_hi)<br class="">printVarSummary(vAvgTime_lo)<br class=""><br class="">zAvgTime_hi = dim_avg_n_Wrap(zDJF_NAO_hi,0)<br class="">zAvgTime_lo = dim_avg_n_Wrap(zDJF_NAO_lo,0)<br class=""><br class="">printVarSummary(zAvgTime_hi)<br class="">printVarSummary(zAvgTime_lo)<br class=""><br class="">; dirty way to copy metadata over first.<br class="">diff_u = uAvgTime_hi;<br class="">diff_v = vAvgTime_hi;<br class="">diff_z = zAvgTime_hi;<br class=""><br class="">diff_u = uAvgTime_hi - uAvgTime_lo<br class="">diff_v = vAvgTime_hi - vAvgTime_lo<br class="">diff_z = zAvgTime_hi - zAvgTime_lo<br class=""><br class="">printVarSummary(diff_u)<br class="">printVarSummary(diff_v)<br class="">printVarSummary(diff_z)<br class=""><br class="">;************************************************<br class="">; create plot<br class="">;************************************************<br class="">wks = gsn_open_wks("ps","Panel_NAO_z_500") <span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span>; open a ps file<br class="">plot = new(3,graphic)<span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span>; create a plot array<br class=""><br class="">;---- set common resources for all plots<br class="">res <span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span>= True<br class="">res@gsnDraw<span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span>= False <span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span>; dont draw<br class="">res@gsnFrame<span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span>= False<span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span>; dont advance frame<br class="">res@cnInfoLabelOn<span class="Apple-tab-span" style="white-space: pre;"> </span>= False<span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span>; trn off cn info label<br class=""></div><div class="">res@gsnAddCyclic = False<span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span>; has to do with wrapping the longitude at 0/360<br class="">;************************************************<br class="">; Choose a subregion<br class="">;************************************************<br class="">res@mpMaxLatF = 90 <span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span>;maximum latitude<br class="">res@mpMinLatF = 45 <span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span>;minimum latitude<br class="">res@mpMaxLonF = 357.5 <span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span>;<span class="Apple-tab-span" style="white-space: pre;"> </span>;maximum longitude<br class="">res@mpMinLonF = 270 <span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span>;minimum longitude</div><div class=""><br class=""></div><div class="">;************************************************<br class="">; Map Options<br class="">;************************************************</div><div class="">res@mpOutlineBoundarySets = "National" <span class="Apple-tab-span" style="white-space: pre;"> </span>; turn on country boundaries<br class="">res@mpGeophysicalLineColor = "Navy" <span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span>; color of cont. outlines<br class="">res@mpGeophysicalLineThicknessF = 1.5 <span class="Apple-tab-span" style="white-space: pre;"> </span>; thickness of outlines<br class=""><br class="">;***********************************************<br class="">; ----wind vector plot<br class="">;***********************************************<br class="">vcres = res<br class="">vcres@vcRefAnnoOrthogonalPosF = -1.0 <span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span>; move ref vector up<br class="">vcres@vcRefMagnitudeF = 10.0 <span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span>; define vector ref mag<br class="">vcres@vcRefLengthF = 0.045 <span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span>; define length of vec ref<br class="">vcres@vcGlyphStyle = "CurlyVector" <span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span>; turn on curly vectors<br class="">vcres@vcMinDistanceF = 0.017<br class="">vcres@mpFillOn = False <span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span>; turn off gray fill<br class=""><br class="">;vcres@gsnLeftString = "DJF High NAO"<br class="">winds_hi = gsn_csm_vector(wks,uAvgTime_hi,vAvgTime_hi,vcres)<br class="">;vcres@gsnLeftString = "DJF Low NAO"<br class="">winds_lo = gsn_csm_vector(wks,uAvgTime_lo,vAvgTime_lo,vcres)<br class="">;vcres@gsnLeftString = "Difference of High - Low"<br class="">winds_diff = gsn_csm_vector(wks, diff_u, diff_v,vcres)<br class="">;************************************************<br class="">;---- geopotential height filled contour plot<br class="">;***********************************************<br class="">zfres = res<br class="">zfres@cnFillOn = True<br class="">;zfres@cnLevelSelectionMode = "ExplicitLevels<br class="">;zfres@cnLevels = ispan(-20,90,5)<br class="">zfres@lbLabelFontHeightF = 0.015<br class="">zfres@lbOrientation = "Vertical"<br class="">zfres@pmLabelBarOrthogonalPosF = -0.005<br class=""><br class="">plot(0) = gsn_csm_contour_map_ce(wks,zAvgTime_hi,zfres)<br class="">plot(1) = gsn_csm_contour_map_ce(wks,zAvgTime_lo,zfres)<br class="">plot(2) = gsn_csm_contour_map_ce(wks,diff_z,zfres)<br class=""><br class="">overlay(plot(0),winds_hi)<br class="">overlay(plot(1),winds_lo)<br class="">overlay(plot(2),winds_diff)<br class=""><br class="">;************************************************<br class="">; create panel<br class="">;************************************************<br class="">resP = True<span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span>; modify the panel plot<br class="">resP@txString = "NAO 500mb"<br class="">gsn_panel(wks,plot,(/3,1/),resP)<span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span>; now draw as one plot;<br class=""><br class=""><br class="">end</div><div class=""><br class=""></div><div class=""><br class=""><div><blockquote type="cite" class=""><div class="">On Jul 29, 2015, at 1:57 PM, <a href="mailto:mberdahl@envsci.rutgers.edu" class="">mberdahl@envsci.rutgers.edu</a> wrote:</div><br class="Apple-interchange-newline"><div class=""><span style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">Hi Alex,</span><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">Great - thanks for the suggestions! That does help, but does not totally</span><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">solve my problem. Now I get three panel plots of the correct region and</span><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">all on one page which is great. But only the filled contours show. There</span><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">are no overlaid vectors or map of coastal boundaries. Any suggestions</span><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">would be greatly appreciated.</span><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">Thanks so much!</span><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">Mira</span><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><blockquote type="cite" style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class="">Hi Mira,<br class=""><br class="">I think the issue here is that you set a bunch of resources as res but<br class="">then set vres to True instead of res. I think switching it to vres = res<br class="">will give you the maps you were expecting. I see you also set res = True<br class="">twice, might be helpful to remove the double.<br class=""><br class="">Hope that helps,<br class="">-Alex<br class=""><br class=""><br class=""><blockquote type="cite" class="">On Jul 27, 2015, at 9:39 PM, <a href="mailto:mberdahl@envsci.rutgers.edu" class="">mberdahl@envsci.rutgers.edu</a> wrote:<br class=""><br class="">Hi Alex,<br class="">Thanks for the suggestion. I tried this and now I do not have the same<br class="">error message any more, so that is good. However, I am now just getting<br class="">3<br class="">blank plots of the whole world, even though I'd like plots with overlaid<br class="">geopotential height and winds, for just a region over Greenland. These<br class="">are<br class="">followed by a fourth figure that has filled contours of the correct<br class="">region, but no map or wind vectors. I'll attach the figures as an<br class="">attachment here. They are not panel plots as I want, but each figure is<br class="">on its own page. I'll copy the updated code below.<br class="">Any thoughts?<br class="">Thanks!<br class="">Mira<br class=""><br class=""><br class=""><br class=""><br class="">;*****************************************<br class="">; ; plot average winds overlayed on geopotental filled contours for the<br class="">years that are extreme<br class="">; the years are found with my matlab script findExtremeYrs.m<br class="">; this particlar script does the years that correspond to high<br class="">correlation<br class="">with DJF, SE precip and NAO<br class="">;*****************************************<br class="">; the original data goes from 1948 January to April 2015.<br class="">; I will cut it out to 2014 December, so we have something divisible by<br class="">12<br class="">so we can do seasonal averages...<br class=""><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="">load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"<br class="">;************************************************<br class="">begin<br class="">;************************************************<br class="">; read in netCDF file s<br class="">;************************************************<br class="">a = addfile("uwnd.mon.mean.alllevels.nc","r") ; u winds<br class="">b = addfile("vwnd.mon.mean.alllevels.nc","r") ; v winds<br class=""><br class="">c = addfile("../Geopotential/hgt.mon.mean.nc","r") ; geopotential<br class="">heights.<br class=""><br class="">;************************************************<br class="">; read in zonal [u] and meridional [v] winds (July)<br class="">;************************************************<br class=""><br class="">u = a->uwnd(0:803,{500},{45:90},{270:357.5})<br class="">v = b->vwnd(0:803,{500},{45:90},{270:357.5}) ; Get u, v, time (1),level<br class="">(1000hpa),latitude(-90:90) and longitude(0:360) data.<br class="">z = c->hgt(0:803,{500},{45:90},{270:357.5}) ; get geopotenial<br class="">heights...<br class=""><br class="">printVarSummary(u)<br class="">printVarSummary(v)<br class="">printVarSummary(z)<br class=""><br class="">; Calculate the seasonal averages.<br class="">uDJF = month_to_season(u, "DJF")<br class="">vDJF = month_to_season(v, "DJF")<br class="">zDJF = month_to_season(z, "DJF")<br class=""><br class="">printVarSummary(uDJF)<br class="">printVarSummary(vDJF)<br class="">printVarSummary(zDJF)<br class=""><br class="">; from the matlab script i wrote: findExtremeYrs, i pulled out the<br class="">extreme<br class="">years (> or < 1std) that i want to average and plot here.<br class=""><br class="">; for ans = 4 (NAO)<br class="">; yearList_hi = 1973 1975 1983 1989 1995<br class="">2000 2007 2012<br class="">; yearList_lo = 1963 1964 1965 1969 1977<br class="">1979 1996 1997 2010 2011<br class=""><br class=""><br class="">; this data starts at 1948 (this is index 0), so 1953=5, 1963=15 etc.<br class=""><br class="">uDJF_NAO_hi = uDJF((/25,27,35,41,47,52,59,64/),:,:)<br class="">uDJF_NAO_lo = uDJF((/15,16,17,21,29,31,48,49,62,63/),:,:)<br class=""><br class="">vDJF_NAO_hi = vDJF((/25,27,35,41,47,52,59,64/),:,:)<br class="">vDJF_NAO_lo = vDJF((/15,16,17,21,29,31,48,49,62,63/),:,:)<br class=""><br class="">zDJF_NAO_hi = zDJF((/25,27,35,41,47,52,59,64/),:,:)<br class="">zDJF_NAO_lo = zDJF((/15,16,17,21,29,31,48,49,62,63/),:,:)<br class=""><br class=""><br class="">uAvgTime_hi = dim_avg_n_Wrap(uDJF_NAO_hi,0)<br class="">uAvgTime_lo = dim_avg_n_Wrap(uDJF_NAO_lo,0)<br class=""><br class="">printVarSummary(uAvgTime_hi)<br class="">printVarSummary(uAvgTime_lo)<br class=""><br class="">vAvgTime_hi = dim_avg_n_Wrap(vDJF_NAO_hi,0)<br class="">vAvgTime_lo = dim_avg_n_Wrap(vDJF_NAO_lo,0)<br class=""><br class="">printVarSummary(vAvgTime_hi)<br class="">printVarSummary(vAvgTime_lo)<br class=""><br class="">zAvgTime_hi = dim_avg_n_Wrap(zDJF_NAO_hi,0)<br class="">zAvgTime_lo = dim_avg_n_Wrap(zDJF_NAO_lo,0)<br class=""><br class="">printVarSummary(zAvgTime_hi)<br class="">printVarSummary(zAvgTime_lo)<br class=""><br class="">; dirty way to copy metadata over first.<br class="">diff_u = uAvgTime_hi;<br class="">diff_v = vAvgTime_hi;<br class="">diff_z = zAvgTime_hi;<br class=""><br class="">diff_u = uAvgTime_hi - uAvgTime_lo<br class="">diff_v = vAvgTime_hi - vAvgTime_lo<br class="">diff_z = zAvgTime_hi - zAvgTime_lo<br class=""><br class="">printVarSummary(diff_u)<br class="">printVarSummary(diff_v)<br class="">printVarSummary(diff_z)<br class=""><br class="">;************************************************<br class="">; create plot<br class="">;************************************************<br class="">wks = gsn_open_wks("ps","Panel_NAO_z_500")<span class="Apple-converted-space"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span>; open a ps file<br class="">plot = new(3,graphic)<span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span>; create a plot array<br class=""><br class="">;---- set common resources for all plots<br class="">res<span class="Apple-converted-space"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span>= True<br class="">res@gsnDraw<span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span>= False<span class="Apple-converted-space"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span>; dont draw<br class="">res@gsnFrame<span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span>= False<span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span>; dont advance frame<br class="">res@cnInfoLabelOn<span class="Apple-tab-span" style="white-space: pre;"> </span>= False<span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span>; trn off cn info label<br class="">res = True<span class="Apple-converted-space"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span>; plot mods desired<br class="">res@gsnAddCyclic = False<span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span>; has to do with wrapping the longitude at<br class="">0/360<br class="">;************************************************<br class="">; Choose a subregion<br class="">;************************************************<br class="">res@mpMaxLatF = 90<span class="Apple-converted-space"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span>;maximum latitude<br class="">res@mpMinLatF = 45<span class="Apple-converted-space"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span>;minimum latitude<br class="">res@mpMaxLonF = 357.5<span class="Apple-converted-space"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span>;<span class="Apple-tab-span" style="white-space: pre;"> </span>;maximum longitude<br class="">res@mpMinLonF = 270<span class="Apple-converted-space"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span>;minimum longitude<br class=""><br class="">;***********************************************<br class="">; ----wind vector plot<br class="">;***********************************************<br class="">vcres = True<br class="">vcres@vcRefAnnoOrthogonalPosF = -1.0<span class="Apple-converted-space"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span>; move ref vector up<br class="">vcres@vcRefMagnitudeF = 10.0<span class="Apple-converted-space"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span>; define vector ref mag<br class="">vcres@vcRefLengthF = 0.045<span class="Apple-converted-space"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span>; define length of vec ref<br class="">vcres@vcGlyphStyle = "CurlyVector"<span class="Apple-converted-space"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span>; turn on curly vectors<br class="">vcres@vcMinDistanceF = 0.017<br class="">vcres@mpFillOn = False<span class="Apple-converted-space"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span>; turn off gray fill<br class="">vcres@mpOutlineBoundarySets = "National"<span class="Apple-converted-space"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span>; turn on country boundaries<br class="">vcres@mpGeophysicalLineColor = "Navy"<span class="Apple-converted-space"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span>; color of cont. outlines<br class="">vcres@mpGeophysicalLineThicknessF = 1.5<span class="Apple-converted-space"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span>; thickness of outlines<br class=""><br class=""><br class="">;vcres@gsnLeftString = "DJF High NAO"<br class="">winds_hi = gsn_csm_vector_map_ce(wks,uAvgTime_hi,vAvgTime_hi,vcres)<br class="">;vcres@gsnLeftString = "DJF Low NAO"<br class="">winds_lo = gsn_csm_vector_map_ce(wks,uAvgTime_lo,vAvgTime_lo,vcres)<br class="">;vcres@gsnLeftString = "Difference of High - Low"<br class="">winds_diff = gsn_csm_vector_map_ce(wks, diff_u, diff_v,vcres)<br class="">;************************************************<br class="">;---- geopotential height filled contour plot<br class="">;***********************************************<br class="">zfres = res<br class="">zfres@cnFillOn = True<br class="">;zfres@cnLevelSelectionMode = "ExplicitLevels<br class="">;zfres@cnLevels = ispan(-20,90,5)<br class="">zfres@lbLabelFontHeightF = 0.015<br class="">zfres@lbOrientation = "Vertical"<br class="">zfres@pmLabelBarOrthogonalPosF = -0.005<br class=""><br class=""><br class="">contour_zf_hi = gsn_csm_contour(wks,zAvgTime_hi,zfres)<br class="">contour_zf_lo = gsn_csm_contour(wks,zAvgTime_lo,zfres)<br class="">contour_zf_diff = gsn_csm_contour(wks,diff_z,zfres)<br class=""><br class="">plot(0) = gsn_csm_contour(wks,zAvgTime_hi,zfres)<br class="">plot(1) = gsn_csm_contour(wks,zAvgTime_lo,zfres)<br class="">plot(2) = gsn_csm_contour(wks,diff_z,zfres)<br class=""><br class="">overlay(contour_zf_hi,winds_hi)<br class="">overlay(contour_zf_lo,winds_lo)<br class="">overlay(contour_zf_diff,winds_diff)<br class=""><br class=""><br class="">;************************************************<br class="">; create panel<br class="">;************************************************<br class="">resP = True<span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span>; modify the panel plot<br class="">resP@txString = "NAO 500mb"<br class="">gsn_panel(wks,plot,(/3,1/),resP)<span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span>; now draw as one plot;<br class=""><br class=""><br class="">end<br class=""><br class=""><br class=""><br class=""><br class=""><br class=""><br class=""><br class=""><br class=""><br class=""><br class=""><br class=""><br class=""><blockquote type="cite" class="">Mira,<br class=""><br class="">Try this:<br class=""><br class=""><br class="">plot(0) = gsn_csm_contour(wks,zAvgTime_hi,zfres)<br class="">plot(1) = gsn_csm_contour(wks,zAvgTime_lo,zfres)<br class="">plot(2) = gsn_csm_contour(wks,diff_z,zfres)<br class=""><br class="">overlay(contour_zf_hi,winds_hi)<br class="">overlay(contour_zf_lo,winds_lo)<br class="">overlay(contour_diff,winds_diff)<br class=""><br class="">I don't believe overlay returns anything and you are trying to save it<br class="">as<br class="">a<br class="">graphic. The code above would replace the lines in your script.<br class=""><br class="">Hope that helps,<br class=""><br class="">-Alex<br class=""><br class=""><br class=""><br class=""><br class="">On Sunday, July 26, 2015, <<a href="mailto:mberdahl@envsci.rutgers.edu" class="">mberdahl@envsci.rutgers.edu</a><br class=""><javascript:_e(%7B%7D,'cvml','<a href="mailto:mberdahl@envsci.rutgers.edu" class="">mberdahl@envsci.rutgers.edu</a><br class=""><<a href="mailto:mberdahl@envsci.rutgers.edu" class="">mailto:mberdahl@envsci.rutgers.edu</a>>');>> wrote:<br class=""><br class=""><blockquote type="cite" class="">Hi all,<br class=""><br class="">I'm trying to adapt a working code I had that was plotting 3 panel<br class="">plots<br class="">of wind vectors in a certain region. Now, all I want to do is add an<br class="">overlay of geopotential height (filled contours) underneath each of<br class="">these<br class="">plots, but the way I'm approaching it is wrong. I've tried<br class="">simplifying<br class="">it<br class="">to do this without panels, just a single plot, but I still don't get<br class="">what<br class="">I need. The error message I am getting now is:<br class=""><br class="">fatal:syntax error: line 153 in file plotWinds_z_NAO_level.ncl before<br class="">or<br class="">near overlay<br class="">plot(0)= overlay<br class="">---------------^<br class=""><br class=""><br class="">My code is copied below. Any suggestions are greatly appreciated.<br class="">Thanks,<br class="">Mira<br class=""><br class="">;*****************************************<br class="">; ; plot average winds overlayed on geopotental filled contours for<br class="">the<br class="">years that are extreme<br class="">; the years are found with my matlab script findExtremeYrs.m<br class="">; this particlar script does the years that correspond to high<br class="">correlation<br class="">with DJF, SE precip and NAO<br class="">;*****************************************<br class="">; the original data goes from 1948 January to April 2015.<br class="">; I will cut it out to 2014 December, so we have something divisible<br class="">by<br class="">12<br class="">so we can do seasonal averages...<br class=""><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="">load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"<br class="">;************************************************<br class="">begin<br class="">;************************************************<br class="">; read in netCDF file s<br class="">;************************************************<br class="">a = addfile("uwnd.mon.mean.alllevels.nc","r") ; u winds<br class="">b = addfile("vwnd.mon.mean.alllevels.nc","r") ; v winds<br class=""><br class="">c = addfile("../Geopotential/hgt.mon.mean.nc","r") ; geopotential<br class="">heights.<br class=""><br class="">;************************************************<br class="">; read in zonal [u] and meridional [v] winds (July)<br class="">;************************************************<br class=""><br class="">u = a->uwnd(0:803,{500},{45:90},{270:357.5})<br class="">v = b->vwnd(0:803,{500},{45:90},{270:357.5}) ; Get u, v, time<br class="">(1),level<br class="">(1000hpa),latitude(-90:90) and longitude(0:360) data.<br class="">z = c->hgt(0:803,{500},{45:90},{270:357.5}) ; get geopotenial<br class="">heights...<br class=""><br class="">printVarSummary(u)<br class="">printVarSummary(v)<br class="">printVarSummary(z)<br class=""><br class="">; Calculate the seasonal averages.<br class="">uDJF = month_to_season(u, "DJF")<br class="">vDJF = month_to_season(v, "DJF")<br class="">zDJF = month_to_season(z, "DJF")<br class=""><br class="">printVarSummary(uDJF)<br class="">printVarSummary(vDJF)<br class="">printVarSummary(zDJF)<br class=""><br class="">; from the matlab script i wrote: findExtremeYrs, i pulled out the<br class="">extreme<br class="">years (> or < 1std) that i want to average and plot here.<br class=""><br class="">; for ans = 4 (NAO)<br class="">; yearList_hi = 1973 1975 1983 1989 1995<br class="">2000 2007 2012<br class="">; yearList_lo = 1963 1964 1965 1969 1977<br class="">1979 1996 1997 2010 2011<br class=""><br class=""><br class="">; this data starts at 1948 (this is index 0), so 1953=5, 1963=15 etc.<br class=""><br class="">uDJF_NAO_hi = uDJF((/25,27,35,41,47,52,59,64/),:,:)<br class="">uDJF_NAO_lo = uDJF((/15,16,17,21,29,31,48,49,62,63/),:,:)<br class=""><br class="">vDJF_NAO_hi = vDJF((/25,27,35,41,47,52,59,64/),:,:)<br class="">vDJF_NAO_lo = vDJF((/15,16,17,21,29,31,48,49,62,63/),:,:)<br class=""><br class="">zDJF_NAO_hi = zDJF((/25,27,35,41,47,52,59,64/),:,:)<br class="">zDJF_NAO_lo = zDJF((/15,16,17,21,29,31,48,49,62,63/),:,:)<br class=""><br class=""><br class="">uAvgTime_hi = dim_avg_n_Wrap(uDJF_NAO_hi,0)<br class="">uAvgTime_lo = dim_avg_n_Wrap(uDJF_NAO_lo,0)<br class=""><br class="">printVarSummary(uAvgTime_hi)<br class="">printVarSummary(uAvgTime_lo)<br class=""><br class="">vAvgTime_hi = dim_avg_n_Wrap(vDJF_NAO_hi,0)<br class="">vAvgTime_lo = dim_avg_n_Wrap(vDJF_NAO_lo,0)<br class=""><br class="">printVarSummary(vAvgTime_hi)<br class="">printVarSummary(vAvgTime_lo)<br class=""><br class="">zAvgTime_hi = dim_avg_n_Wrap(zDJF_NAO_hi,0)<br class="">zAvgTime_lo = dim_avg_n_Wrap(zDJF_NAO_lo,0)<br class=""><br class="">printVarSummary(zAvgTime_hi)<br class="">printVarSummary(zAvgTime_lo)<br class=""><br class="">; dirty way to copy metadata over first.<br class="">diff_u = uAvgTime_hi;<br class="">diff_v = vAvgTime_hi;<br class="">diff_z = zAvgTime_hi;<br class=""><br class="">diff_u = uAvgTime_hi - uAvgTime_lo<br class="">diff_v = vAvgTime_hi - vAvgTime_lo<br class="">diff_z = zAvgTime_hi - zAvgTime_lo<br class=""><br class="">printVarSummary(diff_u)<br class="">printVarSummary(diff_v)<br class="">printVarSummary(diff_z)<br class=""><br class="">;************************************************<br class="">; create plot<br class="">;************************************************<br class="">wks = gsn_open_wks("ps","Panel_NAO_z_500") ; open a ps<br class="">file<br class="">plot = new(3,graphic) ; create a plot array<br class=""><br class="">;---- set common resources for all plots<br class="">res = True<br class="">res@gsnDraw = False ; dont draw<br class="">res@gsnFrame = False ; dont advance frame<br class="">res@cnInfoLabelOn = False ; trn off cn info<br class="">label<br class="">res = True ; plot mods desired<br class="">res@gsnAddCyclic = False ; has to do with<br class="">wrapping<br class="">the longitude at 0/360<br class="">;************************************************<br class="">; Choose a subregion<br class="">;************************************************<br class="">res@mpMaxLatF = 90 ;maximum latitude<br class="">res@mpMinLatF = 45 ;minimum latitude<br class="">res@mpMaxLonF = 357.5 ; ;maximum longitude<br class="">res@mpMinLonF = 270 ;minimum longitude<br class=""><br class="">;***********************************************<br class="">; ----wind vector plot<br class="">;***********************************************<br class="">vcres = True<br class="">vcres@vcRefAnnoOrthogonalPosF = -1.0 ; move ref vector up<br class="">vcres@vcRefMagnitudeF = 10.0 ; define vector ref<br class="">mag<br class="">vcres@vcRefLengthF = 0.045 ; define length of vec<br class="">ref<br class="">vcres@vcGlyphStyle = "CurlyVector" ; turn on curly<br class="">vectors<br class="">vcres@vcMinDistanceF = 0.017<br class="">vcres@mpFillOn = False ; turn off gray fill<br class="">vcres@mpOutlineBoundarySets = "National" ; turn on country<br class="">boundaries<br class="">vcres@mpGeophysicalLineColor = "Navy" ; color of cont.<br class="">outlines<br class="">vcres@mpGeophysicalLineThicknessF = 1.5 ; thickness of<br class="">outlines<br class=""><br class=""><br class="">;vcres@gsnLeftString = "DJF High NAO"<br class="">winds_hi = gsn_csm_vector_map_ce(wks,uAvgTime_hi,vAvgTime_hi,vcres)<br class="">;vcres@gsnLeftString = "DJF Low NAO"<br class="">winds_lo = gsn_csm_vector_map_ce(wks,uAvgTime_lo,vAvgTime_lo,vcres)<br class="">;vcres@gsnLeftString = "Difference of High - Low"<br class="">winds_diff = gsn_csm_vector_map_ce(wks, diff_u, diff_v,vcres)<br class="">;************************************************<br class="">;---- geopotential height filled contour plot<br class="">;***********************************************<br class="">zfres = res<br class="">zfres@cnFillOn = True<br class="">;zfres@cnLevelSelectionMode = "ExplicitLevels<br class="">;zfres@cnLevels = ispan(-20,90,5)<br class="">zfres@lbLabelFontHeightF = 0.015<br class="">zfres@lbOrientation = "Vertical"<br class="">zfres@pmLabelBarOrthogonalPosF = -0.005<br class=""><br class=""><br class="">contour_zf_hi = gsn_csm_contour(wks,zAvgTime_hi,zfres)<br class="">contour_zf_lo = gsn_csm_contour(wks,zAvgTime_lo,zfres)<br class="">contour_zf_diff = gsn_csm_contour(wks,diff_z,zfres)<br class=""><br class="">plot(0) = overlay(contour_zf_hi,winds_hi)<br class="">plot(1) = overlay(contour_zf_lo,winds_lo)<br class="">plot(2) = overlay(contour_diff,winds_diff)<br class=""><br class=""><br class="">;************************************************<br class="">; create panel<br class="">;************************************************<br class="">resP = True ; modify the panel<br class="">plot<br class="">resP@txString = "NAO 500mb"<br class="">gsn_panel(wks,plot,(/3,1/),resP) ; now draw as one<br class="">plot;<br class=""><br class=""><br class="">end<br class=""><br class=""><br class="">_______________________________________________<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=""><br class=""></blockquote><br class=""></blockquote><Panel_NAO_z_500.pdf></blockquote></blockquote></div></blockquote></div><br class=""></div></body></html>