<div dir="ltr"><div class="gmail_default" style="font-family:times new roman,serif">Thank you Adam Phillips. I am able to plot in terms of (/5.e-5,1.e-4,1.e-3,1.e-2,.1,1,10/). But I am unable to normalize (P/Pmax). <br><br></div><div class="gmail_default" style="font-family:times new roman,serif">Regards<br></div><div class="gmail_default" style="font-family:times new roman,serif">Subrahmanyam Sappati<br></div><div class="gmail_default" style="font-family:times new roman,serif"><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Aug 21, 2015 at 1:25 AM, Adam Phillips <span dir="ltr"><<a href="mailto:asphilli@ucar.edu" target="_blank">asphilli@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><div><div><div><div><div>Hi Subrahmanyam,<br></div>Yes, I think you just need to set your contour levels appropriately.<br><br></div>Try using these resources:<br>res@cnLevelSelectionMode = "ExplicitLevels"<br></div>res@cnLevels = (/5.e-5,1.e-4,1.e-3,1.e-2,.1,1,10/)<br></div>; You can add levels to the above line as needed<br><br></div>That should do it. If that doesn't help or if you have further issues please respond to the ncl-talk email list..<br></div>Adam<br><div><div><br><div><br><br><br><br></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="h5">On Thu, Aug 20, 2015 at 12:55 AM, Subrahmanyam Sappati <span dir="ltr"><<a href="mailto:subrahmanyamsappati@gmail.com" target="_blank">subrahmanyamsappati@gmail.com</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="h5"><div dir="ltr"><div class="gmail_default" style="font-family:times new roman,serif">Dear All,<br><br>I am able to plot Joint PDF of two sets (data set 1, data set 2) with following ncl code (Refer fig - Fig-10.png). Is it possible to plot normalized Joint PDF in logarithmic scale (10^(n)) ? (please refer Ref.png)<br></div><br><div class="gmail_default" style="font-family:times new roman,serif">Here is my ncl code.<br><br>;-----------------------------------------------------------------------------<br>; Plot Phase Space Diagram<br>;----------------------------------------------------------------------------<br>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl" <br>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl" <br>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl" <br>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/diagnostics_cam.ncl" <br><br>begin<br><br>dir="./"<br><br>;o1=asciiread(dir+"2nd-NELP-Ow199.dat",(/103389/),"float")<br>;o2=asciiread(dir+"2nd-Diff-N-ELP-Hw-Ow.dat",(/103389/),"float")<br>o2=asciiread(dir+"cos-angle.dat",(/103389/),"float")<br>o1=asciiread(dir+"PTC-N31-Hw-Ow.dat",(/103389/),"float")<br> printVarSummary(o1)<br><br><br> opt = True<br> opt@bin_nice = True<br><br> pdf2obs = pdfxy(o1,o2, 0, 0, opt) ; 25 is the default <br>printVarSummary(pdf2obs)<br>; print(pdf2obs)<br><br>;************************************************<br>; create plots<br>;************************************************<br>; wks = gsn_open_wks("pdf","FIG-10") ; open a ps file<br> wks = gsn_open_wks("png","FIG-10") ; open a ps file<br>; gsn_define_colormap(wks,"prcp_2")<br> gsn_define_colormap(wks,"amwg256")<br><br><br> res = True<br> res@gsnDraw = True <br> res@gsnFrame = True <br> res@cnInfoLabelOn = False<br> res@cnFillOn = True ; Turn on color<br> res@cnLinesOn = False <br> res@cnLineLabelsOn = False <br> res@lbLabelAutoStride = True <br> res@gsnSpreadColors = True<br> res@lbOrientation = "vertical" ; vertical label bar<br><br> res@trYMinF = 0.6 ; set minimum Y-axis value<br> res@trYMaxF = 1.0 ; set maximum Y-axis value<br> res@trXMinF = -0.8 ; set minimum X-axis value<br> res@trXMaxF = 1.8 ; set maximum X-axis value<br><br><br> res@tmXMajorGrid = True ; implement x grid<br> res@tmXMajorGridThicknessF = 0.5 ; 2.0 is default<br> res@tmXMajorGridLineDashPattern = 2 ; select short dash lines<br><br> res@tmYMajorGrid = True ; implement y grid<br> res@tmYMajorGridThicknessF = 0.5 ; 2.0 is default<br> res@tmYMajorGridLineDashPattern = 2 ; select short dash lines<br><br> res@cnLevelSelectionMode = "ManualLevels" ; manually set the contour levels with the following 3 resources<br> res@cnMinLevelValF = 0.001 ; set the minimum contour level<br> res@cnMaxLevelValF = 6.0 ; set the maximum contour leve<br> res@cnLevelSpacingF = 0.04 ; set the interval between contours<br><br> res@tiXAxisString = "bond Length (A)"<br> res@tiYAxisString = "tragectory"<br><br><br> plot = gsn_csm_contour (wks,pdf2obs, res)<br><br><br>end<span><font color="#888888"><br><br></font></span></div><span><font color="#888888"><br>-- <br><div><div dir="ltr"><br><span style="font-family:times new roman,serif"><span style="color:rgb(0,0,153)">Subrahmanyam.Sappati</span><br><span style="color:rgb(102,0,204)">C/o Dr. Prasenjit Ghosh</span><br><span style="color:rgb(102,0,204)">IISER PUNE<br>Pin code : 411 008</span><br><br><span style="color:rgb(204,0,0)">Phone : <a href="tel:%2B91%2094%2003%2059%2035%2018" value="+919403593518" target="_blank">+91 94 03 59 35 18</a></span></span><br><br></div></div>
</font></span></div>
<br></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><br>
<br></blockquote></div><span class="HOEnZb"><font color="#888888"><br><br clear="all"><br>-- <br><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div><div><span><font color="#888888">Adam Phillips <br></font></span></div><span><font color="#888888">Associate Scientist, </font></span><span><font color="#888888">Climate and Global Dynamics Laboratory, NCAR<br></font></span></div></div><div><span><font color="#888888"><a href="http://www.cgd.ucar.edu/staff/asphilli/" target="_blank">www.cgd.ucar.edu/staff/asphilli/</a> </font></span><span><font color="#888888">303-497-1726 </font></span></div><span><font color="#888888"></font></span><div><div><span><font color="#888888"><br></font></span><div><span><font color="#888888"><a href="http://www.cgd.ucar.edu/staff/asphilli" target="_blank"></a></font></span></div></div></div></div></div></div></div></div></div></div></div>
</font></span></div>
</blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature"><div dir="ltr"><br><span style="font-family:times new roman,serif"><span style="color:rgb(0,0,153)">Subrahmanyam.Sappati</span><br><span style="color:rgb(102,0,204)">C/o Dr. Prasenjit Ghosh</span><br><span style="color:rgb(102,0,204)">IISER PUNE<br>Pin code : 411 008</span><br><br><span style="color:rgb(204,0,0)">Phone : +91 94 03 59 35 18</span></span><br><br></div></div>
</div>