S6302019

"S DEFlNfTiONS { parametęr (fefińWonś }

| Hs=18 ((m) odległość od stropu warstwy nieprzepuszczalne)' do swobodnego I te=2.3o-6 !(m/s) współczynnik filtracji | ćett=l 5*60 !($) delta t I deł*=8.5 ?<m) óetiax n=0.35 ! porowatość .

I INITiAt VALUES

ECJUATłONS {PDE‘s, one for each yariabłe > h: df(łi>=Hs*K^łde®(deb(*n)*dMgr3d<h)) {one possibility)

!h: div(0rad^j)}“O I

IIICONSTRAINTS {Integra! constraints}

11 BOUNOARIES {The domaln definition}

I REGION 1    {For each materia! region}

! 8TART(-10,10) {Walk the domain boundaiy} value(h>=10 LINĘ TO (-4.10) naturai(h)=o linę to (-4,6.5) value(h)=y linę TO (-4,6)

natural(h)=0 fineto (-4,-0.5) to (-3.98,-0.5) TO (-3.98,0) naturalno linę (0,0)

natural(h)=0 linę to (O, -6.5) to (-10. -6.5) to close

jf! [{value(h)=10 UNE TO (-4,10)

I j nalural(h)=0 linę to (-4,-0 5) TO {-3,98,-0.5) TO (-3.98.C0 Hj natural(h)=0 linę (0,0)

H natural(h)=0 linę to (0. -6.5) to (-10, -6.5) to close ■ ^}

li TIME 0 TO 60*3600 {tf time dependent}

11MONITOR3 {show progress)

i

U |

|^:j for 1=1 by delt to 4*delt by 3600 to 60*3600 fContour(h) painted

J? PLOTS {save result displays}

for t=1 by delt to 4*dełt by 12*3600 to 60*3600