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',