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 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 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):
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):
#--------------
# 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]')