// Synchronize J evolution with the rotor cycle

if (srate > 500) then 
   t1YXinept = trunc(t1YXinept*1.0e-6*srate + 0.5)/(1.0e-6*srate)  
   t2YXinept = trunc(t2YXinept*1.0e-6*srate + 0.5)/(1.0e-6*srate) 
endif   












