[ncl-talk] Joint probability distribution in Logarithmic scale

[Dennis Shea]

You must provide better descriptions of the issue.**Why** can't you normalize?

pdf2obs = pdf2obs/max(pdf2obs)

On Thu, Aug 20, 2015 at 11:57 PM, Subrahmanyam Sappati
<subrahmanyamsappati at gmail.com> wrote:
> 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).
>
> Regards
> Subrahmanyam Sappati
>
>
> On Fri, Aug 21, 2015 at 1:25 AM, Adam Phillips <asphilli at ucar.edu> wrote:
>>
>> 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
>>
>
>
>
> --
>
> 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
>