Hallo Dimitris, is there any reason why GWATX and GWATY are computed in the main routine? These fields are only used in the dynamics part (dynsolver, ostres), and there's a debug statement in growth.F. Can I move the assignment of GWATX/Y into dynsolver? Martin