*
* Build 3 tet meshes at 1000m, 500m and 250m resolution
* for SCEC strike-slip benchmark problem.
*
* Carl Gable
* gable@lanl.gov
*
*
* Define variables that will be used in point distribution.
* These define the min/max extents of the box.
*
define / XMIN / 0.0
define / XMAX / 24000.
define / YMIN / 0.0
define / YMAX / 24000.
define / ZMIN / -24000.
define / ZMAX / 0.0
*
* Build mesh with 1000m point spacing
*
define / CMO_NAME / cmo1_tet
define / OUTPUT_table /  tet_1000m_res.table
define / OUTPUT_gmv /    tet_1000m_res.gmv
define / OUTPUT_avs /    tet_1000m_res.inp
define / OUTPUT_zone /   tet_1000m_res
cmo / create / CMO_NAME / / / tet
createpts / xyz / 25 25 25 / XMIN YMIN ZMIN / XMAX YMAX ZMAX / 1 1 1
cmo / setatt / CMO_NAME / imt / 1 0 0 / 1
cmo / setatt / CMO_NAME / itp / 1 0 0 / 0
connect / noadd
resetpts / itp
infile / identify_fault_nodes.mlgi
cmo / delete / CMO_NAME
*
* Build mesh with 500m point spacing
*
define / CMO_NAME / cmo2_tet
define / OUTPUT_table /  tet_0500m_res.table
define / OUTPUT_gmv /    tet_0500m_res.gmv
define / OUTPUT_avs /    tet_0500m_res.inp
define / OUTPUT_zone /   tet_0500m_res
cmo / create / CMO_NAME / / / tet
createpts / xyz / 49 49 49 / XMIN YMIN ZMIN / XMAX YMAX ZMAX / 1 1 1
cmo / setatt / CMO_NAME / imt / 1 0 0 / 1
cmo / setatt / CMO_NAME / itp / 1 0 0 / 0
connect / noadd
resetpts / itp
infile / identify_fault_nodes.mlgi
cmo / delete / CMO_NAME
*
* Build mesh with 250m point spacing
*
define / CMO_NAME / cmo3_tet
define / OUTPUT_table /  tet_0250m_res.table
define / OUTPUT_gmv /    tet_0250m_res.gmv
define / OUTPUT_avs /    tet_0250m_res.inp
define / OUTPUT_zone /   tet_0250m_res
cmo / create / CMO_NAME / / / tet
createpts / xyz / 97 97 97 / XMIN YMIN ZMIN / XMAX YMAX ZMAX / 1 1 1
cmo / setatt / CMO_NAME / imt / 1 0 0 / 1
cmo / setatt / CMO_NAME / itp / 1 0 0 / 0
connect / noadd
resetpts / itp
infile / identify_fault_nodes.mlgi
cmo / delete / CMO_NAME
*
*
*
finish