Why check for normality of data in a sample? from what I understand the assumption of normality (that must be assumed if one wants to use parametric tests) refers to the sampling distribution of the mean. It does not imply that the distriubtion is normal in the population itself nor a in a given sample. With the CLT we can assume that with a large enough sample this assumption is met:

given random and independent samples of N observations each, the distribution of sample means
approaches normality as the size of N increases, regardless of the shape of the population
distribution

Therefore (from my understanding):
If we have a large sample ( I know the N > 30 is debatable, but say the N = 500 simple random sample) OR the variable of interest form population that we are sampling from is known to be normal, then the sampling distribution has to be normal and therefore we can use parametric tests. For example I stumbled across this in one of the scientific journals:

With large enough sample sizes (> 30 or 40), the violation of the normality assumption should not cause major problems (4); this implies that we can use parametric procedures even when the data are not normally distributed (8). If we have samples consisting of hundreds of observations, we can ignore the distribution of the data (3). According to the central limit theorem, (a) if the sample data are approximately normal then the sampling distribution too will be normal; (b) in large samples (> 30 or 40), the sampling distribution tends to be normal, regardless of the shape of the data (2, 8); and (c) means of random samples from any distribution will themselves have normal distribution (3).
Link to the article

Going back to my question, if the normality assumption DOES NOT imply that the variable of interest in a sample is normaly distributed (only that the sampling distriubtion of the mean is normal). why do people plot histograms to see if the data in a sample is normally distributed? How is this relevant for the assumption of normality?
 A: Some of the trouble is that the central limit theorem is that the convergence is a limit, not realized in finite sample sizes...
...except that it is realized in finite sample sizes if the data have a normal distribution. In fact, the distribution of the sample mean is normal if the sample size is just $1$.
Therefore, if you have normal data, you immediately have a normal sampling distribution.
Aother reason is that we sometimes do inference about parameters other than means. Another parameter we might want to test is the variance. Let's do a simulation to test the equality of variances for two equal exponential distributions.
library(ggplot2)
set.seed(2021)
B <- 1000
N <- 300
mean_ps <- var_ps <- rep(NA, B)
for (i in 1:B){
  
  x <- rexp(N, 1)
  y <- rexp(N, 1)
  mean_ps[i] <- t.test(x, y)$p.value
  var_ps[i] <- var.test(x, y)$p.value
}
s <- seq(0, 1, 0.001)
d1 <- data.frame(
  p = s,
  Quantile = ecdf(mean_ps)(s),
  Test = "Mean"
)
d2 <- data.frame(
  p = s,
  Quantile = ecdf(var_ps)(s),
  Test = "Variance"
)
d <- rbind(d1, d2)
ggplot(d, aes(x = p, y = Quantile, col = Test)) +
  geom_point() +
  geom_abline(slope = 1, intercept = 0) + 
  theme_bw()


The F-test of variance equality rejects the true null hypothesis much more often that the t-test of mean equality ($35\%$ rejections at $\alpha = 5\%$), which behaves approximately correctly (uniform distribution of p-values).
Finally, it's easy to check the data. To check the distribution of sample means, we need multiple sample means. There are ways of approximating such a distribution (e.g., bootstrap), but we don't have that distribution "for free" when we collect the data.
A: To quote a source:

Said another way, CLT is a statistical theory stating that given a sufficiently large sample size from a population with a finite level of variance, the mean of all samples from the same population will be approximately equal to the mean of the population. Furthermore, all the samples will follow an approximate normal distribution pattern, with all variances being approximately equal to the variance of the population, divided by each sample's size.

So, per your question: "why do people plot histograms to see if the data in a sample is normally distributed?", one possible answer is that we are actually checking for "an approximate normal distribution pattern", per the source above.
