Changeset 5832 for anuga_work/development/anuga_1d/comp_flux_ext.c
 Timestamp:
 Oct 11, 2008, 12:38:18 PM (16 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

anuga_work/development/anuga_1d/comp_flux_ext.c
r5831 r5832 188 188 // Update timestep based on edge i and possibly neighbour n 189 189 if (max_speed > epsilon) { 190 191 192 timestep = min(timestep, 0.5*cfl*areas[k]/max_speed); //Here, CFL=1.0 is assumed. ?????????????????????????????????????????????193 194 timestep = min(timestep, 0.5*cfl*areas[n]/max_speed); //Here, CFL=1.0 is assumed. ?????????????????????????????????????????????195 190 // Original CFL calculation 191 192 timestep = min(timestep, 0.5*cfl*areas[k]/max_speed); 193 if (n>=0) { 194 timestep = min(timestep, 0.5*cfl*areas[n]/max_speed); 195 } 196 196 } 197 197 } // End edge i (and neighbour n) … … 263 263 // Call underlying flux computation routine and update 264 264 // the explicit update arrays 265 timestep = _compute_fluxes_ext( 266 cfl, 265 timestep = _compute_fluxes_ext(cfl, 267 266 timestep, 268 267 epsilon, … … 329 328 g = get_python_double(domain,"g"); 330 329 h0 = get_python_double(domain,"h0"); 331 cfl = get_python_double(domain,"cfl");330 cfl = get_python_double(domain,"CFL"); 332 331 333 332 … … 358 357 // the explicit update arrays 359 358 timestep = _compute_fluxes_ext( 360 361 362 363 364 365 366 367 368 369 370 371 372 373 374 375 376 377 359 cfl, 360 timestep, 361 epsilon, 362 g, 363 h0, 364 (long*) neighbours > data, 365 (long*) neighbour_vertices > data, 366 (double*) normals > data, 367 (double*) areas > data, 368 (double*) stage_vertex_values > data, 369 (double*) xmom_vertex_values > data, 370 (double*) bed_vertex_values > data, 371 (double*) stage_boundary_values > data, 372 (double*) xmom_boundary_values > data, 373 (double*) stage_explicit_update > data, 374 (double*) xmom_explicit_update > data, 375 number_of_elements, 376 (double*) max_speed_array > data); 378 377 379 378
Note: See TracChangeset
for help on using the changeset viewer.