util1_jan/read_array_m.m

   1 function a=read_array_m(f,num1,num2);\r
   2 % a=read_array_m(f,num1,num2)\r
   3 % read array produced by matching call write_array_m\r
   4 % in module_fr_sfire_util.F\r
   5 %\r
   6 % a=read_array_m(f)\r
   7 % read array from file f\r
   8 \r
   9 % Jan Mandel, 2008\r
  10 \r
  11 if exist('num1','var'),\r
  12     if ~exist('num2','var'),\r
  13         num2=-1;\r
  14     end\r
  15     f=[file_name(f,num1,num2),'.txt'];\r
  16 end\r
  17 b=load(f);\r
  18 l=length(b);\r
  19 its=b(1);\r
  20 ite=b(2);\r
  21 jts=b(3);\r
  22 jte=b(4);\r
  23 kts=b(5);\r
  24 kte=b(6);\r
  25 m=(ite-its+1);\r
  26 n=(jte-jts+1);\r
  27 o=(kte-kts+1);\r
  28 s=m*n*o+6;\r
  29 fprintf(1,'matrix size %i:%i %i:%i %i:%i from file %s\n',its,ite,jts,jte,kts,kte,f)\r
  30 if l~=s,\r
  31     error(sprintf('incorrect file length, should be %i',s))\r
  32 end\r
  33 a=reshape(b(7:s),[m,n,o]);\r
  34 end\r