<div dir="ltr"><div class="gmail_default" style="font-size:small">Xianoi,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I've already replied offline about the conda issues, and I'm hoping Kevin can address these.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Meanwhile, with regard to your python script, I think you are on the right track. However, I don't think you want to transpose the data before you plot it. If it is dimensioned latitude x longitude, with latitude being the leftmost dimension, then you don't need to transpose it. If it is dimensioned longitude by latitude, then you need to transpose it.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">To show the figure on the screen (we call it an "X11 window"), then change "wks_type" to be "x11" instead of "ps". You need to click on the X11 window with your left mouse button if you want the plot to go away.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Also, you don't want to set these two options:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small"><div style="font-size:12.8px"> resources.nglDraw = False</div><div style="font-size:12.8px"> resources.nglFrame = False</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">because these settings would cause your plot to not be drawn at all. Either remove those lines or comment them out. They are both True by default.</div></div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">And yes, once you have all the parameters (we call them "resources") set up, you should be able to use the same ones to plot other data. The tiMainString is probably the only thing you'll need to change.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">However, if you decide you don't like the contour levels that PyNGL chooses for you, then you might want to change these for your plots. You can look at the cn05p.py for an example of setting contour levels using a min value, max value, and a spacing (for equally spaced contours), or ctnccl.py for example of explicitly setting contour levels to an unequally-spaced array of values.</div>
<div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Finally, if your data is taking a long time to plot because it is a large array, then you can try setting:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small"><div style="font-size:12.8px"> resources.cnFillMode = "RasterFill" # default is AreaFill</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">Raster contouring can be noticeably faster than the default "AreaFill", which is smooth contouring. You can also set:</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px"> resources.cnRasterSmoothingOn = True</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">which will smooth your raster contours, but still should be faster than AreaFill.</div></div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I will attach a Python script that plots the "topo" variable from old and partial ISCCP file that I have. The data on this file is one-dimensional, along with the lat/lon arrays. Here's what a dump of the file looks like:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default"><div class="gmail_default"><font face="monospace, monospace"> file global attributes:</font></div><div class="gmail_default"><font face="monospace, monospace"> description : ISSCP reduced grid</font></div><div class="gmail_default"><font face="monospace, monospace"> dimensions:</font></div><div class="gmail_default"><font face="monospace, monospace"> nodes = 6596</font></div><div class="gmail_default"><font face="monospace, monospace"> variables:</font></div><div class="gmail_default"><font face="monospace, monospace"> float lon ( nodes )</font></div><div class="gmail_default"><font face="monospace, monospace"> units : degrees east</font></div><div class="gmail_default"><font face="monospace, monospace"> long_name : longitude</font></div><div class="gmail_default"><font face="monospace, monospace"><br></font></div><div class="gmail_default"><font face="monospace, monospace"> float lat ( nodes )</font></div><div class="gmail_default"><font face="monospace, monospace"> units : degrees north</font></div><div class="gmail_default"><font face="monospace, monospace"> long_name : latitude</font></div><div class="gmail_default"><font face="monospace, monospace"><br></font></div><div class="gmail_default"><font face="monospace, monospace"> float topo ( nodes )</font></div><div class="gmail_default"><font face="monospace, monospace"> units : meters above sea level (I think)</font></div><div class="gmail_default"><font face="monospace, monospace"> long_name : topology</font></div><div class="gmail_default"><br></div></div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I'll include the script and the second and third images just in case. I use "png" as the output format.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">--Mary</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Oct 24, 2016 at 6:50 AM, Xiaoni Wang <span dir="ltr"><<a href="mailto:xiaoni.wang@obspm.fr" target="_blank">xiaoni.wang@obspm.fr</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word">Dear colleague,<div> I tried to install Pyngl on my computer in order to draw 2D figures with ISCCP grids. Before this, I had used condo to install python2 and other python tools in /obs/xxwang/anaconda2/.</div><div><br></div><div> 1) To install Pyngl, I followed the suggestion on your website:</div><div> <font color="#831100">conda create --name pyn_test python=2</font></div><div><font color="#831100"> source activate pyn_test</font></div><div><font color="#831100"> conda install --channel ncar --channel khallock pynio pyngl</font></div><div><br></div><div> Then I modified my .bashrc by add the new path : export PATH="/obs/xxw/anaconda2/bin:<font color="#7c2a00">/<wbr>obs/xxw/anaconda2/envs/pyn_<wbr>test/bin</font>:$PATH »</div><div> I can see the environments:</div><div> pyn_test * /obs/xxw/anaconda2/envs/pyn_<wbr>test</div><div> root /obs/xxw/anaconda2</div><div><br></div><div> Then I tried to import Nio, but I have an error: ImportError: No module named Nio</div><div><br></div><div> Could you help me to figure out why ? </div><div><br></div><div> 2) I have another question about using the Pyngl for plotting data on ISCCP grids. Maybe this question would be (too) simple….</div><div><br></div><div> I have my data on ISCCP equal area grids (0.25 degree in latitude, and varying espace in longitude, with more grids near equator and less grids near pole). My data are global data. In my test:</div><div> ISCCP latitude : lat_isccp, 1D, (256915,) ),</div><div> ISCCP longitude: lon_isccp, (1D, (256915,)) </div><div> data: emis (1D , (256915,)). </div><div><br></div><div> I would like to plot emis data on 2D map. I do not know if I should define <span style="white-space:pre-wrap">resources.mpProjection ? </span><span style="white-space:pre-wrap">how to show figure on screen, or h</span><span style="white-space:pre-wrap">ow to save this figure ? or If I have several parameters to plot, I guess that I can use these configurations for all the figures (except </span>resources.tiMainString )<span style="white-space:pre-wrap">? </span></div><div> In the below is my example code for plotting one parameter ‘emis’ . I adapted from your SEAM example.</div><div><br></div><div> wks_type = "ps"</div><div> wks = Ngl.open_wks(wks_type, "isccp")</div><div> resources = Ngl.Resources()</div><div> resources.sfXArray = lon_isccp</div><div> resources.sfYArray = lat_isccp</div><div> resources.nglSpreadColorStart = 176</div><div> resources.nglSpreadColorEnd = 2</div><div> resources.cnFillOn = True</div><div> resources.cnFillMode = "AreaFill"</div><div> resources.cnFillPalette = "BlAqGrYeOrReVi200"</div><div> resources.cnLinesOn = False</div><div> resources.cnLineLabelsOn = False</div><div> resources.lbBoxLinesOn = False</div><div> resources.lbLabelFontHeightF = 0.015</div><div> resources.tiMainString = "Emis"</div><div> resources.nglDraw = False</div><div> resources.nglFrame = False</div><div><br></div><div> zt = np.transpose(emis)</div><div> contour = Ngl.contour(wks,zt,resources)</div><div> <span style="white-space:pre-wrap">Ngl.end()</span></div><div><span style="white-space:pre-wrap"><br></span></div><div><span style="white-space:pre-wrap"> Thank you very much for your help in advance !!</span></div><div><span style="white-space:pre-wrap"><br></span></div><div><span style="white-space:pre-wrap"><br></span></div><div><span style="white-space:pre-wrap">Best regards</span></div><div><span style="white-space:pre-wrap"><br></span></div><div><span style="white-space:pre-wrap">Xiaoni W</span></div><div><br></div><div><br></div><div> </div><div> </div><div> </div></div><br>______________________________<wbr>_________________<br>
pyngl-talk mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/pyngl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/pyngl-talk</a><br>
<br></blockquote></div><br></div>