Nowadays, the main difficulties encountered in cryosurgery are the exact position of freezing contour, the depth and extent of tissue freezing, and temperature distribution in tissue. In this study, for the first time, superficial cryosurgery treatment has been modeled by finite element method using macros. In this modeling, the tissue is treated as nonideal materials freezing over a temperature range, and the thermophysical properties of which are temperature dependent. The enthalpy method is applied to solve this highly nonlinear problem. After solving this problem for superficial cryoprobe, good agreement was seen in penetration depth of freezing in tissue and the extent of freezing on tissue between our results and the results of other studies. This problem was solved for a new probe by applying cooling and heating loads simultaneously. According to calculated results, for the new probe, the depth of freezing and the quantity of freezing degrease to 1mmand 20%, respectively.