cellarea.m

来自「大气模型计算」· M 代码 · 共 28 行

% CELLAREA Calculate area of TOMS cell. Input is refl matrix
% indices.
%
% A = cellarea(i,j,sy,sx)
%
% i,j are the indices of the sy x sx REFL matrix.
%
% $Id: cellarea.m,v 1.7 2006/05/17 14:39:17 danji Exp $

function A = cellarea(i,j,sy,sx);

CONST.EMR = 6371.01e3;
CONST.d2r = pi/180;

[theta,phi] = idx2rad(i,j,sy,sx);

% Grid size
dphi = (180/sy)*CONST.d2r;
dtheta = (360/sx)*CONST.d2r;

% Diagonal points
phimax = phi + dphi/2;
phimin = phi - dphi/2;

A = CONST.EMR^2*dtheta*(cos(phimin)-cos(phimax));

return