Modeling issues of infinite dimensional systems is studied in this paper. Although the modeling problem has been solved to some extent, use of decomposition techniques still pose several difficulties. A prime one of this is the amount of data to be processed. Method of snapshots integrated with POD is a remedy. The second difficulty is the fact that the decomposition followed by a projection yields an autonomous set of finite dimensional ODEs that is not useful for developing a concise understanding of the input operator of the system. A numerical approach to handle this issue is presented in this paper. As the example, we study 2D heat flow problem. The results obtained confirm the theoretical claims of the paper and emphasize that the technique presented here is not only applicable to infinite dimensional linear systems but also to nonlinear ones