<div dir="ltr">Thanks a lot! <br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Apr 29, 2019 at 10:01 PM Dennis Shea <<a href="mailto:shea@ucar.edu">shea@ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>NCL is <span style="color:rgb(255,0,0)"><b>not</b></span> a statistics package and ncl-talk is not an expert on assorted distributions. <br></div><div>Please ask a statistician questions about the Weibull distribution. <br></div><div>Other tools [<b>R</b>, Matlab, Python] may be a better option. <br></div><div>For sure, <b>R</b> has multiple functions that may be useful. <br></div><div>-------<br></div><div dir="ltr"><a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/weibull.shtml" target="_blank"><b>http://www.ncl.ucar.edu/Document/Functions/Built-in/weibull.shtml</b></a></div><div dir="ltr"><br></div><div>Please <b>carefully</b> read the documentation.</div><div>Example one shows a comparison of NCL's and R's result.</div>---<br><div><br></div><div>TRMM: [time | 733] x [latitude | 120] x [longitude | 100]</div><div><br></div><div><div>From experience, daily TRMM has <b>*MANY*</b> 0 values.</div></div><div>You <b>may</b> want to sett all 0 value to _FillValue</div><div><br></div><div>TRMM = <b>where</b>(<b>ismissing</b>(TRMM),TRMM@_FillValue, TRMM)<br></div><div><br></div><div>; Weibull at each grid point<br></div><div><pre><pre><pre><span style="font-family:arial,helvetica,sans-serif">wTRMM = <b>weibull</b>(TRMM, False, 0) ; wTRMM(<b>2</b>,120,100)<br>wTRMM@long_name = "Weibull: Grid point: Shape, Scale"<br>copy_VarCoords(TRMM(0,:,:), wTRMM(0,:,:)) </span><b><br></b></pre><pre><span style="font-family:arial,helvetica,sans-serif"><b>printVarSummary</b>(wTRMM)<br><b>printMinMax</b>(wTRMM,0)<br><b>print</b>("---")<br></span><br><span style="font-family:arial,helvetica,sans-serif">; Weibull over entire region</span><pre><pre><span style="font-family:arial,helvetica,sans-serif">WTRMM = <b>weibull</b>(<b>ndtooned</b>(TRMM), False, 0) ; WTRMM(<b>2</b>)<br>WTRMM@long_name = "Weibull: All points: Shape, Scale"<br>printVarSummary(WTRMM)<br>print(WTRMM,0)<br>print("---")<br></span></pre></pre><span style="font-family:arial,helvetica,sans-serif">==============================</span><span style="font-family:arial,helvetica,sans-serif">===========<br><br>opt = True
opt@confidence = 0.90</span><pre><pre><span style="font-family:arial,helvetica,sans-serif">wTRMM = <b>weibull</b>(TRMM, opt, 0 ) ; wTRMM(6,120,100)<br>wTRMM@long_name = "Weibull: Grid point: Shape, Scale; ShapeLow, ShapeHigh, ScaleLow, ScaleHigh"<br>copy_VarCoords(TRMM(0,:,:), wTRMM(0,:,:))<b><br>printVarSummary</b>(wTRMM)<b><br>printMinMax</b>(wTRMM,0)<b><br>print</b>("---")</span><br> <b><br></b><br><span style="font-family:arial,helvetica,sans-serif">WTRMM = <b>weibull</b>(<b>ndtooned</b>(TRMM), opt, 0) ; WTRMM(<b>6</b>)<br>WTRMM@long_name = "Weibull: All points: Shape, Scale; ShapeLow, ShapeHigh, ScaleLow, ScaleHigh"<br>print(WTRMM,0)<br>print("---")<br><br>==============================</span><span style="font-family:arial,helvetica,sans-serif">============<br></span></pre><pre><span style="font-family:arial,helvetica,sans-serif"><b>You</b> will have to calculate the CDF and PDF. NCL does not have functions to do this.<br><br></span></pre><pre><span style="font-family:arial,helvetica,sans-serif">Good Luck</span><br></pre></pre></pre></pre></pre></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, Apr 27, 2019 at 9:19 PM Md. Jalal Uddin <<a href="mailto:dmjalal90@gmail.com" target="_blank">dmjalal90@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr">Hi all,<div><br></div><div>I want to calculate PDF and CDF with Weibull distribution using shape 5 and scale 2 parameters. I found one function in NCL that calculate the extreme value (<a href="https://www.ncl.ucar.edu/Document/Functions/Ext_val/extval_weibull.shtml" target="_blank">https://www.ncl.ucar.edu/Document/Functions/Ext_val/extval_weibull.shtml</a> ).</div><div><br></div><div>However, I don't want to calculate the extreme value. I would like to produce a similar figure in the attached figure (adapted from Prat and Nelson 2013a). </div><div><br></div><div>Any suggestions would be appreciated.</div><div><br></div><div><div>Variable: TRMM</div><div>Type: float</div><div>Total Size: 35184000 bytes</div><div> 8796000 values</div><div>Number of Dimensions: 3</div><div>Dimensions and sizes: [time | 733] x [latitude | 120] x [longitude | 100]</div><div>Coordinates:</div><div> time: [ 0..131112]</div><div> latitude: [0.125..29.875]</div><div> longitude: [75.125..99.875]</div><div>Number Of Attributes: 7</div><div> long_name : precipitation (mm/hr)</div><div> _FillValue : -9999.9</div><div> missing_value : -9999.9</div></div><div><br></div><div><br></div><div>Best regards,</div><div>Jalal<br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail-m_-148493187503154460gmail-m_-4906245794854826201gmail-m_-4913048755411867974gmail-m_8551023608905312279gmail-m_3238675581061288564gmail_signature"><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><font size="2" color="#000000"><b>Md. Jalal Uddin</b></font></div><div dir="ltr">MSc in Applied Meteorology (English Language)<br>Nanjing University of Information, Science and Technology, China</div><div dir="ltr"><span style="letter-spacing:0.2px">B.Sc. in Disaster Management (Hons.) </span><br></div><div dir="ltr"><span style="letter-spacing:0.2px">Patuakhali Science and Technology University, Bangladesh.</span></div><div dir="ltr">Cell: +8613260859092, <span style="letter-spacing:0.2px">+8801792052662</span></div><div dir="ltr"><font size="2">Web: <a href="http://www.dmjalal90.weebly.com" target="_blank">www.dmjalal90.weebly.com</a> <br>Facebook: jalal.hossen.39 <br>LinkedIn: <a href="https://bd.linkedin.com/in/md-jalal-uddin-80a026b0" target="_blank">https://bd.linkedin.com/in/md-jalal-uddin-80a026b0</a> <br>Twitter: dmjalal90 <br>Skype: dmjalal90</font><br></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></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>
</blockquote></div>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><font color="#000000" size="2"><b>Md. Jalal Uddin</b></font></div><div dir="ltr">MSc in Applied Meteorology (English Language)<br>Nanjing University of Information, Science and Technology, China</div><div dir="ltr"><span style="letter-spacing:0.2px">B.Sc. in Disaster Management (Hons.) </span><br></div><div dir="ltr"><span style="letter-spacing:0.2px">Patuakhali Science and Technology University, Bangladesh.</span></div><div dir="ltr">Cell: +8613260859092, <span style="letter-spacing:0.2px">+8801792052662</span></div><div dir="ltr"><font size="2">Web: <a href="http://www.dmjalal90.weebly.com" target="_blank">www.dmjalal90.weebly.com</a> <br>Facebook: jalal.hossen.39 <br>LinkedIn: <a href="https://bd.linkedin.com/in/md-jalal-uddin-80a026b0" target="_blank">https://bd.linkedin.com/in/md-jalal-uddin-80a026b0</a> <br>Twitter: dmjalal90 <br>Skype: dmjalal90</font><br></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div>