elif opt_dict['a'] is not None:
Amps = [opt_dict['a']]
else:
Amps = [PPU.params['Amplitude']]
for iq, Q in enumerate(Qs):
for ik, K in enumerate(Ks):
dirname = "Q%1.2fK%1.2f" % (Q, K)
print "Working in {} directory".format(dirname)
fzs, lvec, nDim = GU.load_scal_field(dirname + '/OutFz',
data_format=data_format)
dfs = PPU.Fz2df(fzs,
dz=dz,
k0=PPU.params['kCantilever'],
f0=PPU.params['f0Cantilever'],
n=Amp / dz)
# print "TYT", fzs.shape
for p in options.points:
xmin = float(p[0].split('x')[0])
ymin = float(p[0].split('x')[1])
zmin = float(p[0].split('x')[2])
xmax = float(p[1].split('x')[0])
ymax = float(p[1].split('x')[1])
zmax = float(p[1].split('x')[2])
npoints = float(p[2])
print opt_dict['disp']
if opt_dict['disp']:
'./OutFzTip_base', data_format=data_format)
except:
print("error: ", sys.exc_info())
print("cannot load : ",
'./OutFzTip_base.' + data_format)
for iA, Amp in enumerate(Amps):
for iT, TbQ in enumerate(TbQs):
if (TbQ == 0.0):
AmpStr = "/Amp%2.2f" % Amp
print("Amp= ", AmpStr)
dirNameAmp = dirname + AmpStr
if not os.path.exists(dirNameAmp):
os.makedirs(dirNameAmp)
dfs = PPU.Fz2df(fzs,
dz=dz,
k0=PPU.params['kCantilever'],
f0=PPU.params['f0Cantilever'],
n=int(Amp / dz))
else:
AmpStr = "/Amp%2.2f_qTip%2.2f" % (Amp, TbQ)
print("Amp= ", AmpStr)
dirNameAmp = dirname + AmpStr
if not os.path.exists(dirNameAmp):
os.makedirs(dirNameAmp)
dfs = PPU.Fz2df(fzs + TbQ * fzt,
dz=dz,
k0=PPU.params['kCantilever'],
f0=PPU.params['f0Cantilever'],
n=int(Amp / dz))
if opt_dict['save_df']:
GU.save_scal_field(dirNameAmp + '/df',
