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