Discussions are carried out on the vertical discretization of current atmospheric models.It is pointed out that there exist problems in the integration of the hydrostatic equation and the computation of vertical advec...Discussions are carried out on the vertical discretization of current atmospheric models.It is pointed out that there exist problems in the integration of the hydrostatic equation and the computation of vertical advection,vertical diffusion and so on.Then some possible ways for solving or alleviating them are suggested.Finally,the choice of vertical coordinate and basis functions is discussed.展开更多
Two schemes for vertical discretization of the model are proposed,one with equal △lnσ and the other in terms of Tschebyscheff polynomials.It is proved that in adiabatic and inviscid cases if the meteorological eleme...Two schemes for vertical discretization of the model are proposed,one with equal △lnσ and the other in terms of Tschebyscheff polynomials.It is proved that in adiabatic and inviscid cases if the meteorological elements and related physical quantities are continuous in time and in the horizontal,the total energy and total mass are conserved within a high approximation respectively, and there is a correct conversion between total kinetic and total potential energy.Numerical computations show that the schemes both have high accuracy.For example,in integrating the hydrostatic equation the computational errors of geopotential height resulting from the schemes are much less than those resulting from EC79 in a-coordinate.展开更多
Through analysis and numerical computation of ECMWF's discrete scheme of hydrostatic equation (Baede et al. 1979),it has been found that in the case of equal △σ there exist systematic errors in the scheme.The er...Through analysis and numerical computation of ECMWF's discrete scheme of hydrostatic equation (Baede et al. 1979),it has been found that in the case of equal △σ there exist systematic errors in the scheme.The error E caused by taking the arithmetic mean of the geopotential heights of two adjacent half a-levels as the geopotential height of the cor- responding integer a-level,increases with height and has an unacceptable maximum in the vicinity of the top of the at- mosphere;however,the errors caused by the temperature treatment are generally small.On the other hand,if an uneven △σ-scheme in which the levels in the upper and lower atmosphere are denser than those in the middle atmosphere,is adopted,then E can be much reduced.However,if the resolution of the original equal Art-scheme doubles,then E can only be found to be much reduced in the troposphere and that in the vicinity of the atmospheric top is still unacceptable. Two numerical schemes for improvement have been presented.Of them one is the same as the ECMWF's scheme, but with equal △lnσ,and the other is to integrate the equation by the use of Tschebyscheff polynomials T and to adopt the zeros of TN as the atmospheric levels where N is the total number of levels.The results show that with both schemes the computational errors can be much reduced,especially the latter,due to the fact that the errors of three statistical types are generally less than the root mean square error of the geopotential heights reported in TEMP.展开更多
文摘Discussions are carried out on the vertical discretization of current atmospheric models.It is pointed out that there exist problems in the integration of the hydrostatic equation and the computation of vertical advection,vertical diffusion and so on.Then some possible ways for solving or alleviating them are suggested.Finally,the choice of vertical coordinate and basis functions is discussed.
文摘Two schemes for vertical discretization of the model are proposed,one with equal △lnσ and the other in terms of Tschebyscheff polynomials.It is proved that in adiabatic and inviscid cases if the meteorological elements and related physical quantities are continuous in time and in the horizontal,the total energy and total mass are conserved within a high approximation respectively, and there is a correct conversion between total kinetic and total potential energy.Numerical computations show that the schemes both have high accuracy.For example,in integrating the hydrostatic equation the computational errors of geopotential height resulting from the schemes are much less than those resulting from EC79 in a-coordinate.
文摘Through analysis and numerical computation of ECMWF's discrete scheme of hydrostatic equation (Baede et al. 1979),it has been found that in the case of equal △σ there exist systematic errors in the scheme.The error E caused by taking the arithmetic mean of the geopotential heights of two adjacent half a-levels as the geopotential height of the cor- responding integer a-level,increases with height and has an unacceptable maximum in the vicinity of the top of the at- mosphere;however,the errors caused by the temperature treatment are generally small.On the other hand,if an uneven △σ-scheme in which the levels in the upper and lower atmosphere are denser than those in the middle atmosphere,is adopted,then E can be much reduced.However,if the resolution of the original equal Art-scheme doubles,then E can only be found to be much reduced in the troposphere and that in the vicinity of the atmospheric top is still unacceptable. Two numerical schemes for improvement have been presented.Of them one is the same as the ECMWF's scheme, but with equal △lnσ,and the other is to integrate the equation by the use of Tschebyscheff polynomials T and to adopt the zeros of TN as the atmospheric levels where N is the total number of levels.The results show that with both schemes the computational errors can be much reduced,especially the latter,due to the fact that the errors of three statistical types are generally less than the root mean square error of the geopotential heights reported in TEMP.
