Albert's method ported to IDL

works.
Added parapeak, will tweak next time