grid1.number_of_cells)
# Age of the plants is randomnly generated by the constructor of vegca
# component. Seedlings of shrubs and trees are determined by their age
# within the constructor as well. Hence the random vegetation type
# field is called for types GRASS, SHRUB, TREE and BARE only. -SN 10Mar15
grid['cell']['VegetationType'] = np.arange(0, 6)
# Create radiation, soil moisture and Vegetation objects
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)
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
##########
# component. Seedlings of shrubs and trees are determined by their age
# within the constructor as well. Hence the random vegetation type
# field is called for types GRASS, SHRUB, TREE and BARE only. -SN 10Mar15
# Dummy grid for PET
grid1 = rmg(5, 4, 5)
grid1['node']['Elevation'] = 1700. * np.ones(grid1.number_of_nodes)
grid1['cell']['VegetationType'] = np.arange(0, 6)
# Create radiation, soil moisture and Vegetation objects
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
##########
