fid = fopen('gpcp.2000jan.satl.dat','r','b');
[A,count]=fread(fid,[8184,30],'float');
B=reshape(A,[132,62,30]);