[ncl-talk] Joint probability distribution in Logarithmic scale

[Adam Phillips]

Hi Subrahmanyam,
Yes, I think you just need to set your contour levels appropriately.

Try using these resources:
res at cnLevelSelectionMode = "ExplicitLevels"
res at cnLevels = (/5.e-5,1.e-4,1.e-3,1.e-2,.1,1,10/)
; You can add levels to the above line as needed

That should do it. If that doesn't help or if you have further issues
please respond to the ncl-talk email list..
Adam






On Thu, Aug 20, 2015 at 12:55 AM, Subrahmanyam Sappati <
subrahmanyamsappati at gmail.com> wrote:

> Dear All,
>
> 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)
>
> Here is my ncl code.
>
>
> ​;-----------------------------------------------------------------------------
> ; Plot Phase Space Diagram
>
> ;----------------------------------------------------------------------------
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/diagnostics_cam.ncl"
>
> begin
>
> dir="./"
>
> ;o1=asciiread(dir+"2nd-NELP-Ow199.dat",(/103389/),"float")
> ;o2=asciiread(dir+"2nd-Diff-N-ELP-Hw-Ow.dat",(/103389/),"float")
> o2=asciiread(dir+"cos-angle.dat",(/103389/),"float")
> o1=asciiread(dir+"PTC-N31-Hw-Ow.dat",(/103389/),"float")
>  printVarSummary(o1)
>
>
>   opt          = True
>   opt at bin_nice = True
>
>   pdf2obs   = pdfxy(o1,o2, 0, 0, opt)         ; 25 is the default
> printVarSummary(pdf2obs)
> ;  print(pdf2obs)
>
> ;************************************************
> ; create plots
> ;************************************************
> ;  wks = gsn_open_wks("pdf","FIG-10")               ; open a ps file
>   wks = gsn_open_wks("png","FIG-10")               ; open a ps file
> ;  gsn_define_colormap(wks,"prcp_2")
>   gsn_define_colormap(wks,"amwg256")
>
>
>   res      = True
>   res at gsnDraw       = True
>   res at gsnFrame      = True
>   res at cnInfoLabelOn = False
>   res at cnFillOn          = True    ; Turn on color
>   res at cnLinesOn         = False
>   res at cnLineLabelsOn    = False
>   res at lbLabelAutoStride = True
>   res at gsnSpreadColors   = True
>   res at lbOrientation =   "vertical"          ; vertical label bar
>
>   res at trYMinF = 0.6             ; set minimum Y-axis value
>   res at trYMaxF = 1.0             ; set maximum Y-axis value
>   res at trXMinF = -0.8             ; set minimum X-axis value
>   res at trXMaxF = 1.8            ; set maximum X-axis value
>
>
>   res at tmXMajorGrid                = True          ; implement x grid
>   res at tmXMajorGridThicknessF      = 0.5           ; 2.0 is default
>   res at tmXMajorGridLineDashPattern = 2             ; select short dash
> lines
>
>   res at tmYMajorGrid                = True          ; implement y grid
>   res at tmYMajorGridThicknessF      = 0.5           ; 2.0 is default
>   res at tmYMajorGridLineDashPattern = 2             ; select short dash
> lines
>
>   res at cnLevelSelectionMode = "ManualLevels"     ; manually set the
> contour levels with the following 3 resources
>   res at cnMinLevelValF  = 0.001                     ; set the minimum
> contour level
>   res at cnMaxLevelValF  = 6.0                     ; set the maximum contour
> leve
>   res at cnLevelSpacingF = 0.04                     ; set the interval
> between contours
>
>   res at tiXAxisString   = "bond Length (A)"
>   res at tiYAxisString   = "tragectory"
>
>
>   plot = gsn_csm_contour (wks,pdf2obs, res)
>
>
> end
>
>
> --
>
> Subrahmanyam.Sappati
> C/o Dr. Prasenjit Ghosh
> IISER PUNE
> Pin code : 411 008
>
> Phone : +91 94 03 59 35 18
>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>


-- 
Adam Phillips
Associate Scientist,  Climate and Global Dynamics Laboratory, NCAR
www.cgd.ucar.edu/staff/asphilli/   303-497-1726

<http://www.cgd.ucar.edu/staff/asphilli>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20150820/d31c74cf/attachment.html