def opt_fun(par): try: # distribute the parameters SpatialVarFun.Function(par, kub=SpatialVarFun.Kub, klb=SpatialVarFun.Klb, Maskingum=SpatialVarFun.Maskingum) self.Parameters = SpatialVarFun.Par3d #run the model Wrapper.FW1(self) # calculate performance of the model try: # error = self.OF(self.QGauges, self.qout, self.quz_routed, self.qlz_translated,*[self.GaugesTable]) error = self.OF(self.QGauges, self.qout, *[self.GaugesTable]) except TypeError: # if no of inputs less than what the function needs assert False, "the objective function you have entered needs more inputs please enter then in a list as *args" # print error if printError != 0: print(error) print(par) fail = 0 except: error = np.nan fail = 1 return error, [], fail
def RunFW1(self): """ ======================================================================= RunDistwithLake(PrecPath, Evap_Path, TempPath, DemPath, FlowAccPath, FlowDPath, ParPath, p2) ======================================================================= this function runs the conceptual distributed hydrological model Inputs: ---------- 1-Paths: 1-PrecPath: [String] path to the Folder contains precipitation rasters 2-Evap_Path: [String] path to the Folder contains Evapotranspiration rasters 3-TempPath: [String] path to the Folder contains Temperature rasters 4-FlowAccPath: [String] path to the Flow Accumulation raster of the catchment (it should include the raster name and extension) 5-FlowDPath: [String] path to the Flow Direction raster of the catchment (it should include the raster name and extension) 7-ParPath: [String] path to the Folder contains parameters rasters of the catchment 8-p2: [List] list of unoptimized parameters p2[0] = tfac, 1 for hourly, 0.25 for 15 min time step and 24 for daily time step p2[1] = catchment area in km2 Outputs: ---------- 1- st: [4D array] state variables 2- q_out: [1D array] calculated Discharge at the outlet of the catchment 3- q_uz: [3D array] Distributed discharge for each cell Example: ---------- PrecPath = prec_path="meteodata/4000/calib/prec" Evap_Path = evap_path="meteodata/4000/calib/evap" TempPath = temp_path="meteodata/4000/calib/temp" DemPath = "GIS/4000/dem4000.tif" FlowAccPath = "GIS/4000/acc4000.tif" FlowDPath = "GIS/4000/fd4000.tif" ParPath = "meteodata/4000/parameters" p2=[1, 227.31] st, q_out, q_uz_routed = RunModel(PrecPath,Evap_Path,TempPath,DemPath, FlowAccPath,FlowDPath,ParPath,p2) """ # input data validation # input dimensions assert np.shape(self.Prec)[0] == self.rows and np.shape( self.ET )[0] == self.rows and np.shape(self.Temp)[0] == self.rows and np.shape( self.Parameters )[0] == self.rows, "all input data should have the same number of rows" assert np.shape(self.Prec)[1] == self.cols and np.shape( self.ET )[1] == self.cols and np.shape(self.Temp)[1] == self.cols and np.shape( self.Parameters )[1] == self.cols, "all input data should have the same number of columns" assert np.shape(self.Prec)[2] == np.shape(self.ET)[2] and np.shape( self.Temp )[2], "all meteorological input data should have the same length" #run the model Wrapper.FW1(self) print("Model Run has finished")