;************************************************* ; vector_3.ncl ;************************************************ load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl" ;************************************************ begin ;************************************************ ; read in netCDF file ;************************************************ a = addfile("../uwnd.mon.mean.nc","r") b = addfile("../vwnd.mon.mean.nc","r") ;************************************************ ; read in zonal [u] and meridional [v] winds (July) ;************************************************ do t = 0,2 do z = 1,3 if(z.eq.1) then nivel="1000 hPa" end if if(z.eq.2) then nivel="925 hPa" end if if(z.eq.3) then nivel="850 hPa" end if u = short2flt(a->uwnd(t,z,:,:)) v = short2flt(b->vwnd(t,z,:,:)) ;************************************************ ; create plot ;************************************************ wks = gsn_open_wks("ps","vetor") ; open a ps file vcres = True ; plot mods desired vcres@tiMainString = "Vetor Vento" vcres@gsnLeftString = t vcres@gsnCenterString = nivel vcres@gsnRightString = "m/s" vcres@pmTickMarkDisplayMode = "Always" ; Turn on map tickmarks vcres@mpDataSetName = "Earth..4" vcres@mpDataBaseVersion = "MediumRes" vcres@mpOutlineOn = True vcres@mpOutlineSpecifiers = (/"Brazil:states"/) vcres@mpFillOn = True ; Turn on map fill vcres@mpFillBoundarySets = "National" vcres@mpFillAreaSpecifiers = (/"Brazil:states","China:states","India:states"/) vcres@gsnMaximize = True ; Ajusta o grafico a janela vcres@gsnMajorLonSpacing = 10 ; escala do eixo x, variacao vcres@gsnMajorLatSpacing = 10 ; escala do eixo x, variacao vcres@gsnFrame = False ; so we can draw time stamp vcres@vcRefAnnoOrthogonalPosF = -1.0 ; move ref vector up vcres@vcRefMagnitudeF = 10.0 ; define vector ref mag vcres@vcRefLengthF = 0.045 ; define length of vec ref vcres@vcGlyphStyle = "CurlyVector" ; turn on curly vectors vcres@vcMinDistanceF = 0.017 vcres@mpMinLonF = -100. ; longitude oeste vcres@mpMaxLonF = -20. ; longitude leste vcres@mpMinLatF = -60. ; latitude sul vcres@mpMaxLatF = 20. ; latitude norte vector = gsn_csm_vector_map_ce(wks,u,v,vcres) frame (wks) end do end do end