* Use lagrit executable
* lagrit < 02_grid_colors.in
* reads grid.inp created from gridder program

define GRID_NAME grid.inp

define AVS_NAME1 brick_clr.inp
define GMV_NAME1 brick_clr.gmv
define AVS_NAME2 brick_rmtop.inp
define GMV_NAME2 brick_rmtop.gmv

define MINX 740487.
define MINY 313988.
define MINZ -180.
define MAXX 747387.
define MAXY 321888.
define MAXZ 20.

define TOP 4

*
*cmo / create / cmo1 / / / hex
*cmo / select / cmo1
*rzbrick / xyz / NX NY NZ / MINX MINY MINZ / MAXX,MAXY,MAXZ 

read avs GRID_NAME cmo1
  cmo setatt cmo1 imt 1
  resetpts/itp

* color material using the surface geology 
infile 02infile_geom.in 
setpts/no_interface
settets/geometry
cmo printatt cmo1 itetclr minmax
cmo printatt cmo1 imt minmax

dump avs AVS_NAME1 cmo1
dump gmv GMV_NAME1 cmo1


* Remove top above water table 
rmmat/ TOP /element
rmpoint element 
rmpoint compress
resetpts itp

cmo printatt cmo1 itetclr minmax
cmo printatt cmo1 imt1 minmax
resetpts/cell_color/4
cmo printatt cmo1 imt1 minmax

dump avs AVS_NAME2 cmo1
dump gmv GMV_NAME2 cmo1

finish