#
# Create some node attributes to hold elevation and isopach
#
cmo / addatt / MONAME / z_aq_bot / vdouble / scalar / nnodes
cmo / addatt / MONAME / z_aq_top / vdouble / scalar / nnodes
cmo / addatt / MONAME / z_cu_bot / vdouble / scalar / nnodes
cmo / addatt / MONAME / z_cu_top / vdouble / scalar / nnodes
cmo / addatt / MONAME / z_aq_isopach / vdouble / scalar / nnodes
cmo / addatt / MONAME / z_cu_isopach / vdouble / scalar / nnodes
#
# Copy zic surface elevations into node attributes
#
cmo / copyatt / MONAME / M1 / z_aq_bot / zic
cmo / copyatt / MONAME / M2 / z_aq_top / zic
cmo / copyatt / MONAME / M2 / z_cu_bot / zic
cmo / copyatt / MONAME / M3 / z_cu_top / zic
#
# Compute isopach as z_top - z_bottom
#
math / sub / MONAME / z_aq_isopach / 1 0 0 / &
             MONAME / z_aq_top / MONAME / z_aq_bot
math / sub / MONAME / z_cu_isopach / 1 0 0 / &
             MONAME / z_cu_top / MONAME / z_cu_bot
#
# Output some statistics to the screen
#
cmo / printatt / MONAME / z_aq_isopach / minmax
cmo / printatt / MONAME / z_cu_isopach / minmax

finish