<div dir="ltr">Dear Grads Users<div><br></div><div>I am trying to build contingency table for model verification using grads. Currently my model and observation are both at exactly same resolution. So I have written the following script. But it is not giving me the count of Yes/Yes, No/No etc correctly. How am I supposed to build contingency table using grads if someone can please help.</div>

<div><br></div><div>Following is the script. </div><div><br></div><div><div>&#39;reinit&#39;</div><div><br></div><div>&#39;open model.ctl&#39;</div><div>&#39;set t 3&#39;</div><div>&#39;set y 57 208&#39;</div><div>&#39;set x 72 228&#39;</div>
<div>&#39;define a = rain&#39;</div><div>&#39;close 1&#39;</div><div><br></div><div>&#39;open observation.ctl&#39;</div><div>&#39;set t 3&#39;</div><div>&#39;set y 57 208&#39;</div><div>&#39;set x 72 228&#39;</div><div>&#39;define b = prec&#39;</div>
<div><br></div><div>n1=0</div><div>n2=0</div><div>n3=0</div><div>n4=0</div><div>j=72<br></div><div>while (j&lt;=228)</div><div>&#39;set y &#39;j</div><div>i=57</div><div>while (i&lt;=208)</div><div>&#39;set x &#39;i</div>
<div>if(a&gt;0&amp;b&gt;0)</div><div>n1=n1+1</div><div>else</div><div>if(a&lt;=0&amp;b&lt;=0)</div><div>n2=n2+1</div><div>else</div><div>if(a&lt;=0&amp;b&gt;0)</div><div>n3=n3+1</div><div>else</div><div>if(a&gt;0&amp;b&lt;=0)</div>
<div>n4=n4+1</div><div>endif</div><div>endif</div><div>endif</div><div>endif</div><div>i=i+1</div><div>endwhile</div><div>j=j+1</div><div>endwhile</div><div>say n1</div><div>say n2</div><div>say n3</div><div>say n4</div></div>
<div>&#39;close 1&#39;</div><div><br></div><div>Or is anyother method available for making contingency tables for forecast verification with grads?</div><div><br></div><div>Please let me know.</div><div><br></div><div>Regards</div>
<div>Amulya Chevuturi</div><div><br></div><div><br></div></div>