Skip to content
GitLab
Menu
Projects
Groups
Snippets
Loading...
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in
Toggle navigation
Menu
Open sidebar
mpim-sw
cdo
Commits
cb64b952
Commit
cb64b952
authored
Sep 25, 2017
by
Uwe Schulzweida
Browse files
cleanup.
parent
416eaa59
Changes
1
Hide whitespace changes
Inline
Side-by-side
src/afterburnerlib.cc
View file @
cb64b952
...
...
@@ -1665,8 +1665,7 @@ void after_processML(struct Control *globs, struct Variable *vars)
if
(
vars
[
code
].
fourier
==
NULL
)
{
fieldSize
=
vars
[
code
].
hlev
*
globs
->
DimFC
;
vars
[
code
].
fourier
=
alloc_dp
(
fieldSize
,
FieldName
(
code
,
"fourier"
));
vars
[
code
].
fourier
=
alloc_dp
(
fieldSize
,
FieldName
(
code
,
"fourier"
));
sp2fc
(
vars
[
code
].
spectral
,
vars
[
code
].
fourier
,
globs
->
poli
,
vars
[
code
].
hlev
,
globs
->
Latitudes
,
globs
->
Fouriers
,
globs
->
Truncation
);
}
...
...
@@ -1766,8 +1765,7 @@ void after_processML(struct Control *globs, struct Variable *vars)
if
(
vars
[
code
].
hybrid
==
NULL
)
{
fieldSize
=
globs
->
DimGP
*
vars
[
code
].
hlev
;
vars
[
code
].
hybrid
=
alloc_dp
(
fieldSize
,
FieldName
(
code
,
"hybrid"
));
vars
[
code
].
hybrid
=
alloc_dp
(
fieldSize
,
FieldName
(
code
,
"hybrid"
));
after_FC2GP
(
vars
[
code
].
fourier
,
vars
[
code
].
hybrid
,
globs
->
Latitudes
,
globs
->
Longitudes
,
vars
[
code
].
hlev
,
globs
->
Fouriers
);
}
...
...
@@ -1795,7 +1793,7 @@ void after_processML(struct Control *globs, struct Variable *vars)
if
(
globs
->
Orography
==
NULL
)
{
globs
->
Orography
=
alloc_dp
(
globs
->
DimGP
,
"Orography"
);
globs
->
Orography
=
alloc_dp
(
globs
->
DimGP
,
"Orography"
);
if
(
vars
[
GEOPOTENTIAL
].
hybrid
)
after_copy_array
(
globs
->
Orography
,
vars
[
GEOPOTENTIAL
].
hybrid
,
globs
->
DimGP
);
else
...
...
@@ -2221,11 +2219,9 @@ void after_processML(struct Control *globs, struct Variable *vars)
vars
[
VELOPOT
].
needed
||
vars
[
STREAM
].
needed
)
{
if
(
vars
[
DIVERGENCE
].
spectral
==
NULL
)
vars
[
DIVERGENCE
].
spectral
=
alloc_dp
(
globs
->
DimSP
*
globs
->
NumLevelRequest
,
"vars[DIVERGENCE].spectral"
);
vars
[
DIVERGENCE
].
spectral
=
alloc_dp
(
globs
->
DimSP
*
globs
->
NumLevelRequest
,
"vars[DIVERGENCE].spectral"
);
if
(
vars
[
VORTICITY
].
spectral
==
NULL
)
vars
[
VORTICITY
].
spectral
=
alloc_dp
(
globs
->
DimSP
*
globs
->
NumLevelRequest
,
"vars[VORTICITY].spectral"
);
vars
[
VORTICITY
].
spectral
=
alloc_dp
(
globs
->
DimSP
*
globs
->
NumLevelRequest
,
"vars[VORTICITY].spectral"
);
if
(
(
vars
[
U_WIND
].
fourier
==
0
||
vars
[
V_WIND
].
fourier
==
0
)
&&
globs
->
NumLevelRequest
)
Error
(
"uwind or vwind missing!"
);
uv2dv
(
vars
[
U_WIND
].
fourier
,
vars
[
V_WIND
].
fourier
,
...
...
Write
Preview
Markdown
is supported
0%
Try again
or
attach a new file
.
Attach a file
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment