qgread.m

集合卡尔曼滤波(EnKF) 数据同化方法可以避免了EKF 中协方差演变方程预报过程中出现的计算不准确和关于协方差矩阵的大量数据的存储问题,最主要的是可以有效的控制估计误差方差的增长,改善预报的效果。

% function [psi, q, psimean, qmean, t] = qgread(fname)
%
% Reads PSI and Q fields from the QG model output file. Also reads time array
% for the dumps in the output file and calculates mean fields.
%
% @param fname - output file name
% @return psi - 3D array with PSI fields (nr, 129, 129)
% @return q - 3D array with Q fields (nr, 129, 129)
% @return psimean - mean PSI field (129, 129)
% @return qmean - mean Q field (129, 129)
% @return t - array with time values for the dumps in the output file

% File:           qread.m
%
% Created:        31/08/2007
%
% Last modified:  08/02/2008
%
% Author:         Pavel Sakov
%                 CSIRO Marine and Atmospheric Research
%                 NERSC
% Purpose:        Reads PSI and Q fields from the QG model output file. Also
%                 reads time array for the dumps in the output file and
%                 calculates mean fields.
%
% Description:
%
% Revisions:

%% Copyright (C) 2008 Pavel Sakov
%% 
%% This file is part of EnKF-Matlab. EnKF-Matlab is a free software. See 
%% LICENSE for details.

function [psi, q, psimean, qmean, t] = qgread(fname)

    nc = netcdf(fname, 'nowrite');
    psi = nc{'psi'}(:, :, :);
    q = nc{'q'}(:, :, :);
    psimean = squeeze(mean(psi));
    qmean = squeeze(mean(q));
    t = nc{'t'}(:);

  return