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)