*
*  Read in romanian_surf*.inp and define geometry
*

* this is the box around the problem
* make it a little bigger than your data
* the cmo name must match the main input file
cmo / select / cmo1
surface / box / reflect / box / &
  740480. 313980. -200. /  747389. 321889. 50.

*
* read each surface and create an interface
* surface 3 is the water table
*
read / avs / romanian_surf001.inp / cmotmp / 1 1 0
 surface/s1/intrface/sheet/cmotmp
 cmo/release/cmotmp
 cmo/select/cmo1
read / avs / romanian_surf002.inp / cmotmp / 1 1 0
 surface/s2/intrface/sheet/cmotmp
 cmo/release/cmotmp
 cmo/select/cmo1
read / avs / romanian_surf003.inp / cmotmp / 1 1 0
 surface/s3/intrface/sheet/cmotmp
 cmo/release/cmotmp
 cmo/select/cmo1

*
* create regions and material regions for geometry
* mat points are less and equal to their upper surface
*
region/ F / le box and ( lt s1 and lt s2 and le s3 ) 
region/ E / le box and ( ge s1 and lt s2 and le s3 ) 
region/ D / le box and ( gt s1 and ge s2 and le s3 ) 
region/ abvh2o / le box and ( gt s3 ) 
* no materials at interfaces
mregion/ F / le box and ( lt s1 and lt s2 and lt s3 )
mregion/ E / le box and ( gt s1 and lt s2 and lt s3 ) 
mregion/ D / le box and ( gt s1 and gt s2 and lt s3 )
mregion/ abvh2o / le box and ( gt s3 )

finish