- rh1 and rh2 are initialisations of N1 and N3
- we use rather init_N1 and init_N3 params of discrete-time