
<html dir="ltr">

<head>

<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">

<style id="owaParaStyle" type="text/css">P {margin-top:0;margin-bottom:0;}</style>

</head>

<body ocsi="0" fpstyle="1">

<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">Hi all<br>

<br>

I am trying to determine the significance of precipitation anomalies at each gridpoint over the Indian Ocean.<br>

I have calculated the pdfs at each gridpoint, in the script attached below and would like to display the spatial map of percentages from the pdfs over the Indian ocean to see how extreme the anomalies are over the chosen domain from 1000-1849AD.<br>

If anyone has advice on how to do this it would be much appreciated.<br>

<br>

Kindest Regards<br>

Melissa<br>

<br>

<br>

Code:<br>

;----------------------------------------------------------------------<br>

; pdf_1.ncl<br>

;----------------------------------------------------------------------<br>

load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl&quot;&nbsp;&nbsp; <br>

load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl&quot;&nbsp;&nbsp;&nbsp; <br>

load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl&quot; <br>

load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl&quot;<br>

begin<br>

<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>

&nbsp;<br>

h=addfile(&quot;/mnt/geog/ml382/McGill_Project/pr/Scen2/PDF/pdf_ens_correct.nc&quot;, &quot;r&quot;)&nbsp;

<br>

<br>

&nbsp;<br>

&nbsp; latS&nbsp;&nbsp; =&nbsp; -70.<br>

&nbsp; latN&nbsp;&nbsp; =&nbsp; 50. <br>

&nbsp; lonL&nbsp;&nbsp; =&nbsp; 0.<br>

&nbsp; lonR&nbsp;&nbsp; =&nbsp; 150.<br>

<br>

&nbsp;<br>

&nbsp; pr&nbsp;&nbsp;&nbsp; = h-&gt;pr<br>

&nbsp; lat&nbsp;&nbsp; = h-&gt;lat<br>

&nbsp; lon&nbsp;&nbsp; = h-&gt;lon<br>

<br>

nlat&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = dimsizes(lat)<br>

mlon&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = dimsizes(lon)<br>

<br>

&nbsp; opt&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = True<br>

&nbsp; opt@bin_nice = True<br>

&nbsp; <br>

&nbsp;&nbsp; do nl = 0, nlat-1<br>

&nbsp;&nbsp;&nbsp; do ml = 0, mlon-1<br>

<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nbins = 25<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; pdf_gridpt = pdfx( pr(:,nl,ml),nbins,opt )<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; printVarSummary(pdf_gridpt)<br>

<br>

&nbsp;&nbsp;&nbsp; <br>

printVarSummary(pdf_gridpt)<br>

print(pdf_gridpt)<br>

<br>

<br>

<br>

;create plot<br>

<br>

&nbsp; res&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = True&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; plot mods desired<br>

<br>

&nbsp; res@cnFillOn&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = True&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; turn on color for contours<br>

&nbsp; res@cnLinesOn&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = False&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; turn off contour lines<br>

&nbsp; res@cnLineLabelsOn&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = False&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; turn off contour line labels<br>

&nbsp; res@gsnScalarContour&nbsp;&nbsp;&nbsp;&nbsp; = True&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; contours desired<br>

&nbsp; ;res@cnFillPalette&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = cmap(17:,:)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; subset the color map<br>

<br>

&nbsp;<br>

&nbsp; res@lbOrientation&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = &quot;Vertical&quot;&nbsp;&nbsp;&nbsp;&nbsp; ; vertical label bar<br>

&nbsp; ;res@pmLabelBarOrthogonalPosF = -0.01&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; move label bar closer<br>

<br>

; note, when doing a subregion, NCL determines the range of the data from<br>

; the full domain. If you wish to just consider the domain you are plotting,<br>

; you must manually set those levels.<br>

<br>

&nbsp; ;res@cnLevelSelectionMode = &quot;ManualLevels&quot;&nbsp;&nbsp;&nbsp;&nbsp; ; set manual contour levels<br>

&nbsp; ;res@cnMinLevelValF&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 24.0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; set min contour level<br>

&nbsp; ;res@cnMaxLevelValF&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 29&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; set max contour level<br>

&nbsp; ;res@cnLevelSpacingF&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0.10&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; set contour spacing<br>

<br>

&nbsp; res@vcRefMagnitudeF&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1.0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; define vector ref mag<br>

&nbsp; res@vcRefLengthF&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0.025&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; define length of vec ref<br>

&nbsp; ;res@vcRefAnnoOrthogonalPosF&nbsp;&nbsp; = -1.0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; move ref vector<br>

&nbsp; res@vcRefAnnoArrowLineColor&nbsp;&nbsp; = &quot;black&quot;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; change ref vector color<br>

