<div dir="ltr"><div>I hesitate to answer this because I do not want to end up doing the programming. <b><br></b></div><div><b>You must do the work.</b></div><div><b><br></b></div><div>To get you started:<b><br></b></div><div><br></div><div>  u(nens,time,lev,lat,lon), v(nens,time,lev,lat,lon)</div><div><br></div><div>  pHigh   = 250              <br></div><div>  pLow    = 850<br></div><div>  ushear = u(:,:,{pHigh,:,:)-u(:,:,{pLow},:,:)</div><div>  vshear = v(:,:,{pHigh},:,:)-v(:,:,{pLow},:,:)</div><div>  ushear@long_name = "USHEAR: ("+pHigh+"-"+pLow+")"<br></div><div>  vshear@long_name = "VSHEAR: ("+pHigh+"-"+pLow+")"</div><div>  ushear@units           = u@units</div><div>  vshear@units           = v@units</div><div>  copy_VarCoords(u(:,:,0,:,:), ushear)</div><div>  copy_VarCoords(v(:,:,0,:,:), vshear)</div><div><br></div><div>  printVarSummary(ushear)                          ; (nens,time,lat,lon)<br></div><div>  printVarSummary(vshear)</div><div><br></div><div>  shear = sqrt(ushear^2 + vshear^2)<br>  shear@long_name = "SHEAR" <br>  shear@units     = ushear@units<br>  copy_varCoords(ushear, shear)<br>  printVarSummary(shear) <br></div><div><br></div><div>----</div><div>Perhaps:  <br></div><div><br></div><div>[1] clat(nens,ntim) and clon(nens,ntim)  ; storm center for each ensemble member at each time.</div><div>[2] Use <a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/gc_latlon.shtml" target="_blank"><b>gc_latlon</b></a> to calculate the distance from the storm center(s)</div><div>[3] Calculate areal means according to distance. <br></div><div><br></div><div>----<br></div><div>Another example that might be useful:<br></div><div><a href="http://www.ncl.ucar.edu/Applications/regrid.shtml" target="_blank"><b>http://www.ncl.ucar.edu/Applications/regrid.shtml</b></a></div><div>See Example: regrid_15.ncl</div><div><br></div><div>Download: <a href="http://www.ncl.ucar.edu/Applications/Scripts/grid2geocircle.ncl" target="_blank"><b>http://www.ncl.ucar.edu/Applications/Scripts/grid2geocircle.ncl</b></a></div><div><b>*Carefully* </b>read the documentation <br></div><div>
</div><div><pre><br></pre></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Nov 29, 2019 at 9:44 AM Goriparthi Pavani via ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><div style="font-family:arial,helvetica,sans-serif;font-size:12pt;color:rgb(0,0,153)"><div>Hi,</div><div>I want to plot magnitude of shear for each ensemble member by taking average of the difference in wind velocities between the 850 hPa and 250hPa levels at every grid point that is between 200 and 500km from the surface center for a tropical cyclone.</div><div>Domain used is moving nest with 3km resolution.<br></div><div><br></div><div>Could you please help me with the ncl script to do the area average of wind shear magnitude between 200 and 500km from surface center.</div><div><br></div><div>Regards</div><div>Pavani<br></div></div></div>_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a></blockquote></div>