Arne, The problem is not really the loop range but that dyU is not initialized in the first column (look at ini_spherical_grid.F). For this reason we supply recip_dyU(...) so that we canavoid dividing by zero in the overlaps. So replace /dyU with *recip_dyU and that should fix it also. A.