def Initialize_( data,grid,grid1, elevation ):
## Plant types are defined as following:
# GRASS = 0; SHRUB = 1; TREE = 2; BARE = 3;
# SHRUBSEEDLING = 4; TREESEEDLING = 5
## Initialize random plant type field
grid['cell']['VegetationType'] = np.random.choice([0,1,2,3,0,2], \
grid.number_of_cells)
## Assign plant type for representative ecohydrologic simulations
grid1['cell']['VegetationType'] = np.arange(0,6)
grid1['node']['Elevation'] = 1700. * np.ones(grid1.number_of_nodes)
grid['node']['Elevation'] = elevation
PD_D = PrecipitationDistribution(mean_storm = data['mean_storm_dry'], \
mean_interstorm = data['mean_interstorm_dry'],
mean_storm_depth = data['mean_storm_depth_dry'])
PD_W = PrecipitationDistribution(mean_storm = data['mean_storm_wet'], \
mean_interstorm = data['mean_interstorm_wet'],
mean_storm_depth = data['mean_storm_depth_wet'])
Rad = Radiation( grid )
Rad_PET = Radiation( grid1 )
PET_Tree = PotentialEvapotranspiration( grid1, method = data['PET_method'], \
MeanTmaxF = data['MeanTmaxF_tree'],
DeltaD = data['DeltaD'] )
PET_Shrub = PotentialEvapotranspiration( grid1, method = data['PET_method'], \
MeanTmaxF = data['MeanTmaxF_shrub'],
DeltaD = data['DeltaD'] )
PET_Grass = PotentialEvapotranspiration( grid1, method = data['PET_method'], \
MeanTmaxF = data['MeanTmaxF_grass'],
DeltaD = data['DeltaD'] )
SM = SoilMoisture( grid, data ) # Soil Moisture object
VEG = Vegetation( grid, data ) # Vegetation object
vegca = VegCA( grid, data ) # Cellular automaton object
## Initializing inputs for Soil Moisture object
grid['cell']['LiveLeafAreaIndex'] = 1.6 * np.ones( grid.number_of_cells )
SM._SO = 0.59 * np.ones(grid.number_of_cells) # Initializing Soil Moisture
return PD_D, PD_W, Rad, Rad_PET, PET_Tree, PET_Shrub, PET_Grass, SM, \
VEG, vegca
mean_storm_depth = data['mean_storm_depth_dry'])
PD_W = PrecipitationDistribution(mean_storm = data['mean_storm_wet'], \
mean_interstorm = data['mean_interstorm_wet'],
mean_storm_depth = data['mean_storm_depth_wet'])
Rad = Radiation( grid )
PET_Tree = PotentialEvapotranspiration( grid1, method = data['PET_method'], \
MeanTmaxF = data['MeanTmaxF_tree'],
DeltaD = data['DeltaD'] )
PET_Shrub = PotentialEvapotranspiration( grid1, method = data['PET_method'], \
MeanTmaxF = data['MeanTmaxF_shrub'],
DeltaD = data['DeltaD'] )
PET_Grass = PotentialEvapotranspiration( grid1, method = data['PET_method'], \
MeanTmaxF = data['MeanTmaxF_grass'],
DeltaD = data['DeltaD'] )
SM = SoilMoisture( grid, data ) # Soil Moisture object
VEG = Vegetation( grid, data ) # Vegetation object
vegca = VegCA( grid1, data ) # Cellular automaton object
##########
n = 6600 # Defining number of storms the model will be run
##########
## Create arrays to store modeled data
P = np.empty(n) # Record precipitation
Tb = np.empty(n) # Record inter storm duration
Tr = np.empty(n) # Record storm duration
Time = np.empty(n) # To record time elapsed from the start of simulation
CumWaterStress = np.empty([n/50, grid1.number_of_cells]) # Cum Water Stress
VegType = np.empty([n/50, grid1.number_of_cells],dtype = int)
mean_storm_depth = data['mean_storm_depth_dry'])
PD_W = PrecipitationDistribution(mean_storm = data['mean_storm_wet'], \
mean_interstorm = data['mean_interstorm_wet'],
mean_storm_depth = data['mean_storm_depth_wet'])
Rad = Radiation(grid)
PET_Tree = PotentialEvapotranspiration( grid1, method = data['PET_method'], \
MeanTmaxF = data['MeanTmaxF_tree'],
DeltaD = data['DeltaD'] )
PET_Shrub = PotentialEvapotranspiration( grid1, method = data['PET_method'], \
MeanTmaxF = data['MeanTmaxF_shrub'],
DeltaD = data['DeltaD'] )
PET_Grass = PotentialEvapotranspiration( grid1, method = data['PET_method'], \
MeanTmaxF = data['MeanTmaxF_grass'],
DeltaD = data['DeltaD'] )
SM = SoilMoisture(grid, data) # Soil Moisture object
VEG = Vegetation(grid, data) # Vegetation object
vegca = VegCA(grid1, data) # Cellular automaton object
##########
n = data['n_short'] # Defining number of storms the model will be run
##########
## Create arrays to store modeled data
P = np.empty(n) # Record precipitation
Tb = np.empty(n) # Record inter storm duration
Tr = np.empty(n) # Record storm duration
Time = np.empty(n) # To record time elapsed from the start of simulation
CumWaterStress = np.empty([n / 50, grid1.number_of_cells]) # Cum Water Stress
CumWS = np.empty([n / 50, grid.number_of_cells
mean_storm_depth = data['mean_storm_depth_dry'])
PD_W = PrecipitationDistribution(mean_storm = data['mean_storm_wet'], \
mean_interstorm = data['mean_interstorm_wet'],
mean_storm_depth = data['mean_storm_depth_wet'])
Rad = Radiation( grid )
Rad_PET = Radiation( grid1 )
PET_Tree = PotentialEvapotranspiration( grid1, method = data['PET_method'], \
MeanTmaxF = data['MeanTmaxF_tree'],
DeltaD = data['DeltaD'] )
PET_Shrub = PotentialEvapotranspiration( grid1, method = data['PET_method'], \
MeanTmaxF = data['MeanTmaxF_shrub'],
DeltaD = data['DeltaD'] )
PET_Grass = PotentialEvapotranspiration( grid1, method = data['PET_method'], \
MeanTmaxF = data['MeanTmaxF_grass'],
DeltaD = data['DeltaD'] )
SM = SoilMoisture( grid, data ) # Soil Moisture object
VEG = Vegetation( grid, data ) # Vegetation object
vegca = VegCA( grid, data ) # Cellular automaton object
##########
n = data['n_long_DEM'] # Defining number of storms the model will be run
##########
## Create arrays to store modeled data
P = np.empty(n) # Record precipitation
Tb = np.empty(n) # Record inter storm duration
Tr = np.empty(n) # Record storm duration
Time = np.empty(n) # To record time elapsed from the start of simulation
CumWaterStress = np.empty([n/55, grid.number_of_cells]) # Cum Water Stress
VegType = np.empty([n/55, grid.number_of_cells],dtype=int)
mean_storm_depth = data['mean_storm_depth_dry'])
PD_W = PrecipitationDistribution(mean_storm = data['mean_storm_wet'], \
mean_interstorm = data['mean_interstorm_wet'],
mean_storm_depth = data['mean_storm_depth_wet'])
Rad = Radiation(grid)
Rad_PET = Radiation(grid1)
PET_Tree = PotentialEvapotranspiration( grid1, method = data['PET_method'], \
MeanTmaxF = data['MeanTmaxF_tree'],
DeltaD = data['DeltaD'] )
PET_Shrub = PotentialEvapotranspiration( grid1, method = data['PET_method'], \
MeanTmaxF = data['MeanTmaxF_shrub'],
DeltaD = data['DeltaD'] )
PET_Grass = PotentialEvapotranspiration( grid1, method = data['PET_method'], \
MeanTmaxF = data['MeanTmaxF_grass'],
DeltaD = data['DeltaD'] )
SM = SoilMoisture(grid, data) # Soil Moisture object
VEG = Vegetation(grid, data) # Vegetation object
vegca = VegCA(grid, data) # Cellular automaton object
##########
n = data['n_long_DEM'] # Defining number of storms the model will be run
##########
## Create arrays to store modeled data
P = np.empty(n) # Record precipitation
Tb = np.empty(n) # Record inter storm duration
Tr = np.empty(n) # Record storm duration
Time = np.empty(n) # To record time elapsed from the start of simulation
CumWaterStress = np.empty([n / 55, grid.number_of_cells]) # Cum Water Stress
VegType = np.empty([n / 55, grid.number_of_cells], dtype=int)