
292
J.
S.
FREDERIKSEN
mensional basic Northern Hemisphere flow in a multilevel primitive equa-
tion model.
4.1.
Model Details
This study was camed out using the a-coordinate spherical primitive
equation spectral model developed by Bourke (1974) and Bourke
et
al.
(1977), and used by Frederiksen
(
198
la,b) for studying life cycles of nonlin-
ear baroclinic waves growing on initially zonally averaged flows. The equa-
tions for vorticity, divergence, thermodynamics, and continuity are de-
scribed in Section 2 of Frederiksen
(
198
1
a). Since initial three-dimensional
basic states are considered constant, forcing functions are included on the
right-hand sides of
Eqs.
(2.la)-(2.ld) of Frederiksen (1981a). The model
also includes topography,
so
the term
-OW.
is added to the right-hand side
of the divergence equation, where
a.
is the geopotential height ofthe North-
em
Hemisphere topography.
A five-level version of the model with five equispaced
a
levels of
a
=
0.1,
0.3,0.5,0.7, and 0.9 and a rhomboidal truncation scheme with
a
truncation
number
J
=
21 are used. The model includes biharmonic diffusion and
surface drag. For the basic state we take the monthly averaged flow field for
January 1978, which has zonal variations that are qualitatively quite similar
to those of the average winter field(shown in Fig.
1)
in the upper troposphere.
The streamfunction fields for the basic state are shown in Fig.
1
of Frederik-
sen and Pun
(
1985). The forcing functions mentioned above are chosen in
such a way that, in the absence of drag and diffusion, the flow field for
January 1978 is an
exact
steady-state solution of the equations of motion.
However, if the basic state is perturbed, the perturbation will extract energy
from it and grow.
For the small-amplitude perturbation to the basic state, the fastest growing
linear instability mode for the January 1978 basic state obtained from a
three-dimensional instability calculation in a quasi-geostrophic model (Fre-
deriksen, 1983a) is used. The disturbance amplitude was chosen to
be
suffi-
ciently large,
so
that time-stepping errors are comparatively small, but small
enough
so
that the different interesting stages of evolution can be studied.
The disturbance kinetic energy per unit mass is
6.03
X
m2 sec2, and
the disturbance streamfunction at
a
=
0.7 is shown in Fig. 7. We see that the
disturbance is a monopole cyclogenesis mode with largest amplitudes in the
storm tracks in the Atlantic and Pacific Oceans. That is, it is similar to Mode
1
of
Case
1
in Fig. 2, but
is
somewhat more localized.
The model was integrated for a
period
of
10
days from the perturbed initial
conditions and the evolution of the disturbance was examined. Here the