Subprogram move_marker_on_characteristic (marker, dIperp, Blocal, coeff)

Change the kick in vpar to stay exactly on characteristics.

Separate the estimate of the total power absorption used in the E-field normalisation, which should be set to zero at every renormalisation of the E-field, and the time-cumlative diagnostics.

What should the bounding bosx be used for, only checking if points are realistic, or should it also be used for the diagnostic grid.

Add 2D grid as part of input-namelist

Rebuild input to be come in one single structure, e.g. RFOF_struct

Write alternative resonance condition routines that old checks if the particle has crossed the resonance (no extrapolation to the upcoming resonace)

Test if the available techniques are quick enough?

Feed suggested

Read namelist in contructor_RFOF_cumlative_diagnostics

Subprogram RFOF_diagnostics::contructor_RFOF_cumlative_diagnostics (diagno, nfreq, nnphi)

RZ-boundaing box as input

rho-theta diagnostics - box interpolation for speed

rho diagnostics - take grid from ASCOT

Type RFOF_Efield_update

output E-field normalisation to file

Subprogram RFOF_Efield_update::RFOF_merge_absorption_diagno_from_mpi_nodes (diagno, ierr)

Handle errors!!

Subprogram RFOF_Efield_update::RFOF_merge_power_absorption_diagno_from_mpi_nodes (pabs, pabs2, count, ierr)

Handle errors!!

Subprogram RFOF_kick::quasilinear_RF_kick_steinbrecher_integrator (marker, mem, Blocal, RFlocal, RFglobal, diagno, MPI_node_Id)

Improve the Monte Carlo stepping scheme to something adaptive. Introduce a check to see that the RF step has given a physical orbit.

What to do with the diffusion and drift coefficients going into the diagnostics?

Type RFOF_main

The normalisation of the radial electric is still not implemented. This will require communication between MPI-nodes. A suggestion for treating this could be for the user (the person coupling to their Orbit code) to implement a wrapper for the nessecary MPI routines.

Subprogram RFOF_main::RFOF_master (time, time_previous_step, RFglobal, marker, resonanceMemory, diagno, MPI_nod_Id, interaction_failed_particle_overshot_resonance, time_at_resonance)

MPI nod Id as input

Subprogram RFOF_parameters::initialise_RFOF_parameters (Rmin_bound, Rmax_bound, Zmin_bound, Zmax_bound)

Put input parameters in namelist.

Subprogram RFOF_parameters::simplify__static_resonance_position_during_RF_kick

Yet to be implemented in RFOF_kick

Subprogram RFOF_parameters::width_of_rf_resonance_layer

Check this approximation - probably need better model.

Subprogram RFOF_random_numbers::gaussian_random_number ()

STILL TO BE IMPLEMENTED!

Subprogram RFOF_resonance_condition::close_enough_to_resonance (t_res, R_res, z_res, time, time_previous_step, marker)

Improve condition for what particle is "resonant". At the moment the condition measures the distance to the resonance and compares it to an number which has to be specified; there is still no physical justification for why this number has been choosen.

Type RFOF_types::RFOF_state

test this works with the structs from ASCOT!

Type RFOF_waves

Implement interpolation of wave field on grid. (1) Build wrapper module for Diercx and derived type including spline coefficients (2) Use spline-type in rf_wave_global for each quantity you wish to spline (3) Add in constructor the generation of the spline coefficients

Subprogram wrap_move_marker_on_characteristic (Iperp, state, errorFlag)

Make this into routine which is not wrapped

Subprogram wrap_RF_diffusion_and_move_marker (Iperp, state, diffusion, errorFlag)

Optimize, the present version is a quick wrapper including already wrapped routines.