def run1(): config = readInput(fname) tm0 = time.time() tsdat = run(config) write(tsdat, config, outname) print('TurbSim exited normally, runtime was %g seconds' % (time.time() - tm0))
def run2(): tsinput = ptsin.read(fname) tsinput['URef'] = 16.0 tsr = ptsm.cfg2tsrun(tsinput) tsr.cohere = pyts.cohereModels.main.nwtc() tsr.stress = pyts.stressModels.main.uniform(0,0,0) rho = 0.5 tmax = 50.0 tsr.phase = Rinker(rho, np.pi) ### pgraf turned it off for testing! cg = tsr.grid tsr.grid = tsGrid(center=cg.center, ny=cg.n_y, nz=cg.n_z, height=cg.height, width=cg.width, time_sec=tmax, dt=cg.dt) tm0 = time.time() tsdata = tsr() ## actually runs turbsim ptsm.write(tsdata, tsinput, fname=outname+"2") print('TurbSim exited normally, runtime was %g seconds' % (time.time() - tm0))
def run2(): tsinput = ptsin.read(fname) tsinput['URef'] = 16.0 tsr = ptsm.cfg2tsrun(tsinput) tsr.cohere = pyts.cohereModels.main.nwtc() tsr.stress = pyts.stressModels.main.uniform(0, 0, 0) rho = 0.5 tmax = 50.0 tsr.phase = Rinker(rho, np.pi) ### pgraf turned it off for testing! cg = tsr.grid tsr.grid = tsGrid(center=cg.center, ny=cg.n_y, nz=cg.n_z, height=cg.height, width=cg.width, time_sec=tmax, dt=cg.dt) tm0 = time.time() tsdata = tsr() ## actually runs turbsim ptsm.write(tsdata, tsinput, fname=outname + "2") print('TurbSim exited normally, runtime was %g seconds' % (time.time() - tm0))
tsr.randgen.seed(tsr.RandSeed) else: tsr.randgen.seed(rs) ## this does nothing! np.random.seed(rs) ### this does the trick! ###tsr = self.add_phase_dist(tsr, rho, tmax) tsr.cohere = pyts.cohereModels.main.nwtc() tsr.stress = pyts.stressModels.main.uniform(0,0,0) tsr.phase = Rinker(rho, np.pi) cg = tsr.grid tsr.grid = tsGrid(center=cg.center, ny=cg.n_y, nz=cg.n_z, height=cg.height, width=cg.width, time_sec=tmax, dt=cg.dt) tsdata = tsr() ## actually runs turbsim ptsm.write(tsdata, tsinput, fname=tsoutname) # here we provide the means to link turbsim to fast: self.tswind_file = tsoutname self.tswind_dir = run_dir if __name__ == "__main__": case_inputs = {} case_inputs[("TMax")] = {'vals':[10.], 'group':0} case_inputs[("Vhub")] = {'vals':[10., 11., 12.], 'group':1} case_inputs[("Rho")] = {'vals':[.2, .25, .3], 'group':1} case_inputs[("RandSeed1")] = {'vals':[123,234], 'group':2} case_list, case_names = CaseGen_General(case_inputs, dir_matrix='', namebase='pyTurbsim_testing') filedict = {
def run1(): config = readInput(fname) tm0 = time.time() tsdat = run(config) write(tsdat, config, outname) print('TurbSim exited normally, runtime was %g seconds' % (time.time() - tm0))
def execute(self): case = self.inputs print("CASE", case.fst_params) ws=case.fst_params['Vhub'] rho = case.fst_params['Rho'] #case.fst_params['rho'] ####### TODO: how does this get here? rs = case.fst_params['RandSeed1'] if 'RandSeed1' in case.fst_params else None tmax = 2 ## should not be hard default ## if ('TMax' in case.fst_params): ## Note, this gets set via "AnalTime" in input files--FAST peculiarity ? ## tmax = case.fst_params['TMax'] # run TurbSim to generate the wind: # for now, turbsim params we mess with are possibly: TMax, RandomSeed, Tmax. These should generate # new runs, otherwise we should just use wind file we already have # for now, just differentiate by wind speed ts_case_name = "TurbSim-Vhub%.4f" % ws if rs != None: ts_case_name = "%s-Rseed%d" % (ts_case_name, rs) run_dir = os.path.join(self.basedir, ts_case_name) if (not os.path.exists(run_dir)): os.mkdir(run_dir) self._logger.info("running TurbSim in %s " % run_dir) print("running TurbSim in " , run_dir) # self.rawts.run_dir = run_dir tsdict = dict({"URef": ws, "AnalysisTime":tmax, "UsableTime":tmax}.items() + case.fst_params.items()) # self.rawts.set_dict(tsdict) print(case.fst_params) tsoutname = self.ts_file.replace("inp", "wnd") #self.rawts.ts_file.replace("inp", "wnd") tsoutname = os.path.join(run_dir, tsoutname) tssumname = tsoutname.replace("wnd", "sum") reuse_run = False if (False and os.path.isfile(tsoutname) and os.path.isfile(tssumname)): # maybe there's an old results we can use: while (not reuse_run): ln = file(tssumname).readlines() if (ln != None and len(ln) > 0): ln1 = ln[-1] # check last line2 lines (sometimes Turbsim inexplicably writes a final blank line!) ln1 = ln1.split(".") ln2 = ln[-2] # check last line2 lines (sometimes Turbsim inexplicably writes a final blank line!) ln2 = ln2.split(".") if ((len(ln1) > 0 and ln1[0] == "Processing complete") or (len(ln2) > 0 and ln2[0] == "Processing complete")): print("re-using previous TurbSim output %s for ws = %f" % (tsoutname, ws)) reuse_run = True if (not reuse_run): time.sleep(2) print("waiting for ", tsoutname) self._logger.info("waiting for %s" % tsoutname) self._logger.info("DONE waiting for %s" % tsoutname) if (not reuse_run): tsinput = ptsin.read(os.path.join(self.ts_dir, self.ts_file)) tsr = ptsm.cfg2tsrun(tsinput) tsr = self.add_phase_dist(tsr, rho, tmax) # tsdata=ptsm.run(tsinput) tsdata = tsr() ## actually runs turbsim dphi_prob = tsr.phase.delta_phi_prob ptsm.write(tsdata, tsinput, fname=tsoutname) #fout = file("uhub.out", "w") #hubdat = tsdata.uturb[0, self.mid, self.mid, :] #for i in range(len(hubdat)): # fout.write("%d %e %e\n" % ( i, tsdata.time[i], hubdat[i])) #fout.close() print("dphi prob is ", dphi_prob) fout = file(os.path.join(run_dir, "delta_phis_prob.out"), "w") fout.write("%e\n" % dphi_prob) fout.close() # tsr = self.create_tsr(rho, ws, tmax, self.dt) # tsdat = tsr() # print("writing data to ", tsdat, tsoutname) # tsdat.write_bladed(tsoutname) ### check for errors!!?? ### # here we link turbsim -> fast self.tswind_file = tsoutname self.tswind_dir = run_dir
class runTurbSimpy(Component): """ a component to run TurbSim. will try to make it aware of whether the wind file already exists""" inputs = Instance(IECRunCaseBaseVT, iotype='in') tswind_file = Str(iotype='out') tswind_dir = Str(iotype='out') def create_tsr(self, rho, U, tmax, dt): tsr = api.tsrun() tsr.grid = api.tsGrid(center=10, ny=self.ny, nz=self.nz, height=10, width=10, time_sec=tmax, dt=dt) tsr.prof = api.profModels.pl(U, 90) tsr.spec = api.specModels.tidal(self.ustar, 10) #tsr.cohere = api.cohereModels.nwtc() tsr.cohere = pyts.cohereModels.main.none() tsr.stress = api.stressModels.uniform(0.0, 0.0, 0.0) tsr.cohere = pyts.cohereModels.main.none() tsr.stress = pyts.stressModels.main.uniform(0, 0, 0) from pyts.phaseModels.main import Rinker, Uniform tsr.phase = Rinker(rho, self.mu) #tsr.phase = Uniform() # tsr.stress=np.zeros(tsr.grid.shape,dtype='float32') self.tsr = tsr return tsr def add_phase_dist(self, tsr, rho, tmax): # tsr.cohere = pyts.cohereModels.main.none() tsr.cohere = pyts.cohereModels.main.nwtc() tsr.stress = pyts.stressModels.main.uniform(0, 0, 0) tsr.phase = Rinker(rho, self.mu) cg = tsr.grid tsr.grid = tsGrid(center=cg.center, ny=cg.n_y, nz=cg.n_z, height=cg.height, width=cg.width, time_sec=tmax, dt=cg.dt) return tsr def __init__(self, filedict): super(runTurbSimpy, self).__init__() # self.rawts.ts_exe = filedict['ts_exe'] self.ts_dir = filedict['ts_dir'] self.ts_file = filedict['ts_file'] self.ustar = .8 ## self.U = 17. self.ny = 15 self.nz = 15 self.dt = 0.05 # these (should/have in past) come from TurbSim template file self.mid = self.nz / 2 ## self.rho = 0.9999 self.mu = np.pi ### np.random.seed(1) ## probably don't want this in this context # self.rawts.run_name = self.run_name self.basedir = os.path.join(os.getcwd(), "allts_runs") if 'run_dir' in filedict: self.basedir = os.path.join(os.getcwd(), filedict['run_dir']) if (not os.path.exists(self.basedir)): os.mkdir(self.basedir) def execute(self): case = self.inputs print("CASE", case.fst_params) ws = case.fst_params['Vhub'] rho = case.fst_params[ 'Rho'] #case.fst_params['rho'] ####### TODO: how does this get here? rs = case.fst_params[ 'RandSeed1'] if 'RandSeed1' in case.fst_params else None tmax = 2 ## should not be hard default ## if ( 'TMax' in case.fst_params ): ## Note, this gets set via "AnalTime" in input files--FAST peculiarity ? ## tmax = case.fst_params['TMax'] # run TurbSim to generate the wind: # for now, turbsim params we mess with are possibly: TMax, RandomSeed, Tmax. These should generate # new runs, otherwise we should just use wind file we already have # for now, just differentiate by wind speed ts_case_name = "TurbSim-Vhub%.4f" % ws if rs != None: ts_case_name = "%s-Rseed%d" % (ts_case_name, rs) run_dir = os.path.join(self.basedir, ts_case_name) if (not os.path.exists(run_dir)): os.mkdir(run_dir) self._logger.info("running TurbSim in %s " % run_dir) print("running TurbSim in ", run_dir) # self.rawts.run_dir = run_dir tsdict = dict({ "URef": ws, "AnalysisTime": tmax, "UsableTime": tmax }.items() + case.fst_params.items()) # self.rawts.set_dict(tsdict) print(case.fst_params) tsoutname = self.ts_file.replace( "inp", "wnd") #self.rawts.ts_file.replace("inp", "wnd") tsoutname = os.path.join(run_dir, tsoutname) tssumname = tsoutname.replace("wnd", "sum") reuse_run = False if (False and os.path.isfile(tsoutname) and os.path.isfile(tssumname)): # maybe there's an old results we can use: while (not reuse_run): ln = file(tssumname).readlines() if (ln != None and len(ln) > 0): ln1 = ln[ -1] # check last line2 lines (sometimes Turbsim inexplicably writes a final blank line!) ln1 = ln1.split(".") ln2 = ln[ -2] # check last line2 lines (sometimes Turbsim inexplicably writes a final blank line!) ln2 = ln2.split(".") if ((len(ln1) > 0 and ln1[0] == "Processing complete") or (len(ln2) > 0 and ln2[0] == "Processing complete")): print( "re-using previous TurbSim output %s for ws = %f" % (tsoutname, ws)) reuse_run = True if (not reuse_run): time.sleep(2) print("waiting for ", tsoutname) self._logger.info("waiting for %s" % tsoutname) self._logger.info("DONE waiting for %s" % tsoutname) if (not reuse_run): tsinput = ptsin.read(os.path.join(self.ts_dir, self.ts_file)) tsr = ptsm.cfg2tsrun(tsinput) tsr = self.add_phase_dist(tsr, rho, tmax) # tsdata=ptsm.run(tsinput) tsdata = tsr() ## actually runs turbsim dphi_prob = tsr.phase.delta_phi_prob ptsm.write(tsdata, tsinput, fname=tsoutname) #fout = file("uhub.out", "w") #hubdat = tsdata.uturb[0, self.mid, self.mid, :] #for i in range(len(hubdat)): # fout.write("%d %e %e\n" % ( i, tsdata.time[i], hubdat[i])) #fout.close() print("dphi prob is ", dphi_prob) fout = file(os.path.join(run_dir, "delta_phis_prob.out"), "w") fout.write("%e\n" % dphi_prob) fout.close() # tsr = self.create_tsr(rho, ws, tmax, self.dt) # tsdat = tsr() # print("writing data to ", tsdat, tsoutname) # tsdat.write_bladed(tsoutname) ### check for errors!!?? ### # here we link turbsim -> fast self.tswind_file = tsoutname self.tswind_dir = run_dir
#!/usr/bin/python """ This is the pyTurbSim 'executable script', which utilizes the :mod:`pyts.runConfig` package. """ import sys from pyts.runInput.main import readInput, run, write import time if len(sys.argv) > 1: fname = sys.argv[1] else: fname = 'HydroTurbSim.inp' config = readInput(fname) tm0 = time.time() tsdat = run(config) write(tsdat, config, fname) print 'TurbSim exited normally, runtime was %g seconds' % (time.time() - tm0)