params.h

基于格子Boltzmann方法开源可视化软件的源代码 具有很高的实用价值。对学习LBM方法及其软件开发非常游泳

  {
  fprintf( o, "rho_A[1]             %s\n", "--"                          );
  fprintf( o, "rho_B[1]             %s\n", "--"                          );
  }
  else
  {
    printf(
      "read_params() -- "
      "Unhandled case "
      "NUM_FLUID_COMPONENTS = %d .  "
      "Exiting!\n", 
      NUM_FLUID_COMPONENTS);
    exit(1);
  }
#if INAMURO_SIGMA_COMPONENT
  fprintf( o, "rho_sigma            %f\n", lattice->param.rho_sigma      );
  fprintf( o, "rho_sigma_in         %f\n", lattice->param.rho_sigma_in   );
  fprintf( o, "rho_sigma_out        %f\n", lattice->param.rho_sigma_out  );
  fprintf( o, "u_sigma              %f\n", lattice->param.u_sigma        );
  fprintf( o, "u_sigma_in           %f\n", lattice->param.u_sigma_in     );
  fprintf( o, "u_sigma_out          %f\n", lattice->param.u_sigma_out    );
  fprintf( o, "sigma_start          %d\n", lattice->param.sigma_start    );
  fprintf( o, "sigma_stop           %d\n", lattice->param.sigma_stop     );
  fprintf( o, "sigma_btc_rate       %d\n", lattice->param.sigma_btc_rate );
  fprintf( o, "sigma_btc_spot       %d\n", lattice->param.sigma_btc_spot );
#else /* !( INAMURO_SIGMA_COMPONENT) */
  fprintf( o, "rho_sigma            %s\n", "--"                          );
  fprintf( o, "rho_sigma_in         %s\n", "--"                          );
  fprintf( o, "rho_sigma_out        %s\n", "--"                          );
  fprintf( o, "u_sigma              %s\n", "--"                          );
  fprintf( o, "u_sigma_in           %s\n", "--"                          );
  fprintf( o, "u_sigma_out          %s\n", "--"                          );
  fprintf( o, "sigma_start          %s\n", "--"                          );
  fprintf( o, "sigma_stop           %s\n", "--"                          );
  fprintf( o, "sigma_btc_rate       %s\n", "--"                          );
  fprintf( o, "sigma_btc_spot       %s\n", "--"                          );
#endif /* INAMURO_SIGMA_COMPONENT */
  fprintf( o, "rho_in               %f\n", lattice->param.rho_in         );
  fprintf( o, "rho_out              %f\n", lattice->param.rho_out        );
  fprintf( o, "ux_in                %f\n", lattice->param.ux_in          );
  fprintf( o, "ux_out               %f\n", lattice->param.ux_out         );
  fprintf( o, "uy_in                %f\n", lattice->param.uy_in          );
  fprintf( o, "uy_out               %f\n", lattice->param.uy_out         );
  fprintf( o, "G               %f\n", lattice->param.G         );
  fprintf( o, "Gads[0]      %f\n", lattice->param.Gads[0]);
  if( NUM_FLUID_COMPONENTS==2)
  {
  fprintf( o, "Gads[1]      %f\n", lattice->param.Gads[1]);
  }
  else if( NUM_FLUID_COMPONENTS==1)
  {
  fprintf( o, "Gads[1]      %s\n", "--"                          );
  }
  else
  {
    printf(
      "read_params() -- "
      "Unhandled case "
      "NUM_FLUID_COMPONENTS = %d .  "
      "Exiting!\n", 
      NUM_FLUID_COMPONENTS);
    exit(1);
  }
  fprintf( o, "periodic_x[0]        %d\n", lattice->periodic_x[0]        );
  fprintf( o, "periodic_y[0]        %d\n", lattice->periodic_y[0]        );
  if( NUM_FLUID_COMPONENTS==2)
  {
  fprintf( o, "periodic_x[1]        %d\n", lattice->periodic_x[1]        );
  fprintf( o, "periodic_y[1]        %d\n", lattice->periodic_y[1]        );
  }
  else if( NUM_FLUID_COMPONENTS==1)
  {
  fprintf( o, "periodic_x[1]        %s\n", "--"                          );
  fprintf( o, "periodic_y[1]        %s\n", "--"                          );
  }
  else
  {
    printf(
      "read_params() -- "
      "Unhandled case "
      "NUM_FLUID_COMPONENTS = %d .  "
      "Exiting!\n", 
      NUM_FLUID_COMPONENTS);
    exit(1);
  }
  fprintf( o, "ns                   %d\n", lattice->param.ns_flag        );
  fprintf( o, "ns                   %f\n", lattice->param.ns             );
  fprintf( o, "slice_x              %d\n", lattice->param.slice_x        );
  fprintf( o, "slice_y              %d\n", lattice->param.slice_y        );
  fprintf( o, "ic_poisseuille       %d\n", lattice->param.ic_poisseuille );
  fprintf( o, "bc_poisseuille       %d\n", lattice->param.bc_poisseuille );
  fprintf( o, "bc_slip_north        %d\n", lattice->param.bc_slip_north  );
