pro arealcomp, metric,mud,arealvec,vbin_size ; procedure computes areal integration of the array metric tmetric=metric sz=size(tmetric) nrpt=sz(1) tmetric(mud)=0. tt1=tmetric*0.0 tt1(*,*)=1.0 tt1(mud)=0 tt2=tmetric*tt1 arealvec=fltarr(nrpt) for i=0, nrpt-1 do begin tt3=total(tt2(i,*))*vbin_size arealvec(i)=tt3 endfor ;read,'enter number ',jnk end