def RequestInformation(self):
from paraview import util
from PVGPpy.read import ubcExtent
# Preview the mesh file and get the mesh extents
ext = ubcExtent(FileName_Mesh)
# Set the mesh extents
util.SetOutputWholeExtent(self, ext)
def RequestInformation():
import datetime
from paraview import util
sys.path.insert(0, r'EMC_SRC_PATH')
import IrisEMC_Paraview_Lib as Lib
Latitude_Begin, Latitude_End, Longitude_Begin, Longitude_End = Lib.get_area(
Area, Latitude_Begin, Latitude_End, Longitude_Begin, Longitude_End)
if not Time_End.strip():
Time_End = datetime.datetime.today().strftime('%Y-%m-%d')
query = Lib.earthquakeQuery % (Time_Begin, Time_End, Magnitude_Begin,
Magnitude_End, Depth_Begin, Depth_End,
Latitude_Begin, Latitude_End,
Longitude_Begin, Longitude_End)
Alternate_FileName = Alternate_FileName.strip()
if len(Alternate_FileName) <= 0:
eq_file = Lib.query2filename(query,
url=Lib.earthquakeKeys[Data_Source])
query = '?'.join([Lib.earthquakeKeys[Data_Source], query])
file_found, address, source = Lib.find_file(
eq_file, loc=r'EMC_EARTHQUAKES_PATH', query=query)
else:
file_found, address, source = Lib.find_file(
Alternate_FileName, loc=r'EMC_EARTHQUAKES_PATH')
if not file_found:
raise Exception(
'earthquake catalog file "' + address +
'" not found! Please provide the full path or UR for the file. Aborting.'
)
(params, lines) = Lib.read_geocsv(address)
num = len(lines)
# ABSOLUTELY NECESSARY FOR THE READER TO WORK:
util.SetOutputWholeExtent(self, [0, num - 1, 0, num - 1, 0, num - 1])
def RequestInformation():
from paraview import util
sys.path.insert(0, "EMC_SRC_PATH")
import IrisEMC_Paraview_Lib as lib
Latitude_Begin, Latitude_End, Longitude_Begin, Longitude_End = lib.getArea(
Area, Latitude_Begin, Latitude_End, Longitude_Begin, Longitude_End)
query = lib.earthquakeQuery % (Start_Time, Magnitude_Begin, Magnitude_End,
Depth_Begin, Depth_End, Latitude_Begin,
Latitude_End, Longitude_Begin,
Longitude_End, Max_Count)
if len(Alternate_FileName) <= 0:
eqFile = lib.query2fileName(query, url=lib.earthquakeKeys[Data_Source])
query = '?'.join([lib.earthquakeKeys[Data_Source], query])
fileFound, address, source = lib.findFile(eqFile,
loc='EMC_EARTHQUAKES_PATH',
query=query)
else:
fileFound, address, source = lib.findFile(Alternate_FileName,
loc='EMC_EARTHQUAKES_PATH')
if not fileFound:
raise Exception(
'earthquake catalog file "' + address +
'" not found! Please provide the full path or UR for the file. Aborting.'
)
(params, lines) = lib.readGcsv(address)
num = len(lines)
# ABSOLUTELY NECESSARY FOR THE READER TO WORK:
util.SetOutputWholeExtent(self, [0, num - 1, 0, num - 1, 0, num - 1])
def RequestInformation():
from paraview import util
import csv
if FileName_Mesh == 'absolute path':
raise Exception('No mesh file selected. Aborting.')
# Read in the mesh and set up structure of output grid from mesh file input
with open(FileName_Mesh) as f:
reader = csv.reader(f, delimiter=Delimiter_Field)
h = reader.next()
# TODO: ignore header lines if start with '!'
# Ignore header lines if start with '!'
while h[0][0] == '!':
h = reader.next()
# Number of CELLS in each axial direction
n1, n2, n3 = int(h[0]), int(h[1]), int(h[2])
# change to number of POINTS in each axial direction
n1, n2, n3 = n1 + 1, n2 + 1, n3 + 1
# ABSOLUTELY NECESSARY FOR THE FILTER TO WORK:
util.SetOutputWholeExtent(self, [0, n1 - 1, 0, n2 - 1, 0, n3 - 1])
f.close()
def RequestInformation():
from paraview import util
from PVGPpy.filt import refoldidx
# Setup the ImageData
idx = refoldidx(SEPlib=SEPlib, swapXY=Transpose_XY)
nx, ny, nz = extent[idx[0]], extent[idx[1]], extent[idx[2]]
# ABSOLUTELY NECESSARY FOR THE FILTER TO WORK:
util.SetOutputWholeExtent(self, [0, nx - 1, 0, ny - 1, 0, nz - 1])
def RequestInformation(self):
from paraview import util
from PVGPpy.read import setOutputTimesteps, readPVGPGridExtents
# This is necessary to set time steps
setOutputTimesteps(self, FileNames[0], dt=Time_Step)
# NOTE: if using time series, they all must have the same extents
ext = readPVGPGridExtents(FileNames[0])
util.SetOutputWholeExtent(self, ext)
def RequestInformation(self, request, inInfo, outInfo):
logger.debug("Requesting information...")
# For the `vtkUniformGrid` output we need to provide extents
# so that it gets rendered at all.
N = self.num_points_per_dim
util.SetOutputWholeExtent(self, [0, N - 1, 0, N - 1, 0, N - 1])
logger.debug(f"Information object: {outInfo.GetInformationObject(0)}")
return 1
def RequestInformation():
from paraview import util
from PVGPpy.read import ubcExtent
if FileName_Mesh == 'absolute path':
raise Exception('No mesh file selected. Aborting.')
# Preview the mesh file and get the mesh extents
ext = ubcExtent(FileName_Mesh)
# Set the mesh extents
util.SetOutputWholeExtent(self, ext)
def RequestInformation():
sys.path.insert(0, r'EMC_SRC_PATH')
import IrisEMC_Paraview_Lib as lib
from paraview import util
from os.path import splitext
import IrisEMC_Paraview_Utils as utils
import IrisEMC_Paraview_Param as param
depth_begin = float(Depth_begin)
depth_end = float(Depth_end)
File_name = File_name.strip()
ext = None
if File_name in list(param.filesDict.values()) or not (
File_name.lower().endswith(param.filesExtDict['ssl'].lower())
or File_name.lower().endswith(param.filesExtDict['geo'].lower())):
if utils.support_nc():
ext = param.filesExtDict['ssl']
else:
ext = param.filesExtDict['geo']
fileFound, address, source = lib.find_file(File_name,
loc=r'EMC_MODELS_PATH',
ext=ext)
if not fileFound:
raise Exception('model file "' + address + '" not found! Aborting.')
Latitude_begin, Latitude_end, Longitude_begin, Longitude_end = lib.get_area(
Area, Latitude_begin, Latitude_end, Longitude_begin, Longitude_end)
this_filename, extension = splitext(address)
if extension.lower() == '.nc':
nx, ny, nz = lib.read_netcdf_model(address,
Latitude_variable,
Longitude_variable,
Depth_variable,
(Latitude_begin, Longitude_begin),
(Latitude_end, Longitude_end),
depth_begin,
depth_end,
Vertical_Scaling,
inc=Sampling,
extent=True)
else:
try:
nx, ny, nz = lib.read_geocsv_model_3d(
address, (Latitude_begin, Longitude_begin),
(Latitude_end, Longitude_end),
depth_begin,
depth_end,
Vertical_Scaling,
inc=Sampling,
extent=True)
except Exception:
raise Exception('cannot recognize model file "' + address +
'"! Aborting.')
# ABSOLUTELY NECESSARY FOR THE READER TO WORK:
util.SetOutputWholeExtent(self, [0, nx, 0, ny, 0, nz])
def RequestInformation():
from paraview import util
# Setup the ImageData
# Cartesian: d1=x, d2=y, d3=z
nx, ny, nz = extent[0], extent[1], extent[2]
if SEPlib:
# SEPlib: d1=z, d2=x, d3=y
nz, nx, ny = extent[0], extent[1], extent[2]
# ABSOLUTELY NECESSARY FOR THE FILTER TO WORK:
util.SetOutputWholeExtent(self, [0, nx - 1, 0, ny - 1, 0, nz - 1])
def RequestInformation():
from paraview import util
from PVGPpy.read import sgemsExtent, setOutputTimesteps, getTimeStepFileIndex
# This is necessary to set time steps
setOutputTimesteps(self, FileNames, dt=Time_Step)
# Only grab extent for first file... requires all to have same extent
ext = sgemsExtent(FileNames[0],
deli=Delimiter_Field,
useTab=Use_tab_delimiter)
util.SetOutputWholeExtent(self, ext)
def RequestInformation():
from paraview import util
import csv
with open(FileName) as f:
reader = csv.reader(f, delimiter=Delimiter_Field)
h = reader.next()
n1, n2, n3 = int(h[0]), int(h[1]), int(h[2])
# ABSOLUTELY NECESSARY FOR THE FILTER TO WORK:
util.SetOutputWholeExtent(self, [0, n1 - 1, 0, n2 - 1, 0, n3 - 1])
f.close()
def RequestInformation():
from os.path import splitext
sys.path.insert(0, r'EMC_SRC_PATH')
import IrisEMC_Paraview_Lib as lib
from paraview import util
import IrisEMC_Paraview_Utils as utils
import IrisEMC_Paraview_Param as param
Latitude_Begin, Latitude_End, Longitude_Begin, Longitude_End = lib.get_area(
Area, Latitude_Begin, Latitude_End, Longitude_Begin, Longitude_End)
if len(Alternate_FileName.strip()) > 0:
FileName = Alternate_FileName.strip()
Label = ' '.join(['SLAB', lib.file_name(Alternate_FileName).strip()])
else:
FileName = lib.usgsSlabKeys[Slab]
Label = ' '.join(['USGS Slab 1.0 -', lib.usgsSlabValues[Slab].strip()])
FileName = FileName.strip()
ext = None
if FileName in list(param.usgsSlabDict.keys()):
if utils.support_nc():
ext = param.usgsSlabExtDict['ssl']
else:
ext = param.usgsSlabExtDict['geo']
fileFound, address, source = lib.find_file(FileName,
loc=r'EMC_SLABS_PATH',
ext=ext)
if not fileFound:
raise Exception('model file "' + address + '" not found! Aborting.')
this_filename, extension = splitext(address)
if extension.lower() == '.grd':
nx, ny, nz = lib.read_slab_file(address,
(Latitude_Begin, Longitude_Begin),
(Latitude_End, Longitude_End),
inc=Sampling,
extent=True)
elif extension.lower() == '.csv':
nx, ny, nz = lib.read_geocsv_model_2d(
address, (Latitude_Begin, Longitude_Begin),
(Latitude_End, Longitude_End),
Sampling,
0,
extent=True)
else:
raise Exception('cannot recognize model file "' + address +
'"! Aborting.')
# ABSOLUTELY NECESSARY FOR THE READER TO WORK:
util.SetOutputWholeExtent(self, [0, nx, 0, ny, 0, nz])
def RequestInformation():
sys.path.insert(0, "EMC_SRC_PATH")
import IrisEMC_Paraview_Lib as lib
from paraview import util
fileFound,address,source = lib.findFile(FileName,loc='EMC_MODELS_PATH')
if not fileFound:
raise Exception('model file "'+address+'" not found! Aborting.')
Latitude_Begin,Latitude_End,Longitude_Begin,Longitude_End = lib.getArea(Area,Latitude_Begin,Latitude_End,Longitude_Begin,Longitude_End)
nx,ny,nz = lib.find2DnetCDFExtent(address,Latitude_Variable,Longitude_Variable,(Latitude_Begin,Longitude_Begin),(Latitude_End,Longitude_End),inc=Sampling)
# ABSOLUTELY NECESSARY FOR THE READER TO WORK:
util.SetOutputWholeExtent(self, [0,nx, 0,ny, 0,nz+1])
def RequestInformation():
sys.path.insert(0, "EMC_SRC_PATH")
import IrisEMC_Paraview_Lib as lib
from paraview import util
fileFound, address, source = lib.findFile(lib.etopo5File,
loc='EMC_MODELS_PATH')
Latitude_Begin, Latitude_End, Longitude_Begin, Longitude_End = lib.getArea(
Area, Latitude_Begin, Latitude_End, Longitude_Begin, Longitude_End)
nx, ny, nz = lib.findTopoExtent(address, (Latitude_Begin, Longitude_Begin),
(Latitude_End, Longitude_End), Sampling)
# ABSOLUTELY NECESSARY FOR THE READER TO WORK:
util.SetOutputWholeExtent(self, [0, nx - 1, 0, ny - 1, 0, nz - 1])
def RequestInformation():
sys.path.insert(0, r'EMC_SRC_PATH')
from os.path import splitext
import IrisEMC_Paraview_Lib as lib
from paraview import util
import IrisEMC_Paraview_Utils as utils
import IrisEMC_Paraview_Param as param
if len(Alternate_FileName.strip()) > 0:
file_name = Alternate_FileName.strip()
else:
file_name = lib.topoKeys[TopoFile]
file_name = file_name.strip()
ext = None
if file_name in list(param.topoDict.keys()):
if utils.support_nc():
ext = param.topoExtDict['ssl']
else:
ext = param.topoExtDict['geo']
fileFound, address, source = lib.find_file(file_name,
loc=r'EMC_MODELS_PATH',
ext=ext)
Latitude_Begin, Latitude_End, Longitude_Begin, Longitude_End = lib.get_area(
Area, Latitude_Begin, Latitude_End, Longitude_Begin, Longitude_End)
this_filename, extension = splitext(address)
if extension.lower() in ['.nc', '.grd']:
nx, ny, nz = lib.read_netcdf_topo_file(
address, (Latitude_Begin, Longitude_Begin),
(Latitude_End, Longitude_End),
Sampling,
Roughness,
lon_var=Longitude_Variable,
lat_var=Latitude_Variable,
elev_var=Elevation_Variable,
extent=True)
else:
try:
nx, ny, nz = lib.read_geocsv_model_2d(
address, (Latitude_Begin, Longitude_Begin),
(Latitude_End, Longitude_End),
Sampling,
Roughness,
extent=True)
except Exception:
raise Exception('cannot recognize model file "' + address +
'"! Aborting.')
# ABSOLUTELY NECESSARY FOR THE READER TO WORK:
util.SetOutputWholeExtent(self, [0, nx, 0, ny, 0, nz])
def RequestInformation():
from paraview import util
sys.path.insert(0, "EMC_SRC_PATH")
import IrisEMC_Paraview_Lib as lib
if len(File_name.strip()) <= 0:
fileFound = False
address = ''
else:
fileFound, address, source = lib.find_file(File_name.strip(),
loc=r'EMC_VOLCANOES_PATH')
if not fileFound:
raise Exception(
'volcano file "' + address +
'" not found! Please provide the full path or UR for the file. Aborting.'
)
(params, lines) = lib.read_geocsv(address)
num = len(lines)
# ABSOLUTELY NECESSARY FOR THE READER TO WORK:
util.SetOutputWholeExtent(self, [0, num - 1, 0, num - 1, 0, num - 1])
def RequestInformation():
from paraview import util
sys.path.insert(0, "EMC_SRC_PATH")
import IrisEMC_Paraview_Lib as lib
query = lib.volcanoLocationsQuery
if len(Alternate_FileName) <= 0:
volcanoFile = lib.query2fileName(query)
query = '?'.join([lib.volcanoLocationsKeys[Data_Source], query])
fileFound, address, source = lib.findFile(volcanoFile,
loc='EMC_VOLCANOES_PATH',
query=query)
else:
fileFound, address, source = lib.findFile(Alternate_FileName,
loc='EMC_VOLCANOES_PATH')
if not fileFound:
raise Exception(
'volcano file "' + address +
'" not found! Please provide the full path or UR for the file. Aborting.'
)
(params, lines) = lib.readGcsv(address)
num = len(lines)
# ABSOLUTELY NECESSARY FOR THE READER TO WORK:
util.SetOutputWholeExtent(self, [0, num - 1, 0, num - 1, 0, num - 1])
def RequestInformation():
sys.path.insert(0, "EMC_SRC_PATH")
import IrisEMC_Paraview_Lib as lib
from paraview import util
Latitude_Begin, Latitude_End, Longitude_Begin, Longitude_End = lib.getArea(
Area, Latitude_Begin, Latitude_End, Longitude_Begin, Longitude_End)
if len(Alternate_FileName.strip()) > 0:
FileName = Alternate_FileName
Label = ' '.join(['SLAB', lib.fileName(Alternate_FileName).strip()])
else:
FileName = lib.usgsSlabKeys[Slab]
Label = ' '.join(['USGS Slab 1.0 -', lib.usgsSlabValues[Slab].strip()])
fileFound, address, source = lib.findFile(FileName, loc='EMC_SLABS_PATH')
if not fileFound:
raise Exception('model file "' + address + '" not found! Aborting.')
nx, ny, nz = lib.getSlabExtent(address, (Latitude_Begin, Longitude_Begin),
(Latitude_End, Longitude_End),
inc=Sampling)
# ABSOLUTELY NECESSARY FOR THE READER TO WORK:
util.SetOutputWholeExtent(self, [0, nx - 1, 0, ny - 1, 0, nz - 1])
def RequestInformation():
from paraview import util
pdi = self.GetInput()
util.SetOutputWholeExtent(
self, [0, NumberOfSpaceBins, 0, NumberOfVeloBins, 0, 0])
def RequestInformation():
from paraview import util
# ABSOLUTELY NECESSARY FOR THE FILTER TO WORK:
nx, ny, nz = extent[0], extent[1], extent[2]
util.SetOutputWholeExtent(self, [0, nx - 1, 0, ny - 1, 0, nz - 1])
