if o == '-A': areamap = a.strip() if o == '-c': configfile = a.strip() if o == '-O': T0 = datetime.strptime(a, '%Y-%m-%d %H:%M:%S') global pointer dw_CreateDwRun(dwdir) config = ConfigParser.SafeConfigParser() config.optionxform = str config.read(caseId + "/" + configfile) timestepsecs = int( configget(config, "model", "timestepsecs", str(timestepsecs))) if fewsrun: timeSteps = wflow_adapt.getTimeStepsfromRuninfo( runinfoFile, timestepsecs) T0 = wflow_adapt.getStartTimefromRuninfo(runinfoFile) print timeSteps #: we need one delwaq calculation timesteps less than hydrology # timeSteps = timeSteps # need one more hydrological timestep as dw timestep firstTimeStep = 0 logger = pcrut.setlogger(dwdir + "/debug/wflow_delwaq.log", "wflow_delwaq") #caseid = "default_hbv" logger.info("T0 of run: " + str(T0)) boundids = len(sourcesMap) # extra number of exchnages for all bounds #Number of exchnages is elements minus number of outflows!! # Get subcatchment data
print "WARN: -P Does not work at the moment" if o == '-p': exec "multdynapars =" + a print "WARN: -p Does not work at the moment" if o == '-C': caseName = a if o == '-R': runId = a if o == '-c': configfile = a if o == '-s': timestepsecs = int(a) if o == '-T': _lastTimeStep=int(a) if o == '-S': _firstTimeStep=int(a) if o == '-h': usage() if o == '-f': NoOverWrite = 1 if fewsrun: _lastTimeStep = wflow_adapt.getTimeStepsfromRuninfo(runinfoFile) * 86400/timestepsecs _firstTimeStep = 1 myModel = WflowModel(wflow_cloneMap, caseName,runId,configfile) dynModelFw = wf_DynamicFramework(myModel, _lastTimeStep,firstTimestep=_firstTimeStep) dynModelFw.createRunId() for o, a in opts: if o == '-X': configset(myModel.config,'model','OverWriteInit','1',overwrite=True) if o == '-I': configset(myModel.config,'model','reinit','1',overwrite=True) if o == '-i': configset(myModel.config,'model','intbl',a,overwrite=True) if o == '-s': configset(myModel.config,'model','timestepsecs',a,overwrite=True) if o == '-x': configset(myModel.config,'model','sCatch',a,overwrite=True) if o == '-c': configset(myModel.config,'model','configfile', a,overwrite=True) if o == '-M': configset(myModel.config,'model','MassWasting',"1",overwrite=True) if o == '-h': usage()
global pointer dw_CreateDwRun(dwdir) config = ConfigParser.SafeConfigParser() config.optionxform = str config.read(caseId + "/" + configfile) timestepsecs = int(configget(config,"model","timestepsecs",str(timestepsecs))) if fewsrun: timeSteps = wflow_adapt.getTimeStepsfromRuninfo(runinfoFile,timestepsecs) T0 = wflow_adapt.getStartTimefromRuninfo(runinfoFile) print timeSteps #: we need one delwaq calculation timesteps less than hydrology # timeSteps = timeSteps # need one more hydrological timestep as dw timestep firstTimeStep = 0 logger = pcrut.setlogger(dwdir + "/debug/wflow_delwaq.log","wflow_delwaq") #caseid = "default_hbv" logger.info("T0 of run: " + str(T0)) boundids = len(sourcesMap) # extra number of exchnages for all bounds #Number of exchnages is elements minus number of outflows!!