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