<div dir="ltr"><div><div><div><div><div><div><div>The contour plot routines expect (y_axis,x_axis) order<br><br></div><div>  print(RF_fnl&amp;nlev)    ; look ... does this have a long_name or units attribute?<br><br></div><div><br></div><div>  work = RF_fnl(nlev|:,nscan|:,npix|:)  ; convenience for illustration<br></div><div>  printVarSummary(work)                    ;    [nlev | 80] x   [nscan | 9247] x [npix | 49] <br><br></div><div>  np       = ??   ; your choice<br></div><div>  work2 = work(:,:,np)<br></div><div>  printVarSummary(work2)                   ;  [nlev | 80] x   [nscan | 9247] <br><br></div><div>  res@tmYLFormat = &quot;f&quot;                       ; No unnecessary zeroes<br></div><div>  res@tiYAxisString = &quot;level (km)&quot;        ; if not specified, long_name will be used<br></div>  plot = gsn_csm_contour(wks,work2, res)<br><br></div>You will still get index values along the x-axis (longitude axis)<br><br>---<br></div>You can plot the lat/lon at a specified &#39;npix&#39; index usinig lat2d/lon2d<br><br></div>See: <a href="http://www.ncl.ucar.edu/Applications/narr.shtml">http://www.ncl.ucar.edu/Applications/narr.shtml</a><br></div>Example 6<br><br></div>It is not for satellite data but the principle is the same.<br></div><div><br></div><div><br>===<br></div><div>This is *not* a &#39;pure&#39; zonal plot. Fot that you must interpolate to a user specified grid.<br></div><div>Use ESMF regridding.<br><br></div><div>Good luck<br></div><br><div><div><div>        <br><div><div><div><div><div><br></div></div></div></div></div></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Jun 17, 2016 at 5:40 AM, Geeta Geeta <span dir="ltr">&lt;<a href="mailto:geetag54@yahoo.com" target="_blank">geetag54@yahoo.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div style="color:#000;background-color:#fff;font-family:HelveticaNeue,Helvetica Neue,Helvetica,Arial,Lucida Grande,Sans-Serif;font-size:16px"><div style="outline:none 0px;display:table;width:844px;padding-top:12px;padding-left:0px;font-family:&#39;Helvetica Neue&#39;,&#39;Segoe UI&#39;,Helvetica,Arial,&#39;Lucida Grande&#39;,sans-serif;font-size:13px"><div style="display:table-cell;width:auto;word-wrap:break-word;word-break:break-word"><div><div><div style="font-family:HelveticaNeue,&#39;Helvetica Neue&#39;,Helvetica,Arial,&#39;Lucida Grande&#39;,sans-serif;font-size:16px"><div dir="ltr">I have swath data and I am extracting the rain rate from it This is 3 dimensional.</div><div><br></div><div>Variable: RF_fnl</div><div>Type: double</div><div>Total Size: 289985920 bytes</div><div>            36248240 values</div><div>Number of Dimensions: 3</div><div>Dimensions and sizes:   [nscan | 9247] x [npix | 49] x [nlev | 80]</div><div>Coordinates:</div><div>            nlev: [20..0.25]</div><div>Number Of Attributes: 3</div><div>  lat2d :       &lt;ARRAY of 453103 elements&gt;</div><div>  lon2d :       &lt;ARRAY of 453103 elements&gt;</div><div>  _FillValue :  -9999</div><div dir="ltr"><br></div><div dir="ltr">The 3rd dim  has 80 levels that corresponds to data from 20km to 0.25km at 250m resolution. </div><div dir="ltr"><br></div><div dir="ltr">In the vertical cross section (Ht vs Longitude Plot) plot for the region of interest (attatched). </div><div dir="ltr">My problem s how do I convert these nlevels from 0-79 as shown in the X axis to 0.25 to 20 km.</div><div dir="ltr"><br></div><div dir="ltr">also how to change 0-120 on the Y axis to the longitude values</div></div></div></div></div></div><span class="HOEnZb"><font color="#888888"><div style="margin-top:30px;font-family:&#39;Helvetica Neue&#39;,&#39;Segoe UI&#39;,Helvetica,Arial,&#39;Lucida Grande&#39;,sans-serif;font-size:13px"><div><div><div style="width:844px;border-top-width:0px;border-top-color:rgb(236,236,236)"><div style="padding-bottom:0px"><ul><li title="Vert_Plot.png (54.3 KB)" style="background-image:url(https://thumbp10-sg3.mail.yahoo.com/tn?sid=2786232916&amp;mid=AFHFCmoAAHFwV2PetQrF0O3xFaY&amp;midoffset=2_0_0_1_594916&amp;partid=2&amp;f=1901&amp;fid=Inbox&amp;ymreqid=bb465376-2dd3-db24-01e3-6d0141010000&amp;m=ThumbnailService)"><div dir="ltr"><br></div></li></ul></div></div></div></div></div><div></div><div> </div><div>Geeta.</div></font></span></div></div><br>_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">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><br>
<br></blockquote></div><br></div>