- Timestamp:
- Sep 4, 2007, 5:36:01 PM (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_core/source/anuga/shallow_water/shallow_water_domain.py
r4687 r4701 173 173 self.forcing_terms.append(gravity) 174 174 175 # Stored output175 # Stored output 176 176 self.store = True 177 177 self.format = 'sww' … … 179 179 self.minimum_storable_height = minimum_storable_height 180 180 self.quantities_to_be_stored = ['stage','xmomentum','ymomentum'] 181 181 182 182 183 … … 234 235 235 236 #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. 239 241 self.minimum_allowed_height = minimum_allowed_height 240 242 self.H0 = minimum_allowed_height … … 262 264 self.points_file_block_line_size = points_file_block_line_size 263 265 266 264 267 def set_quantities_to_be_stored(self, q): 265 268 """Specify which quantities will be stored in the sww file. … … 286 289 q = [q] # Turn argument into a list 287 290 288 # Check correcness291 # Check correcness 289 292 for quantity_name in q: 290 293 msg = 'Quantity %s is not a valid conserved quantity'\ … … 930 933 #MH090605 if second order, 931 934 #perform the extrapolation and limiting on 932 #all of the conserved quantitie 935 #all of the conserved quantities 933 936 934 937 if (domain._order_ == 1): … … 944 947 for name in domain.conserved_quantities: 945 948 Q = domain.quantities[name] 949 946 950 if domain._order_ == 1: 947 951 Q.extrapolate_first_order() 948 952 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 949 965 Q.extrapolate_second_order() 950 966 Q.limit()
Note: See TracChangeset
for help on using the changeset viewer.