I figured out the problem with my fortran code reading the lake output. I should have specified double precision (real*8). Now it seems to work OK. Does anyone know the difference between "hi" and "hice", from the lake output? I previously emailed this question but did not receive a reply. Why is ice thickness? Thanks, Mike