Increasing number of massive globular clusters (GCs) in the Milky Way are now turned out to host multiple stellar populations having different heavy element abundances enriched by supernovae. Recent observations have further shown that [CNO/Fe] is also enhanced in metal-rich subpopulations in most of these GCs, including ω Cen and M22 (Marino et al. 2011, 2012). In order to reflect this in our population modeling, we have expanded the parameter space of Y 2 isochrones and horizontal-branch (HB) evolutionary tracks to include the cases of normal and enhanced nitrogen abundances ([N/Fe] = 0.0, 0.8, and 1.6). The observed variations in the total CNO content were reproduced by interpolating these nitrogen enhanced stellar models. Our test simulations with varying N and O abundances show that, once the total CNO sum ([CNO/Fe]) is held constant, both N and O have almost identical effects on the HR diagram (see Fig. 1).