c cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc c subroutine nsumav(wzs,tim,aout,sdout) c implicit none include 'itg.par' real aout(mz),wa(nez,mz),sdout(mz),wzs(nez,mz,nz),tim, * av,sdv integer i,m,n include 'itg.cmn' c do 10 m=1,mz do 10 i=1,ne wa(i,m)=0.0 10 continue c do 20 n=1,nd do 20 m=1,md do 20 i=1,ne wa(i,m)=wa(i,m)+wzs(i,m,n) 20 continue c do 30 m=1,md call timeav(wa(1,m),tim,av,sdv) aout(m)=av sdout(m)=sdv 30 continue c return end