function phi=dtophi(dmm); % dtophi(dmm) provides a phi size, given % grain size in mms. % phi=-1*log(dmm)/log(2);