BUFFER IN(32,f,b,L1*m1*2,ios)
           call istat(32,f,uuu,L1,M1,K,IOS)