<html><body><div style="color:#000; background-color:#fff; font-family:HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, Sans-Serif;font-size:16px"><div id="yui_3_16_0_1_1432615248546_17553"><span>Hi</span></div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr"><span id="yui_3_16_0_1_1432615248546_17582">I have to plot 850hpa Vector wind for a regional domain (0-40N/50-100E). My dataset is FNL reanalysis. I tried following scripts. (png file is attached)</span></div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style="">;</div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style="">load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"</div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style="">load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"</div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style="">load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"</div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style="">; Open the file----------</div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style="">begin</div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style="">a = addfile("~/Mesomodel/gfs/test.nc","r")</div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style="">;printVarSummary(a)</div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style="">;</div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style="">; -----Read data------</div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style="">Uwind = a->UGRD_P0_L100_GLL0(20,:,:); Read the variable to memory</div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style="">Vwind = a->VGRD_P0_L100_GLL0(20,:,:); Read the variable to memory</div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style="">Uwind@_FillValue = -999</div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style="">Vwind@_FillValue = -999</div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style="">printVarSummary(Uwind)</div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style="">printVarSummary(Vwind)</div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style="">;</div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style="">; ----Create Plot-----------------</div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style="">wks = gsn_open_wks("png", "wind850")</div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style="">res = True</div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style="">res@tiMainString = "WIND 850hPa "</div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style="">;pltres = True</div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style="">;</div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style="">; Map Resources</div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style="">mpres = True</div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style="">res@mpFillOn =False</div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style="">res@mpOutlineBoundarySets ="National"</div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style="">res@mpGeophysicalLineColor ="Black"</div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style="">res@mpMinLatF = 0.0</div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style="">res@mpMaxLatF = 40.0</div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style="">res@mpMinLonF = 50.0</div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style="">res@mpMaxLonF = 100.0</div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style="">res@mpDataBaseVersion ="MediumRes"</div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style="">;plot = gsn_csm_map_ce(wks,res) ;Draws a map</div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style="">;plot = gsn_csm_vector_map(wks,Uwind,Vwind,True)</div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style="">;plot = gsn_csm_map_ce(wks,True)</div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style="">;</div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style=""> </div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr"></div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style="">;vecres =True</div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style="">res@vcMinFracLengthF =1.0 ; Length of min vector</div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style="">res@vcRefMagnitudeF =3.0 ; Makes Vectors larger</div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style="">res@vcRefLengthF =0.045 ; ref vec Lengthlarger</div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style="">res@vcGlyphStyle ="CurlyVector" ; ref vec Lengthlarger</div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style="">res@vcMinDistanceF =0.01 ; ref vec Lengthlarger</div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style="">res@vcMonoLineArrowColor =False ; ref vec Lengthlarger</div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style=""> </div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style="">plot = gsn_csm_vector_map(wks,Uwind,Vwind,True)</div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style="">end</div><div id="yui_3_16_0_1_1432615248546_17553" dir="ltr" class="" style=""> </div><div dir="ltr" class="" style="" id="yui_3_16_0_1_1432615248546_17703">With this script, the problem is it is NOT Plotting the data for the regional domain. </div><div></div><div id="yui_3_16_0_1_1432615248546_17522" dir="ltr"><br></div><div id="yui_3_16_0_1_1432615248546_17522" dir="ltr"><br></div><div id="yui_3_16_0_1_1432615248546_17522" dir="ltr">A little change in the above script (changes given below)</div><div id="yui_3_16_0_1_1432615248546_17522" dir="ltr" class="" style="">plot = gsn_csm_map_ce(wks,res) ;Draws a map</div><div id="yui_3_16_0_1_1432615248546_17522" dir="ltr"><br></div><div id="yui_3_16_0_1_1432615248546_17522" dir="ltr">draws a regional map BUT without data. </div><div id="yui_3_16_0_1_1432615248546_17522" dir="ltr"><br></div><div id="yui_3_16_0_1_1432615248546_17522" dir="ltr">Both figures attatched. </div><div id="yui_3_16_0_1_1432615248546_17522" dir="ltr"><br></div><div id="yui_3_16_0_1_1432615248546_17522" dir="ltr">Any suggestions............ </div><div id="yui_3_16_0_1_1432615248546_17521"><div id="yui_3_16_0_1_1432615248546_17520">Geeta. </div></div></div></body></html>