# now get the forcings that wflow expects # The bmi is such that you can get the input variables and the output variables. However, the # input variable list also contains the in/out variables. So to # get the input only we subtract the two lists. invars = LA_model.get_input_var_names() outvars = LA_model.get_output_var_names() inputmstacks = list(set(invars) - set(outvars)) # In[]: Investigate start time, end time and time step of both models print "WFlow:" LA_dt = LA_model.get_time_step() # LA_start = LA_model.get_start_time() timeutc = adapter.getStartTimefromRuninfo("inmaps/runinfo.xml") print timeutc LA_start = calendar.timegm(timeutc.timetuple()) timeutc = adapter.getEndTimefromRuninfo("inmaps/runinfo.xml") LA_end = calendar.timegm(timeutc.timetuple()) # LA_end = LA_model.get_end_time() print LA_dt print timeutc print LA_start print LA_end print "RTC-Tools" RTC_dt = RTC_model.get_time_step() RTC_start = RTC_model.get_start_time() RTC_end = RTC_model.get_end_time()
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
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 logger.info("Reading basemaps")
# now get the forcings that wflow expects # The bmi is such that you can get the input variables and the output variables. However, the # input variable list also contains the in/out variables. So to # get the input only we subtract the two lists. invars = LA_model.get_input_var_names() outvars = LA_model.get_output_var_names() inputmstacks = list(set(invars) - set(outvars)) # In[]: Investigate start time, end time and time step of both models print("WFlow:") LA_dt = LA_model.get_time_step() # LA_start = LA_model.get_start_time() timeutc = adapter.getStartTimefromRuninfo("inmaps/runinfo.xml") print(timeutc) LA_start = calendar.timegm(timeutc.timetuple()) timeutc = adapter.getEndTimefromRuninfo("inmaps/runinfo.xml") LA_end = calendar.timegm(timeutc.timetuple()) # LA_end = LA_model.get_end_time() print(LA_dt) print(timeutc) print(LA_start) print(LA_end) print("RTC-Tools") RTC_dt = RTC_model.get_time_step() RTC_start = RTC_model.get_start_time() RTC_end = RTC_model.get_end_time()