def interpolate_edot(frame, xtrench):
# domain bounds in km
fi.xmin = xtrench + left3
fi.xmax = xtrench + right3
fi.zmin = down3
fi.zmax = up3
# resolution in km
fi.dx = dx
fi.dz = dz
xx, zz, edot = fi.interpolate(frame, 'srII')
return xx, zz, edot
def interpolate_chamber(frame, xtrench):
# domain bounds in km
fi.xmin = xtrench + left
fi.xmax = xtrench + right
fi.zmin = down
fi.zmax = up
# resolution in km
fi.dx = dx
fi.dz = dz
xx, zz, C = fi.interpolate(frame, 'chamber')
return xx, zz, C
def interpolate_dens(frame, xtrench):
# domain bounds in km
fi.xmin = xtrench + left3
fi.xmax = xtrench + right3
fi.zmin = down3
fi.zmax = up3
# resolution in km
fi.dx = dx
fi.dz = dz
xx, zz, dens = fi.interpolate(frame, 'density')
return xx, zz, dens
def interpolate_phase(frame, xtrench, kinds):
if (kinds == 1):
# domain bounds in km
fi.xmin = xtrench + left
fi.xmax = xtrench + right
fi.zmin = down
fi.zmax = up
# resolution in km
fi.dx = dx
fi.dz = dz
xx, zz, ph = fi.interpolate(frame, 'phase')
return xx, zz, ph
def interpolate_sII(frame, xtrench):
# domain bounds in km
fi.xmin = xtrench + left
fi.xmax = xtrench + right
fi.zmin = down
fi.zmax = up
# resolution in km
fi.dx = dx
fi.dz = dz
sx, sz, s = fi.interpolate(frame, 'sII')
return sx, sz, s
