'reinit'

**Calculate the mixed reanalysis average**
'open indian.ctl'
'define average =           avgrain'
'close 1'

'open cfsr.ctl'
'define average = average + avgrain'
'close 1'

'open merra.ctl'
'define average = average + avgrain'

'define average = average/3'
'close 1'


**Standard deviation**
'open indian.ctl'
'open cfsr.ctl'
'open merra.ctl'

count = 1

'define avgsum = 0'

WHILE (count <= 5)
  'define avgsum = avgsum + pow((saavg1'count'.1 - average),2) + pow((saavg2'count'.1 - average),2) + pow((saavg3'count'.1 - average),2)+ pow((saavg1'count'.2 - average),2) + pow((saavg2'count'.2 - average),2) + pow((saavg3'count'.2 - average),2)+ pow((saavg1'count'.3 - average),2) + pow((saavg2'count'.3 - average),2) + pow((saavg3'count'.3 - average),2)'

  count = count + 1

ENDWHILE **count

'define stddev = sqrt(avgsum/(45-1))'

'set gxout fwrite'
'set fwrite evapsource.grd'
'set x 1 192'
'set y 1 94'

'd average'
'd stddev'

'disable fwrite'
'set gxout contour'