read / gmv / bm5_2d_uniform.gmv / cmo_all
*
* Refine around the tip of the fault x = 2000. , y = -1600.
*
pset/prefine/geom/rtp/1 0 0/0.,0.,0./300.,360.,360./2000.,-1600.,0.

refine/rivara///edge/pset get prefine/100.///inclusive
resetpts / itp
recon 0
pset/prefine/geom/rtp/1 0 0/0.,0.,0./350.,360.,360./2000.,-1600.,0.
smooth / position / esug / pset get prefine
recon 0
pset/prefine/geom/rtp/1 0 0/0.,0.,0./350.,360.,360./2000.,-1600.,0.
smooth / position / esug / pset get prefine

dump / gmv / output_rivara100.gmv / cmo_all

pset/prefine/geom/rtp/1 0 0/0.,0.,0./150.,360.,360./2000.,-1600.,0.

refine/rivara///edge/pset get prefine/60.///inclusive
resetpts / itp
recon 0
pset/prefine/geom/rtp/1 0 0/0.,0.,0./170.,360.,360./2000.,-1600.,0.
smooth / position / esug / pset get prefine
recon 0
pset/prefine/geom/rtp/1 0 0/0.,0.,0./170.,360.,360./2000.,-1600.,0.
smooth / position / esug / pset get prefine
recon 0
dump / gmv / output_rivara060.gmv / cmo_all

pset/prefine/geom/rtp/1 0 0/0.,0.,0./ 90.,360.,360./2000.,-1600.,0.
refine/rivara///edge/pset get prefine/40.///inclusive
resetpts / itp
recon 0
pset/prefine/geom/rtp/1 0 0/0.,0.,0./100.,360.,360./2000.,-1600.,0.
smooth / position / esug / pset get prefine
recon 0
pset/prefine/geom/rtp/1 0 0/0.,0.,0./100.,360.,360./2000.,-1600.,0.
smooth / position / esug / pset get prefine
recon 1
dump / gmv / output_rivara040.gmv / cmo_all

finish