Sponge layers are no longer necessary as of version 2.3 to stabilize the model. If
need arises they can be turn on by defining the preprocessor macro SPONGE
.
In the standard code v2.2, artificial ocean boundary is used and no sponge
functions are turned on. To deal with the boundary instability a sponge layer can be
at latitudes higher than (this cutoff value
is set by the local variable
y0
in subroutine parinit
).
There are five options of sponge layer functions.
The first function spngh0
acts to relax land surface temperature to climatology.
The sponge function spngh1
reduces heat and moisture flux and
spngh2
increase
diffusion in the boundary layer. spngh3
increases momentum damping,
while spngh4
reduces excessive precipitation in high latitudes.
In the version 2.1 release, only spngh0
,spngh1
and spngh3
are active; sponge layer terms spngh2
and spngh4
are commented
out. The functional descriptions of the sponge layer are set in
subroutine parinit
.
Sponge layer functions are defined as:
where is the latitude (in degrees) at the jth grid point.
In Figure 3.1, this is the location where the
u
(and T
) variables are located.