Earth
Lodowce rzeźbią doliny
Wypełnij rzeczną dolinę w kształcie V lodowcem grubym na sześćset metrów i patrz, jak mija sto tysięcy lat. Lód płynie najszybciej w środku powierzchni, a wcale przy łożu; ścinanie, które ciągnie po łożu, jest największe na zboczach V i znika w jego wąskim dnie, więc ślizganie i abrazja ścinają zbocza, a dolina poszerza się ku U – stare V zostaje narysowane dla porównania. Czysty lód żłobi dużo mniej niż lód pełen skał; bez ślizgania nie żłobi nic; pozwól lodowcowi nacierać i cofać się w cyklach. Każdy metr sześcienny utracony przez łoże trafia do moreny. W kanale półkolistym wynik to zamknięte rozwiązanie Poissona.

Co się dzieje
Przekrój doliny szerokiej na 2400 m na siatce 61 na 40. Prędkość lodu newtonowskiego w dół doliny (10¹³ Pa·s, spadek 0,05) rozwiązuje eta grad²u = −rho g sin alpha z u = 0 na łożu i du/dz = 0 na powierzchni; nieregularne łoże obsługują różnice Shortley-Weller, a rozwiązanie nadrelaksacja od poprzedniego. Ścinanie na łożu napędza ślizganie proporcjonalnie – deklarowane uproszczenie, ślizganie nie wchodzi do rozwiązania – a ślizganie ściera łoże (Hallet), skalowane przez gruz w lodzie i wytrzymałość skały. Każdy metr sześcienny trafia do moreny. Kroki 1 tysiąca lat; sekunda to tysiąc lat przy skali czasu 1.
Co wypróbować
Patrz, jak V poszerza się ku U, i odczytaj wykładnik profilu; porównaj lód czysty i brudny; zatrzymaj ślizganie; pozwól lodowcowi nacierać i cofać się w cyklach; sprawdź półkole ze wzorem zamkniętym.
Model i jego ograniczenia
Przekrój bez przepływu wzdłuż doliny, lód newtonowski, ślizganie proporcjonalne do obliczonego ścinania, jedno prawo abrazji; detrakcja i pękanie skał nie są modelowane. Profile należą do tej wyidealizowanej doliny.
Równania modułu
Równania η·∇²u = −ρ·g·sin α, u = 0 on the bed, ∂u/∂z = 0 at the surface τ_b = η·|∇u| at the bed, u_b = c_s·τ_b dz_b/dt = −K·D/S·u_b·√(1 + (∂z_b/∂y)²) semicircle: u = ρ·g·sin α·(R² − r²)/(4η), τ_b = ρ·g·sin α·R/2


