       { gamma_s.alpha() = -0.5*p.tau[j]*p.mu[j];
         invert_invgamma(gamma_s, gamma_n, ntilde[j]);
         dg::blas1::axpby(p.a[j]*p.mu[j], ntilde[j], 0., chi);
-        dg::blas1::plus( chi, p.a[j]*p.mu[j]); 
+        dg::blas1::plus( chi, p.a[j]*p.mu[j]);
         dg::blas1::pointwiseDot(chi, binv, chi);
         dg::blas1::pointwiseDot(chi, binv, chi);
         pol.set_chi(chi);

  err_out = nc_def_grp( ncid_out, "err_etime", &err_etime_wgrp_id);
  for (unsigned i = 0 ; i < num_err_etime; i++)
  { err_in = nc_inq_varid(ncid_in, err_etime[i].data(), &err_etime_rval_id[i]);
    err_out = nc_def_var( err_etime_wgrp_id, err_etime[i].data(), NC_DOUBLE, 1, &dim_et_id, &err_etime_wval_id[i]);
    err_out = nc_def_var_deflate( err_etime_wgrp_id, err_etime_wval_id[i], shuffle_flag, compress_flag, compress_level);
  }
  err_out = nc_enddef(ncid_out);
  ////////remap data////////
  //err_etime
  double transfer_etime[num_etime];
  static size_t count_etime[] = {num_etime};
  static size_t start_etime[] = {0};
  for (unsigned i = 0; i < num_err_etime; i++)
  { err_in = nc_get_vara_double(ncid_in, err_etime_rval_id[i], start_etime, count_etime, transfer_etime);
    err_out = nc_put_vara_double(err_etime_wgrp_id, err_etime_wval_id[i], start_etime, count_etime, transfer_etime);
  }
  //timestepping
  int add_to_spatial_idx = 0;
  for (unsigned i = 0; i <= p.maxout; i++)

  reader.parse(input, js, false);
  const imp::Parameters p(js);
  p.display(std::cout);
  // dimensions
  //.nc parameter: etime not found in .json file
  int etdim_r_id, xdim_r_id, ydim_r_id, tdim_r_id;
  size_t num_etime, num_x, num_y, num_time;
  err_in = nc_inq_dimid(ncid_in, "energy_time", &etdim_r_id);
  err_in = nc_inq_dimlen(ncid_in, etdim_r_id, &num_etime);
  err_in = nc_inq_dimid(ncid_in, "x", &xdim_r_id);
  err_in = nc_inq_dimlen(ncid_in, xdim_r_id, &num_x);
  err_in = nc_inq_dimid(ncid_in, "y", &ydim_r_id);
  err_in = nc_inq_dimlen(ncid_in, ydim_r_id, &num_y);
  err_in = nc_inq_dimid(ncid_in, "time", &tdim_r_id);
  err_in = nc_inq_dimlen(ncid_in, tdim_r_id, &num_time);
  ////////compose data////////
  const size_t num_species = 3, num_fields = 2;
  const size_t num_spatial = 14;
