Today

more code study to realize what inputs are needed:

• ipart: particle initial data. Can be left empty and just add particles. [need x, y, velocities, water depth?, area, smoothing length, flag (default 1)]
• ibottom: can also add each, need: x, y, bottom elevation, smoothing length[ (dx+dy)/2 * coef[=1.2] ], manning coefficient (friction-related, see manning's formula on wikipedia), area (dx*dy)
• ipartv; virtual particles, responsible for closed boundaries: x,y, normalized vector into domain [x,y], spacing
• iopenbc: in/outflow, x/y-start, length, unit normal x/y, entering/exiting x-y plane flag???, spacing, steady, time, water depth, vel in time..
• area_lim; refinement parameters. Specify a grid in xy and specify for each cell maximum area allowed before splitting (scaled by max overall area?) default set to 1.E18
  (basically, this is obsolete: if no refinement, can just skip reading this file)

Rsesults so far:

• Parameters v1 demoed, and works also in solver.
• Plan for next part: Implement structures for Bottom and OpenBC particles, and add manually using

1. addParticle(...) methods
2. possibly SWE-SPHysics input files for defaults and testing.

ALSO:

• output from sphysics originally actually is surface elevation (not water depth!)
• ...=dw + h_t
• ... :arrow_right: could output this myself as well to compare
• ...anyway, h_t is calculated from the elevation variably thingy of the bottom particles, HENCE, the bottom particle variably thingy is an elevation-type variable (increases with upwards position)