; get some continuum data
 cont
 filein,'peaks.fits'
 get,index=1
 ; setup the gauss fit
 gregion,[60,80]
 ngauss,[1]
 gmaxiter,500
 gparamvalues, 0, 0, [400000.,70.,100.]
 ; find and show the fit
 gauss
 gshow