script for plotting fnat, irms and lrms