def read_input_files(self): rti = self.rti #------------------------------------------------------- # All grids are assumed to have data type of Float32. #------------------------------------------------------- # NB! Green-Ampt and Smith-Parlange currently only # support ONE layer (n_layers == 1). #------------------------------------------------------- for k in xrange(self.n_layers): Ks = model_input.read_next(self.Ks_unit[k], self.Ks_type[k], rti) if (Ks != None): self.Ks[k] = Ks Ki = model_input.read_next(self.Ki_unit[k], self.Ki_type[k], rti) if (Ki != None): self.Ki[k] = Ki qs = model_input.read_next(self.qs_unit[k], self.qs_type[k], rti) if (qs != None): self.qs[k] = qs qi = model_input.read_next(self.qi_unit[k], self.qi_type[k], rti) if (qi != None): self.qi[k] = qi G = model_input.read_next(self.G_unit[k], self.G_type[k], rti) if (G != None): self.G[k] = G
def read_input_files(self): rti = self.rti #------------------------------------------------------- # All grids are assumed to have data type of Float32. #------------------------------------------------------- # NB! Green-Ampt and Smith-Parlange currently only # support ONE layer (n_layers == 1). #------------------------------------------------------- for k in range(self.n_layers): Ks = model_input.read_next(self.Ks_unit[k], self.Ks_type[k], rti) if (Ks is not None): self.Ks[k] = Ks Ki = model_input.read_next(self.Ki_unit[k], self.Ki_type[k], rti) if (Ki is not None): self.Ki[k] = Ki qs = model_input.read_next(self.qs_unit[k], self.qs_type[k], rti) if (qs is not None): self.qs[k] = qs qi = model_input.read_next(self.qi_unit[k], self.qi_type[k], rti) if (qi is not None): self.qi[k] = qi G = model_input.read_next(self.G_unit[k], self.G_type[k], rti) if (G is not None): self.G[k] = G
def read_input_files(self): rti = self.rti #------------------------------------------------------- # All grids are assumed to have data type of Float32. #------------------------------------------------------- # NB! Beven, Green-Ampt & Smith-Parlange currently only # support ONE layer (n_layers == 1). #------------------------------------------------------- for k in xrange(self.n_layers): Ks = model_input.read_next(self.Ks_unit[k], self.Ks_type[k], rti) if (Ks is not None): self.Ks[k] = Ks qs = model_input.read_next(self.qs_unit[k], self.qs_type[k], rti) if (qs is not None): self.qs[k] = qs qi = model_input.read_next(self.qi_unit[k], self.qi_type[k], rti) if (qi is not None): self.qi[k] = qi f = model_input.read_next(self.f_unit[k], self.f_type[k], rti) if (f is not None): self.f[k] = f C = model_input.read_next(self.C_unit[k], self.C_type[k], rti) if (C is not None): self.C[k] = C
def read_input_files(self): rti = self.rti #------------------------------------------------------- # All grids are assumed to have a data type of Float32. #------------------------------------------------------- alpha = model_input.read_next(self.alpha_unit, self.alpha_type, rti) if (alpha != None): self.alpha = alpha K_soil = model_input.read_next(self.K_soil_unit, self.K_soil_type, rti) if (K_soil != None): self.K_soil = K_soil soil_x = model_input.read_next(self.soil_x_unit, self.soil_x_type, rti) if (soil_x != None): self.soil_x = soil_x T_soil_x = model_input.read_next(self.T_soil_x_unit, self.T_soil_x_type, rti) if (T_soil_x != None): self.T_soil_x = T_soil_x
def read_input_files(self): rti = self.rti #------------------------------------------------------- # All grids are assumed to have a data type of Float32. #------------------------------------------------------- ## alpha = model_input.read_next(self.alpha_unit, self.alpha_type, rti) ## if (alpha is not None): self.alpha = alpha K_soil = model_input.read_next(self.K_soil_unit, self.K_soil_type, rti) if (K_soil is not None): self.K_soil = K_soil soil_x = model_input.read_next(self.soil_x_unit, self.soil_x_type, rti) if (soil_x is not None): self.soil_x = soil_x T_soil_x = model_input.read_next(self.T_soil_x_unit, self.T_soil_x_type, rti) if (T_soil_x is not None): self.T_soil_x = T_soil_x
def read_input_files(self): rti = self.rti #------------------------------------------------------- # All grids are assumed to have a data type of Float32. #------------------------------------------------------- c0 = model_input.read_next(self.c0_unit, self.c0_type, rti) if (c0 is not None): self.c0 = c0 T0 = model_input.read_next(self.T0_unit, self.T0_type, rti) if (T0 is not None): self.T0 = T0 rho_snow = model_input.read_next(self.rho_snow_unit, self.rho_snow_type, rti) if (rho_snow is not None): self.rho_snow = rho_snow h0_snow = model_input.read_next(self.h0_snow_unit, self.h0_snow_type, rti) if (h0_snow is not None): self.h0_snow = h0_snow h0_swe = model_input.read_next(self.h0_swe_unit, self.h0_swe_type, rti) if (h0_swe is not None): self.h0_swe = h0_swe
def read_input_files(self): rti = self.rti #-------------------------------------------------------- # All grids are assumed to have a data type of Float32. #-------------------------------------------------------- Cp_snow = model_input.read_next(self.Cp_snow_unit, self.Cp_snow_type, rti) if (Cp_snow != None): self.Cp_snow = Cp_snow rho_snow = model_input.read_next(self.rho_snow_unit, self.rho_snow_type, rti) if (rho_snow != None): self.rho_snow = rho_snow T0 = model_input.read_next(self.T0_unit, self.T0_type, rti) if (T0 != None): self.T0 = T0 h0_snow = model_input.read_next(self.h0_snow_unit, self.h0_snow_type, rti) if (h0_snow != None): self.h0_snow = h0_snow h0_swe = model_input.read_next(self.h0_swe_unit, self.h0_swe_type, rti) if (h0_swe != None): self.h0_swe = h0_swe
def read_input_files(self): rti = self.rti #------------------------------------------------------- # All grids are assumed to have data type of Float32. #------------------------------------------------------- # This method works for Green-Ampt and Smith-Parlange # but must be overridden for Richards 1D. #------------------------------------------------------- # NB! Green-Ampt and Smith-Parlange currently only # support ONE layer (n_layers == 1). #------------------------------------------------------- for k in xrange(self.n_layers): Ks = model_input.read_next(self.Ks_unit[k], self.Ks_type[k], rti) if (Ks is not None): self.Ks[k] = Ks Ki = model_input.read_next(self.Ki_unit[k], self.Ki_type[k], rti) if (Ki is not None): self.Ki[k] = Ki qs = model_input.read_next(self.qs_unit[k], self.qs_type[k], rti) if (qs is not None): self.qs[k] = qs qi = model_input.read_next(self.qi_unit[k], self.qi_type[k], rti) if (qi is not None): self.qi[k] = qi G = model_input.read_next(self.G_unit[k], self.G_type[k], rti) if (G is not None): self.G[k] = G gam = model_input.read_next(self.gam_unit[k], self.gam_type[k], rti) if (gam is not None): self.gam[k] = gam
def read_input_files(self): rti = self.rti #------------------------------------------------------- # All grids are assumed to have data type of Float32. #------------------------------------------------------- # This method works for Green-Ampt and Smith-Parlange # but must be overridden for Richards 1D. #------------------------------------------------------- # NB! Green-Ampt and Smith-Parlange currently only # support ONE layer (n_layers == 1). #------------------------------------------------------- for k in range(self.n_layers): Ks = model_input.read_next(self.Ks_unit[k], self.Ks_type[k], rti) if (Ks is not None): self.Ks_list[k] = Ks Ki = model_input.read_next(self.Ki_unit[k], self.Ki_type[k], rti) if (Ki is not None): self.Ki_list[k] = Ki qs = model_input.read_next(self.qs_unit[k], self.qs_type[k], rti) if (qs is not None): self.qs_list[k] = qs qi = model_input.read_next(self.qi_unit[k], self.qi_type[k], rti) if (qi is not None): self.qi_list[k] = qi G = model_input.read_next(self.G_unit[k], self.G_type[k], rti) if (G is not None): self.G_list[k] = G gam = model_input.read_next(self.gam_unit[k], self.gam_type[k], rti) if (gam is not None): self.gam_list[k] = gam
def read_input_files(self): #-------------- # For testing #-------------- ## print 'self.ET_type =', self.ET_type ## print 'self.ET_file =', self.ET_file rti = self.rti #------------------------------------------------------- # All grids are assumed to have a data type of Float32. #------------------------------------------------------- ET = model_input.read_next(self.ET_unit, self.ET_type, rti, factor=self.mmph_to_mps) if (ET != None): self.ET = ET print 'min(ET) =', ET.min() * self.mps_to_mmph, ' [mmph]' print 'min(ET) =', ET.max() * self.mps_to_mmph, ' [mmph]'
def read_input_files(self): #-------------- # For testing #-------------- ## print 'self.ET_type =', self.ET_type ## print 'self.ET_file =', self.ET_file rti = self.rti #------------------------------------------------------- # All grids are assumed to have a data type of Float32. #------------------------------------------------------- ET = model_input.read_next(self.ET_unit, self.ET_type, rti, factor=self.mmph_to_mps) if (ET is not None): self.ET = ET print('min(ET) =', ET.min() * self.mps_to_mmph, ' [mmph]') print('min(ET) =', ET.max() * self.mps_to_mmph, ' [mmph]')