Beispiel #1
0
        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
Beispiel #2
0
         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()
Beispiel #3
0
        
    

   
    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!!