mat2asc.m

来自「麻省理工学院的人工智能工具箱,很珍贵,希望对大家有用!」· M 代码 · 共 36 行

%*******Contenido de dmn****************************
dmn.lvl=[1000,850];
dmn.tim=[00,06,12,18];
dmn.par=[130];
%*******Contenido de dmn****************************
indObj=1:5538;
[YY,MM,DD]=datevec((indObj)+datenum(1979,1,1)-1);

Camino='./mat/';
Camino2='./asc/';

nP=size(dmn.par,2);
nH=size(dmn.tim,2);
nN=size(dmn.lvl,2);

%diamax=fix(dmn.tim(nH)/24);
for P=1:nP
   for H=1:nH
      hora=dmn.tim(H);
      for N=1:nN
         eval(['load ' Camino 'tim' num2str(hora,'%02d') '/par' num2str(dmn.par(P),'%03d') '_' num2str(dmn.lvl(N),'%04d') '.mat' ' campo']);  

         disp([Camino 'tim' num2str(hora,'%02d') '/par' num2str(dmn.par(P),'%03d') '_' num2str(dmn.lvl(N),'%04d')]);
         
         fid=fopen([Camino2 'tim' num2str(hora,'%02d') '/par' num2str(dmn.par(P),'%03d') '_' num2str(dmn.lvl(N),'%04d') '.asc'],'wb');
         for i=1:size(campo,1)
            fprintf(fid,'%04d%02d%02d',YY(i),MM(i),DD(i));
            fprintf(fid,' %8g',campo(i,:));
            fprintf(fid,'\n');
         end
         fclose(fid);
         
      end
   end
end