qgread.m

来自「集合卡尔曼滤波(EnKF) 数据同化方法可以避免了EKF 中协方差演变方程预报过」· M 代码 · 共 45 行

% 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