def run_attm(self):
""" Program sequence """
#====================================================
# Initialization Process
#====================================================
print '==================='
print ' Initializing ATTM'
print '==================='
read_control.read_control(self)
initialize.initialize(self)
read_layers.read_layers(self)
model_domain.model_domain(self)
create_attm_cohort_arrays.create_attm_cohort_arrays(self)
if self.Simulation_area.lower() == 'barrow':
initial_cohort_population.barrow_initial_cohort_population(self)
initial_cohort_check.barrow_initial_cohort_check(self)
cohort_present.barrow_cohort_present(self)
elif self.Simulation_area.lower() == 'tanana':
initial_cohort_population.tanana_initial_cohort_population(self)
initial_cohort_check.tanana_initial_cohort_check(self)
cohort_present.tanana_cohort_present(self)
#=======================================
# READ MET Data & Calculate Degree Days
#=======================================
initialize.Met(self)
#++++++++++++++++++++++++++++++++++++++++++++++
# ========================================
# INITIALIZE BARROW COHORT PROPERTIES
# ========================================
#++++++++++++++++++++++++++++++++++++++++++++++
if self.Simulation_area.lower() == 'barrow':
print '=================================== '
print ' Initializing Lake & Pond Properties'
print '===================================='
initialize.LakePond(self)
set_lake_pond_depth.set_lake_pond_depth(self)
set_lake_ice_depth_constant.set_lake_ice_depth_constant(self)
set_ice_thickness_array.set_ice_thickness_array(self)
climate_expansion_arrays.set_climate_expansion_arrays(self)
set_pond_growth_array.set_pond_growth_array(self)
print '====================================='
print ' Initializing Terrestrial Properties'
print '====================================='
initialize.Terrestrial_Barrow(self)
read_ice_content.read_ice_content(self)
read_drainage_efficiency.read_drainage_efficiency(self)
read_initial_ALD.read_initial_ALD(self)
set_ALD_constant.set_ALD_constant(self)
set_ALD_array.set_ALD_array(self)
set_protective_layer.set_protective_layer(self)
set_initial_cumulative_probability.set_initial_cumulative_probability(self)
# Initializing Terrestrial Cohort Properties
initialize.Wet_NPG(self)
initialize.Wet_LCP(self)
initialize.Wet_CLC(self)
initialize.Wet_FCP(self)
initialize.Wet_HCP(self)
# Other needed information [in the future]
initial_cohort_age.initial_cohort_age(self)
elif self.Simulation_area.lower() == 'tanana':
print '======================================'
print ' Initializing Terrestrial Properties '
print '======================================'
initialize.Terrestrial_Tanana(self)
print '=================================================='
print ' Starting the MAIN LOOP '
print '=================================================='
initialize.run(self)
for time in range(0, self.stop):
if time == 0:
if self.Simulation_area.lower() == 'barrow':
cohorts.initial_barrow(self)
elif self.Simulation_area.lower() == 'tanana':
cohorts.initial_tanana(self)
print ' at time step: ', time
# ++++++++++++++++++++++++++++++++++++++
# Check for significant climatic event
# ++++++++++++++++++++++++++++++++++++++
check_climate_event.check_climate_event(self)
# ----------------------------------------------------------
# Looping over elements
# ----------------------------------------------------------
for element in range(0, self.ATTM_nrows * self.ATTM_ncols):
# ----------------------------------------------------
# Define the total fractional area of cohorts for
# each element
# ----------------------------------------------------
cohort_start = cohort_check.cohort_start(self, element, time)
# ----------------------------------------------------
# Expand/Infill lake & ponds by prescribed rates
# ----------------------------------------------------
lake_pond_expansion.lake_pond_expansion(self, element)
lake_pond_expansion.pond_infill(self, element, time)
# ----------------------------------------------------------
# Set active layer depth
# ---------------------------------------------------------
active_layer_depth.active_layer_depth(self, time, element)
# ----------------------------------
# Cycle through terrestrial cohorts
# ----------------------------------
check_Wet_NPG.check_Wet_NPG(self, element, time)
check_Wet_LCP.check_Wet_LCP(self, element, time)
check_Wet_CLC.check_Wet_CLC(self, element, time)
check_Wet_FCP.check_Wet_FCP(self, element, time)
check_Wet_HCP.check_Wet_HCP(self, element, time)
# ----------------------------------
# Set pond/lake ice thickness depth
# ----------------------------------
ice_thickness.ice_thickness(self, time, element)
# ------------------------------
# Cycle through ponds and lakes
# ------------------------------
check_Ponds.check_Ponds(self, element, time)
check_Lakes.check_Lakes(self, element, time)
# -------------------------------------------------
# Cohort Fraction Check (mass balance of cohorts)
# -------------------------------------------------
cohort_check.cohort_check(self, element, time, cohort_start)
if time == self.stop-1:
if self.Simulation_area.lower() == 'barrow':
cohorts.final_barrow(self)
elif self.Simulation_area.lower() == 'tanana':
cohorts.final_tanana(self)
# ========================================================================
# END MAIN LOOP
# ========================================================================
# ========================================================================
# OUTPUT RESULTS (if requested)
# ========================================================================
# - - - - - - - - -
# Fractional Areas
# - - - - - - - - -
Output_cohorts_by_year.Output_cohorts_by_year(self, time)
# - - - - - - - - - - - - -
# Dominant Fractional Area
# - - - - - - - - - - - - -
Output_cohorts_by_year.dominant_cohort(self) # Terrestrial_Control
Output_cohorts_by_year.dominant_fractional_plot(self, time) # Terrestrial_Control
# =================================
# OUTPUT ANIMATIONS (if requested)
# =================================
# - - - - - - - - - - - - - - -
# Fractional Area of Cohorts
# - - - - - - - - - - - - - - - -
Output_cohorts_by_year.write_Fractions_avi(self)
Output_cohorts_by_year.write_Dominant_Cohort_avi(self) # Terrestrial_Control
# -------------------
# Simulation End Time
# -------------------
clock.finish(self)
#===========================
# Output Simulation Results
#===========================
if self.results_onscreen.lower() == 'yes':
results.on_screen(self)
if self.archive_simulation.lower() == 'yes':
results.on_file(self)
# ================
# Archive Results
# ================
if self.archive_simulation.lower() == 'yes':
#----------------------------------------------------------------------------------------------------------
# Create the tarfile
#----------------------------------------------------------------------------------------------------------
self.archive_file =tarfile.open(self.control['Run_dir']+self.Output_directory+str('/Archive/')+ \
self.archive_time+str('_')+self.simulation_name+".tar.gz", mode='w:gz')
#----------------------------------------------------------------------------------------------------------
archive.read_archive(self)
archive.archive(self)
print '----------------------------------------'
print ' Simulation Complete '
print '----------------------------------------'
def run_barrow(self, time):
for time in range(0, self.stop):
if time == 0:
cohorts.initial_barrow(self)
print ' at time step: ', time
# ++++++++++++++++++++++++++++++++++++++
# Check for significant climatic event
# ++++++++++++++++++++++++++++++++++++++
check_climate_event.check_climate_event(self)
# ----------------------------------------------------------
# Looping over elements
# ----------------------------------------------------------
for element in range(0, self.ATTM_nrows * self.ATTM_ncols):
# ----------------------------------------------------
# Define the total fractional area of cohorts for
# each element
# ----------------------------------------------------
cohort_start = cohort_check.cohort_start(self, element, time)
# ----------------------------------------------------
# Expand/Infill lake & ponds by prescribed rates
# ----------------------------------------------------
lake_pond_expansion.lake_pond_expansion(self, element)
lake_pond_expansion.pond_infill(self, element, time)
# ----------------------------------------------------------
# Set active layer depth
# ---------------------------------------------------------
active_layer_depth.active_layer_depth(self, time, element)
# ----------------------------------
# Cycle through terrestrial cohorts
# ----------------------------------
check_Wet_NPG.check_Wet_NPG(self, element, time)
check_Wet_LCP.check_Wet_LCP(self, element, time)
check_Wet_CLC.check_Wet_CLC(self, element, time)
check_Wet_FCP.check_Wet_FCP(self, element, time)
check_Wet_HCP.check_Wet_HCP(self, element, time)
# ----------------------------------
# Set pond/lake ice thickness depth
# ----------------------------------
ice_thickness.ice_thickness(self, time, element)
# ------------------------------
# Cycle through ponds and lakes
# ------------------------------
check_Ponds.check_Ponds(self, element, time)
check_Lakes.check_Lakes(self, element, time)
# -------------------------------------------------
# Cohort Fraction Check (mass balance of cohorts)
# -------------------------------------------------
cohort_check.cohort_check(self, element, time, cohort_start)
if time == self.stop -1:
cohorts.final_barrow(self)
# ========================================================================
# END MAIN LOOP
# ========================================================================
# ========================================================================
# OUTPUT RESULTS (if requested)
# ========================================================================
# - - - - - - - - -
# Fractional Areas
# - - - - - - - - -
Output_cohorts_by_year.Output_cohorts_by_year(self, time)
# - - - - - - - - - - - - -
# Dominant Fractional Area
# - - - - - - - - - - - - -
Output_cohorts_by_year.dominant_cohort(self) # Terrestrial_Control
Output_cohorts_by_year.dominant_fractional_plot(self, time) # Terrestrial_Control
# =================================
# OUTPUT ANIMATIONS (if requested)
# =================================
# - - - - - - - - - - - - - - -
# Fractional Area of Cohorts
# - - - - - - - - - - - - - - - -
Output_cohorts_by_year.write_Fractions_avi(self)
Output_cohorts_by_year.write_Dominant_Cohort_avi(self) # Terrestrial_Control
def run_barrow(self, time):
for time in range(0, self.stop):
if time == 0:
cohorts.initial_barrow(self)
print ' at time step: ', time
# ++++++++++++++++++++++++++++++++++++++
# Check for significant climatic event
# ++++++++++++++++++++++++++++++++++++++
check_climate_event.check_climate_event(self)
# ----------------------------------------------------------
# Looping over elements
# ----------------------------------------------------------
for element in range(0, self.ATTM_nrows * self.ATTM_ncols):
# ----------------------------------------------------
# Define the total fractional area of cohorts for
# each element
# ----------------------------------------------------
cohort_start = cohort_check.cohort_start(self, element, time)
# ----------------------------------------------------
# Expand/Infill lake & ponds by prescribed rates
# ----------------------------------------------------
lake_pond_expansion.lake_pond_expansion(self, element)
lake_pond_expansion.pond_infill(self, element, time)
# ----------------------------------------------------------
# Set active layer depth
# ---------------------------------------------------------
active_layer_depth.active_layer_depth(self, time, element)
# ----------------------------------
# Cycle through terrestrial cohorts
# ----------------------------------
check_Meadow_WT.check_Meadow_WT_Y(self, element, time)
check_Meadow_WT.check_Meadow_WT_M(self, element, time)
check_Meadow_WT.check_Meadow_WT_O(self, element, time)
check_LCP_WT.check_LCP_WT_Y(self, element, time)
check_LCP_WT.check_LCP_WT_M(self, element, time)
check_LCP_WT.check_LCP_WT_O(self, element, time)
check_CLC_WT.check_CLC_WT_Y(self, element, time)
check_CLC_WT.check_CLC_WT_M(self, element, time)
check_CLC_WT.check_CLC_WT_O(self, element, time)
check_FCP_WT.check_FCP_WT_Y(self, element, time)
check_FCP_WT.check_FCP_WT_M(self, element, time)
check_FCP_WT.check_FCP_WT_O(self, element, time)
check_HCP_WT.check_HCP_WT_Y(self, element, time)
check_HCP_WT.check_HCP_WT_M(self, element, time)
check_HCP_WT.check_HCP_WT_O(self, element, time)
#=====================================================
# NOTE: 17 Oct 2016. The following are place holders
# until we figure out if/how these sets of cohorts
# can transition into other cohorts.
#-----------------------------------------------------
#check_CoastalWaters_WT.check_CoastalWaters_WT_O(self, element, time)
#check_DrainedSlope_WT.check_DrainedSlope_WT_Y(self, element, time)
#check_DrainedSlope_WT.check_DrainedSlope_WT_M(self, element, time)
#check_DrainedSlope_WT.check_DrainedSlope_WT_O(self, element, time)
#check_NoData_WT.check_NoData_WT_O(self, element, time)
#check_SandDunes_WT.check_SandDunes_WT_Y(self, element, time)
#check_SandDunes_WT.check_SandDunes_WT_M(self, element, time)
#check_SandDunes_WT.check_SandDunes_WT_O(self, element, time)
#check_SaturatedBarrens_WT.check_SaturatedBarrens_WT_Y(self, element, time)
#check_SaturatedBarrens_WT.check_SaturatedBarrens_WT_M(self, element, time)
#check_SaturatedBarrens_WT.check_SaturatedBarrens_WT_O(self, element, time)
#check_Shrubs_WT.check_Shrubs_WT_O(self, element, time)
#check_Urban_WT.check_Urban_WT(self, element, time)
#====================================================
# Note: 17 Oct 2016. The following checks are pretty much obsolete
# at this point. Will clean up once everything is working well.
# ----------------------------------------------------
# check_Wet_NPG.check_Wet_NPG(self, element, time)
# check_Wet_LCP.check_Wet_LCP(self, element, time)
# check_Wet_CLC.check_Wet_CLC(self, element, time)
# check_Wet_FCP.check_Wet_FCP(self, element, time)
# check_Wet_HCP.check_Wet_HCP(self, element, time)
#=====================================================
# ----------------------------------
# Set pond/lake ice thickness depth
# ----------------------------------
ice_thickness.ice_thickness(self, time, element)
# ------------------------------
# Cycle through ponds and lakes
# ------------------------------
check_Ponds_WT.check_Ponds_WT_Y(self, element, time)
check_Ponds_WT.check_Ponds_WT_M(self, element, time)
check_Ponds_WT.check_Ponds_WT_O(self, element, time)
check_Lakes_WT.check_LargeLakes_WT_Y(self, element, time)
check_Lakes_WT.check_LargeLakes_WT_M(self, element, time)
check_Lakes_WT.check_LargeLakes_WT_O(self, element, time)
check_Lakes_WT.check_MediumLakes_WT_Y(self, element, time)
check_Lakes_WT.check_MediumLakes_WT_M(self, element, time)
check_Lakes_WT.check_MediumLakes_WT_O(self, element, time)
check_Lakes_WT.check_SmallLakes_WT_Y(self, element, time)
check_Lakes_WT.check_SmallLakes_WT_M(self, element, time)
check_Lakes_WT.check_SmallLakes_WT_O(self, element, time)
# -------------------------------------------------
# Cohort Fraction Check (mass balance of cohorts)
# -------------------------------------------------
cohort_check.cohort_check(self, element, time, cohort_start)
if time == self.stop - 1:
cohorts.final_barrow(self)
# ========================================================================
# END MAIN LOOP
# ========================================================================
# ========================================================================
# OUTPUT RESULTS (if requested)
# ========================================================================
# - - - - - - - - -
# Fractional Areas
# - - - - - - - - -
Output_cohorts_by_year.Output_cohorts_by_year(self, time)
# - - - - - - - - - - - - -
# Dominant Fractional Area
# - - - - - - - - - - - - -
Output_cohorts_by_year.dominant_cohort(self) # Terrestrial_Control
Output_cohorts_by_year.dominant_fractional_plot(
self, time) # Terrestrial_Control
# =================================
# OUTPUT ANIMATIONS (if requested)
# =================================
# - - - - - - - - - - - - - - -
# Fractional Area of Cohorts
# - - - - - - - - - - - - - - - -
Output_cohorts_by_year.write_Fractions_avi(self)
Output_cohorts_by_year.write_Dominant_Cohort_avi(
self) # Terrestrial_Control