&nbsp; res@vcRefAnnoArrowUseVecColor = False&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; don't use vec color for ref<br>

<br>

&nbsp; res@vcGlyphStyle&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = &quot;CurlyVector&quot;&nbsp;&nbsp;&nbsp;&nbsp; ; turn on curley vectors<br>

&nbsp; res@vcLineArrowColor&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = &quot;black&quot;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; change vector color<br>

&nbsp; res@vcLineArrowThicknessF&nbsp;&nbsp; = 2.0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; change vector thickness<br>

&nbsp; res@vcVectorDrawOrder&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = &quot;PostDraw&quot;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; draw vectors last<br>

&nbsp; res@vcMinDistanceF = 0.02<br>

&nbsp;<br>

&nbsp;&nbsp; ;res@gsnCenterString = &quot;Difference in MME Wind Anomaly for Scen1 for the Non-Monsoon Season&quot;<br>

&nbsp;&nbsp; res@gsnLeftString&nbsp;&nbsp; = &quot;&quot;<br>

&nbsp;&nbsp; res@gsnMainString&nbsp;&nbsp; = &quot;&quot; <br>

&nbsp;&nbsp; res@gsnRightString&nbsp; = &quot;&quot;<br>

<br>

&nbsp; res@mpMinLonF&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; lonL&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; select a subregion<br>

&nbsp; res@mpMaxLonF&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; lonR<br>

&nbsp; res@mpMinLatF&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; latS<br>

&nbsp; res@mpMaxLatF&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; latN<br>

<br>

<br>

&nbsp;pres&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = True<br>

<br>

&nbsp; pres@mpFillOn = False<br>

&nbsp; pres@cnFillOn&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = True&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; color on<br>

&nbsp; pres@cnLinesOn&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = False&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; turn off contour lines<br>

&nbsp; pres@gsnScalarContour&nbsp;&nbsp;&nbsp;&nbsp; = True&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; vectors over contours<br>

&nbsp; pres@gsnSpreadColors&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = True&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; use full colormap<br>

&nbsp; ;res@gsnSpreadColorEnd&nbsp;&nbsp;&nbsp; = -4&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; last color used<br>

&nbsp; pres@gsnAddCyclic&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = False<br>

&nbsp; ;pres@gsnDraw&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = False&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; do not draw the plot<br>

&nbsp; ;pres@gsnFrame&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = False&nbsp; <br>

&nbsp; pres@lbLabelBarOn&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = True<br>

<br>

pres@gsnCenterString = &quot;Difference in MME Wind Anomaly for Scen1 for the Non-Monsoon Season&quot;<br>

&nbsp;&nbsp; pres@gsnLeftString&nbsp;&nbsp; = &quot;&quot;<br>

&nbsp;&nbsp; pres@gsnMainString&nbsp;&nbsp; = &quot;&quot; <br>

&nbsp;&nbsp; pres@gsnRightString&nbsp; = &quot;&quot;<br>

<br>

&nbsp; pres@mpMinLonF&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; lonL&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; select a subregion<br>

&nbsp; pres@mpMaxLonF&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; lonR<br>

&nbsp; pres@mpMinLatF&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; latS<br>

&nbsp; pres@mpMaxLatF&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; latN<br>

<br>

<br>

&nbsp;pres@cnLevelSelectionMode= &quot;ManualLevels&quot;&nbsp; ; manual levels<br>

&nbsp; pres@cnMinLevelValF&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 0.2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; min level<br>

&nbsp; pres@cnMaxLevelValF&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 2.8&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; max level<br>

&nbsp; pres@cnLevelSpacingF&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 0.2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; contour spacing<br>

<br>

&nbsp; <br>

&nbsp; wks&nbsp; = gsn_open_wks (&quot;X11&quot;,&quot;PDF_mounttambora&quot;)<br>

&nbsp; gsn_define_colormap(wks,&quot;WhiteBlueGreenYellowRed&quot;)&nbsp; ; choose colormap<br>

<br>

&nbsp; prres&nbsp; = True<br>

&nbsp; prres@gsnCenterString = &quot;Precipitation PDF&quot;<br>

<br>

&nbsp; <br>

&nbsp;plot = gsn_csm_xy (wks, pdf_gridpt@bin_center, pdf_gridpt, prres)<br>

&nbsp;plot1= gsn_csm_contour_map_ce(wks, pdf_gridpt, res)<br>

&nbsp;<br>

<br>

&nbsp;&nbsp;&nbsp;&nbsp; end do<br>

&nbsp;&nbsp;&nbsp; end do<br>

end<br>

<br>

</div>

</body>

</html>