Ignore:
Timestamp:
Sep 4, 2007, 5:36:01 PM (16 years ago)
Author:
ole
Message:

Work towards ticket:192

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_core/source/anuga/shallow_water/shallow_water_domain.py

    r4687 r4701  
    173173        self.forcing_terms.append(gravity)
    174174
    175         #Stored output
     175        # Stored output
    176176        self.store = True
    177177        self.format = 'sww'
     
    179179        self.minimum_storable_height = minimum_storable_height
    180180        self.quantities_to_be_stored = ['stage','xmomentum','ymomentum']
     181       
    181182               
    182183
     
    234235
    235236        #FIXME (Ole): rename H0 to minimum_allowed_height_in_flux_computation
    236         #rename tight_slope_limiters to tight_slope_limiters.
    237         #Maybe use histogram to identify isolated extreme speeds and deal with them adaptively
    238         #similarly to how we used to use 1 order steps to recover.
     237
     238        #FIXME (Ole): Maybe use histogram to identify isolated extreme speeds
     239        #and deal with them adaptively similarly to how we used to use 1 order
     240        #steps to recover.
    239241        self.minimum_allowed_height = minimum_allowed_height
    240242        self.H0 = minimum_allowed_height   
     
    262264        self.points_file_block_line_size = points_file_block_line_size
    263265       
     266       
    264267    def set_quantities_to_be_stored(self, q):
    265268        """Specify which quantities will be stored in the sww file.
     
    286289            q = [q] # Turn argument into a list
    287290
    288         #Check correcness
     291        # Check correcness
    289292        for quantity_name in q:
    290293            msg = 'Quantity %s is not a valid conserved quantity'\
     
    930933        #MH090605 if second order,
    931934        #perform the extrapolation and limiting on
    932         #all of the conserved quantitie
     935        #all of the conserved quantities
    933936
    934937        if (domain._order_ == 1):
     
    944947        for name in domain.conserved_quantities:
    945948            Q = domain.quantities[name]
     949
    946950            if domain._order_ == 1:
    947951                Q.extrapolate_first_order()
    948952            elif domain._order_ == 2:
     953
     954                # Experiment
     955                #if name == 'stage':
     956                #    #print name, 'second'
     957                #    Q.extrapolate_second_order()
     958                #    Q.limit()
     959                #else:
     960                #    #print name, 'first'               
     961                #    Q.extrapolate_first_order()
     962                #    #Q.extrapolate_second_order()
     963                #    #Q.limit()               
     964               
    949965                Q.extrapolate_second_order()
    950966                Q.limit()
Note: See TracChangeset for help on using the changeset viewer.