wrfv2_fire/phys/next_mat.m

   1    function mat=next_mat(f)\r
   2     %mat=reshape(a(k:k+mm*nn),mm,nn);\r
   3     [sz,count]=fscanf(f,'%g',[1 2]);\r
   4     switch count\r
   5         case 0,\r
   6             mat=[];\r
   7         case 1\r
   8             error('not enough terms in size')\r
   9         otherwise\r
  10         %fprintf('reading matrix size %g %g\n',sz)\r
  11         if any(any(sz~=round(sz) | sz < 0)),\r
  12             error('bad size')\r
  13         end\r
  14         [mat,count]=fscanf(f,'%g',sz);\r
  15         if(count<sz(1)*sz(2)),\r
  16             error('not enough terms in matrix')\r
  17         end\r
  18     end\r
  19     end\r