Mesh tally function of monte Carlo method can give a detailed and intensive calculation of flux distribution in specific volumes. To realize such function in JMCT mesh tally function are designed and realized. It supports non-uniform mesh in three kinds of orthogonal geometry (xyz of rectangular coordinates, rθz of cylindrical coordinates, and rθφ of spherical coordinates). Algorithm for rectangular coordinates is discussed. Calculation on DAYAWAN reactor core pin-by-pin model, Venus benchmark model and a 1-D ITER model verifies preliminarily correctness of JMCT mesh tally. Furthermore, U-array benchmark model is used to test serial performance of JMCT mesh tally. Both JMCT and MCNP5 use same xyz mesh grids and run under same condition. It shows that JMCT takes less time consuming and has higher performance dealing with xyz geometry.