We develop a two-dimensional model for numerical simulation of thermal
convection in the Martian lower atmosphere.
Its dynamical core is based on that of the cumulus ensemble model
developed by Nakajima (1994), but
the surface boundary layer, the radiative transfer and the grand
surface thermal model are newly implemented in the model.
The source code of the model can be seen Nakajika and Odaka (2000).
Detail mathematical expressions of the model and its finite
differenced form are shown in appendix A
and appendix B.
