In MEMO, initialization is performed with suitable diagnostic methods: a mass-consistent initial wind field is formulated using an objective analysis model and scalar fields are initialized using appropriate interpolating techniques (Kunz, R., 1991). Data needed to apply the diagnostic methods may be derived either from observations or from larger scale simulations.
Suitable boundary conditions have to be imposed for the wind velocity components ,  and , the potential temperature  and pressure  at all boundaries. At open boundaries, wave reflection and deformation may be minimized by the use of so-called radiation conditions (Orlanski 1976). 
According to the experience gained so far with the model MEMO, neglecting large scale environmental information might result in instabilities in case of simulations over longer time periods.
For the nonhydrostatic part of the mesoscale pressure perturbation, homogeneous Neumann boundary conditions are used at lateral boundaries. With these conditions, the wind velocity component perpendicular to the boundary remains unaffected by the  pressure change.
At the upper boundary, Neumann boundary conditions are imposed for the horizontal velocity components and the potential temperature. To ensure non-reflectivity, a radiative condition is used for the hydrostatic part of the mesoscale pressure perturbation 
 at that boundary. Hence, vertically propagating internal gravity waves are allowed to leave the computational domain (Klemp and Durran 1983). For the nonhydrostatic part of the mesoscale pressure perturbation, homogeneous staggered Dirichlet conditions are imposed. Being justified by the fact that nonhydrostatic effects are negligible at large heights, this condition is necessary, if singularity of the elliptic pressure equation is to be avoided in view of the Neumann boundary conditions at all other boundaries.
The lower boundary coincides with the ground (or, more precisely, a height above ground corresponding to its aerodynamic roughness). For the non-hydrostatic part of the mesoscale pressure perturbation, inhomogeneous Neumann conditions are imposed at that boundary. All other conditions at the lower boundary follow from the assumption that the –Obukhov similarity theory is valid.
The one way interactive nesting facility is possible within MEMO. Thus, successive simulations on grids of increasing resolution are possible. During these simulations, the results of the application to a coarse grid are used as boundary conditions for the application to the finer grid (Kunz and Moussiopoulos, 1995).