**Script to obtain the spatial RMSE of monthly
**climatology rainfall from aphrodite and ccam 
**for comparison 07Sept2011
n=1
while(n<10)
'sdfopen /home/mar/Model/analysis/aphro/monthly/rraphro_20phil.0'n'.nc'
'define aph=precip'
'close 1'
'sdfopen /home/mar/Model/analysis/a1b/monthly_climat/rr/rra1b_1971to00.0'n'.nc'
'define dif=aph-rnd24'
'define sq=dif*dif'
'define am=amean(sq,lon=116.82,lon=126.72,lat=4.86,lat=20.92)'
'define rmse=sqrt(am)'
'd rmse'
rain=subwrd(result,4)
output=rain
res=write(rmsea1b.csv,output,append)
'undefine aph'
'undefine dif'
'undefine sq'
'undefine am'
'undefine rmse'
'close 1'
n=n+1
endwhile
************************************************************
j=10
while(j<13)
'sdfopen /home/mar/Model/analysis/aphro/monthly/rraphro_20phil.'j'.nc'
'define aph=precip'
'close 1'
'sdfopen /home/mar/Model/analysis/a1b/monthly_climat/rr/rra1b_1971to00.'j'.nc'
'define dif=aph-rnd24'
'define sq=dif*dif'
'define am=amean(sq,lon=116.82,lon=126.72,lat=4.86,lat=20.92)'
'define rmse=sqrt(am)'
'd rmse'
rain=subwrd(result,4)
output=rain
res=write(rmsea1b.csv,output,append)
'undefine aph'
'undefine dif'
'undefine sq'
'undefine am'
'undefine rmse'
'close 1'
j=j+1
endwhile
'quit'