function main(args)
'reinit'
ctl=subwrd(args,1)
out_PREC=subwrd(args,2)
out_DESVIO=subwrd(args,3)

'open  ' ctl

'media=ave(PREC,t=1,t=4839)'

'sigma=sqrt(ave(pow(prec-media,2),t=1,t=4839))'

'Pmed=media'
'Desvio=sigma'

'set sdfwrite -3dt ' out_PREC
'sdfwrite  Pmed'


'set sdfwrite -3dt ' out_DESVIO
'sdfwrite Desvio'


'quit'
return