<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; color: rgb(0, 0, 0); font-size: 14px; font-family: Calibri, sans-serif; ">
<div>Sunmin</div>
<div><br>
</div>
<div>You're on the right lines trying to mask out values which would cause division by zero, but your mask() statement is incorrect.</div>
<div><br>
</div>
<div>Since the denominator is (1-prob^2), you need to remove values of prob=1 (NOT prob = 0 as you have done). You can also use the where() function instead of mask(), i.e.</div>
<div><br>
</div>
<div>> probt = where(prob.ne.1., <span style="font-family: Times; line-height: 16px; ">prob*sqrt((Nr-2)/(1-prob^2)), prob@_FillValue)</span></div>
<div><br>
</div>
<div>Hope that helps,</div>
<div><br>
</div>
<div>Will</div>
<div><br>
</div>
<div><br>
</div>
<span id="OLK_SRC_BODY_SECTION">
<div style="font-family:Calibri; font-size:11pt; text-align:left; color:black; BORDER-BOTTOM: medium none; BORDER-LEFT: medium none; PADDING-BOTTOM: 0in; PADDING-LEFT: 0in; PADDING-RIGHT: 0in; BORDER-TOP: #b5c4df 1pt solid; BORDER-RIGHT: medium none; PADDING-TOP: 3pt">
<span style="font-weight:bold">From: </span>Sunmin Park <<a href="mailto:mireiyue@gmail.com">mireiyue@gmail.com</a>><br>
<span style="font-weight:bold">Date: </span>Friday, 16 January 2015 9:03 AM<br>
<span style="font-weight:bold">To: </span>"<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>" <<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>><br>
<span style="font-weight:bold">Subject: </span>[ncl-talk] fatal:divide: Division by 0, Can't continue<br>
</div>
<div><br>
</div>
<span style="color: rgb(0, 0, 0); font-family: Times; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: 16px; orphans: 2; text-align: -webkit-auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; display: inline !important; float: none; ">probt
= prob1*sqrt((Nr-2)/(1-prob1^2))</span></span>
<p style="font-size:10pt; line-height:10pt; calibri"><br>
<br>
University of Tasmania Electronic Communications Policy (December, 2014). <br>
This email is confidential, and is for the intended recipient only. Access, disclosure, copying, distribution, or reliance on any of it by anyone outside the intended recipient organisation is prohibited and may be a criminal offence. Please delete if obtained
in error and email confirmation to the sender. The views expressed in this email are not necessarily the views of the University of Tasmania, unless clearly intended otherwise.
</p>
</body>
</html>