Matt, this sounds like an input file format problem. That is, hydrogThetaFile expects formatted FORTRAN records and you provide it with direct access binary input, or vice versa. The zeros that appear in T.0000000000.data are the land mask. Cheers, Dimitris