load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl" ;---Generate dummy u, v, w arrays N = 20 PI = 3.14159 I = ispan(0,N-1,1) u = 10. * cos(conform_dims((/N,N/),(2.0 * PI / N) * I, 0)) v = 10. * cos(conform_dims((/N,N/),(2.0 * PI / N) * I, 1)) w = u^2+v^2 wks = gsn_open_wks("x11","overlay1d") vcres = True ; Vector resources vcres@vcLineArrowThicknessF = 3.0 ; 3x as thick vcres@tiMainString = "Vector plot" vector = gsn_csm_vector(wks,u,v,vcres) cnres = True ; Contour resources cnres@cnFillOn = True ; Turn on color fill cnres@lbOrientation = "Vertical" ; Make labelbar vertical cnres@tiMainString = "Contour plot" contour = gsn_csm_contour(wks,w,cnres) overlay(contour,vector) ; Overlay vector plot on contour plot draw(contour) ; This will draw both plots frame(wks)