!
! gaussclumps.INIT
TASK\FILE       "input file                   " INP_IMAGE$ = cube.gdf
TASK\REAL       "rms [K]                      " RMS$ = 1.
TASK\CHARACTER  "beam FWHM [m] or [s]         " FWHM_BEAM$ = 13S
TASK\REAL       "initial guess [beam FWHM]    " FWHM_START$ = 1.5
TASK\REAL       "velocity resolution [km/s]   " VELO_RES$ = .5
TASK\REAL       "initial guess [VELO_RES]     " VELO_START$ =  2.
TASK\REAL       "stiffness S0[3]              " STIFFNESS$[3] = 1. 1. 1.
TASK\REAL       "threshhold [K]               " THRESHHOLD$ = 1.
TASK\INTEGER    "max number of clumps         " MAXCLUMP$ = 5000
TASK\REAL       "contrast [fraction of maximum]  " CONTRAST$ = -10.
TASK\REAL       "minimum integral [fration of total]  " MININTEGRAL$ = 0.
TASK\REAL       "aperture FWHM [initial guess]" APERTURE_FWHM$ = 3.
TASK\REAL       "aperture cutoff [initial guess]       " APERTURE_LMTS$ = 6.
TASK\LOGICAL    "fit DV/DR [y/n]              " FIT_DVDR$ = .TRUE.
TASK\LOGICAL    "fit center position [y/n]    " FIT_CENTERPOS$ = .TRUE.
TASK\LOGICAL    "fit amplitude [y/n]          " FIT_AMPLITUDE$ = .TRUE.
TASK\INTEGER    "maximum number of iterations " MAXIT$ = 100
TASK\REAL       "epsilon                      " EPSILON$ =  0.01
TASK\REAL       "flambda                      " FLAMBDA$ =  0.001
TASK\LOGICAL    "update output images         " UPD_IMAGE$ = .false.
TASK\FILE       "RESIDUAL IMAGE               " RES_IMAGE$ = cube_res.gdf
TASK\FILE       "FIT IMAGE                    " FIT_IMAGE$ = cube_fit.gdf
TASK\GO
