def run_pgrtrans(self, **kwargs):
# do nothing if output data already exist
if len(np.shape(pgrtrans.ivals)) > 0:
print('ERROR in run_pgrtrans: data already exist!')
return
# set arguments
self.inputs = grtrans_inputs(**kwargs)
# call pgrtrans routine with arguments:
self.print_input_list()
pgrtrans.grtrans_main(
self.inputs.standard, self.inputs.mumin, self.inputs.mumax,
self.inputs.nmu, self.inputs.phi0, self.inputs.spin,
self.inputs.uout, self.inputs.uin, self.inputs.rcut,
self.inputs.nrotype, self.inputs.gridvals, self.inputs.nn,
self.inputs.i1, self.inputs.i2, self.inputs.fname, self.inputs.dt,
self.inputs.nt, self.inputs.nload, self.inputs.nmdot,
self.inputs.mdotmin, self.inputs.mdotmax, self.inputs.ename,
self.inputs.mbh, self.inputs.nfreq, self.inputs.fmin,
self.inputs.fmax, self.inputs.muval, self.inputs.gmin,
self.inputs.gmax, self.inputs.p1, self.inputs.p2,
self.inputs.jetalpha, self.inputs.stype, self.inputs.use_geokerr,
self.inputs.nvals, self.inputs.iname, self.inputs.cflag,
self.inputs.extra, self.inputs.debug, self.inputs.outfile,
self.inputs.fdfile, self.inputs.fhfile, self.inputs.fgfile,
self.inputs.fsim, self.inputs.fnt, self.inputs.findf,
self.inputs.fnfiles, self.inputs.fjonfix, self.inputs.pnw,
self.inputs.pnfreq_tab, self.inputs.pnr, self.inputs.foffset,
self.inputs.fdindf, self.inputs.fmagcrit, self.inputs.hrspot,
self.inputs.hr0spot, self.inputs.hn0spot, self.inputs.ntscl,
self.inputs.nrscl, self.inputs.pwmin, self.inputs.pwmax,
self.inputs.pfmin, self.inputs.pfmax, self.inputs.prmax,
self.inputs.psigt, self.inputs.pfcol, self.inputs.tmdot,
self.inputs.snscl, self.inputs.snnthscl, self.inputs.snnthp,
self.inputs.sbeta, self.inputs.sbl06, self.inputs.snp,
self.inputs.snt, self.inputs.srin, self.inputs.srout,
self.inputs.sthin, self.inputs.sthout, self.inputs.sphiin,
self.inputs.sphiout, self.inputs.fscalefac, self.inputs.sigcut,
self.inputs.epcoefindx, self.inputs.epotherargs, self.inputs.nep)
# read output
self.ivals = pgrtrans.ivals.copy()
self.ab = pgrtrans.ab.copy()
self.nu = pgrtrans.freqs.copy()
self.nx = int(self.inputs.nn[0])
self.ny = int(self.inputs.nn[1])
self.nvals = int(self.inputs.nvals)
self.del_pgrtrans_data()
def run_pgrtrans(self,**kwargs):
# do nothing if output data already exist
if len(np.shape(pgrtrans.ivals)) > 0:
print 'ERROR in run_pgrtrans: data already exist!'
return
# set arguments
self.inputs=grtrans_inputs(**kwargs)
# call pgrtrans routine with arguments:
self.print_input_list()
pgrtrans.grtrans_main(self.inputs.standard,self.inputs.mumin,self.inputs.mumax,self.inputs.nmu,self.inputs.phi0,self.inputs.spin,self.inputs.uout,self.inputs.uin,self.inputs.rcut,self.inputs.nrotype,self.inputs.gridvals,self.inputs.nn,self.inputs.i1,self.inputs.i2,self.inputs.fname,self.inputs.dt,self.inputs.nt,self.inputs.nload,self.inputs.nmdot,self.inputs.mdotmin,self.inputs.mdotmax,self.inputs.ename,self.inputs.mbh,self.inputs.nfreq,self.inputs.fmin,self.inputs.fmax,self.inputs.muval,self.inputs.gmin,self.inputs.gmax,self.inputs.p1,self.inputs.p2,self.inputs.jetalpha,self.inputs.stype,self.inputs.use_geokerr,self.inputs.nvals,self.inputs.iname,self.inputs.cflag,self.inputs.extra,self.inputs.debug,self.inputs.outfile,self.inputs.fdfile,self.inputs.fhfile,self.inputs.fgfile,self.inputs.fsim,self.inputs.fnt,self.inputs.findf,self.inputs.fnfiles,self.inputs.fjonfix,self.inputs.pnw,self.inputs.pnfreq_tab,self.inputs.pnr,self.inputs.foffset,self.inputs.fdindf,self.inputs.fmagcrit,self.inputs.hrspot,self.inputs.hr0spot,self.inputs.hn0spot,self.inputs.ntscl,self.inputs.nrscl,self.inputs.pwmin,self.inputs.pwmax,self.inputs.pfmin,self.inputs.pfmax,self.inputs.prmax,self.inputs.psigt,self.inputs.pfcol,self.inputs.tmdot,self.inputs.snscl,self.inputs.snnthscl,self.inputs.snnthp,self.inputs.sbeta,self.inputs.sbl06,self.inputs.snp,self.inputs.snt,self.inputs.srin,self.inputs.srout,self.inputs.sthin,self.inputs.sthout,self.inputs.sphiin,self.inputs.sphiout,self.inputs.epcoefindx,self.inputs.epotherargs,self.inputs.nep)
# read output
self.ivals = pgrtrans.ivals.copy()
self.ab = pgrtrans.ab.copy()
self.nu = pgrtrans.freqs.copy()
self.nx = self.inputs.nn[0]
self.ny = self.inputs.nn[1]
self.del_pgrtrans_data()
fmagcrit = int(0)
hrspot = float(0)
hr0spot = float(0)
hn0spot = float(0)
ntscl = float(0)
nrscl = float(0)
pwmin = float(0)
pwmax = float(0)
pfmin = float(0)
pfmax = float(0)
prmax = float(0)
psigt = float(0)
pfcol = float(0)
tmdot = float(0)
snscl = float(0)
snnthscl = float(0)
snnthp = float(0)
sbeta = float(0)
pgrtrans.grtrans_main(standard, mumin, mumax, nmu, phi0, spin, uout, uin, rcut,
nrotype, gridvals, nn, fname, dt, nt, nload, nmdot,
mdotmin, mdotmax, ename, mbh, nfreq, fmin, fmax, muval,
gmin, gmax, p1, p2, jetalpha, stype, use_geokerr, nvals,
iname, cflag, extra, ofile, fdfile, fhfile, fgfile, fsim,
fnt, findf, fnfiles, fjonfix, pnw, pnfreq_tab, pnr,
foffset, fdindf, fmagcrit, hrspot, hr0spot, hn0spot,
ntscl, nrscl, pwmin, pwmax, pfmin, pfmax, prmax, psigt,
pfcol, tmdot, snscl, snnthscl, snnthp, sbeta)
pgrtrans.del_pgrtrans_data()
pgrtrans.grtrans_main(
standard,
mumin,
mumax,
nmu,
phi0,
spin,
uout,
uin,
rcut,
nrotype,
gridvals,
nn,
fname,
dt,
nt,
nload,
nmdot,
mdotmin,
mdotmax,
ename,
mbh,
nfreq,
fmin,
fmax,
muval,
gmin,
gmax,
p1,
p2,
jetalpha,
stype,
use_geokerr,
nvals,
iname,
cflag,
extra,
ofile,
fdfile,
fhfile,
fgfile,
fsim,
fnt,
findf,
fnfiles,
fjonfix,
pnw,
pnfreq_tab,
pnr,
foffset,
fdindf,
fmagcrit,
hrspot,
hr0spot,
hn0spot,
ntscl,
nrscl,
pwmin,
pwmax,
pfmin,
pfmax,
prmax,
psigt,
pfcol,
tmdot,
snscl,
snnthscl,
snnthp,
sbeta,
)
