S6303313

180**5 1!

Ć0ORDINATES cattesian2 (coordinale system. 1Ó.2D.30. etc}

3383. |

1487 I 7139 | 1 391 e-4 | 4.386e-4 || 4.7029-7 1 6.1180-7 3

j;oqtours=8 erriim=i e-3

Hs=16;S|t(m) odległość od stropu warstwy nieprzepuszczalnej

k=2.3e-6 KnVs) współczynnik ONrtcjt

delt=60!(s) delta i

dete=0 5 !(m) delta*

n=0.35! porowatość

l

do swobodnego zwierciadła wody

INITIAL VALUĘS h=10

I EOUATłONS    {PDE‘s, one fot each rariabfe }

4 h: n*dt(h)=HsVdrv(grad(h)) (one possibility}

|!h: div(grad(h))=0

! C0NSTRA1NTS {Integral constraints}

BOUNDARtES {The domain definition)

REGION 1    {For each materia! region )

START(-10,10) {Walk the domain boundaiy} value(h)=10 UNETO(-4,10)

«naturał(h>=0 linę to (>4,6.5) value(h)=y linę TO (-4,6) natura!(h)=0 lin* to (-4.-0.5) to (-3.98.--0.5) TO (-3.98,0) natural(h)=0 ime (0,0)

‘1    natura!(fi)=0 linę to (0, *6.5) to (-10. -6.5) to close

I (value(h)=10 UNETO(-4,1Q)

I naturai(h)=0 linę to (-4,-0.5) TO (-3.98,-0.5) TO (-3.98,0)