#if INAMURO_SIGMA_COMPONENT
  fprintf( o, "bc_sigma_slip        %d\n", lattice->param.bc_sigma_slip  );
  fprintf( o, "bc_sigma_walls       %d\n", lattice->param.bc_sigma_walls );
#else /* !( INAMURO_SIGMA_COMPONENT) */
  fprintf( o, "bc_sigma_slip        %s\n", "--"                          );
  fprintf( o, "bc_sigma_walls       %s\n", "--"                          );
#endif /* INAMURO_SIGMA_COMPONENT */
  fprintf( o, "pressure_n_in[0]     %d\n", lattice->param.pressure_n_in[0]  );
  fprintf( o, "pressure_s_in[0]     %d\n", lattice->param.pressure_s_in[0]  );
  fprintf( o, "pressure_n_out[0]    %d\n", lattice->param.pressure_n_out[0] );
  fprintf( o, "pressure_s_out[0]    %d\n", lattice->param.pressure_s_out[0] );
  fprintf( o, "velocity_n_in[0]     %d\n", lattice->param.velocity_n_in[0]  );
  fprintf( o, "velocity_s_in[0]     %d\n", lattice->param.velocity_s_in[0]  );
  fprintf( o, "velocity_n_out[0]    %d\n", lattice->param.velocity_n_out[0] );
  fprintf( o, "velocity_s_out[0]    %d\n", lattice->param.velocity_s_out[0] );
  fprintf( o, "pressure_e_in[0]     %d\n", lattice->param.pressure_e_in[0]  );
  fprintf( o, "pressure_w_in[0]     %d\n", lattice->param.pressure_w_in[0]  );
  fprintf( o, "pressure_e_out[0]    %d\n", lattice->param.pressure_e_out[0] );
  fprintf( o, "pressure_w_out[0]    %d\n", lattice->param.pressure_w_out[0] );
  fprintf( o, "velocity_e_in[0]     %d\n", lattice->param.velocity_e_in[0]  );
  fprintf( o, "velocity_w_in[0]     %d\n", lattice->param.velocity_w_in[0]  );
  fprintf( o, "velocity_e_out[0]    %d\n", lattice->param.velocity_e_out[0] );
  fprintf( o, "velocity_w_out[0]    %d\n", lattice->param.velocity_w_out[0] );
  if( NUM_FLUID_COMPONENTS==2)
  {
  fprintf( o, "pressure_n_in[1]     %d\n", lattice->param.pressure_n_in[1]  );
  fprintf( o, "pressure_s_in[1]     %d\n", lattice->param.pressure_s_in[1]  );
  fprintf( o, "pressure_n_out[1]    %d\n", lattice->param.pressure_n_out[1] );
  fprintf( o, "pressure_s_out[1]    %d\n", lattice->param.pressure_s_out[1] );
  fprintf( o, "velocity_n_in[1]     %d\n", lattice->param.velocity_n_in[1]  );
  fprintf( o, "velocity_s_in[1]     %d\n", lattice->param.velocity_s_in[1]  );
  fprintf( o, "velocity_n_out[1]    %d\n", lattice->param.velocity_n_out[1] );
  fprintf( o, "velocity_s_out[1]    %d\n", lattice->param.velocity_s_out[1] );
  fprintf( o, "pressure_e_in[1]     %d\n", lattice->param.pressure_e_in[1]  );
  fprintf( o, "pressure_w_in[1]     %d\n", lattice->param.pressure_w_in[1]  );
  fprintf( o, "pressure_e_out[1]    %d\n", lattice->param.pressure_e_out[1] );
  fprintf( o, "pressure_w_out[1]    %d\n", lattice->param.pressure_w_out[1] );
  fprintf( o, "velocity_e_in[1]     %d\n", lattice->param.velocity_e_in[1]  );
  fprintf( o, "velocity_w_in[1]     %d\n", lattice->param.velocity_w_in[1]  );
  fprintf( o, "velocity_e_out[1]    %d\n", lattice->param.velocity_e_out[1] );
  fprintf( o, "velocity_w_out[1]    %d\n", lattice->param.velocity_w_out[1] );
  }
  else if( NUM_FLUID_COMPONENTS==1)
  {
  fprintf( o, "pressure_n_in[1]     %s\n", "--"                             );
  fprintf( o, "pressure_s_in[1]     %s\n", "--"                             );
  fprintf( o, "pressure_n_out[1]    %s\n", "--"                             );
  fprintf( o, "pressure_s_out[1]    %s\n", "--"                             );
  fprintf( o, "velocity_n_in[1]     %s\n", "--"                             );
  fprintf( o, "velocity_s_in[1]     %s\n", "--"                             );
  fprintf( o, "velocity_n_out[1]    %s\n", "--"                             );
  fprintf( o, "velocity_s_out[1]    %s\n", "--"                             );
  fprintf( o, "pressure_e_in[1]     %s\n", "--"                             );
  fprintf( o, "pressure_w_in[1]     %s\n", "--"                             );
  fprintf( o, "pressure_e_out[1]    %s\n", "--"                             );
  fprintf( o, "pressure_w_out[1]    %s\n", "--"                             );
  fprintf( o, "velocity_e_in[1]     %s\n", "--"                             );
  fprintf( o, "velocity_w_in[1]     %s\n", "--"                             );
  fprintf( o, "velocity_e_out[1]    %s\n", "--"                             );
  fprintf( o, "velocity_w_out[1]    %s\n", "--"                             );
  }
  else
  {
    printf(
      "read_params() -- "
      "Unhandled case "
      "NUM_FLUID_COMPONENTS = %d .  "
      "Exiting!\n", 
      NUM_FLUID_COMPONENTS);
    exit(1);
  }
  if( INAMURO_SIGMA_COMPONENT)
  {
  fprintf( o, "constcon_n_in        %d\n", lattice->param.constcon_n_in  );
  fprintf( o, "constcon_s_in        %d\n", lattice->param.constcon_s_in  );
  fprintf( o, "constcon_n_out       %d\n", lattice->param.constcon_n_out );
  fprintf( o, "constcon_s_out       %d\n", lattice->param.constcon_s_out );
  fprintf( o, "constflx_n_in        %d\n", lattice->param.constflx_n_in  );
  fprintf( o, "constflx_s_in        %d\n", lattice->param.constflx_s_in  );
  fprintf( o, "constflx_n_out       %d\n", lattice->param.constflx_n_out );
  fprintf( o, "constflx_s_out       %d\n", lattice->param.constflx_s_out );
  fprintf( o, "constcon_n_in        %d\n", lattice->param.constcon_e_in  );
  fprintf( o, "constcon_s_in        %d\n", lattice->param.constcon_w_in  );
  fprintf( o, "constcon_n_out       %d\n", lattice->param.constcon_e_out );
  fprintf( o, "constcon_s_out       %d\n", lattice->param.constcon_w_out );
  fprintf( o, "constflx_n_in        %d\n", lattice->param.constflx_e_in  );
  fprintf( o, "constflx_s_in        %d\n", lattice->param.constflx_w_in  );
  fprintf( o, "constflx_n_out       %d\n", lattice->param.constflx_e_out );
  fprintf( o, "constflx_s_out       %d\n", lattice->param.constflx_w_out );
  fprintf( o, "zeroconcgrad_n       %d\n", lattice->param.zeroconcgrad_n );
  fprintf( o, "zeroconcgrad_s       %d\n", lattice->param.zeroconcgrad_s );
  fprintf( o, "zeroconcgrad_n       %d\n", lattice->param.zeroconcgrad_e );
  fprintf( o, "zeroconcgrad_s       %d\n", lattice->param.zeroconcgrad_w );
  fprintf( o, "zeroconcgrad_full    %d\n", lattice->param.zeroconcgrad_full );
  }
  else
  {
  fprintf( o, "constcon_n_in        %s\n", "--"                             );
  fprintf( o, "constcon_s_in        %s\n", "--"                             );
  fprintf( o, "constcon_n_out       %s\n", "--"                             );
  fprintf( o, "constcon_s_out       %s\n", "--"                             );
  fprintf( o, "constflx_n_in        %s\n", "--"                             );
  fprintf( o, "constflx_s_in        %s\n", "--"                             );
  fprintf( o, "constflx_n_out       %s\n", "--"                             );
  fprintf( o, "constflx_s_out       %s\n", "--"                             );
  fprintf( o, "constcon_e_in        %s\n", "--"                             );
  fprintf( o, "constcon_w_in        %s\n", "--"                             );
  fprintf( o, "constcon_e_out       %s\n", "--"                             );
  fprintf( o, "constcon_w_out       %s\n", "--"                             );
  fprintf( o, "constflx_e_in        %s\n", "--"                             );
  fprintf( o, "constflx_w_in        %s\n", "--"                             );
  fprintf( o, "constflx_e_out       %s\n", "--"                             );
  fprintf( o, "constflx_w_out       %s\n", "--"                             );
  fprintf( o, "zeroconcgrad_n       %s\n", "--"                             );
  fprintf( o, "zeroconcgrad_s       %s\n", "--"                             );
  fprintf( o, "zeroconcgrad_e       %s\n", "--"                             );
  fprintf( o, "zeroconcgrad_w       %s\n", "--"                             );
  fprintf( o, "zeroconcgrad_full    %s\n", "--"                             );
  }

  fprintf( o, "plot_scale_dynamic   %d\n", lattice->param.plot_scale_dynamic);
  fprintf( o, "use_colormap         %d\n", lattice->param.use_colormap      );
  fprintf( o, "initial_condition    %d\n", lattice->param.initial_condition );
  fprintf( o, "x0                   %f\n", lattice->param.x0                );
  fprintf( o, "y0                   %f\n", lattice->param.y0                );
  fprintf( o, "r0                   %f\n", lattice->param.r0                );
  fprintf( o, "cut                  %f\n", lattice->param.cut               );
  fprintf( o, "x1                   %f\n", lattice->param.x1                );
  fprintf( o, "x2                   %f\n", lattice->param.x2                );
  fprintf( o, "y1                   %f\n", lattice->param.y1                );
  fprintf( o, "y2                   %f\n", lattice->param.y2                );
  fprintf( o, "rel_x1               %f\n", lattice->param.rel_x1            );
  fprintf( o, "rel_x2               %f\n", lattice->param.rel_x2            );
  fprintf( o, "rel_y1               %f\n", lattice->param.rel_y1            );
  fprintf( o, "rel_y2               %f\n", lattice->param.rel_y2            );
  fprintf( o, "dump_rho             %d\n", lattice->param.dump_rho          );
  fprintf( o, "dump_u               %d\n", lattice->param.dump_u            );
  fprintf( o, "dump_force           %d\n", lattice->param.dump_force        );
  fprintf( o, "dump_vor             %d\n", lattice->param.dump_vor          );

  fclose(o);

#if VERBOSITY_LEVEL >= 1
  printf("lbio.c: dump_params() -- Wrote file \"params.dat\".\n");
#endif /* VERBOSITY_LEVEL >= 1 */

} /* void dump_params( struct lattice_struct *lattice) */