The function zroot_watterstrom calls the parser

Since this function is used for curvilinear coordinates, the parser is called for each grid point.