Changeset 6459
- Timestamp:
- Mar 5, 2009, 9:45:32 AM (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/australia_ph2/dampier/run_model.py
r6445 r6459 141 141 # Prepare time boundary 142 142 TB.prepare_timeboundary(project.boundary_csv) 143 f = file_function(project.boundary_csv[:-4] + '.tms' ,domain)143 f = file_function(project.boundary_csv[:-4] + '.tms') 144 144 145 145 Br = Reflective_boundary(domain) 146 146 Bt = Transmissive_stage_zero_momentum_boundary(domain) 147 147 Bd = Dirichlet_boundary([project.tide, 0, 0]) 148 Bf = Field_boundary(project.event_sts+'.sts',149 domain, mean_stage=project.tide,150 time_thinning=1,151 default_boundary=Bd,152 boundary_polygon=bounding_polygon_sts,153 use_cache=True,154 verbose=True)155 Tb = Time_boundary(domain,f)156 148 157 domain.set_boundary({'back': Br, 149 if project.wave == 'Bf': 150 Bf = Field_boundary(project.event_sts+'.sts', 151 domain, mean_stage=project.tide, 152 time_thinning=1, 153 default_boundary=Bd, 154 boundary_polygon=bounding_polygon_sts, 155 use_cache=True, 156 verbose=True) 157 domain.set_boundary({'back': Br, 158 158 'side': Bd, 159 'ocean': Tb}) 159 'ocean': Bf}) 160 161 elif project.wave == 'Tb': 162 Tb = Time_boundary(domain,f,default_boundary=Bd ) 163 164 domain.set_boundary({'back': Br, 165 'side': Bd, 166 'ocean': Tb}) 167 else: 168 print 'No wave specified in project script (Bf or Tb)' 169 160 170 161 171 #-------------------------------------------------------------------------------
Note: See TracChangeset
for help on using the changeset viewer.