SAR calculation confuses me.

[h3lg3]

Hey!
Following my bug reports #221 and #223, I am still confused about the SAR calculation and wonder if it needs to be fixed. Generally I observed, that SAR changes throughout the sequence, even though the sequence is highly repetitive. I am not an expert on SAR calculations, so please correct me if I am wrong.
In the first plot you can see the evolution of SAR for the write_epi.py example sequence if Nx and Ny are changed to 256 and n_slices is changed to 40. The second plot shows a zoomed version of the original SAR and time vectors together with its interpolated values to 1 second increments.
Isn't the interpolation completely off? I wonder if you could clarify the issue.
Thanks a lot
Helge

https://github.com/imr-framework/pypulseq/blob/dev/examples/scripts/write_epi.py

SAR calculation for whole sequence

Actual SAR and interpolated SAR