[gradsusr] Model Verification

[Amulya Chevuturi]

Dear Grads Users,

Recently I asked a question about correction in a script building
contingency tables (Asnwered by Dr. Eric Altshuler). I am attaching the
script at the end of the mail.

But this script is for when the model and observation temporal and spatial
extent are same. Even the model and observations resolution is same for
this script to run. Can anyone tell me how to build 2*2 contingency tables
for when the model and observations datasets are at different resolutions
and extent.

With regards,
Amulya Chevuturi
School of Environmental Sciences
Jawaharlal Nehru University
New Delhi, India


******Script for contingency table*************
'reinit'

'open model.ctl'
'set t 3'
'set y 57 208'
'set x 72 228'
'define a = rain'
'close 1'

'open observation.ctl'
'set t 3'
'set y 57 208'
'set x 72 228'
'define b = prec'

n1=0
n2=0
n3=0
n4=0
j=72
while (j<=228)
'set y 'j
i=57
while (i<=208)
'set x 'i
'd a'
c=subwrd(result,4)
'd b'
d=subwrd(result,4)
if(c>0&d>0)
n1=n1+1
else
if(c<=0&d<=0)
n2=n2+1
else
if(c<=0&d>0)
n3=n3+1
else
if(c>0&d<=0)
n4=n4+1
endif
endif
endif
endif
i=i+1
endwhile
j=j+1
endwhile
say n1
say n2
say n3
say n4
'close 1'
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20140501/8b2ceb51/attachment.html