<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>
</head>
<body dir="ltr">
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
In the example below, I have 20 values of precipitation in inches.</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Each value is the highest daily precipitation for one of 20 years.</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
I want to estimate the 1-day/100-year GEV precipitation level.</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
For some reason, my older script included multiplying the shape</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
parameter by -1. That gave a GEV value of 65757", definitely wrong.</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
If I get rid of the *-1 part, the GEV value becomes more reasonable,</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
1.73", but too low compared to the 20 annual max values. So I have</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
two questions.</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<ol>
<li>Is there something wrong with the script below?<span></span></li><li>Is there ever a need to multiply any of the parameters by -1? I</li><li style="display:block">did it once before but have no recollection of why I did it.</li></ol>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Thanks, Michael</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
yearmax=(/1.83,1.69,3.24,1.99,1.69,1.79,2.78,1.70,4.85,2.35,3.37,1.72,3.96,2.35,2.29,3.33,3.97,3.53,3.11,1.69/)
<div>param=extval_mlegev(yearmax,0,False)</div>
<div>param_mle_location=param(0)</div>
<div>param_mle_scale=param(1)</div>
<div>param_mle_shape=param(2)*-1.</div>
<div>coeff=-1.*log(1.-1./100.)</div>
<div>gev=param_mle_scale/param_mle_shape*[coeff^(-1.*param_mle_shape)-1.]+param_mle_location</div>
<span>print("param_mle_location="+param_mle_location+" param_mle_scale="+param_mle_scale+" param_mle_shape="+param_mle_shape+" coeff="+coeff+" gev= "+gev)</span><br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span><br>
</span></div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span style="background-color: rgb(255, 255, 0);">param_mle_location=1.65715 param_mle_scale=0.216018 param_mle_shape=2.98224 coeff=0.0100503 gev= 65757.5</span><span><br>
</span></div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span><br>
</span></div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span><span style="background-color:rgb(255, 255, 255);display:inline !important">earmax=(/1.83,1.69,3.24,1.99,1.69,1.79,2.78,1.70,4.85,2.35,3.37,1.72,3.96,2.35,2.29,3.33,3.97,3.53,3.11,1.69/)</span>
<div style="margin:0px;background-color:rgb(255, 255, 255)">param=extval_mlegev(yearmax,0,False)</div>
<div style="margin:0px;background-color:rgb(255, 255, 255)">param_mle_location=param(0)</div>
<div style="margin:0px;background-color:rgb(255, 255, 255)">param_mle_scale=param(1)</div>
<div style="margin:0px;background-color:rgb(255, 255, 255)">param_mle_shape=param(2)*-1.</div>
<div style="margin:0px;background-color:rgb(255, 255, 255)">coeff=-1.*log(1.-1./100.)</div>
<div style="margin:0px;background-color:rgb(255, 255, 255)">gev=param_mle_scale/param_mle_shape*[coeff^(-1.*param_mle_shape)-1.]+param_mle_location</div>
<span style="margin:0px;background-color:rgb(255, 255, 255)">print("param_mle_location="+param_mle_location+" param_mle_scale="+param_mle_scale+" param_mle_shape="+param_mle_shape+" coeff="+coeff+" gev= "+gev)</span><br>
</span></div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span><span style="margin:0px;background-color:rgb(255, 255, 255)"><br>
</span></span></div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span><span style="margin: 0px; background-color: rgb(255, 255, 0);">param_mle_location=1.65715 param_mle_scale=0.216018 param_mle_shape=-2.98224 coeff=0.0100503 gev= 1.72958</span><span style="margin:0px;background-color:rgb(255, 255, 255)"><br>
</span></span></div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span><span style="margin:0px;background-color:rgb(255, 255, 255)"><br>
</span></span></div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span><br>
</span></div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span><br>
</span></div>
<div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div id="Signature">
<div>
<div id="divtagdefaultwrapper" style="font-size:12pt; color:#000000; background-color:#FFFFFF; font-family:Calibri,Arial,Helvetica,sans-serif">
Michael Notaro<br>
<div>Associate Director<br>
</div>
<div>Nelson Institute Center for Climatic Research<br>
</div>
<div>University of Wisconsin-Madison<br>
</div>
<div>Phone: (608) 261-1503<br>
</div>
<div>Email: mnotaro@wisc.edu<br>
</div>
</div>
</div>
</div>
</div>
</body>
</html>