This is a fresh repository for DoPhot ready to go as v4.0 to github.  Models for an extended pseudogaussian 'extpgauss', and a sersic are not yet ready to be implemented in full and so are not included with 4.0, but their remnants do remain in the code as a guide for others to make their own models if they so choose.
copyright added in LEGAL document
LEGAL changed to LICENSE
Updates to manual made.  Typos and such corrected.  Round 1.
parameter MAX_PERF added as an optional maximum number of perfect stars to use for computing average shape values for stars.  Only the brightest MAX_PERF stars will be used.  Default is set to 2000.
extended pseudogauss model branch created with extpgauss model added back in for debugging.  all debugging outputs are in working_data which are NOT backed up
extpgauss model now fits beta4 in log space and empirically subtracts warmstart stars.  However fit still models large objects not well fit in initial pass as many smaller objects.  This is a problem.
extpgauss now fits for beta4 and beta6, but not in log space because that didn't seem to be the fitting problem.  The log space option remains in the code (onefit, twofit, extpgauss), but commented out.  The problem of one object being fit as many still remains.

I am switching to master and having only the first NPERF_MAX objects be averaged into the star shape parameters in the hopes that having only the higher signal to noise objects determine the shape will help.
fitting has been moved back to log spcae as negative beta4 was becoming an issue.  commit done in preparation for upgrade to master.  type 9 will be checked with pgauss model.
Objects which fail on first fit of 7+ parameter model now attempt a PGAUSS model before declaring object of type 9.  Struct model and object model_.which_model switches from 0 to 1 for objects using the pgauss model rather than the default.
output files now indicate if model used is one chosen (0) or pgauss (1)
Covariance matrix problem fixed.  Covariance matrix is ocrrect ofr beta4&6 fitting.
Shape now fits the pgauss model before other models and only if pgauss converges does it use the pgauss parameters as starting parameters for the other model fits.  More objects converge on the other fits with this ordering.
Shape now only checks for the expanded model convergence once, on the first shape analysis of a given object, rather than on every shape pass.  An array in model_struct, tested, has been added to this end.  1 implies that the object has been tested for 7+ parameter model convergence. Additionally, redundancies in shape have been cleaned up and print statements have been added so the code is cleaner and the log file is easier to read.
shape commented and cleaned of obsolete TWOFPAR.  also addlims now checks for EXTPGAUSS model beta6 coming from fit param.
comments added.  no substantial changes.  time to merge
dophot now excludes isearch function on first pass rather than extra including shape and improve fits on first pass.  This allows fow paravg and varipar of warmstar parameters befoe isearch fits.  Additionally.  parupd now updates a specified number of parameters, the number of parameters updated is the last integer passed, thus allowing for the first shape fit of an object to inform the sky, intensity, and x,y values of the first fit for stars.  these values are superior to the crudestat values.  However, neither the informed average values for stars before isearch nor the improved first guess from shape before improve correct the overfit stars in the extpgauss model.  the overfitting is very robust.  lastly, the shadow parameter of the second star of a type three object is now updated in shape as well as the first star.  the lack of similarity in the update of these two objects seemed incorrect.
addstar now takes ints rather than pointers to ints in preparation for empirical and type 2 v type 1 toggle addstar wrapper
toggle between empirical, analytic, star and galaxy object addition has ben given its own routine add_analytic_or_empirical_obj which wraps addstar.  makes shape and and bestab cleaner.  improve needs to be updated still because the empirical switch is more complicated there.
missing routines from previous commit added