TIR.NEW can be read by, open(10,file='TIR.TBL',status='old') do 100 n=1,1300 read(10,1000) BT, (tbl(n,nband),nband=1,4) 100 continue 1000 format(e12.7,1x,e12.7,1x,e12.7,1x,e12.7,1x,e12.7) close(10) In the LUT, the radiance value (RD) are difined as, RD[W m^-2 sr^-1]=L1BDNvalue*gain(nband)*10. And cange RD to BT by the LUT.(we use linear interpolation between the values of lines.) Details of the algorithm are written in the following reforence. http://www.eorc.nasda.go.jp/ADEOS/UsersGuide/Doc_and_Paper/final_repo.html ----------------------------------------------------------- The Version 4 SST coefficients are , MCSST = -29.760851 + 1.1126003 T11 + 4.2436047 (T11-T12)ave - 0.6637208 (T11-T10) + 0.6855296 (T11-T12)ave *(m-1) - 0.3704848 (T11-T10) * (m-1) T10,T11,T12 : brightness temperatures of Band 10,11,12. (T11-T12)ave : 20pixel x 20line spatial moving-average centered each grid (for noise reduction). Brigitness temparatures are converted from L1B "bandwidth-integrated" radiance using the table TIR.NEW. See below. Sakaida, F., M. Moriyama, H. Murakami, H. Oaku, Y. Mitomi, A Mukaida and H. Kawamura, The sea surface temperature product algorithm of the Ocean Color and Temperature Scanner (OCTS) and its accuracy, J. Oceanogr., 54, 437-442, 1998.