<div dir="ltr"><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">You said you are getting a bad plot, but you need to provide more information than this. Please attach the PNG images so we can see what they look like. Even better, can you provide me with your data files? You can use our ftp, if the files are not too large:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default"><a href="http://www.ncl.ucar.edu/ftp_files.shtml">http://www.ncl.ucar.edu/ftp_files.shtml</a><br></div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I'm not sure exactly what you are trying to do, since you are creating two completely different plots, both with vectors:</div><div class="gmail_default" style="font-size:small"><blockquote style="margin:0px 0px 0px 40px;border:none;padding:0px"><div class="gmail_default" style="font-size:small"><br></div>plot = gsn_csm_pres_hgt_vector(wks,U(0,:,:,{170}),W(0,:,:,{170}),wscale(0,:,:,{170}),res ) <br>plot_ov = gsn_csm_vector_map(wks,ud(0,1,:,:),vd(0,1,:,:),res)</blockquote></div><br><br>The gsn_csm_pres_hgt_vector function is what draws the Hovmueller vector plot. Since you didn't set gsnDraw and gsnFrame to False for this plot, the first PNG image should be this plot. Does it look correct at all?<br><br>Then with your second plot, you are plotting vectors over a map. This should be the second PNG image. Does this look correct?<div><br></div><div><div class="gmail_default" style="font-size:small">In order for the gsn_csm_pres_hgt_vector plot function to work, the first data array must be the one you want to contour, and the second two data arrays must be the U/V vector arrays. Did you double-check that this is what you have?</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Please read the documentation at:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default"><a href="http://www.ncl.ucar.edu/Document/Graphics/Interfaces/gsn_csm_pres_hgt_vector.shtml">http://www.ncl.ucar.edu/Document/Graphics/Interfaces/gsn_csm_pres_hgt_vector.shtml</a><br></div><div class="gmail_default"><br></div><div class="gmail_default">to make sure you are calling this function correctly. You can also see a couple of examples at:</div><div class="gmail_default"><br></div><div class="gmail_default">h_lat_7.ncl</div><div class="gmail_default"><a href="http://www.ncl.ucar.edu/Applications/height_lat.shtml#ex7">http://www.ncl.ucar.edu/Applications/height_lat.shtml#ex7</a><br></div><div class="gmail_default"><br></div><div class="gmail_default">vector_5.ncl</div><div class="gmail_default"><a href="http://www.ncl.ucar.edu/Applications/vector.shtml#ex5">http://www.ncl.ucar.edu/Applications/vector.shtml#ex5</a><br></div><div class="gmail_default"><br></div><div class="gmail_default">--Mary</div><div class="gmail_default"><br></div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Jan 10, 2018 at 11:30 PM, Ankita Sharma <span dir="ltr"><<a href="mailto:ankita@iiserb.ac.in" target="_blank">ankita@iiserb.ac.in</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_default" style="font-family:comic sans ms,sans-serif;color:#073763">Dear Mary,</div><div class="gmail_default" style="font-family:comic sans ms,sans-serif;color:#073763"><br></div><div class="gmail_default" style="font-family:comic sans ms,sans-serif;color:#073763">Thank you for your time and reply.</div><div class="gmail_default" style="font-family:comic sans ms,sans-serif;color:#073763">As I have said earlier, I was trying to plot hovmoller vector (pressure/height vs longitude) since last several weeks but still not able to do that. I just want to ask shall I give up NCL and plot the diagram for the same in GRADS or ArcGIS?</div><div class="gmail_default" style="font-family:comic sans ms,sans-serif;color:#073763"><br></div><div class="gmail_default" style="font-family:comic sans ms,sans-serif;color:#073763">Thanking you</div><div><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Jan 10, 2018 at 8:00 PM, Mary Haley <span dir="ltr"><<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div style="font-size:small"><br>As I said before, you are trying to overlay a pressure/latitude plot on a lat/lon plot.</div><div style="font-size:small"><br></div><div style="font-size:small">You cannot call overlay on two plots that have different coordinate systems. The Hovmueller plot is a plot that is pressure versus latitude, at a single specified longitude. You are then trying to overlay this on a map plot that is latitude versus longitude. You simply cannot do this, as NCL has no idea how to transform data in pressure / latitude space to latitude / longitude space.</div><div style="font-size:small"><br></div><div style="font-size:small">Do you really need to overlay the Hovmueller plot on a map? Why not just draw the Hovmueller plot as-is?</div><span class="m_-5722018072170113727HOEnZb"><font color="#888888"><div style="font-size:small"><br></div><div style="font-size:small">--Mary</div><div style="font-size:small"><br></div></font></span></div><div class="m_-5722018072170113727HOEnZb"><div class="m_-5722018072170113727h5"><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Jan 9, 2018 at 4:01 AM, Ankita Sharma <span dir="ltr"><<a href="mailto:ankita@iiserb.ac.in" target="_blank">ankita@iiserb.ac.in</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div style="font-family:comic sans ms,sans-serif;color:#073763">I am getting a bad plot, and I am not able to plot the vector winds in Hovmoller plot.</div></div><div class="m_-5722018072170113727m_4959151146809643928HOEnZb"><div class="m_-5722018072170113727m_4959151146809643928h5"><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Jan 8, 2018 at 11:29 PM, Mary Haley <span dir="ltr"><<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div style="font-size:small">Dear Ankita,</div><div style="font-size:small"><br></div><div style="font-size:small">In general, it is hard for us to simply look at a script and determine what the problem is without information about what is wrong. That is, are you getting an error? Are you getting a bad plot?</div><div style="font-size:small"><br></div><div style="font-size:small">However, I see you are trying to overlay a pressure/hgt plot, which is pressure versus latitude, and overlay this on a map (lat/lon) plot, which doesn't make sense. These are two different coordinate systems and it simply won't work.</div><div style="font-size:small"><br></div><div style="font-size:small">Please explain in more detail about what you are trying to plot. It would also help if you can provide the data.</div><div style="font-size:small"><br></div><div style="font-size:small">Thanks,</div><div style="font-size:small"><br></div><div style="font-size:small">--Mary</div><div style="font-size:small"><br></div><div style="font-size:small"><br></div><div style="font-size:small"><br></div></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="m_-5722018072170113727m_4959151146809643928m_-690408185290703186h5">On Sat, Jan 6, 2018 at 12:11 AM, Ankita Sharma via ncl-talk <span dir="ltr"><<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="m_-5722018072170113727m_4959151146809643928m_-690408185290703186h5"><div dir="ltr"><div style="font-family:"comic sans ms",sans-serif"><font color="#000000">Hello,</font></div><div style="font-family:"comic sans ms",sans-serif"><font color="#000000"><br></font></div><div style="font-family:"comic sans ms",sans-serif"><font color="#000000"> I have seen scripts regarding my issues and also tried to correct my scripts using concepts of that scripts. I am new to NCL and have tried every possible measure to plot Hovmueller diagram but still, I'm stuck into it.</font><span style="color:rgb(7,55,99)"></span></div><div style="font-family:"comic sans ms",sans-serif"><font color="#000000">It's my humble request to NCL committee that please have a look at my script and help me with my issues in plotting vector Hovmueller diagram.</font></div><div style="font-family:"comic sans ms",sans-serif"><font color="#000000">I am sending you the script upon which I have been working from last one month.</font></div><div style="font-family:"comic sans ms",sans-serif"><font color="#000000"><br></font></div><div style="font-family:"comic sans ms",sans-serif"><span style="color:rgb(0,0,255)"><br></span></div><div style="font-family:"comic sans ms",sans-serif"><span style="color:rgb(0,0,255)">begin</span><br></div><div><font face="comic sans ms, sans-serif" color="#0000ff"><div><br></div><div>;---File handling</div><div> <span style="white-space:pre-wrap"> </span></div><div> <span style="white-space:pre-wrap"> </span>fn = "<a href="http://omega.mon.mean.nc" target="_blank">omega.mon.mean.nc</a>" ; define filename</div><div> <span style="white-space:pre-wrap"> </span>omega = addfile(fn,"r")<span style="white-space:pre-wrap"> </span> <span style="white-space:pre-wrap"> </span>; open netcdf file</div><div><span style="white-space:pre-wrap"> </span>fn = "<a href="http://uwnd.mon.mean.nc" target="_blank">uwnd.mon.mean.nc</a>" ; define filename</div><div> <span style="white-space:pre-wrap"> </span>uwnd = addfile(fn,"r") <span style="white-space:pre-wrap"> </span>; open netcdf file</div><div><br></div><div>;---Read needed variables from file</div><div> </div><div> <span style="white-space:pre-wrap"> </span>U = uwnd->uwnd</div><div> <span style="white-space:pre-wrap"> </span>W = omega->omega</div><div> <span style="white-space:pre-wrap"> </span>P0mb = 1000.</div><div><span style="white-space:pre-wrap"> </span>hyam = uwnd->uwnd <span style="white-space:pre-wrap"> </span>; get a coefficiants</div><div> <span style="white-space:pre-wrap"> </span>hybm = omega->omega <span style="white-space:pre-wrap"> </span>; get b coefficiants</div><div><br></div><div>;---Define other arguments required</div><div> interp = 2 </div><div> pnew = (/ 1000,900,850,800,750,700,550,5<wbr>00,250,200,150,100/)</div><div> pnew@units = "mb" </div><div>;</div><div>; Omega is significantly smaller than u, so we will scale it so that some vertical motion is visible</div><div>;</div><div> <span style="white-space:pre-wrap"> </span></div><div><span style="white-space:pre-wrap"> </span>wAve = avg(W(0,:,:,{170})) <span style="white-space:pre-wrap"> </span>; used for scaling</div><div> <span style="white-space:pre-wrap"> </span>uAve = avg(U(0,:,:,{170}))</div><div> <span style="white-space:pre-wrap"> </span>scale = fabs(uAve/wAve)</div><div> <span style="white-space:pre-wrap"> </span>wscale = W*scale <span style="white-space:pre-wrap"> </span>; now scale</div><div><br></div><div> <span style="white-space:pre-wrap"> </span>copy_VarCoords(W, wscale) <span style="white-space:pre-wrap"> </span>; copy coordinate variables</div><div><br></div><div><br></div><div>;*****************************<wbr>********************</div><div>; calculate divergence: Use Wrap to include meta data</div><div>;*****************************<wbr>********************</div><div><br></div><div><span style="white-space:pre-wrap"> </span>f = addfile ("<a href="http://uwnd.mon.mean.nc" target="_blank">uwnd.mon.mean.nc</a>", "r")</div><div> <span style="white-space:pre-wrap"> </span>u = f->uwnd</div><div><span style="white-space:pre-wrap"> </span>f = addfile ("<a href="http://vwnd.mon.mean.nc" target="_blank">vwnd.mon.mean.nc</a>", "r")</div><div> <span style="white-space:pre-wrap"> </span>v = f->vwnd</div><div><br></div><div> <span style="white-space:pre-wrap"> </span>div = uv2dvG_Wrap(u,v) ; u,v ==> divergence</div><div><br></div><div>;*****************************<wbr>********************</div><div>; calculate divergent wind components </div><div>;*****************************<wbr>******************** </div><div> <span style="white-space:pre-wrap"> </span>ud = new ( dimsizes(u), typeof(u), u@_FillValue )</div><div><span style="white-space:pre-wrap"> </span>vd = new ( dimsizes(v), typeof(v), v@_FillValue )</div><div><br></div><div> <span style="white-space:pre-wrap"> </span>dv2uvg(div,ud,vd) ; div ==> divergent wind components</div><div><br></div><div> <span style="white-space:pre-wrap"> </span>copy_VarCoords(u, ud ) </div><div> <span style="white-space:pre-wrap"> </span>copy_VarCoords(u, vd ) </div><div> <span style="white-space:pre-wrap"> </span>ud@long_name = "Zonal Divergent Wind"</div><div> <span style="white-space:pre-wrap"> </span>ud@units = u@units</div><div> <span style="white-space:pre-wrap"> </span>vd@long_name = "Meridional Divergent Wind"</div><div> <span style="white-space:pre-wrap"> </span>vd@units = v@units </div><div><br></div><div>;---Create plot</div><div> </div><div><span style="white-space:pre-wrap"> </span>wks = gsn_open_wks ("png", "Vector" ) <span style="white-space:pre-wrap"> </span>; send graphics to PNG file</div><div><br></div><div> <span style="white-space:pre-wrap"> </span>res = True <span style="white-space:pre-wrap"> </span>; plot mods desired</div><div> <span style="white-space:pre-wrap"> </span>res@tiMainString = "Vector Hovmoller" <span style="white-space:pre-wrap"> </span>; title</div><div><br></div><div> <span style="white-space:pre-wrap"> </span>res@cnLineLabelsOn = True <span style="white-space:pre-wrap"> </span>; turn off line labels</div><div> <span style="white-space:pre-wrap"> </span>res@cnFillOn = True <span style="white-space:pre-wrap"> </span>; turn on color fill</div><div> <span style="white-space:pre-wrap"> </span>res@cnFillPalette = "GMT_polar" <span style="white-space:pre-wrap"> </span>; choose color map</div><div><br></div><div> <span style="white-space:pre-wrap"> </span>res@lbLabelStride = 2 <span style="white-space:pre-wrap"> </span>; every other color</div><div><br></div><div> <span style="white-space:pre-wrap"> </span>res@vcRefMagnitudeF = 3.0 <span style="white-space:pre-wrap"> </span>; define vector ref mag</div><div> <span style="white-space:pre-wrap"> </span>res@vcRefLengthF = 0.045 <span style="white-space:pre-wrap"> </span>; define length of vec ref</div><div> <span style="white-space:pre-wrap"> </span>res@vcGlyphStyle = "CurlyVector" <span style="white-space:pre-wrap"> </span>; turn on curly vectors</div><div> <span style="white-space:pre-wrap"> </span>res@vcMinDistanceF = 0.01 <span style="white-space:pre-wrap"> </span>; thin out vectors</div><div> <span style="white-space:pre-wrap"> </span>res@vcMapDirection = True</div><div><br></div><div><span style="white-space:pre-wrap"> </span>res = True</div><div> <span style="white-space:pre-wrap"> </span>res@vcRefMagnitudeF = 4. ; make vectors larger</div><div> <span style="white-space:pre-wrap"> </span>res@vcRefLengthF = 0.050 ; reference vector length</div><div> <span style="white-space:pre-wrap"> </span>res@vcGlyphStyle = "CurlyVector" ; turn on curly vectors</div><div> <span style="white-space:pre-wrap"> </span>res@vcMinDistanceF = 0.012 ; thin the vectors</div><div> </div><div><br></div><div><br></div><div>;---Draw plot from pole to pole at 180E</div><div> </div><div><span style="white-space:pre-wrap"> </span>plot = gsn_csm_pres_hgt_vector(wks,U(<wbr>0,:,:,{170}),W(0,:,:,{170}),ws<wbr>cale(0,:,:,{170}),res ) </div><div><span style="white-space:pre-wrap"> </span>plot_ov = gsn_csm_vector_map(wks,ud(0,1,<wbr>:,:),vd(0,1,:,:),res)</div><div><br></div><div><span style="white-space:pre-wrap"> </span>overlay(plot,plot_ov) ; overlay the U-wind plot on the hovmoller plot</div><div> <span style="white-space:pre-wrap"> </span>draw(plot) ; draw the temperature plot (with the U-wind plot overlaid)</div><div> <span style="white-space:pre-wrap"> </span>frame(wks) ; advance the frame</div><div><br></div><div>end</div><div><br></div><div><br></div><div><div style="color:rgb(34,34,34)"><font color="#000000">Thank you,</font></div><div style="color:rgb(34,34,34)"><font color="#000000">Ankita Sharma</font></div></div></font></div></div>
<br></div></div>______________________________<wbr>_________________<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/mailma<wbr>n/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>
</blockquote></div><br></div>
</div></div></blockquote></div><br></div>
</div></div></blockquote></div><br></div></div></div></div>
</blockquote></div><br></div>