subroutine nsumavp(wzs,aout,sdout,net,fave) c implicit none include 'itg.par' include 'itg.cmn' c include 'post.cmn' real aout(mz),wa(fno*nez,mz),sdout(mz), * wzs(fno*nez,mz,nz),fave integer net c Local vars: integer i,m,n real sdv,smax,smin,av c do 10 m=1,mz do 10 i=1,net wa(i,m)=0.0 10 continue c do 20 n=1,nd do 20 m=1,md do 20 i=1,net wa(i,m)=wa(i,m)+wzs(i,m,n) 20 continue c do 30 m=1,md call timeavp(wa(1,m),av,sdv,net,fave,smin,smax) aout(m)=av sdout(m)=sdv 30 continue c return end