$*************************    GENERAL   ***************************************
PROJ 'GW' 'T1'
SET NAUT
SET LEVEL 4.5
SET MAXERR 1
SET DEPMIN 0.01
MODE NONST ONED 
$************************     GRIDS     ***************************************
CGRID REG -1450 0 0 1750 0 1750 0 
INPGRID BOT REG -1450 0 0 1750 0 1 0   $x0 y0 theta nx-1 ny-1 dx dy  
VERT 4
READ BOTTOM -1 'bathy.dat' 3 0 FREE

INPGRID PORO REG -1450 0 0 1750 0 1 0
INPGRID PSIZ REG -1450 0 0 1750 0 1 0
INPGRID HSTRUC REG  -1450 0 0 1750 0 1 0

READINP PORO 1 'poro.dat' 3 0 FREE
READINP PSIZ 1 'psize.dat' 3 0 FREE
READINP HSTRUC 1 'hstru.dat' 3 0 FREE


$***********************   BOUNDARIES  ****************************************
INIT ZERO  
BOUN SIDE W BTYPE WEAK SMOOT 10 SEC ADDBOUNDWAVE CON SPECFILE 'spec.dat'
SPON E 250


$***********************   PHYSICS  *******************************************
BRE 0.6 0.2
FRICTION MANN 0.08
POROsity 4 200 1.1 13
VISC V KEPS

$***********************   NUMERICS  ******************************************
NONHYD BOX 1.0 PREC ILU
DISCRET UPW MOM
TIMEI 0.2 0.6

$***********************   OUTPUTS  *******************************************
SET OUTLEV 1 
BLOCK 'COMPGRID' NOHEAD 'out.mat' LAY 3 TSEC XP DEP BOTL WATL PRESS NHPRSK PRESSK VEL VZ VELK ZK OUTPUT 000000.00 0.25 SEC
COMPUTE 000000.000 0.015 SEC 000100.000 $timeini dtini unit timefin
STOP