<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style id="owaParaStyle" type="text/css">P {margin-top:0;margin-bottom:0;}</style>
</head>
<body ocsi="0" fpstyle="1">
<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">Good day everyone,<br>
<br>
I would like to overlay 12 wind direction vector maps on 12 wind speed maps however I keep getting an error. Can anyone assist?<br>
<br>
Thanking you in advance.<br>
<br>
Regards,<br>
Nkese.<br>
<br>
My script:<br>
<br>
;;create plots<br>
<br>
wks = gsn_open_wks("ps","Caribbeanwinds_monthly6")<br>
plot = new(12,graphic)<br>
plots = new(12,graphic)<br>
<br>
res = True<br>
res@gsnDraw = False ; don't draw<br>
res@gsnFrame = False ; don't advance frame<br>
res@cnFillOn = True ; turn on color<br>
res@gsnMaximize = True <br>
<br>
cmap = read_colormap_file("BkBlAqGrYeOrReViWh200")<br>
res@cnFillPalette = cmap(25:170,:) ; set color map<br>
<br>
res@cnLinesOn = False ; no contour lines<br>
res@lbLabelBarOn = False<br>
<br>
res@gsnLeftString = "Speed" ; change left string<br>
res@gsnRightString = uwnd_practice@units ; assign right string<br>
<br>
res@mpDataSetName = "Earth..4"<br>
res@mpDataBaseVersion = "MediumRes"<br>
res@mpOutlineOn = True<br>
<br>
res@mpGeophysicalLineThicknessF = 2<br>
res@mpNationalLineThicknessF = 2<br>
res@mpFillDrawOrder = "PostDraw"<br>
res@mpFillOn = False ; no map fill<br>
<br>
res@mpMinLatF = 0<br>
res@mpMaxLatF = 30<br>
res@mpMinLonF = -110<br>
res@mpMaxLonF = -30<br>
<br>
vecres = True ; vector only resources<br>
vecres@gsnDraw = False ; don't draw<br>
vecres@gsnFrame = False ; don't advance frame<br>
vecres@vcGlyphStyle = "LineArrow" ; curly vectors<br>
vecres@vcRefMagnitudeF = 10 ; define vector ref mag<br>
vecres@vcRefLengthF = 0.045 ; define length of vec ref<br>
<br>
vecres@vcLineArrowThicknessF = 1<br>
vecres@vcMinDistanceF = 0.05<br>
<br>
vecres@gsnRightString = " " ; turn off right string<br>
vecres@gsnLeftString = " " ; turn off left string<br>
vecres@tiXAxisString = " " ; turn off axis label <br>
vecres@vcRefAnnoOrthogonalPosF = -1.0 ; move ref vector into plot<br>
<br>
months = (/"Jan","Feb","Mar","Apr","May","Jun","Jul","Aug",\<br>
"Sep","Oct","Nov","Dec"/)<br>
<br>
resP = True ; panel only resources<br>
resP@gsnPanelLabelBar = True<br>
resP@gsnMaximize = True ; maximize plots <br>
<br>
<br>
do nmo=0,11<br>
res@gsnCenterString=months(nmo)<br>
vecres@gsnCenterString=months(nmo)<br>
plot(nmo) = gsn_csm_contour_map_ce(wks,wspd_monthly(nmo,:,:),res)<br>
plots(nmo) = gsn_csm_vector_map(wks,uwnd_monthly(nmo,:,:),vwnd_monthly(nmo,:,:),vecres)<br>
overlay(plot(nmo),plots(nmo))<br>
plot = plot(nmo)<br>
end do <br>
<br>
<br>
gsn_panel(wks,plot,(/3,4/),resP)<br>
<br>
<br>
Error:<br>
fatal:NhlAddOverlay: plot class mapPlotClass cannot be overlay plot member<br>
fatal:NhlAddOverlay: plot class mapPlotClass cannot be overlay plot member<br>
fatal:NhlAddOverlay: plot class mapPlotClass cannot be overlay plot member<br>
fatal:NhlAddOverlay: plot class mapPlotClass cannot be overlay plot member<br>
fatal:NhlAddOverlay: plot class mapPlotClass cannot be overlay plot member<br>
fatal:NhlAddOverlay: plot class mapPlotClass cannot be overlay plot member<br>
fatal:NhlAddOverlay: plot class mapPlotClass cannot be overlay plot member<br>
fatal:NhlAddOverlay: plot class mapPlotClass cannot be overlay plot member<br>
fatal:NhlAddOverlay: plot class mapPlotClass cannot be overlay plot member<br>
fatal:NhlAddOverlay: plot class mapPlotClass cannot be overlay plot member<br>
fatal:NhlAddOverlay: plot class mapPlotClass cannot be overlay plot member<br>
fatal:NhlAddOverlay: plot class mapPlotClass cannot be overlay plot member<br>
warning:MapSetTrans: map limits invalid - using maximal area<br>
fatal:MapSetTrans: error initializing map: MAPINT/MDPINT - MAP HAS ZERO AREA<br>
warning:MapSetTrans: map limits invalid - using maximal area<br>
fatal:MapSetTrans: error initializing map: MAPINT/MDPINT - MAP HAS ZERO AREA<br>
warning:MapSetTrans: map limits invalid - using maximal area<br>
fatal:MapSetTrans: error initializing map: MAPINT/MDPINT - MAP HAS ZERO AREA<br>
warning:MapSetTrans: map limits invalid - using maximal area<br>
fatal:MapSetTrans: error initializing map: MAPINT/MDPINT - MAP HAS ZERO AREA<br>
warning:MapSetTrans: map limits invalid - using maximal area<br>
fatal:MapSetTrans: error initializing map: MAPINT/MDPINT - MAP HAS ZERO AREA<br>
warning:MapSetTrans: map limits invalid - using maximal area<br>
fatal:MapSetTrans: error initializing map: MAPINT/MDPINT - MAP HAS ZERO AREA<br>
warning:MapSetTrans: map limits invalid - using maximal area<br>
fatal:MapSetTrans: error initializing map: MAPINT/MDPINT - MAP HAS ZERO AREA<br>
<br>
</div>
CONFIDENTIALITY: This email (including any attachments) may contain confidential, proprietary and/or privileged information. Any duplication, copying, distribution, dissemination, transmission, disclosure or use in any manner of this email (including any attachments)
without the authorisation of the sender is strictly prohibited. If you receive this email (including any attachments) in error, please notify the sender and delete this email (including any attachments) from your system. Thank you.
</body>
</html>