Dear mitgcm, we init temperature for every grid point in a three dimensional arrey in the file hydrogThetaFile='lev_t.bin', why we also set tRef = 15*20 or some other tref? What's the tRef function? Best regards. Brainf brainf at 163.com 2006-08-16