<SYSTEM_TASK:>
Given the following text problem statement, write Python code to implement the functionality described below in problem statement
<END_TASK>
<USER_TASK:>
Problem:
How does one convert a list of Z-scores from the Z-distribution (standard normal distribution, Gaussian distribution) to left-tailed p-values? I have yet to find the magical function in Scipy's stats module to do this, but one must be there.
<ASSISTANT_TASK:>
Python Code:
import numpy as np
import scipy.stats
z_scores = np.array([-3, -2, 0, 2, 2.5])
temp = np.array(z_scores)
p_values = scipy.stats.norm.cdf(temp)
<END_TASK>