I'm trying to adapt a POSDAT routine to account for multiple materials and am having some difficulty understanding how to use the reference temperature(s).

I'm using PREFM(IMAT) as a reference pressure for all the cells, but I can't work out whether I should be doing something similar for TREF. There is a TREFM() array, but it doesn't seem available to POSDAT.

In the nom.inc file, T(I,1) is described as relative temperature in cell ICLMAP(I), relative to TREFM(IMAT). But POSDAT doesn't get TREFM passed.

Anyone know what I should do to get the absolute temperature in a given cell with a given material number?

Hi Steve,

You can find the common blocks available by looking in starin.f (which is in the solve directory of your installation). Do a search for TREFM in starin.f and you will get the common block it is stored in and it's dimensions.

Bye.

Thanks for the info/advice. Yet another undocumented bit of code. I hope it's what I'm after and not going to bite me one day!