Factors in the model equations associated with spherical geometry have
been included in QTCM1 beginning with v2.0. However, in v2.0 through
2.2, these factors appeared only in the dx
terms and were omitted
for instance in
term in the divergence. These have now been corrected in v2.3.
The front-end subroutine for the fatd package have also been made
consistent for spherical geometry.
The term of the momentum
equation is
neglected since it would add computation of several additional terms
(such as
, etc) and is small relative to
especially
in the tropics.
The spherical factors are controlled in subroutine parinit
. Array
cosu
is the cosine factor at each ``whole
gridpoint'' latitude location, i.e.
where the u
variable is defined on the C-grid
(Table 3.1). Array cosv
is the same
at the ``half-grid'' locations, i.e. where the v
variable is
defined on the C-grid. Arrays dxu
, dxv
, dyu
,
and dyv
contain the corresponding grid distances in spherical
coordinates. To remove spherical geometry effects, set both
cosu
and cosv
to 1.