-Yi ' tFLITłłll IUWJ p^fllBITIfWnifWTniWTT
I Hs=i g i(n1) odległość od stropu warstwy nieprzepuszczalnej do swobodnego | |ę=2.3e-6 l(m/s) współczynnik filtracji | delt=15*60 Ke) delta t 1 debe=0.5 !(m) delta x 1 n=0.351 porowatość
( I INfTIAL VALUES
| EOUATIONS {PDE'$, one for each variable}
] h: dt(h>=Hs*k*deltf(de«*n)*dl¥(grad(h)) {ona possibility}
Ilh: div(grad(h))=0
: ICONSTRAINTS {Integra! constraints >
BOUNDARIES (The domain definrtion}
REGION 1 {For each materiał region}
START(-10.10) {Walk the domain boundaiy} value(h)=10 LINĘ TO (-4.10) natural(h)=0 linę to (-4.6.5) vatue(h)=y linę TO (-4.6)
naturai(h)=0 linę to (-4,-0.5) to (-3.98,-0.5) TO (-3.98,0) natural(h)=0 linę (0,0)
natural(h>=0 linę to (0. -6.5) to (-10, -6.5) to close
j] (value(h)=10 LINĘ TO (-4,10)
Ij natuial(h)=0 linę to (-4.-0.5) TO (-3.98,-0.5) TO (-3.98,0) I natural(h)=0 linę (0,0)
I natural(h)=0 linę to (0, -6.5) to (-10, -6.5) to close
I. TIME OTO60*3600 {iflimę dependent}
I MONITORS {show progress)
I for 1=1 by delt to 4*delt by 3600 to 60*3600 Contour(h) painted
I; PLOTS {save result displays)
I fort=1 by delt to 4*dełt by 12*3600 to 60*3600 Contour(h) painted
END