Could the first loop be replaced by a faster way like more matrix operations?

Question

Three are a great amount of calculations in the code below. rlist has about 1000 to 5000 float numbers.

The final goal is to obtain temph001, but I find the calculati