Ignore:
Timestamp:
Oct 11, 2008, 12:38:18 PM (16 years ago)
Author:
steve
Message:
 
File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/development/anuga_1d/comp_flux_ext.c

    r5831 r5832  
    188188                                // Update timestep based on edge i and possibly neighbour n
    189189                                if (max_speed > epsilon) {
    190                         // Original CFL calculation
    191                                        
    192                                         timestep = min(timestep, 0.5*cfl*areas[k]/max_speed); //Here, CFL=1.0 is assumed. ?????????????????????????????????????????????
    193                                         if (n>=0) {
    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                                    }
    196196                                }
    197197            } // End edge i (and neighbour n)
     
    263263  // Call underlying flux computation routine and update
    264264  // the explicit update arrays
    265   timestep = _compute_fluxes_ext(
    266                  cfl,
     265  timestep = _compute_fluxes_ext(cfl,
    267266                                 timestep,
    268267                                 epsilon,
     
    329328    g                 = get_python_double(domain,"g");
    330329    h0                = get_python_double(domain,"h0");
    331         cfl                = get_python_double(domain,"cfl");
     330    cfl               = get_python_double(domain,"CFL");
    332331 
    333332   
     
    358357    // the explicit update arrays
    359358    timestep = _compute_fluxes_ext(
    360                        cfl,
    361                                    timestep,
    362                                    epsilon,
    363                                    g,
    364                                    h0,
    365                                    (long*) neighbours -> data,
    366                                    (long*) neighbour_vertices -> data,
    367                                    (double*) normals -> data,
    368                                    (double*) areas -> data,
    369                                    (double*) stage_vertex_values -> data,
    370                                    (double*) xmom_vertex_values -> data,
    371                                    (double*) bed_vertex_values -> data,
    372                                    (double*) stage_boundary_values -> data,
    373                                    (double*) xmom_boundary_values -> data,
    374                                    (double*) stage_explicit_update -> data,
    375                                    (double*) xmom_explicit_update -> data,
    376                                    number_of_elements,
    377                                    (double*) max_speed_array -> data);
     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);
    378377
    379378
Note: See TracChangeset for help on using the changeset viewer.