4

m ファイルを実行して、2 つの変数を作成していClusますWatts_Map。両方の変数を「.mat」で終わるファイルに保存したいと考えています。変数は両方とも (1152,241,319) として寸法付けされ、1152 は 0.3125 度の増分で 360 度の経度、241 は 30S ～ 30N の緯度として 0.25 度の増分で、319 時間ステップを超えます。コードは最後まで機能し、そこでエラーが発生します。

[Warning: Variable 'Clus' cannot be saved to a MAT-file whose version is older
than 7.3.
To save this variable, use the -v7.3 switch.
Skipping...] 
[Warning: Variable 'Watts_Map' cannot be saved to a MAT-file whose version is
older than 7.3.
To save this variable, use the -v7.3 switch.
Skipping...]

私は Matlab バージョン R2014a を使用しているため、これが最新バージョンであると考えられます。また、同じ正確なコードをより小さな空間ドメイン (ただし 2920 タイムステップ以上) でエラーなしで実行しました。

%  Clear all variables, initialize counter, indexed by timestep
clc;
clear all;

rain = NaN(1152,241,319);
Clus = NaN(1152,241,319);
Areas_Final = NaN(500,319); 
Wattage_Final = NaN(500,319);
Cluster_size = zeros(319,1);
Watts_Map = zeros(1152,241,319);

for year = 2000%:2008;
    Nyear = sprintf('%04d',year);

    %  Call on the files for each year
    filename = (['pr_3hr_GFDL-HIRAM-C360_amip_r1i1p1_' num2str(Nyear) '010100-' num2str(Nyear) '123123_subset_TROPICS.nc']);
    disp(filename)
    disp(year)

    rain_rate = ncread(filename,'pr');

    %  Call on each timestep
    for i = 960:4:2236; % Approx May 1st-Sep 30th

        %  Set extract subset for region, mask land areas, for each
        %  timestep

        disp(i)
        rain = rain_rate(:,:,i);

        % Eliminate bad/no data points
        index_rain = (rain >= (5.4e-05)) & (rain < 1e-2); % 0.2mm/hr is min rain rate

        % Cluster each morning and afternoon matrix
        Clus(:,:,i) = cluster_it(index_rain);

        % Calculate cluster areas
        Areas = cluster_areas_TROPICS(Clus(:,:,i));
        Areas_Final(1:length(Areas),i) = Areas;

        % Calculate cluster wattages
        Wattage = cluster_wattage(Clus(:,:,i),rain);

        Cluster_size(i,1) = max(max(Clus(:,:,i)));

        % Create dummy matricies to populate and use to create the wattage
        % maps
        D = zeros(1152,241);
        E = squeeze(Clus(:,:,i));
        for index = 1:Cluster_size(i);
            D(E == index) = Wattage(index);
        end

        Watts_Map(:,:,i) = D;

        % Clear the dummy matricies
        clear D E

    end

    % Save the output as a .mat file
    file_out = sprintf(num2str(Nyear), year);
    matfile = [file_out '_TROPICS_Watts_Maps_inc_Land_low_rain_threshold.mat'];
    save(matfile, 'Clus', 'Watts_Map');

    % Clear unneeded variables to save memory and prevent overwriting
    clear index_rain rain Areas Wattage Clus file_out filename Areas_Final rain_rate Watts_Map Cluster_size year matfile

end
4

1 に答える 1

9

現在のバージョンでも、デフォルトの形式は v7.3 ではありません。次のように指定する必要があります。

save(matfile, '-v7.3', 'var1', 'var2');

「一般設定」でデフォルトを設定することもできます。

ここに画像の説明を入力

v7.3 は R2008a の時点でデータを圧縮することに注意してください。

-v7.3 MAT

ファイルの圧縮 save 関数の -v7.3 オプションを使用して、サイズが 2 GB を超えるデータ項目を MAT ファイルに格納できます。このオプションは、MATLAB R2006b で導入されました。MATLAB R2008a では、save はこれらの MAT ファイルを圧縮するようになりました。

v7.3 ファイルのファイル サイズが大きくなるのを見てきましたが、それは過去のことであり、おそらく現在は改善されています。

于 2014-12-02T01:20:03.357 に答える