Dear Simon,
if you have a look into the
YUCHKDAT
file of your INT2LM run you will see
that the field “PP” (pressure perturbation) is undefined in the uppermost layer (search for “laf” and/or “lbfd” strings;
there you find the PP values). That is deadly!
The reason certainly is that you used the old reference atmosphere (irefatm=1), which makes no sense for
a model domain height of 30 km, since, for the old reference atmosphere, the reference temperature T0 is decreasing steadily,
and it might become negative (below absolute Zero) if the top of the model domain is too high.
Thus, you have to use the new reference atmosphere, which you get if you set irefatm=2 in Namelist “lmgrid”.
The new ref. atmosphere is much closer to the standard atmosphere.
And there are, in principle, no restrictions with respect to the top height of the model domain.
The only restriction is the top height of the driving large-scale model.
Best regards
Hans-Juergen
Dear Simon,
if you have a look into the YUCHKDAT file of your INT2LM run you will see
that the field “PP” (pressure perturbation) is undefined in the uppermost layer (search for “laf” and/or “lbfd” strings;
there you find the PP values). That is deadly!
The reason certainly is that you used the old reference atmosphere (irefatm=1), which makes no sense for
a model domain height of 30 km, since, for the old reference atmosphere, the reference temperature T0 is decreasing steadily,
and it might become negative (below absolute Zero) if the top of the model domain is too high.
Thus, you have to use the new reference atmosphere, which you get if you set irefatm=2 in Namelist “lmgrid”.
The new ref. atmosphere is much closer to the standard atmosphere.
And there are, in principle, no restrictions with respect to the top height of the model domain.
The only restriction is the top height of the driving large-scale model.
Best regards
Hans-Juergen