Standard deviation of particular location and of area

[Raj Kumar Sharma]

Hi, I have file containing global values of variable for whole year. I would like to compute Average & Standard Deviation for Particular Location and also for defined area.I had written script to compute Std Dev for a particular location, however it is showing error. It will be great if someone can help.
Raj
'reinit''set display color white'clear'set gxout shaded''open pqfull.ctl'ln = 60lt = -10it = 1et = 5
'set lon 'ln'set lat 'lt
'd ave(whu,t='it',t='et')'lin=sublin(result,2)avg = subwrd(lin,4)say 'AAAA 'avgwhile(t<'et')sg = sg + (avg-whu)df =  pow(sg,2)endwhiledt = et-it+1std = df/dt'd sqrt(std)'std = subwrd(result,4)say 'STD DEV = ',stdreturn

Dr Raj Kumar

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090915/9d644665/attachment.html