Index: trunk/Ohana/src/relastro/Makefile =================================================================== --- trunk/Ohana/src/relastro/Makefile (revision 39402) +++ trunk/Ohana/src/relastro/Makefile (revision 39457) @@ -41,4 +41,5 @@ $(SRC)/Shutdown.$(ARCH).o \ $(SRC)/UpdateChips.$(ARCH).o \ +$(SRC)/UpdateStacks.$(ARCH).o \ $(SRC)/UpdateMosaic.$(ARCH).o \ $(SRC)/UpdateObjects.$(ARCH).o \ @@ -47,4 +48,5 @@ $(SRC)/UpdateMeasures.$(ARCH).o \ $(SRC)/GetAstromError.$(ARCH).o \ +$(SRC)/MeasFilterTest.$(ARCH).o \ $(SRC)/args.$(ARCH).o \ $(SRC)/bcatalog.$(ARCH).o \ @@ -120,4 +122,5 @@ $(SRC)/UpdateMeasures.$(ARCH).o \ $(SRC)/GetAstromError.$(ARCH).o \ +$(SRC)/MeasFilterTest.$(ARCH).o \ $(SRC)/args.$(ARCH).o \ $(SRC)/bcatalog.$(ARCH).o \ @@ -195,6 +198,5 @@ $(SRC)/mkpolyterm.$(ARCH).o \ $(SRC)/GetAstromError.$(ARCH).o \ -$(SRC)/ImageOps.$(ARCH).o \ -$(SRC)/MosaicOps.$(ARCH).o \ +$(SRC)/MeasFilterTest.$(ARCH).o \ $(SRC)/extra.$(ARCH).o \ $(SRC)/plotstuff.$(ARCH).o \ @@ -218,6 +220,5 @@ $(SRC)/mkpolyterm.$(ARCH).o \ $(SRC)/GetAstromError.$(ARCH).o \ -$(SRC)/ImageOps.$(ARCH).o \ -$(SRC)/MosaicOps.$(ARCH).o \ +$(SRC)/MeasFilterTest.$(ARCH).o \ $(SRC)/extra.$(ARCH).o \ $(SRC)/plotstuff.$(ARCH).o \ @@ -241,9 +242,15 @@ $(SRC)/mkpolyterm.$(ARCH).o \ $(SRC)/GetAstromError.$(ARCH).o \ -$(SRC)/ImageOps.$(ARCH).o \ -$(SRC)/MosaicOps.$(ARCH).o \ +$(SRC)/MeasFilterTest.$(ARCH).o \ $(SRC)/extra.$(ARCH).o \ $(SRC)/plotstuff.$(ARCH).o \ $(SRC)/fitpoly.$(ARCH).o + +# stupid dependencies: +# UpdateObjects.c:UpdateObjects_Chips needs +# ImageOps.c:MeasFilterTiny forces +# load_images.c:free_astrom_table needs +# select_images.c +# MosaicOps.c $(FITSTK): $(INC)/relastro.h Index: trunk/Ohana/src/relastro/include/relastro.h =================================================================== --- trunk/Ohana/src/relastro/include/relastro.h (revision 39402) +++ trunk/Ohana/src/relastro/include/relastro.h (revision 39457) @@ -88,6 +88,6 @@ off_t *NMEASURE; // allocated Measures per catalog int *index; // lookup table catID -> catalog[i] - int *catIDs; // lookup table catID <- catalog[i] - int maxID; // max catID value to date + unsigned int *catIDs; // lookup table catID <- catalog[i] + unsigned int maxID; // max catID value to date } CatalogSplitter; @@ -259,4 +259,5 @@ int REGION_HOST_ID; int PARALLEL_REGIONS_MANUAL; +char *MANUAL_UNIQUER; int HOST_ID; @@ -416,5 +417,6 @@ int SetSignals (void); -void relastro_client_free (); +void relastro_free (SkyTable *sky, SkyList *skylist); +void relastro_client_free (SkyTable *sky, SkyList *skylist); GSCRegion *find_regions PROTO((Image *image, off_t Nimage, int *Nregions, GSCRegion *fullregion)); @@ -438,4 +440,5 @@ void initImageBins PROTO((Catalog *catalog, int Ncatalog, int FULLINIT)); void initImages PROTO((Image *input, off_t *line_number, off_t N)); +void freeImages PROTO(()); void initMosaicBins PROTO((Catalog *catalog, int Ncatalog)); void initMosaicGrid PROTO((Image *image, off_t Nimage)); @@ -528,4 +531,5 @@ int UpdateSimple (Catalog *catalog, int Ncatalog); int UpdateChips (Catalog *catalog, int Ncatalog, int Nloop); +int UpdateStacks (Catalog *catalog, int Ncatalog); int UpdateMosaic (Catalog *catalog, int Ncatalog); int UpdateMeasures (Catalog *catalog, int Ncatalog); @@ -552,4 +556,5 @@ Mosaic *getmosaics (off_t *N); void initMosaics (Image *image, off_t Nimage); +void freeMosaics (); StarData *getMosaicRaw (Catalog *catalog, int Ncatalog, off_t mos, off_t *Nstars); StarData *getMosaicRef (Catalog *catalog, int Ncatalog, off_t mos, off_t *Nstars); @@ -588,4 +593,5 @@ int initStarMaps (); +void freeStarMaps (); int updateStarMaps(Catalog *catalog); int createStarMapPoints(); @@ -653,5 +659,6 @@ int indexCatalogs (Catalog *catalog, int Ncatalog); -int catID_and_objID_to_seq (int catID, int objID, int *catSeq, off_t *objSeq); +int catID_and_objID_to_seq (unsigned int catID, unsigned int objID, int *catSeq, off_t *objSeq); +void freeCatalogIndexes (int Ncatalog); int markObjects (Catalog *catalog, int Ncatalog); @@ -682,4 +689,5 @@ AstromOffsetTable *get_astrom_table (); void put_astrom_table (AstromOffsetTable *myTable); +void free_astrom_table (); int fit_map (AstromOffsetMap *map, StarData *raw, StarData *ref, int Npts); Index: trunk/Ohana/src/relastro/src/BootstrapOps.c =================================================================== --- trunk/Ohana/src/relastro/src/BootstrapOps.c (revision 39402) +++ trunk/Ohana/src/relastro/src/BootstrapOps.c (revision 39457) @@ -98,4 +98,5 @@ myAbort ("invalid option"); } + free (values); return TRUE; Index: trunk/Ohana/src/relastro/src/BrightCatalog.c =================================================================== --- trunk/Ohana/src/relastro/src/BrightCatalog.c (revision 39402) +++ trunk/Ohana/src/relastro/src/BrightCatalog.c (revision 39457) @@ -358,11 +358,6 @@ ohana_memcheck (1); - my_memdump("start of define_bintable"); gfits_define_bintable_column (&theader, "D", "RA", "ra", "degrees", 1.0, 0.0); - my_memdump("after define_bintable RA"); - gfits_define_bintable_column (&theader, "D", "DEC", "dec", "degrees", 1.0, 0.0); - my_memdump("after define_bintable DEC"); - gfits_define_bintable_column (&theader, "E", "MAG_SYS", "magnitude (sys)", NULL, 1.0, 0.0); gfits_define_bintable_column (&theader, "E", "MAG_CAL", "magnitude (cal)", NULL, 1.0, 0.0); @@ -433,4 +428,8 @@ // add the columns to the output array gfits_set_bintable_column (&theader, &ftable, "RA", R, catalog->Nmeasure); + + fprintf (stderr, "--------------- after set_bintable RA --------------"); + ohana_memdump_file (stderr, TRUE); + gfits_set_bintable_column (&theader, &ftable, "DEC", D, catalog->Nmeasure); gfits_set_bintable_column (&theader, &ftable, "MAG_SYS", M, catalog->Nmeasure); @@ -798,5 +797,7 @@ catalogs->maxID = 0; ALLOCATE (catalogs->index, int, catalogs->maxID + 1); - for (i = 0; i <= catalogs->maxID; i++) catalogs->index[i] = -1; + + unsigned int ID; + for (ID = 0; ID <= catalogs->maxID; ID++) catalogs->index[ID] = -1; catalogs->Nsecfilt = Nsecfilt; @@ -807,7 +808,7 @@ ALLOCATE (catalogs->catalog, Catalog, catalogs->NCATALOG); - ALLOCATE (catalogs->catIDs, int, catalogs->NCATALOG); - ALLOCATE (catalogs->NAVERAGE, off_t, catalogs->NCATALOG); - ALLOCATE (catalogs->NMEASURE, off_t, catalogs->NCATALOG); + ALLOCATE (catalogs->catIDs, unsigned int, catalogs->NCATALOG); + ALLOCATE (catalogs->NAVERAGE, off_t, catalogs->NCATALOG); + ALLOCATE (catalogs->NMEASURE, off_t, catalogs->NCATALOG); for (i = 0; i < catalogs->NCATALOG; i++) { @@ -828,9 +829,20 @@ int BrightCatalogSplitFree (CatalogSplitter *catalogs) { - free (catalogs->catIDs); - free (catalogs->NAVERAGE); - free (catalogs->NMEASURE); - free (catalogs->index); - free (catalogs); +# if (0) + int i; + + for (i = 0; i < catalogs->NCATALOG; i++) { + FREE (catalogs->catalog[i].average); + FREE (catalogs->catalog[i].secfilt); + FREE (catalogs->catalog[i].measureT); + } + FREE (catalogs->catalog); +# endif + + FREE (catalogs->catIDs); + FREE (catalogs->NAVERAGE); + FREE (catalogs->NMEASURE); + FREE (catalogs->index); + FREE (catalogs); return TRUE; } @@ -844,5 +856,5 @@ // find the max value of catID in this BrightCatalog - int catIDmax = 0; + unsigned int catIDmax = 0; for (i = 0; i < bcatalog->Naverage; i++) { catIDmax = MAX(catIDmax, bcatalog->average[i].catID); @@ -850,14 +862,16 @@ // XXX validate the measure ID range here - int maxIDold = catalogs->maxID; + unsigned int maxIDold = catalogs->maxID; catalogs->maxID = MAX (maxIDold, catIDmax); // extend the index array and init REALLOCATE (catalogs->index, int, catalogs->maxID + 1); - for (i = maxIDold + 1; i <= catalogs->maxID; i++) catalogs->index[i] = -1; + + unsigned int id; + for (id = maxIDold + 1; id <= catalogs->maxID; id++) catalogs->index[id] = -1; // identify the new catID values for (i = 0; i < bcatalog->Naverage; i++) { - int catID = bcatalog->average[i].catID; + unsigned int catID = bcatalog->average[i].catID; assert (catID > 0); assert (catID < catalogs->maxID + 1); @@ -881,4 +895,5 @@ if (catalogs->Ncatalog >= catalogs->NCATALOG) { + int oldNCATALOG = catalogs->NCATALOG; catalogs->NCATALOG += 16; // fprintf (stderr, "realloc catalogs->catalog: old: %llx ", (long long) catalogs->catalog); @@ -894,8 +909,8 @@ // fprintf (stderr, "new: %llx - %llx\n", (long long) catalogs->NMEASURE, (long long) (catalogs->NMEASURE + sizeof(off_t)*catalogs->NCATALOG)); - REALLOCATE (catalogs->catIDs, int, catalogs->NCATALOG); + REALLOCATE (catalogs->catIDs, unsigned int, catalogs->NCATALOG); int j; - for (j = catalogs->NCATALOG - 16; j < catalogs->NCATALOG; j++) { + for (j = oldNCATALOG; j < catalogs->NCATALOG; j++) { dvo_catalog_init (&catalogs->catalog[j], TRUE); catalogs->catIDs[j] = 0; @@ -917,5 +932,5 @@ // assign the averages to the corresponding catalog for (i = 0; i < bcatalog->Naverage; i++) { - int ID = bcatalog->average[i].catID; + unsigned int ID = bcatalog->average[i].catID; int Nc = catalogs->index[ID]; assert (Nc > -1); @@ -942,5 +957,5 @@ // XXX what about averef and related links? Do I need them? for (i = 0; i < bcatalog->Nmeasure; i++) { - int ID = bcatalog->measure[i].catID; + unsigned int ID = bcatalog->measure[i].catID; int Nc = catalogs->index[ID]; assert (Nc > -1); Index: trunk/Ohana/src/relastro/src/FitAstromOps.c =================================================================== --- trunk/Ohana/src/relastro/src/FitAstromOps.c (revision 39402) +++ trunk/Ohana/src/relastro/src/FitAstromOps.c (revision 39457) @@ -297,4 +297,5 @@ median = x[(int)(0.5*Npoints)]; } + free (x); return median; Index: trunk/Ohana/src/relastro/src/FitChip.c =================================================================== --- trunk/Ohana/src/relastro/src/FitChip.c (revision 39402) +++ trunk/Ohana/src/relastro/src/FitChip.c (revision 39457) @@ -221,10 +221,11 @@ int i, Ns; - float dL, dM, dLsum, dLsum2, dMsum, dMsum2, *dR; - - ALLOCATE (dR, float, Nstars); + float dL, dM; + + ALLOCATE_PTR (dRvec, double, Nstars); + ALLOCATE_PTR (dLvec, double, Nstars); + ALLOCATE_PTR (dMvec, double, Nstars); Ns = 0; - dLsum = dLsum2 = dMsum = dMsum2 = 0.0; for (i = 0; i < Nstars; i++) { if (raw[i].mask) continue; @@ -235,28 +236,37 @@ dM = raw[i].M - ref[i].M; - dLsum += dL; - dLsum2 += dL*dL; - dMsum += dM; - dMsum2 += dM*dM; - - dR[Ns] = hypot (dL, dM); + dLvec[Ns] = dL; + dMvec[Ns] = dM; + dRvec[Ns] = hypot (dL, dM); Ns++; } - *dLsig = sqrt((dLsum2 / Ns) - SQ(dLsum/Ns)); - *dMsig = sqrt((dMsum2 / Ns) - SQ(dMsum/Ns)); + dsort (dLvec, Ns); + dsort (dMvec, Ns); + + double Slo, Shi; + Slo = VectorFractionInterpolate (dLvec, 0.158655, Ns); + Shi = VectorFractionInterpolate (dLvec, 0.841345, Ns); + *dLsig = (Shi - Slo) / 2.0; + free (dLvec); + + Slo = VectorFractionInterpolate (dMvec, 0.158655, Ns); + Shi = VectorFractionInterpolate (dMvec, 0.841345, Ns); + *dMsig = (Shi - Slo) / 2.0; + free (dMvec); + *nKeep = Ns; if (Ns < 5) { *dRsig = NAN; - free (dR); + free (dRvec); return (FALSE); } // for a 2D Gaussian, 40% of the points are within R = 1 sigma - fsort (dR, Ns); - *dRsig = dR[(int)(0.40*Ns)]; + dsort (dRvec, Ns); + *dRsig = dRvec[(int)(0.40*Ns)]; - free (dR); + free (dRvec); return (TRUE); } Index: trunk/Ohana/src/relastro/src/GetAstromError.c =================================================================== --- trunk/Ohana/src/relastro/src/GetAstromError.c (revision 39402) +++ trunk/Ohana/src/relastro/src/GetAstromError.c (revision 39457) @@ -76,11 +76,11 @@ // weight to ensure the image is tied down to the 2mass frame - if (has2MASS && LoopWeight2MASS) { + if (has2MASS && LoopWeight2MASS && (Nloop >= 0)) { dPtotal = dPtotal / LoopWeight2MASS[Nloop]; } - if (is2MASS && LoopWeight2MASS) { + if (is2MASS && LoopWeight2MASS && (Nloop >= 0)) { dPtotal = dPtotal / LoopWeight2MASS[Nloop]; } - if (isTycho && LoopWeightTycho) { + if (isTycho && LoopWeightTycho && (Nloop >= 0)) { dPtotal = dPtotal / LoopWeightTycho[Nloop]; } Index: trunk/Ohana/src/relastro/src/ImageOps.c =================================================================== --- trunk/Ohana/src/relastro/src/ImageOps.c (revision 39402) +++ trunk/Ohana/src/relastro/src/ImageOps.c (revision 39457) @@ -1,4 +1,3 @@ # include "relastro.h" -int isGPC1chip (int photcode); # define USE_IMAGE_ID 1 @@ -10,16 +9,16 @@ // if we read only a subset of the rows from the Image FITS, LineNumber tells us to which row // each image belongs -static off_t *LineNumber; // match of subset to full image table - -static off_t *N_onImage; // number of measurements on image -static off_t *N_ONIMAGE; // allocated number of measurements on image - -static int *Ncatlist; // catalogs associated with each image -static int *NCATLIST; // catalogs associated with each image -static int **catlist; // catalogs associated with each image - -static IDX_T **MeasureToImage; // link from catalog,measure to image -static IDX_T **ImageToCatalog; // catalog which supplied measurement on image -static IDX_T **ImageToMeasure; // measure reference for measurement on image +static off_t *LineNumber = NULL; // match of subset to full image table + +static off_t *N_onImage = NULL; // number of measurements on image +static off_t *N_ONIMAGE = NULL; // allocated number of measurements on image + +static int *Ncatlist = NULL; // catalogs associated with each image +static int *NCATLIST = NULL; // catalogs associated with each image +static int **catlist = NULL; // catalogs associated with each image + +static IDX_T **MeasureToImage = NULL; // link from catalog,measure to image +static IDX_T **ImageToCatalog = NULL; // catalog which supplied measurement on image +static IDX_T **ImageToMeasure = NULL; // measure reference for measurement on image // if we search by image ID, we sort (imageIDs, imageIdx) by imageIDs to get a sorted @@ -27,9 +26,9 @@ # if USE_IMAGE_ID -static off_t *imageIDs; // list of all image IDs -static off_t *imageIdx; // list of index for image IDs +static off_t *imageIDs = NULL; // list of all image IDs +static off_t *imageIdx = NULL; // list of index for image IDs // as an alternative, we generate imageSeq, which directly maps imageID -> seq -static off_t *imageSeq; // list of index for image IDs +static off_t *imageSeq = NULL; // list of index for image IDs // MAX_ID requires 512M to store the image index @@ -120,4 +119,20 @@ } # endif +} + +void freeImages () { + + FREE (LineNumber); + +# if USE_IMAGE_ID + FREE (imageIDs); + FREE (imageIdx); + FREE (imageSeq); +# else + FREE (start); + FREE (stop); +# endif + + free_astrom_table(); } @@ -197,15 +212,23 @@ for (i = 0; i < Ncatalog; i++) { - free (MeasureToImage[i]); - } - free (MeasureToImage); + FREE (MeasureToImage[i]); + } + FREE (MeasureToImage); for (i = 0; i < Nimage; i++) { - if (ImageToCatalog[i]) { free (ImageToCatalog[i]); } - if (ImageToMeasure[i]) { free (ImageToMeasure[i]); } - } - free (ImageToCatalog); - free (ImageToMeasure); - free (N_onImage); - free (N_ONIMAGE); + FREE (ImageToCatalog[i]); + FREE (ImageToMeasure[i]); + if (catlist) { FREE (catlist[i]); } + } + FREE (Ncatlist); + FREE (NCATLIST); + FREE (catlist); + Ncatlist = NULL; + NCATLIST = NULL; + catlist = NULL; + + FREE (ImageToCatalog); + FREE (ImageToMeasure); + FREE (N_onImage); + FREE (N_ONIMAGE); } @@ -240,6 +263,9 @@ name = GetPhotcodeNamebyCode (image[i].photcode); int showExample = (Nfew < 30); - if (showExample) fprintf (stderr, "image "OFF_T_FMT" (%d, %s) has "OFF_T_FMT" of %d measures (%s, %s) ", i, image[i].imageID, image[i].name, N_onImage[i], image[i].nstar, - ohana_sec_to_date(image[i].tzero), name); + if (showExample) { + char *myDate = ohana_sec_to_date(image[i].tzero); + fprintf (stderr, "image "OFF_T_FMT" (%d, %s) has "OFF_T_FMT" of %d measures (%s, %s) ", i, image[i].imageID, image[i].name, N_onImage[i], image[i].nstar, myDate, name); + free (myDate); + } if (N_onImage[i] < 20) { if (showExample) fprintf (stderr, "*"); @@ -508,4 +534,5 @@ m = ImageToMeasure[im][i]; c = ImageToCatalog[im][i]; + myAssert (c < Ncatalog, "oops"); Measure *measure = &catalog[c].measure[m]; @@ -653,4 +680,5 @@ m = ImageToMeasure[im][i]; c = ImageToCatalog[im][i]; + myAssert (c < Ncatalog, "oops"); MeasureTiny *measure = &catalog[c].measureT[m]; @@ -753,4 +781,5 @@ m = ImageToMeasure[im][i]; c = ImageToCatalog[im][i]; + myAssert (c < Ncatalog, "oops"); MeasureTiny *measure = &catalog[c].measureT[m]; @@ -829,4 +858,5 @@ m = ImageToMeasure[im][i]; c = ImageToCatalog[im][i]; + myAssert (c < Ncatalog, "oops"); // XXX should I use the raw coords or just measure.X,Y -> R,D? @@ -871,4 +901,5 @@ m = ImageToMeasure[im][i]; c = ImageToCatalog[im][i]; + myAssert (c < Ncatalog, "oops"); MeasureTiny *measure = &catalog[c].measureT[m]; @@ -900,4 +931,5 @@ } +# if (0) /** lifted from relphot/StarOps.clean_measures */ void FlagOutliers2D(Catalog *catalog); @@ -1267,2 +1299,3 @@ return TRUE; } +# endif Index: trunk/Ohana/src/relastro/src/ImageTable.c =================================================================== --- trunk/Ohana/src/relastro/src/ImageTable.c (revision 39402) +++ trunk/Ohana/src/relastro/src/ImageTable.c (revision 39457) @@ -33,4 +33,10 @@ } + // gfits_db_free (&db); + gfits_free_header (&db.header); + gfits_free_matrix (&db.matrix); + gfits_free_header (&db.theader); + free (db.filename); + *nimage = Nimage; return image; @@ -61,4 +67,6 @@ dvo_image_unlock (&db); + gfits_db_free (&db); + return TRUE; } Index: trunk/Ohana/src/relastro/src/MeasFilterTest.c =================================================================== --- trunk/Ohana/src/relastro/src/MeasFilterTest.c (revision 39457) +++ trunk/Ohana/src/relastro/src/MeasFilterTest.c (revision 39457) @@ -0,0 +1,369 @@ +# include "relastro.h" + +/** lifted from relphot/StarOps.clean_measures */ +void FlagOutliers2D(Catalog *catalog); + +// operates on Full values (not tiny) +void FlagOutliers (Catalog *catalog) { + + // XXX FlagOutliers is now just using FlagOutliers2D + FlagOutliers2D(catalog); + return; + + int Ndel, Nave; + off_t i, j, k, m, N, Nmax, TOOFEW, Nsecfilt; + double Ns, theta, x, y; + double *R, *D, *dR, *dD; + StatType statsR, statsD; + + Nsecfilt = GetPhotcodeNsecfilt(); + assert(catalog[0].Nsecfilt == Nsecfilt); + + if (VERBOSE2) fprintf (stderr, "marking poor measures\n"); + Nmax = 0; + for (i = 0; i < catalog[0].Naverage; i++) { + Nmax = MAX (Nmax, catalog[0].average[i].Nmeasure); + } + + ALLOCATE (R, double, Nmax); + ALLOCATE (D, double, Nmax); + ALLOCATE (dR, double, Nmax); + ALLOCATE (dD, double, Nmax); + + /* it makes no sense to mark 3-sigma outliers with <5 measurements */ + TOOFEW = MAX (5, SRC_MEAS_TOOFEW); + + Ns = CLIP_THRESH; + Ndel = Nave = 0; + + /* loop over each object in the catalog */ + for (j = 0; j < catalog[0].Naverage; j++) { + + // pointer to this set of measurements + m = catalog[0].average[j].measureOffset; + Measure *measure = &catalog[0].measure[m]; + + /* accumulate list of valid measurements */ + N = 0; + for (k = 0; k < catalog[0].average[j].Nmeasure; k++) { + // skip measurements based on user selected criteria + if (!MeasFilterTest(&measure[k], FALSE)) continue; + R[N] = measure[k].R; + D[N] = measure[k].D; + dR[N] = GetAstromError (&measure[k], ERROR_MODE_RA); + dD[N] = GetAstromError (&measure[k], ERROR_MODE_DEC); + if (isnan(R[N]) || isnan(D[N])) continue; + N++; + } + if (N <= TOOFEW) continue; + + /* 3-sigma clip based on stats of inner 50% */ + initstats ("MEAN"); + liststats (R, dR, N, &statsR); + liststats (D, dD, N, &statsD); + + statsR.sigma = MAX (MIN_ERROR, statsR.sigma); + statsD.sigma = MAX (MIN_ERROR, statsD.sigma); + + /* compare per-object distance to this standard deviation, and flag outliers*/ + N = 0; + for (k = 0; k < catalog[0].average[j].Nmeasure; k++) { + // reset flag on each invocation + measure[k].dbFlags &= ~ID_MEAS_POOR_ASTROM; + + // skip measurements based on user selected criteria + if (!MeasFilterTest(&measure[k], FALSE)) continue; + + x = measure[k].R - statsR.median; + y = measure[k].D - statsD.median; + theta = atan2(y,x); + if ((x*x + y*y) > (SQR(statsR.sigma * Ns * cos(theta)) + + SQR(statsD.sigma * Ns * sin(theta)))) { + measure[k].dbFlags |= ID_MEAS_POOR_ASTROM; + Ndel++; + } + N++; + Nave ++; + } + + // examine results + // relastroVisualPlotOutliers(catalog, catalog[0].average[j].measureOffset, catalog[0].average[j].Nmeasure, statsR, statsD, Ns); + } + + if (VERBOSE) fprintf (stderr, "%d measures marked poor, %d total\n", Ndel, Nave); + free (R); + free (dR); + free (D); + free (dD); +} + + +/** an alternative outlier rejection scheme */ +void FlagOutliers2D (Catalog *catalog) { + + int Ndel, Nave; + off_t i, j, k, m, N, Nmax, TOOFEW, Nsecfilt; + double *index; + double Ns, theta, x, y; + double *R, *D, *dR, *dD, *d2; + StatType statsR, statsD; + + // XXX we are not going to use this for now + return; + + Nsecfilt = GetPhotcodeNsecfilt(); + assert(catalog[0].Nsecfilt == Nsecfilt); + + if (VERBOSE2) fprintf (stderr, "marking poor measures\n"); + Nmax = 0; + for (i = 0; i < catalog[0].Naverage; i++) { + Nmax = MAX (Nmax, catalog[0].average[i].Nmeasure); + } + + ALLOCATE (R, double, Nmax); + ALLOCATE (D, double, Nmax); + ALLOCATE (dR, double, Nmax); + ALLOCATE (dD, double, Nmax); + ALLOCATE (d2, double, Nmax); + ALLOCATE (index, double, Nmax); + + /* it makes no sense to mark 3-sigma outliers with <5 measurements */ + TOOFEW = MAX (5, SRC_MEAS_TOOFEW); + + Ns = CLIP_THRESH; + Ndel = Nave = 0; + + /* loop over each object in the catalog */ + for (j = 0; j < catalog[0].Naverage; j++) { + + // pointer to this set of measurements + m = catalog[0].average[j].measureOffset; + Measure *measure = &catalog[0].measure[m]; + + /* accumulate list of valid measurements */ + N = 0; + for (k = 0; k < catalog[0].average[j].Nmeasure; k++) { + + // reset flag on each invocation + measure[k].dbFlags &= ~ID_MEAS_POOR_ASTROM; + + // skip measurements based on user selected criteria + if (!MeasFilterTest(&measure[k], FALSE)) continue; + R[N] = measure[k].R; + D[N] = measure[k].D; + dR[N] = GetAstromError(&measure[k], ERROR_MODE_RA); + dD[N] = GetAstromError(&measure[k], ERROR_MODE_DEC); + if (isnan(R[N]) || isnan(D[N])) continue; + N++; + } + if (N <= TOOFEW) continue; + + /* calculate mean of all points*/ + initstats ("MEAN"); + liststats (R, dR, N, &statsR); + liststats (D, dD, N, &statsD); + statsR.sigma = MAX (MIN_ERROR, statsR.sigma); + statsD.sigma = MAX (MIN_ERROR, statsD.sigma); + + /* calculate deviations of all points*/ + N = 0; + for (k = 0; k < catalog[0].average[j].Nmeasure; k++) { + // skip bad measurements + if (!MeasFilterTest(&measure[k], FALSE)) continue; + x = measure[k].R - statsR.median; + y = measure[k].D - statsD.median; + theta = atan2(y,x); + d2[N] = (x*x + y*y) / (SQR(statsR.sigma * Ns * cos(theta)) + + SQR(statsD.sigma * Ns * sin(theta))); + index[N] = k; + N++; + } + + // sort d2 + dsortpair(d2, index, N); + N = (N/2 > (N-1)) ? N/2 : N-1; + + // recalculate image center, sigma based on closest 50% of points + for (k = 0; k < N; k++) { + off_t ind = (off_t) index[k]; + R[k] = measure[ind].R; + D[k] = measure[ind].D; + dR[k] = GetAstromError(&measure[ind], ERROR_MODE_RA); + dD[k] = GetAstromError(&measure[ind], ERROR_MODE_DEC); + } + liststats (R, dR, N, &statsR); + liststats (D, dD, N, &statsD); + statsR.sigma = MAX (MIN_ERROR, statsR.sigma); + statsD.sigma = MAX (MIN_ERROR, statsD.sigma); + + // use these new statistics to flag outliers + N = 0; + for (k = 0; k < catalog[0].average[j].Nmeasure; k++) { + //skip bad measurements + if (!MeasFilterTest(&measure[k], FALSE)) continue; + x = measure[k].R - statsR.median; + y = measure[k].D - statsD.median; + theta = atan2(y,x); + d2[N] = (x*x + y*y) / (SQR(statsR.sigma * Ns * cos(theta)) + + SQR(statsD.sigma * Ns * sin(theta))); + if ((d2[N]) > 1) { + measure[k].dbFlags |= ID_MEAS_POOR_ASTROM; + Ndel ++; + } + N++; + Nave++; + } // done rejecting outliers + + // examine results + // relastroVisualPlotOutliers(catalog, catalog[0].average[j].measureOffset, catalog[0].average[j].Nmeasure, statsR, statsD, Ns); + + } // done looping over objects + + if (VERBOSE) fprintf (stderr, "%d measures marked poor, %d total\n", Ndel, Nave); + free (R); + free (dR); + free (D); + free (dD); + free (d2); + free (index); +} + +/** Determine whether a measurement should be included in the analysis, based on supplied filter criteria */ +// we only optionally apply the sigma limit: for object averages, this should not be used +int MeasFilterTestTiny(MeasureTiny *measure, int applySigmaLim) { + int found, k; + long mask; + PhotCode *code; + float mag; + + if (!finite(measure[0].R) || !finite(measure[0].D)) return FALSE; + if (!finite(measure[0].M)) return FALSE; //XXX is this necessary for all relastro tasks? + if (!finite(measure[0].dM)) return FALSE; //XXX is this necessary for all relastro tasks? + + if ((MinBadQF > 0.0) && (isGPC1chip(measure[0].photcode) || isGPC1stack(measure[0].photcode))) { + if (!isfinite(measure[0].psfQF)) { return FALSE; }; + if (measure[0].psfQF < MinBadQF) { return FALSE; }; + } + + /* select measurements by photcode, or equiv photcode, if specified */ + if (NphotcodesKeep > 0) { + found = FALSE; + for (k = 0; (k < NphotcodesKeep) && !found; k++) { + if (photcodesKeep[k][0].code == measure[0].photcode) found = TRUE; + if (photcodesKeep[k][0].code == GetPhotcodeEquivCodebyCode(measure[0].photcode)) found = TRUE; + } + if (!found) return FALSE; + } + + if (NphotcodesSkip > 0) { + found = FALSE; + for (k = 0; (k < NphotcodesSkip) && !found; k++) { + if (photcodesSkip[k][0].code == measure[0].photcode) found = TRUE; + if (photcodesSkip[k][0].code == GetPhotcodeEquivCodebyCode(measure[0].photcode)) found = TRUE; + } + if (found) return FALSE; + } + + /* select measurements by time */ + if (TimeSelect) { + if (measure[0].t < TSTART) return FALSE; + if (measure[0].t > TSTOP) return FALSE; + } + + /* select measurements by quality */ + if (PhotFlagSelect) { + if (PhotFlagBad) { + mask = PhotFlagBad; + } else { + code = GetPhotcodebyCode (measure[0].photcode); + if (!code) return FALSE; + mask = code[0].astromBadMask; + } + if (mask & measure[0].photFlags) return FALSE; + } + + /* select measurements by measurement error */ + // this is a bit convoluted: applySigmaLim is only TRUE when this function is + // called by bcatalog. for UpdateObjects, it is FALSE + if (applySigmaLim && (SIGMA_LIM > 0) && (measure[0].dM > SIGMA_LIM)) { + return FALSE; + } + + /* select measurements by mag limit */ + if (ImagSelect) { + mag = PhotInstTiny (measure, MAG_CLASS_PSF); + if (mag < ImagMin || mag > ImagMax) return FALSE; + } + + return TRUE; +} + +# define SUPER_VERBOSE 0 + +/** Determine whether a measurement should be included in the analysis, based on supplied filter criteria */ +// we only optionally apply the sigma limit: for object averages, this should not be used (should it?) +int MeasFilterTest(Measure *measure, int applySigmaLim) { + int found, k; + long mask; + PhotCode *code; + float mag; + + if (!finite(measure[0].R) || !finite(measure[0].D)) { if (SUPER_VERBOSE) fprintf (stderr, "filter 1\n"); return FALSE; }; + if (!finite(measure[0].M)) { if (SUPER_VERBOSE) fprintf (stderr, "filter 2\n"); return FALSE; }; //XXX is this necessary for all relastro tasks? + if (!finite(measure[0].dM)) { if (SUPER_VERBOSE) fprintf (stderr, "filter 3\n"); return FALSE; }; //XXX is this necessary for all relastro tasks? + + /* select measurements by photcode, or equiv photcode, if specified */ + if (NphotcodesKeep > 0) { + found = FALSE; + for (k = 0; (k < NphotcodesKeep) && !found; k++) { + if (photcodesKeep[k][0].code == measure[0].photcode) found = TRUE; + if (photcodesKeep[k][0].code == GetPhotcodeEquivCodebyCode(measure[0].photcode)) found = TRUE; + } + if (!found) { if (SUPER_VERBOSE) fprintf (stderr, "filter 4\n"); return FALSE; }; + } + + if (NphotcodesSkip > 0) { + found = FALSE; + for (k = 0; (k < NphotcodesSkip) && !found; k++) { + if (photcodesSkip[k][0].code == measure[0].photcode) found = TRUE; + if (photcodesSkip[k][0].code == GetPhotcodeEquivCodebyCode(measure[0].photcode)) found = TRUE; + } + if (found) { if (SUPER_VERBOSE) fprintf (stderr, "filter 5\n"); return FALSE; }; + } + + if ((MinBadQF > 0.0) && isGPC1chip(measure[0].photcode)) { + if (!isfinite(measure[0].psfQF)) { if (SUPER_VERBOSE) fprintf (stderr, "filter 6\n"); return FALSE; }; + if (measure[0].psfQF < MinBadQF) { if (SUPER_VERBOSE) fprintf (stderr, "filter 7\n"); return FALSE; }; + } + + /* select measurements by time */ + if (TimeSelect) { + if (measure[0].t < TSTART) { if (SUPER_VERBOSE) fprintf (stderr, "filter 8\n"); return FALSE; }; + if (measure[0].t > TSTOP) { if (SUPER_VERBOSE) fprintf (stderr, "filter 9\n"); return FALSE; }; + } + + /* select measurements by quality */ + if (PhotFlagSelect) { + if (PhotFlagBad) { + mask = PhotFlagBad; + } else { + code = GetPhotcodebyCode (measure[0].photcode); + if (!code) { if (SUPER_VERBOSE) fprintf (stderr, "filter 10\n"); return FALSE; }; + mask = code[0].astromBadMask; + } + if (mask & measure[0].photFlags) { if (SUPER_VERBOSE) fprintf (stderr, "filter 11\n"); return FALSE; }; + } + + /* select measurements by measurement error */ + if (applySigmaLim && (SIGMA_LIM > 0) && (measure[0].dM > SIGMA_LIM)) { + { if (SUPER_VERBOSE) fprintf (stderr, "filter 12\n"); return FALSE; }; + } + + /* select measurements by mag limit */ + if (ImagSelect) { + mag = PhotInst (measure, MAG_CLASS_PSF); + if (mag < ImagMin || mag > ImagMax) { if (SUPER_VERBOSE) fprintf (stderr, "filter 13\n"); return FALSE; }; + } + + return TRUE; +} Index: trunk/Ohana/src/relastro/src/MosaicOps.c =================================================================== --- trunk/Ohana/src/relastro/src/MosaicOps.c (revision 39402) +++ trunk/Ohana/src/relastro/src/MosaicOps.c (revision 39457) @@ -3,5 +3,5 @@ // array of mosaic definition structures static off_t Nmosaic; -static Mosaic *mosaic; +static Mosaic *mosaic = NULL; // list of all images associated with a mosaic @@ -185,4 +185,21 @@ } +void freeMosaics () { + + off_t i; + + if (!mosaic) return; + + for (i = 0; i < Nmosaic; i++) { + FREE (mosaic_own_images[i]); + } + + FREE (mosaic); + FREE (Nmosaic_own_images); + FREE (Amosaic_own_images); + FREE (mosaic_own_images); + FREE (mosaic_for_images); +} + // return StarData values for detections in the specified image, converting coordinates from the // chip positions: X,Y -> L,M -> P,Q -> R,D @@ -274,4 +291,5 @@ // extend each host image table to include the mosaic 'images' needed by the host int select_mosaics_hostregion (RegionHostTable *regionHosts, Image *image, off_t Nimage) { + OHANA_UNUSED_PARAM(Nimage); int i; @@ -324,6 +342,10 @@ } + free (addMosaic); + host->Nimage += Nadd; } + free (mosaicUsed); + return TRUE; } Index: trunk/Ohana/src/relastro/src/SetSignals.c =================================================================== --- trunk/Ohana/src/relastro/src/SetSignals.c (revision 39402) +++ trunk/Ohana/src/relastro/src/SetSignals.c (revision 39457) @@ -54,4 +54,6 @@ void my_memdump (char *message) { + + return; fprintf (stderr, "%s\n", message); Index: trunk/Ohana/src/relastro/src/StarMaps.c =================================================================== --- trunk/Ohana/src/relastro/src/StarMaps.c (revision 39402) +++ trunk/Ohana/src/relastro/src/StarMaps.c (revision 39457) @@ -25,5 +25,5 @@ } StarMap; -static StarMap *starmap; +static StarMap *starmap = NULL; // static int Nstarmap; -- is it always == Nimages? @@ -55,4 +55,19 @@ } +void freeStarMaps () { + + off_t i, Nimages; + + getimages(&Nimages, NULL); + + for (i = 0; i < Nimages; i++) { + FREE (starmap[i].points); + FREE (starmap[i].stars); + } + + FREE (starmap); + return; +} + int updateStarMaps(Catalog *catalog) { Index: trunk/Ohana/src/relastro/src/UpdateChips.c =================================================================== --- trunk/Ohana/src/relastro/src/UpdateChips.c (revision 39402) +++ trunk/Ohana/src/relastro/src/UpdateChips.c (revision 39457) @@ -152,4 +152,6 @@ } } + + AstromOffsetMapFree (oldCoords.offsetMap); // Apply the modified coords back to the measure.R,D. Note that raw.R,D, ref.L,M, etc @@ -356,4 +358,6 @@ } + AstromOffsetMapFree (oldCoords.offsetMap); + // apply the modified R,D back to the measures setImageRaw (threadinfo->catalog, threadinfo->Ncatalog, i, raw, Nraw, MODE_MOSAIC); @@ -569,4 +573,5 @@ AstromOffsetMapSetOrder (truMap, tmpMap->Nx, tmpMap->Ny, image); AstromOffsetMapCopyData (truMap, tmpMap); + AstromOffsetMapFree (tmpMap); return; } @@ -579,4 +584,5 @@ unlockUpdateChips (); AstromOffsetMapCopyData (image->coords.offsetMap, tmpMap); + AstromOffsetMapFree (tmpMap); return; } Index: trunk/Ohana/src/relastro/src/UpdateMeasures.c =================================================================== --- trunk/Ohana/src/relastro/src/UpdateMeasures.c (revision 39402) +++ trunk/Ohana/src/relastro/src/UpdateMeasures.c (revision 39457) @@ -71,5 +71,5 @@ double dR = 3600.0*fabs(csdec*(oldR - R)); - double dD = 3600.0*fabs(oldD - R); + double dD = 3600.0*fabs(oldD - D); // complain if the new location is far from the old location Index: trunk/Ohana/src/relastro/src/UpdateObjectOffsets.c =================================================================== --- trunk/Ohana/src/relastro/src/UpdateObjectOffsets.c (revision 39402) +++ trunk/Ohana/src/relastro/src/UpdateObjectOffsets.c (revision 39457) @@ -62,5 +62,5 @@ // uses MeasureTiny values, but will also update Measure values when populated // if Galaxy Motion Model correction is active, do not apply at this last stage - UpdateObjects (&catalog, 1, FALSE); + UpdateObjects (&catalog, 1, -1); free_tiny_values(&catalog); @@ -218,5 +218,8 @@ fprintf (stderr, "command: %s\n", command); - if (PARALLEL_MANUAL) continue; + if (PARALLEL_MANUAL) { + free (command); + continue; + } if (PARALLEL_SERIAL) { @@ -236,4 +239,5 @@ group->hosts[i][0].pid = pid; // save for future reference } + free (command); } @@ -338,5 +342,8 @@ fprintf (stderr, "command: %s\n", command); - if (PARALLEL_MANUAL) continue; + if (PARALLEL_MANUAL) { + free (command); + continue; + } if (PARALLEL_SERIAL) { @@ -356,4 +363,5 @@ table->hosts[i].pid = pid; // save for future reference } + free (command); } @@ -366,4 +374,6 @@ } + FreeHostTable (table); + return TRUE; } Index: trunk/Ohana/src/relastro/src/UpdateObjects.c =================================================================== --- trunk/Ohana/src/relastro/src/UpdateObjects.c (revision 39402) +++ trunk/Ohana/src/relastro/src/UpdateObjects.c (revision 39457) @@ -121,4 +121,5 @@ // be called with just MeasureTiny set and Measure == NULL int UpdateObjects_Chips (Average *average, SecFilt *secfilt, MeasureTiny *measure, Measure *measureBig, int Nsecfilt, FitStats *fitStats, int cat, off_t measOff) { + OHANA_UNUSED_PARAM(Nsecfilt); int k; @@ -257,6 +258,7 @@ } + // if we have fitted (and accepted) a parallax model, get the best pm fit and chisq + // given the set of points (mask is respected) if (average[0].flags & ID_OBJ_USE_PAR) { - // get the best pm fit and chisq given the set of points (mask is respected) if (!FitPM_Basic (&fitPM, fitStats->fitdataPM, fitStats->points, fitStats->Npoints)) { average[0].flags |= ID_OBJ_BAD_PM; @@ -551,4 +553,5 @@ // be called with just MeasureTiny set and Measure == NULL int UpdateObjects_Stack (Average *average, SecFilt *secfilt, MeasureTiny *measure, Measure *measureBig, int Nsecfilt, FitStats *fitStats) { + OHANA_UNUSED_PARAM(Nsecfilt); int status; @@ -650,4 +653,5 @@ int UpdateObjects_SelectMeasures (FitStats *fit, Average *average, SecFilt *secfilt, MeasureTiny *measure, Measure *measureBig, int isStack, int *stackEntry) { + OHANA_UNUSED_PARAM(secfilt); // I've already allocated fit->points (and fit->sample) with space for fit->NpointsAlloc entries Index: trunk/Ohana/src/relastro/src/UpdateStacks.c =================================================================== --- trunk/Ohana/src/relastro/src/UpdateStacks.c (revision 39457) +++ trunk/Ohana/src/relastro/src/UpdateStacks.c (revision 39457) @@ -0,0 +1,60 @@ +# include "relastro.h" + +// NOTE: we only measure the systematic floor of the astrometric scatter per stack, no change to the calibration +int UpdateStacks (Catalog *catalog, int Ncatalog) { + + off_t Nimage; + Image *image = getimages (&Nimage, NULL); + + off_t Nskip = 0; + off_t Nmeas = 0; + + // measure the scatter for each stack + for (off_t i = 0; i < Nimage; i++) { + + /* skip all except stack images */ + if (!isGPC1stack(image[i].photcode)) continue; + + /* convert measure coordinates to raw entries */ + off_t Nraw; + StarData *raw = getImageRaw (catalog, Ncatalog, i, &Nraw, MODE_SIMPLE); + if (!raw) { + Nskip ++; + continue; + } + if (Nraw <= IMFIT_TOO_FEW) { + Nskip ++; + free (raw); + continue; + } + + /* convert average coordinates to ref entries */ + off_t Nref; + StarData *ref = getImageRef (catalog, Ncatalog, i, &Nref, MODE_SIMPLE); + if (!ref) { + Nskip ++; + free (raw); + continue; + } + + // note that Nraw & Nref must be equal: if not, we made a programming error in one of these two functions. + assert (Nraw == Nref); + + int Nstat; + float dLsig, dMsig, dRsig; + GetScatterRawRef(&dLsig, &dMsig, &dRsig, &Nstat, raw, ref, Nraw, IMFIT_SYS_SIGMA_LIM); + + // XXX: I need to convert dLsig, dMsig from degrees to pixels + dLsig *= 3600.0 / 0.25; + dMsig *= 3600.0 / 0.25; + + image[i].dXpixSys = dLsig; + image[i].dYpixSys = dMsig; + image[i].nFitAstrom = Nstat; + } + + fprintf (stderr, "UpdateStacks: %d measured, %d skipped\n", (int) Nmeas, (int) Nskip); + return (TRUE); +} + +// XXX 2 hardwired hacks in this file: 1) photcode hardwired for GPC1 stacks, 2) platescale Index: trunk/Ohana/src/relastro/src/args.c =================================================================== --- trunk/Ohana/src/relastro/src/args.c (revision 39402) +++ trunk/Ohana/src/relastro/src/args.c (revision 39457) @@ -51,4 +51,12 @@ remove_argument (N, &argc, argv); RELASTRO_OP = OP_UPDATE_OFFSETS; + } + + // elements needed for parallel regions / parallel images + MANUAL_UNIQUER = NULL; + if ((N = get_argument (argc, argv, "-manual-uniquer"))) { + remove_argument (N, &argc, argv); + MANUAL_UNIQUER = strcreate (argv[N]); + remove_argument (N, &argc, argv); } @@ -372,5 +380,4 @@ snprintf (PHOTCODE_SKIP_LIST, Ntotal, "%s,%s", SuperCOSMOS_SKIP, RawSkip); free (RawSkip); - free (SuperCOSMOS_SKIP); } else { PHOTCODE_SKIP_LIST = RawSkip; @@ -378,4 +385,5 @@ remove_argument (N, &argc, argv); } + free (SuperCOSMOS_SKIP); VERBOSE = VERBOSE2 = FALSE; @@ -562,4 +570,39 @@ if (argc != 1) usage (); return TRUE; +} + +void relastro_free (SkyTable *sky, SkyList *skylist) { + FREE (REGION_FILE); + FREE (IMAGE_TABLE); + FREE (LoopWeight2MASSstr); + FREE (LoopWeightTychostr); + + FREE (PHOTCODE_SKIP_LIST); + FREE (PHOTCODE_KEEP_LIST); + FREE (PHOTCODE_RESET_LIST); + FREE (DCR_RED_COLOR_POS); + FREE (DCR_RED_COLOR_NEG); + FREE (DCR_BLUE_COLOR_POS); + FREE (DCR_BLUE_COLOR_NEG); + + FREE (PHOTCODE_A_LIST); + FREE (PHOTCODE_B_LIST); + FREE (HIGH_SPEED_DIR); + FREE (BCATALOG); + FREE (HOSTDIR); + + // these are set in initialize + FREE(photcodesKeep); + FREE(photcodesSkip); + FREE(photcodesReset); + FREE(photcodesGroupA); + FREE(photcodesGroupB); + + SkyTableFree (sky); + SkyListFree(skylist); + FreePhotcodeTable(); + + ohana_memcheck (VERBOSE); + ohana_memdump (VERBOSE); } @@ -922,5 +965,5 @@ } -void relastro_client_free () { +void relastro_client_free (SkyTable *sky, SkyList *skylist) { FREE (PHOTCODE_SKIP_LIST); FREE (PHOTCODE_KEEP_LIST); @@ -936,4 +979,18 @@ FREE(BCATALOG); FREE (HOSTDIR); + + // these are set in initialize + FREE(photcodesKeep); + FREE(photcodesSkip); + FREE(photcodesReset); + FREE(photcodesGroupA); + FREE(photcodesGroupB); + + SkyTableFree (sky); + SkyListFree(skylist); + FreePhotcodeTable(); + + ohana_memcheck (VERBOSE); + ohana_memdump (VERBOSE); } Index: trunk/Ohana/src/relastro/src/assign_images.c =================================================================== --- trunk/Ohana/src/relastro/src/assign_images.c (revision 39402) +++ trunk/Ohana/src/relastro/src/assign_images.c (revision 39457) @@ -156,4 +156,5 @@ // this is a bit memory expensive : I am making a complete copy of the image table here + // XXX is adding an image, can we just add a pointer? off_t Nsubset = host->Nimage; host->image[Nsubset] = image[j]; @@ -168,4 +169,6 @@ // save the astrometry maps, where they exist + // here we are adding an AstromOffsetTable, but the data are only pointers + // we should NOT free this table with AstromOffsetTableFree() if (image[j].coords.offsetMap) { if (!host->astromTable) { Index: trunk/Ohana/src/relastro/src/bcatalog.c =================================================================== --- trunk/Ohana/src/relastro/src/bcatalog.c (revision 39402) +++ trunk/Ohana/src/relastro/src/bcatalog.c (revision 39457) @@ -184,4 +184,7 @@ offset = catalog[0].average[i].measureOffset + j; + // we do not measure astrometry for the warps, so skip: + if (isGPC1warp(catalog[0].measure[offset].photcode)) continue; + // filter objects based on user supplied criteria, including SIGMA_LIM if (!MeasFilterTest(&catalog[0].measure[offset], TRUE)) { Index: trunk/Ohana/src/relastro/src/dvo_astrom_ops.c =================================================================== --- trunk/Ohana/src/relastro/src/dvo_astrom_ops.c (revision 39402) +++ trunk/Ohana/src/relastro/src/dvo_astrom_ops.c (revision 39457) @@ -2,5 +2,9 @@ /* the Measure carries the instantaneous mean position at the epoch t */ +// average & secfilt no longer used since R and D are now in measure + double getMeanR (MeasureTiny *measure, Average *average, SecFilt *secfilt) { + OHANA_UNUSED_PARAM(average); + OHANA_UNUSED_PARAM(secfilt); double ra; @@ -14,4 +18,6 @@ double getMeanD (MeasureTiny *measure, Average *average, SecFilt *secfilt) { + OHANA_UNUSED_PARAM(average); + OHANA_UNUSED_PARAM(secfilt); double dec; @@ -25,4 +31,6 @@ int setMeanR (double ra_fit, MeasureTiny *measure, Average *average, SecFilt *secfilt) { + OHANA_UNUSED_PARAM(average); + OHANA_UNUSED_PARAM(secfilt); // old: measure[0].dR += (ra_fit - average[0].R) * 3600.0; @@ -34,4 +42,6 @@ int setMeanD (double dec_fit, MeasureTiny *measure, Average *average, SecFilt *secfilt) { + OHANA_UNUSED_PARAM(average); + OHANA_UNUSED_PARAM(secfilt); // old: measure[0].dD += (dec_fit - average[0].D) * 3600.0; @@ -43,4 +53,6 @@ double getMeanR_Big (Measure *measure, Average *average, SecFilt *secfilt) { + OHANA_UNUSED_PARAM(average); + OHANA_UNUSED_PARAM(secfilt); double ra; @@ -54,4 +66,6 @@ double getMeanD_Big (Measure *measure, Average *average, SecFilt *secfilt) { + OHANA_UNUSED_PARAM(average); + OHANA_UNUSED_PARAM(secfilt); double dec; @@ -65,4 +79,6 @@ int setMeanR_Big (double ra_fit, Measure *measure, Average *average, SecFilt *secfilt) { + OHANA_UNUSED_PARAM(average); + OHANA_UNUSED_PARAM(secfilt); // old: measure[0].dR += (ra_fit - average[0].R) * 3600.0; @@ -74,4 +90,6 @@ int setMeanD_Big (double dec_fit, Measure *measure, Average *average, SecFilt *secfilt) { + OHANA_UNUSED_PARAM(average); + OHANA_UNUSED_PARAM(secfilt); // measure[0].dD += (dec_fit - average[0].D) * 3600.0; Index: trunk/Ohana/src/relastro/src/fitobj.c =================================================================== --- trunk/Ohana/src/relastro/src/fitobj.c (revision 39402) +++ trunk/Ohana/src/relastro/src/fitobj.c (revision 39457) @@ -156,4 +156,5 @@ Catalog *mkstar (FitStats *fitStats, double Ro, double Do, double uR, double uD, double plx, int Npoints, int Nbad) { + OHANA_UNUSED_PARAM(fitStats); int i; @@ -239,2 +240,17 @@ return catalog; } + +// these dummy functions are used to avoid including ImageOps, etc +int areImagesMatched () { + return FALSE; +} +float getColorBlue (off_t meas, int cat) { + OHANA_UNUSED_PARAM(meas); + OHANA_UNUSED_PARAM(cat); + return (NAN); +} +float getColorRed (off_t meas, int cat) { + OHANA_UNUSED_PARAM(meas); + OHANA_UNUSED_PARAM(cat); + return (NAN); +} Index: trunk/Ohana/src/relastro/src/fitobj2.c =================================================================== --- trunk/Ohana/src/relastro/src/fitobj2.c (revision 39402) +++ trunk/Ohana/src/relastro/src/fitobj2.c (revision 39457) @@ -173,4 +173,5 @@ Catalog *mkstar (FitStats *fitStats, double Ro, double Do, double uR, double uD, double plx, int Npoints, int Nbad) { + OHANA_UNUSED_PARAM(fitStats); int i; @@ -264,2 +265,17 @@ return catalog; } + +// these dummy functions are used to avoid including ImageOps, etc +int areImagesMatched () { + return FALSE; +} +float getColorBlue (off_t meas, int cat) { + OHANA_UNUSED_PARAM(meas); + OHANA_UNUSED_PARAM(cat); + return (NAN); +} +float getColorRed (off_t meas, int cat) { + OHANA_UNUSED_PARAM(meas); + OHANA_UNUSED_PARAM(cat); + return (NAN); +} Index: trunk/Ohana/src/relastro/src/fitstk.c =================================================================== --- trunk/Ohana/src/relastro/src/fitstk.c (revision 39402) +++ trunk/Ohana/src/relastro/src/fitstk.c (revision 39457) @@ -279,2 +279,17 @@ return catalog; } + +// these dummy functions are used to avoid including ImageOps, etc +int areImagesMatched () { + return FALSE; +} +float getColorBlue (off_t meas, int cat) { + OHANA_UNUSED_PARAM(meas); + OHANA_UNUSED_PARAM(cat); + return (NAN); +} +float getColorRed (off_t meas, int cat) { + OHANA_UNUSED_PARAM(meas); + OHANA_UNUSED_PARAM(cat); + return (NAN); +} Index: trunk/Ohana/src/relastro/src/high_speed_utils.c =================================================================== --- trunk/Ohana/src/relastro/src/high_speed_utils.c (revision 39402) +++ trunk/Ohana/src/relastro/src/high_speed_utils.c (revision 39457) @@ -8,6 +8,4 @@ #define DEFAULT_WHERE_B "((flags & 0x01000000) == 0) && ((flags & 0x04000000) && ((y:nphot > 1) || (z:nphot > 1)) &&((y:err < 0.2) || (z:err < 0.2)))" - -extern char **isolate_elements(int, char **, int *); int print_error(); @@ -76,5 +74,5 @@ // parse it into elements of the where condition char **cstack = NULL; - int Ncstack; + unsigned int Ncstack; cstack = isolate_elements(nWords, words, &Ncstack); free(words); Index: trunk/Ohana/src/relastro/src/indexCatalogs.c =================================================================== --- trunk/Ohana/src/relastro/src/indexCatalogs.c (revision 39402) +++ trunk/Ohana/src/relastro/src/indexCatalogs.c (revision 39457) @@ -1,8 +1,8 @@ # include "relastro.h" -static int catIDmax = 0; -static int *catIDseq = NULL; -static int *objIDmax = NULL; -static int **objIDseq = NULL; +static unsigned int catIDmax = 0; +static int *catIDseq = NULL; +static unsigned int *objIDmax = NULL; +static int **objIDseq = NULL; int indexCatalogs (Catalog *catalog, int Ncatalog) { @@ -19,15 +19,17 @@ ALLOCATE (catIDseq, int, catIDmax + 1); - for (i = 0; i < catIDmax + 1; i++) { - catIDseq[i] = -1; + + unsigned int ID; + for (ID = 0; ID < catIDmax + 1; ID++) { + catIDseq[ID] = -1; } for (i = 0; i < Ncatalog; i++) { - int catID = catalog[i].catID; + unsigned int catID = catalog[i].catID; catIDseq[catID] = i; } - ALLOCATE (objIDmax, int, Ncatalog); - ALLOCATE (objIDseq, int *, Ncatalog); + ALLOCATE (objIDmax, unsigned int, Ncatalog); + ALLOCATE (objIDseq, int *, Ncatalog); for (i = 0; i < Ncatalog; i++) { objIDmax[i] = 0; @@ -37,10 +39,10 @@ ALLOCATE (objIDseq[i], int, objIDmax[i] + 1); - for (j = 0; j < objIDmax[i] + 1; j++) { - objIDseq[i][j] = -1; + for (ID = 0; ID < objIDmax[i] + 1; ID++) { + objIDseq[i][ID] = -1; } for (j = 0; j < catalog[i].Naverage; j++) { - int objID = catalog[i].average[j].objID; + unsigned int objID = catalog[i].average[j].objID; objIDseq[i][objID] = j; } @@ -49,5 +51,16 @@ } -int catID_and_objID_to_seq (int catID, int objID, int *catSeq, off_t *objSeq) { +void freeCatalogIndexes (int Ncatalog) { + int i; + + for (i = 0; i < Ncatalog; i++) { + free (objIDseq[i]); + } + free (objIDmax); + free (objIDseq); + free (catIDseq); +} + +int catID_and_objID_to_seq (unsigned int catID, unsigned int objID, int *catSeq, off_t *objSeq) { if (catID > catIDmax) return FALSE; Index: trunk/Ohana/src/relastro/src/initialize.c =================================================================== --- trunk/Ohana/src/relastro/src/initialize.c (revision 39402) +++ trunk/Ohana/src/relastro/src/initialize.c (revision 39457) @@ -2,4 +2,10 @@ void initialize (int argc, char **argv) { + + photcodesKeep = NULL; + photcodesSkip = NULL; + photcodesReset = NULL; + photcodesGroupA = NULL; + photcodesGroupB = NULL; ConfigInit (&argc, argv); Index: trunk/Ohana/src/relastro/src/launch_region_hosts.c =================================================================== --- trunk/Ohana/src/relastro/src/launch_region_hosts.c (revision 39402) +++ trunk/Ohana/src/relastro/src/launch_region_hosts.c (revision 39457) @@ -161,5 +161,8 @@ fprintf (stderr, "command: %s\n", command); - if (PARALLEL_REGIONS_MANUAL) continue; + if (PARALLEL_REGIONS_MANUAL) { + free (command); + continue; + } // launch the job, then wait for it to be done loading catalogs. force the remote @@ -180,4 +183,6 @@ check_sync_file (syncfile, 1); free (syncfile); + + free (command); } Index: trunk/Ohana/src/relastro/src/load_catalogs.c =================================================================== --- trunk/Ohana/src/relastro/src/load_catalogs.c (revision 39402) +++ trunk/Ohana/src/relastro/src/load_catalogs.c (revision 39457) @@ -142,5 +142,10 @@ int TIME = time(NULL); int PID = getpid(); - snprintf (uniquer, 12, "%05d.%05d", PID, TIME % 100000); + + if (PARALLEL_MANUAL && MANUAL_UNIQUER) { + snprintf (uniquer, 12, "%11s", MANUAL_UNIQUER); + } else { + snprintf (uniquer, 12, "%05d.%05d", PID, TIME % 100000); + } // load the list of hosts @@ -224,5 +229,8 @@ fprintf (stderr, "command: %s\n", command); - if (PARALLEL_MANUAL) continue; + if (PARALLEL_MANUAL) { + free (command); + continue; + } if (PARALLEL_SERIAL) { @@ -242,4 +250,5 @@ table->hosts[i].pid = pid; // save for future reference } + free (command); } @@ -289,6 +298,15 @@ } + FreeHostTable (table); + Catalog *catalog = catalogs->catalog; *Ncatalog = catalogs->Ncatalog; + + // need to free the place-holder catalogs: + for (i = catalogs->Ncatalog; i < catalogs->NCATALOG; i++) { + free (catalogs->catalog[i].average); + free (catalogs->catalog[i].measureT); + free (catalogs->catalog[i].secfilt); + } int Nmeasure = 0; Index: trunk/Ohana/src/relastro/src/load_images.c =================================================================== --- trunk/Ohana/src/relastro/src/load_images.c (revision 39402) +++ trunk/Ohana/src/relastro/src/load_images.c (revision 39457) @@ -27,15 +27,4 @@ MARKTIME("build chip match: %f sec\n", dtime); - char mapfile[DVO_MAX_PATH]; - snprintf (mapfile, DVO_MAX_PATH, "%s/AstroMap.fits", CATDIR); - table = AstromOffsetMapLoad (mapfile, 100000, VERBOSE); - - // assign images.coords.offsetMap -> table->map[i] - if (table) { - AstromOffsetTableMatchChips (image, Nimage, table); - } else { - table = AstromOffsetTableInit (); - } - // select the images which overlap the selected sky regions if (UseAllImages) { @@ -54,4 +43,27 @@ } + /* unlock, if we can (else, unlocked below) and free, if we can */ + int unlockImages = !UPDATE || (RELASTRO_OP == OP_UPDATE_OFFSETS); + if (unlockImages) { + if (subset != image) { + // if we have generated an image subset and we are running UPDATE_OFFSETS, the we can free images here + free (image); + db[0].ftable.buffer = NULL; + BuildChipMatch (subset, Nsubset); + } + dvo_image_unlock (db); + } + + char mapfile[DVO_MAX_PATH]; + snprintf (mapfile, DVO_MAX_PATH, "%s/AstroMap.fits", CATDIR); + table = AstromOffsetMapLoad (mapfile, 100000, VERBOSE); + + // assign images.coords.offsetMap -> table->map[i] + if (table) { + AstromOffsetTableMatchChips (subset, Nsubset, table); + } else { + table = AstromOffsetTableInit (); + } + initImages (subset, LineNumber, Nsubset); MARKTIME(" init images: %f sec\n", dtime); @@ -60,8 +72,4 @@ MARKTIME(" init mosaics: %f sec\n", dtime); - /* unlock, if we can (else, unlocked below) */ - int unlockImages = !UPDATE || (RELASTRO_OP == OP_UPDATE_OFFSETS); - if (unlockImages) dvo_image_unlock (db); - return TRUE; } @@ -102,2 +110,7 @@ } +void free_astrom_table () { + if (!table) return; + AstromOffsetTableFree (table); + free (table); +} Index: trunk/Ohana/src/relastro/src/plot_scatter.c =================================================================== --- trunk/Ohana/src/relastro/src/plot_scatter.c (revision 39402) +++ trunk/Ohana/src/relastro/src/plot_scatter.c (revision 39457) @@ -2,4 +2,6 @@ void plot_scatter (Catalog *catalog, int Ncatalog) { + OHANA_UNUSED_PARAM(catalog); + OHANA_UNUSED_PARAM(Ncatalog); # if (0) Index: trunk/Ohana/src/relastro/src/plotstuff.c =================================================================== --- trunk/Ohana/src/relastro/src/plotstuff.c (revision 39402) +++ trunk/Ohana/src/relastro/src/plotstuff.c (revision 39457) @@ -44,4 +44,5 @@ void JpegPlot (Graphdata *graphmode, int N, char *filename) { + OHANA_UNUSED_PARAM(graphmode); if (Xgraph[N] == 0) return; @@ -52,4 +53,5 @@ void PSPlot (Graphdata *graphmode, int N, char *filename) { + OHANA_UNUSED_PARAM(graphmode); if (Xgraph[N] == 0) return; @@ -79,4 +81,5 @@ void PlotVector (int Npts, double *vect, int mode, int N, char *type) { + OHANA_UNUSED_PARAM(mode); float *values; Index: trunk/Ohana/src/relastro/src/relastro.c =================================================================== --- trunk/Ohana/src/relastro/src/relastro.c (revision 39402) +++ trunk/Ohana/src/relastro/src/relastro.c (revision 39457) @@ -7,7 +7,12 @@ initialize (argc, argv); - SkyTable *sky = SkyTableLoadOptimal (CATDIR, SKY_TABLE, GSCFILE, TRUE, SKY_DEPTH, VERBOSE); - SkyTableSetFilenames (sky, CATDIR, "cpt"); - SkyList *skylist = SkyListByPatch (sky, -1, &UserPatch); + SkyTable *sky = NULL; + SkyList *skylist = NULL; + + if ((RELASTRO_OP != OP_PARALLEL_IMAGES) && (RELASTRO_OP != OP_PARALLEL_REGIONS)) { + sky = SkyTableLoadOptimal (CATDIR, SKY_TABLE, GSCFILE, TRUE, SKY_DEPTH, VERBOSE); + SkyTableSetFilenames (sky, CATDIR, "cpt"); + skylist = SkyListByPatch (sky, -1, &UserPatch); + } switch (RELASTRO_OP) { @@ -15,27 +20,15 @@ /* the object analysis is a separate process iterating over catalogs */ relastro_objects (skylist, 0, NULL); - exit (0); - - case OP_HIGH_SPEED: - /* high-speed is a 2pt cross-correlation process for linking moving objects (high PM) */ - high_speed_catalogs (sky, skylist, 0, NULL); - exit (0); - - case OP_HPM: - hpm_catalogs (sky, skylist, 0, NULL); - exit (0); - - case OP_MERGE_SOURCE: - /* a special method to manually merge unlinked detections of sources togther (not parallel) */ - relastro_merge_source (sky); + relastro_free (sky, skylist); exit (0); case OP_IMAGES: relastro_images (skylist); + relastro_free (sky, skylist); exit (0); - case OP_UPDATE_OFFSETS: + case OP_UPDATE_OFFSETS: { + FITS_DB db; if (!PARALLEL) { - FITS_DB db; set_db (&db); @@ -55,5 +48,14 @@ // iterate over catalogs to make detection coordinates consistant UpdateObjectOffsets (skylist, 0, NULL); + + if (!PARALLEL) { + gfits_db_free (&db); + freeImages(); + freeMosaics (); + } + + relastro_free (sky, skylist); exit (0); + } case OP_PARALLEL_REGIONS: @@ -67,4 +69,18 @@ exit (0); + case OP_HIGH_SPEED: + /* high-speed is a 2pt cross-correlation process for linking moving objects (high PM) */ + high_speed_catalogs (sky, skylist, 0, NULL); + exit (0); + + case OP_HPM: + hpm_catalogs (sky, skylist, 0, NULL); + exit (0); + + case OP_MERGE_SOURCE: + /* a special method to manually merge unlinked detections of sources togther (not parallel) */ + relastro_merge_source (sky); + exit (0); + default: fprintf (stderr, "impossible!\n"); Index: trunk/Ohana/src/relastro/src/relastro_client.c =================================================================== --- trunk/Ohana/src/relastro/src/relastro_client.c (revision 39402) +++ trunk/Ohana/src/relastro/src/relastro_client.c (revision 39457) @@ -55,8 +55,4 @@ } - SkyTableFree (sky); - SkyListFree(skylist); - FreePhotcodeTable(); - int i; for (i = 0; i < Ncatalog; i++) { @@ -66,13 +62,5 @@ BrightCatalogFree(bcatalog); - relastro_client_free (); - FREE(photcodesKeep); - FREE(photcodesSkip); - FREE(photcodesReset); - FREE(photcodesGroupA); - FREE(photcodesGroupB); - - ohana_memcheck (TRUE); - ohana_memdump (TRUE); + relastro_client_free (sky, skylist); break; } @@ -81,16 +69,5 @@ // USAGE: relastro_client -update-objects relastro_objects (skylist, HOST_ID, HOSTDIR); - break; - } - - case OP_HIGH_SPEED: { - // USAGE: relastro_client -high-speed - high_speed_catalogs (sky, skylist, HOST_ID, HOSTDIR); - break; - } - - case OP_HPM: { - // USAGE: relastro_client -high-speed - hpm_catalogs (sky, skylist, HOST_ID, HOSTDIR); + relastro_client_free (sky, skylist); break; } @@ -112,25 +89,24 @@ /* load regions and images based on specified sky patch (default depth) */ load_images (&db, skylist, FALSE, USE_ALL_IMAGES); - - // I can free the database after I have loaded the images... - - // // load the image subset table from the specified location - // off_t Nimage; - // ImageSubset *image = ImageSubsetLoad (IMAGES, &Nimage); - // if (!image) { - // fprintf (stderr, "ERROR loading image subset %s\n", CATDIR); - // exit (2); - // } - // - // // save the available image information in the static array in ImageOps.c - // initImagesSubset (image, NULL, Nimage); - // - // // load the flat-field correction table from CATDIR - // char flatcorrFile[1024]; - // snprintf (flatcorrFile, 1024, "%s/flatcorr.fits", CATDIR); - // FlatCorrectionTable *flatcorr = FlatCorrectionLoad (flatcorrFile, VERBOSE); UpdateObjectOffsets (skylist, HOST_ID, HOSTDIR); + gfits_db_free (&db); + freeImages(); + freeMosaics (); + + relastro_client_free (sky, skylist); + break; + } + + case OP_HIGH_SPEED: { + // USAGE: relastro_client -high-speed + high_speed_catalogs (sky, skylist, HOST_ID, HOSTDIR); + break; + } + + case OP_HPM: { + // USAGE: relastro_client -high-speed + hpm_catalogs (sky, skylist, HOST_ID, HOSTDIR); break; } Index: trunk/Ohana/src/relastro/src/relastro_images.c =================================================================== --- trunk/Ohana/src/relastro/src/relastro_images.c (revision 39402) +++ trunk/Ohana/src/relastro/src/relastro_images.c (revision 39457) @@ -12,4 +12,8 @@ set_db (&db); gfits_db_init (&db); + + // final pass fit mode (use FIT_AVERAGE only for the image loop) + int finalPassMode = FIT_MODE; // start with the globally-defined fit mode + FIT_MODE = FIT_AVERAGE; /* lock and load the image db table */ @@ -77,4 +81,9 @@ MARKTIME("update chips: %f sec\n", dtime); } + + // measure scatter for stacks + UpdateStacks (catalog, Ncatalog); + MARKTIME("UpdateStacks: %f sec\n", dtime); + // create summary plots of the process // relastroVisualSummaryChips(); @@ -93,8 +102,19 @@ } - if (!UPDATE) exit (0); + if (!UPDATE) { + freeStarMaps(); + gfits_db_free (&db); + ohana_memcheck (VERBOSE); + ohana_memdump (VERBOSE); + exit (0); + } // free the image / measurement pointers - freeImageBins (1); + freeImageBins (Ncatalog); + for (i = 0; i < Ncatalog; i++) { + dvo_catalog_free (&catalog[i]); + } + free (catalog); + freeMosaics (); // If we did NOT use all images, then we applied the measured corrections to a subset of @@ -105,4 +125,5 @@ reload_images (&db); } + freeStarMaps(); if (PARALLEL) { @@ -116,4 +137,6 @@ } dvo_image_unlock (&db); + gfits_db_free (&db); + freeImages(); } @@ -123,7 +146,9 @@ // if we do NOT have parallel partitions, we must NOT write them out yet: the act of - // writing out the files byte-swaps the data and makes the values invalide for the + // writing out the files byte-swaps the data and makes the values invalid for the // following calls to UpdateObjectOffsets (which attempt to apply the image values from // the structure in memory) + + FIT_MODE = finalPassMode; // use the user-selected mode for the final pass // iterate over catalogs to make detection coordinates consistant @@ -141,6 +166,8 @@ } dvo_image_unlock (&db); + gfits_db_free (&db); + freeImages(); } - exit (0); + return TRUE; } Index: trunk/Ohana/src/relastro/src/relastro_merge_source.c =================================================================== --- trunk/Ohana/src/relastro/src/relastro_merge_source.c (revision 39402) +++ trunk/Ohana/src/relastro/src/relastro_merge_source.c (revision 39457) @@ -22,10 +22,11 @@ // load data from each region file, only use bright stars for (i = 0; i < sky[0].Nregions; i++) { - if (sky[0].regions[i].index == CAT_ID_SRC) { + myAssert (sky[0].regions[i].index >= 0, "oops"); + if ((unsigned int) sky[0].regions[i].index == CAT_ID_SRC) { catalog_src.filename = sky[0].filename[i]; region_src = &sky[0].regions[i]; } // currently, we only accept dst == src... - if (sky[0].regions[i].index == CAT_ID_DST) { + if ((unsigned int) sky[0].regions[i].index == CAT_ID_DST) { catalog_dst.filename = sky[0].filename[i]; // region_dst = &sky[0].regions[i]; Index: trunk/Ohana/src/relastro/src/relastro_objects.c =================================================================== --- trunk/Ohana/src/relastro/src/relastro_objects.c (revision 39402) +++ trunk/Ohana/src/relastro/src/relastro_objects.c (revision 39457) @@ -84,5 +84,7 @@ populate_tiny_values(&catalog, DVO_TV_MEASURE); - UpdateObjects (&catalog, 1, 0); + // the 3rd argument (-1) is the loop number for re-weighting 2MASS and Tycho. -1 is a special value meaning "ignore" + // we do NOT want to apply special weights to 2MASS and/or Tycho when calculating object final motions. + UpdateObjects (&catalog, 1, -1); free_tiny_values(&catalog); @@ -215,4 +217,5 @@ table->hosts[i].pid = pid; // save for future reference } + free (command); } @@ -224,4 +227,5 @@ HostTableWaitJobsGetIO (table, __FILE__, __LINE__, VERBOSE); } + FreeHostTable (table); return TRUE; Index: trunk/Ohana/src/relastro/src/relastro_parallel_images.c =================================================================== --- trunk/Ohana/src/relastro/src/relastro_parallel_images.c (revision 39402) +++ trunk/Ohana/src/relastro/src/relastro_parallel_images.c (revision 39457) @@ -16,4 +16,8 @@ INITTIME; + + // over-ride user selection here: this program does not set the final object + // astrometry. since it only sets the image parameters it should not use pm or par + FIT_MODE = FIT_AVERAGE; client_logger_init (CATDIR); @@ -141,4 +145,9 @@ LOGRTIME("slurp_meas_pos loop %d on %s, host %d: %f sec\n", i, myHostName, REGION_HOST_ID, dtime); } + + // measure scatter for stacks + UpdateStacks (catalog, Ncatalog); + LOGRTIME("UpdateStacks on %s, host %d: %f sec\n", myHostName, REGION_HOST_ID, dtime); + // create summary plots of the process // relastroVisualSummaryChips(); @@ -176,4 +185,22 @@ LOGRTIME("share image pos loop %d on %s, host %d: %f sec\n", -1, myHostName, REGION_HOST_ID, dtime); + // free the image / measurement pointers + freeImageBins (Ncatalog); + for (i = 0; i < Ncatalog; i++) { + dvo_catalog_free (&catalog[i]); + } + free (catalog); + freeCatalogIndexes(Ncatalog); + + freeMosaics (); + + freeStarMaps(); + + freeImages(); + free (image); + + FreeRegionHostTable (regionHosts); + relastro_free (sky, skylist); + exit (0); } Index: trunk/Ohana/src/relastro/src/relastro_parallel_regions.c =================================================================== --- trunk/Ohana/src/relastro/src/relastro_parallel_regions.c (revision 39402) +++ trunk/Ohana/src/relastro/src/relastro_parallel_regions.c (revision 39457) @@ -83,5 +83,14 @@ } +// testjump: + MARKTIME ("finished relastro -parallel-regions: %f sec total\n", dtime); + + gfits_db_free (&db); + freeImages(); + freeMosaics (); + + FreeRegionHostTable (regionHosts); + relastro_free (sky, skylist); exit (0); Index: trunk/Ohana/src/relastro/src/share_icrf_obj.c =================================================================== --- trunk/Ohana/src/relastro/src/share_icrf_obj.c (revision 39402) +++ trunk/Ohana/src/relastro/src/share_icrf_obj.c (revision 39457) @@ -99,4 +99,5 @@ int cat, ave, meas; ICRFdata (i, &cat, &ave, &meas); + myAssert (cat < Ncatalog, "oops"); Average *average = &catalog[cat].average[ave]; Index: trunk/Ohana/src/relastro/src/share_images_pos.c =================================================================== --- trunk/Ohana/src/relastro/src/share_images_pos.c (revision 39402) +++ trunk/Ohana/src/relastro/src/share_images_pos.c (revision 39457) @@ -141,8 +141,12 @@ AstromOffsetTableFree(table); + free (table); LOGRTIME("image_maps_free host %d loop %d on %s, host %d: %f sec\n", i, nloop, myHostName, REGION_HOST_ID, dtime); } // apply the modified image parameters to my detections + // This loop is skipped because the updated image parameters are applied to + // all measurements in UpdateChips (and those modified measurements are saved + // at the same time) for (i = 0; FALSE && (Ncatalog > 0) && (i < Nimage_pos); i++) { off_t seq = getImageByID (image_pos[i].imageID); Index: trunk/Ohana/src/relastro/src/share_mean_pos.c =================================================================== --- trunk/Ohana/src/relastro/src/share_mean_pos.c (revision 39402) +++ trunk/Ohana/src/relastro/src/share_mean_pos.c (revision 39457) @@ -107,5 +107,6 @@ continue; } - + myAssert (catSeq < Ncatalog, "oops"); + catalog[catSeq].average[objSeq].R = meanpos[i].R; catalog[catSeq].average[objSeq].D = meanpos[i].D; Index: trunk/Ohana/src/relastro/src/share_meas_pos.c =================================================================== --- trunk/Ohana/src/relastro/src/share_meas_pos.c (revision 39402) +++ trunk/Ohana/src/relastro/src/share_meas_pos.c (revision 39457) @@ -58,8 +58,4 @@ off_t i, k; - int Nmeaspos = 0; - MeasPos *measpos = NULL; - ALLOCATE (measpos, MeasPos, 1); - fprintf (stderr, "grabbing meas object pos from other hosts...\n"); @@ -78,40 +74,39 @@ LOGRTIME("meas_load_sync host %d loop %d on %s, host %d: %f sec\n", i, nloop, myHostName, REGION_HOST_ID, dtime); - off_t Nsubset = 0; + off_t Nmeaspos = 0; char *posfile = make_filename (CATDIR, regionHosts->hosts[i].hostname, regionHosts->hosts[i].hostID, "measpos.fits"); - MeasPos *measposSubset = MeasPosLoad (posfile, &Nsubset); + MeasPos *measpos = MeasPosLoad (posfile, &Nmeaspos); free (posfile); LOGRTIME("MeasPosLoad host %d loop %d on %s, host %d: %f sec\n", i, nloop, myHostName, REGION_HOST_ID, dtime); - // merge_meas_pos reallocs measpos and frees the input measposSubset - measpos = merge_meas_pos (measpos, &Nmeaspos, measposSubset, Nsubset); - LOGRTIME("merge_meas_pos host %d loop %d on %s, host %d: %f sec\n", i, nloop, myHostName, REGION_HOST_ID, dtime); + // apply the loaded measurement positions: + for (off_t j = 0; j < Nmeaspos; j++) { + int objID = measpos[j].objID; + int catID = measpos[j].catID; + + // set the meas mag + int catSeq; + off_t objSeq; + if (!catID_and_objID_to_seq (catID, objID, &catSeq, &objSeq)) { + // XXX what should I do if this does not match? + continue; + } + myAssert (catSeq < Ncatalog, "oops"); + + off_t m = catalog[catSeq].average[objSeq].measureOffset; + for (k = 0; k < catalog[catSeq].average[objSeq].Nmeasure; k++, m++) { + if (catalog[catSeq].measureT[m].imageID != measpos[j].imageID) continue; + catalog[catSeq].measureT[m].R = measpos[j].R; + catalog[catSeq].measureT[m].D = measpos[j].D; + if (catalog[catSeq].measure) { + catalog[catSeq].measure[m].R = measpos[j].R; + catalog[catSeq].measure[m].D = measpos[j].D; + } + } + } + free (measpos); + LOGRTIME("match_catID_and_objID/meas loop %d on %s, host %d: %f sec\n", nloop, myHostName, REGION_HOST_ID, dtime); } LOGRTIME("meas_load_done loop %d on %s, host %d: %f sec\n", nloop, myHostName, REGION_HOST_ID, dtime); - - for (i = 0; i < Nmeaspos; i++) { - int objID = measpos[i].objID; - int catID = measpos[i].catID; - - // set the meas mag - int catSeq; - off_t objSeq; - if (!catID_and_objID_to_seq (catID, objID, &catSeq, &objSeq)) { - // XXX what should I do if this does not match? - continue; - } - off_t m = catalog[catSeq].average[objSeq].measureOffset; - for (k = 0; k < catalog[catSeq].average[objSeq].Nmeasure; k++, m++) { - if (catalog[catSeq].measureT[m].imageID != measpos[i].imageID) continue; - catalog[catSeq].measureT[m].R = measpos[i].R; - catalog[catSeq].measureT[m].D = measpos[i].D; - if (catalog[catSeq].measure) { - catalog[catSeq].measure[m].R = measpos[i].R; - catalog[catSeq].measure[m].D = measpos[i].D; - } - } - } - free (measpos); - LOGRTIME("match_catID_and_objID/meas loop %d on %s, host %d: %f sec\n", nloop, myHostName, REGION_HOST_ID, dtime); fprintf (stderr, "DONE grabbing meas object pos from other hosts...\n"); Index: trunk/Ohana/src/relastro/test/HostTable.dat.in =================================================================== --- trunk/Ohana/src/relastro/test/HostTable.dat.in (revision 39457) +++ trunk/Ohana/src/relastro/test/HostTable.dat.in (revision 39457) @@ -0,0 +1,3 @@ +1 HOSTNAME CATDIR.d1 +2 HOSTNAME CATDIR.d2 +3 HOSTNAME CATDIR.d3 Index: trunk/Ohana/src/relastro/test/RegionHosts.dat.in =================================================================== --- trunk/Ohana/src/relastro/test/RegionHosts.dat.in (revision 39457) +++ trunk/Ohana/src/relastro/test/RegionHosts.dat.in (revision 39457) @@ -0,0 +1,2 @@ +1 HOSTNAME 320.0 340.0 -10.0 -0.5 +2 HOSTNAME 320.0 340.0 -0.5 10.0 Index: trunk/Ohana/src/relastro/test/gdb.sh =================================================================== --- trunk/Ohana/src/relastro/test/gdb.sh (revision 39457) +++ trunk/Ohana/src/relastro/test/gdb.sh (revision 39457) @@ -0,0 +1,8 @@ + +define t0 + set $i = 0 + while ($i < Nnomask) + printf "%d : %f %f %f %f %f\n", $i, fitStats->sample[$i].X, fitStats->sample[$i].Y, fitStats->sample[$i].R, fitStats->sample[$i].D, fitStats->sample[$i].T + set $i = $i + 1 + end +end Index: trunk/Ohana/src/relastro/test/run.fakeastro.par =================================================================== --- trunk/Ohana/src/relastro/test/run.fakeastro.par (revision 39457) +++ trunk/Ohana/src/relastro/test/run.fakeastro.par (revision 39457) @@ -0,0 +1,39 @@ +#!/bin/csh -f + +set REALRUN = 1 + +set galaxy = "g0" +set images = "i0" +set region = "330 336 -4 2" +set catdir = "catdir.par" + +rm -rf $catdir.$galaxy $catdir.$images + +if ($REALRUN) then + alias GO 'echo === \!* ===; \!*; if ($status) exit 2' +else + alias GO 'echo === \!* ===' +endif + +alias NO 'echo === skip : \!* ===' + +if (-f template.dat == 0) then + echo "missing template template.dat" + exit 2 +endif + +# make the basic fake data +GO fakeastro -galaxy -region $region -uniform-radec -D CATDIR $catdir.$galaxy -D CATFORMAT PS1_V5 -D SKY_DEPTH 4 -D FAKEASTRO_NSTARS 100000 -D FAKEASTRO_NQSO_ICRF 200 -D FAKEASTRO_NQSO_ZERO 2000 -testing +GO fakeastro -images -input $catdir.$galaxy -output $catdir.$images -input-images template.dat -region $region -D CATFORMAT PS1_V5 -D SKY_DEPTH 3 + +set myHost = `hostname` +set myPath = `pwd` +cat HostTable.dat.in | sed s/HOSTNAME/$myHost/ | sed "s|CATDIR|$myPath/$catdir.$images|" > $catdir.$images/HostTable.dat +# XXX catdir currently needs to include the full path + +GO dvodist -out $catdir.$images + +GO relphot -parallel -averages -region $region -D CATDIR $catdir.$images -update + +GO addstar -parallel -resort -region 320 340 -10 10 -D CATDIR $catdir.$images + Index: trunk/Ohana/src/relastro/test/run.fakeastro.ser =================================================================== --- trunk/Ohana/src/relastro/test/run.fakeastro.ser (revision 39457) +++ trunk/Ohana/src/relastro/test/run.fakeastro.ser (revision 39457) @@ -0,0 +1,37 @@ +#!/bin/csh -f + +set REALRUN = 1 + +set galaxy = "g0" +set images = "i0" +set region = "330 336 -4 2" +set catdir = "catdir.ser" + +rm -rf $catdir.$galaxy $catdir.$images + +if ($REALRUN) then + alias GO 'echo === \!* ===; \!*; if ($status) exit 2' +else + alias GO 'echo === \!* ===' +endif + +alias NO 'echo === skip : \!* ===' + +if (-f template.dat == 0) then + echo "missing template template.dat" + exit 2 +endif + +# make the basic fake data +GO fakeastro -galaxy -region $region -uniform-radec -D CATDIR $catdir.$galaxy -D CATFORMAT PS1_V5 -D SKY_DEPTH 4 -D FAKEASTRO_NSTARS 100000 -D FAKEASTRO_NQSO_ICRF 200 -D FAKEASTRO_NQSO_ZERO 2000 -testing +GO fakeastro -images -input $catdir.$galaxy -output $catdir.$images -input-images template.dat -region $region -D CATFORMAT PS1_V5 -D SKY_DEPTH 3 + + + + + + +GO relphot -averages -region $region -D CATDIR $catdir.$images -update + +GO addstar -resort -region 320 340 -10 10 -D CATDIR $catdir.$images + Index: trunk/Ohana/src/relastro/test/run.relastro.par =================================================================== --- trunk/Ohana/src/relastro/test/run.relastro.par (revision 39457) +++ trunk/Ohana/src/relastro/test/run.relastro.par (revision 39457) @@ -0,0 +1,46 @@ +#!/bin/csh -f + +set REALRUN = 1 + +if ($REALRUN) then + alias GO 'echo === \!* ===; \!*; if ($status) exit 2' +else + alias GO 'echo === \!* ===' +endif + +set catdir = "catdir.par.i0" +set region = "-region 320 340 -10 10" + +set threads = "-threads 4" +set threads = "" + +set options = "" +set options = "$options -v" +set options = "$options -parallel" +set options = "$options -update" +set options = "$options -photcode ICRF_S,ICRF_C,ICRF_X,ICRF_U,ICRF_K" +set options = "$options -D RELASTRO_DPOS_MAX 120" +set options = "$options -max-mean-offset 120" +set options = "$options -D RELASTRO_SRC_MEAS_TOOFEW 3" +set options = "$options $threads" + +if (1) then + set mode = "-images -update-chips" + GO relastro -D CATDIR $catdir $region $options $mode +endif + +if (1) then + set mode = "-update-offsets -pmpar" + GO relastro -D CATDIR $catdir $region $options $mode +endif + +if (1) then + set mode = "-update-objects -pmpar" + GO relastro -D CATDIR $catdir $region $options $mode +endif + +if (1) then + set mode = "-images -update-chips -apply-offsets -pmpar" + GO relastro -D CATDIR $catdir $region $options $mode +endif + Index: trunk/Ohana/src/relastro/test/run.relastro.par.reg =================================================================== --- trunk/Ohana/src/relastro/test/run.relastro.par.reg (revision 39457) +++ trunk/Ohana/src/relastro/test/run.relastro.par.reg (revision 39457) @@ -0,0 +1,37 @@ +#!/bin/csh -f + +set REALRUN = 0 + +if ($REALRUN) then + alias GO 'echo === \!* ===; \!*; if ($status) exit 2' +else + alias GO 'echo === \!* ===' +endif + +set catdir = "`pwd`/catdir.par.i0" + +set myHost = `hostname` +cat RegionHosts.dat.in | sed s/HOSTNAME/$myHost/ > $catdir/RegionHosts.dat + +set region = "-parallel-regions -region-hosts RegionHosts.dat" + +set options = "" +set options = "$options -v" +set options = "$options -parallel" +set options = "$options -update" +set options = "$options -threads 4" +set options = "$options -photcode ICRF_S,ICRF_C,ICRF_X,ICRF_U,ICRF_K" +set options = "$options -D RELASTRO_DPOS_MAX 120" +set options = "$options -max-mean-offset 120" +set options = "$options -D RELASTRO_SRC_MEAS_TOOFEW 3" + +if (0) then + set mode = "-images -update-chips" + GO relastro -D CATDIR $catdir $region $options $mode +endif + +if (1) then + set mode = "-images -update-chips -apply-offsets -pmpar" + GO relastro -D CATDIR $catdir $region $options $mode +endif + Index: trunk/Ohana/src/relastro/test/run.relastro.ser =================================================================== --- trunk/Ohana/src/relastro/test/run.relastro.ser (revision 39457) +++ trunk/Ohana/src/relastro/test/run.relastro.ser (revision 39457) @@ -0,0 +1,44 @@ +#!/bin/csh -f + +set REALRUN = 1 + +if ($REALRUN) then + alias GO 'echo === \!* ===; \!*; if ($status) exit 2' +else + alias GO 'echo === \!* ===' +endif + +set catdir = "catdir.ser.i0" +set region = "-region 320 340 -10 10" + +set threads = "-threads 4" +# set threads = "" + +set options = "" +set options = "$options -v" +set options = "$options -update" +set options = "$options -photcode ICRF_S,ICRF_C,ICRF_X,ICRF_U,ICRF_K" +set options = "$options -D RELASTRO_DPOS_MAX 120" +set options = "$options -max-mean-offset 120" +set options = "$options -D RELASTRO_SRC_MEAS_TOOFEW 3" +set options = "$options $threads" + +if (1) then + set mode = "-images -update-chips" + GO relastro -D CATDIR $catdir $region $options $mode +endif + +if (1) then + set mode = "-update-offsets" + GO relastro -D CATDIR $catdir $region $options $mode +endif + +if (1) then + set mode = "-update-objects -pmpar" + GO relastro -D CATDIR $catdir $region $options $mode +endif + +if (1) then + set mode = "-images -update-chips -apply-offsets -pmpar" + GO relastro -D CATDIR $catdir $region $options $mode +endif Index: trunk/Ohana/src/relastro/test/template.dat =================================================================== --- trunk/Ohana/src/relastro/test/template.dat (revision 39457) +++ trunk/Ohana/src/relastro/test/template.dat (revision 39457) @@ -0,0 +1,722 @@ +SIMPLE = T / BITPIX = 8 / NAXIS = 0 / BSCALE = 1.0000000000 / BZERO = 0.0000000000 / EXTEND = T / NIMAGES = 244 / ZERO_PT = 25.0000000000 / DVO_DBID= '2ef3f3a7e738b596b17bf35a643f296d' / FORMAT = 'PS1_V5 ' / IMAGEID = 245 / END XTENSION= 'BINTABLE' / BITPIX = 8 / NAXIS = 2 / NAXIS1 = 360 / NAXIS2 = 244 / PCOUNT = 0 / GCOUNT = 1 / TFIELDS = 55 / EXTNAME = 'DVO_IMAGE_PS1_V5' / TTYPE1 = 'CRVAL1 ' / coordinate at reference pixel TUNIT1 = ' ' / TFORM1 = 'D ' / TTYPE2 = 'CRVAL2 ' / coordinate at reference pixel TUNIT2 = ' ' / TFORM2 = 'D ' / TTYPE3 = 'CRPIX1 ' / coordinate of reference pixel TUNIT3 = ' ' / TFORM3 = 'E ' / TTYPE4 = 'CRPIX2 ' / coordinate of reference pixel TUNIT4 = ' ' / TFORM4 = 'E ' / TTYPE5 = 'CDELT1 ' / degrees per pixel TUNIT5 = ' ' / TFORM5 = 'E ' / TTYPE6 = 'CDELT2 ' / degrees per pixel TUNIT6 = ' ' / TFORM6 = 'E ' / TTYPE7 = 'PC1_1 ' / rotation matrix TUNIT7 = ' ' / TFORM7 = 'E ' / TTYPE8 = 'PC1_2 ' / rotation matrix TUNIT8 = ' ' / TFORM8 = 'E ' / TTYPE9 = 'PC2_1 ' / rotation matrix TUNIT9 = ' ' / TFORM9 = 'E ' / TTYPE10 = 'PC2_2 ' / rotation matrix TUNIT10 = ' ' / TFORM10 = 'E ' / TTYPE11 = 'POLYTERMS' / higher order warping terms TUNIT11 = ' ' / TFORM11 = '14E ' / TTYPE12 = 'CTYPE ' / coordinate type TUNIT12 = ' ' / TFORM12 = '15A ' / TTYPE13 = 'NPOLYTERMS' / order of polynomial TUNIT13 = ' ' / TFORM13 = 'A ' / TTYPE14 = 'TZERO ' / readout time (row 0) TUNIT14 = ' ' / TFORM14 = 'J ' / TTYPE15 = 'NSTAR ' / number of stars on image TUNIT15 = ' ' / TFORM15 = 'J ' / TTYPE16 = 'SECZ ' / airmass TUNIT16 = 'mag ' / TFORM16 = 'E ' / TTYPE17 = 'NX ' / image width TUNIT17 = ' ' / TFORM17 = 'I ' / TTYPE18 = 'NY ' / image height TUNIT18 = ' ' / TFORM18 = 'I ' / TTYPE19 = 'APMIFIT ' / aperture correction TUNIT19 = 'mag ' / TFORM19 = 'E ' / TTYPE20 = 'DAPMIFIT' / apmifit error TUNIT20 = 'mag ' / TFORM20 = 'E ' / TTYPE21 = 'MCAL ' / calibration mag TUNIT21 = 'mag ' / TFORM21 = 'E ' / TTYPE22 = 'DMCAL ' / error on Mcal TUNIT22 = 'mag ' / TFORM22 = 'E ' / TTYPE23 = 'XM ' / image chisq TUNIT23 = '10*log(value)' / TFORM23 = 'I ' / TTYPE24 = 'PHOTCODE' / identifier for CCD, TUNIT24 = ' ' / TFORM24 = 'I ' / TTYPE25 = 'EXPTIME ' / exposure time TUNIT25 = 'seconds ' / TFORM25 = 'E ' / TTYPE26 = 'ST ' / sidereal time of exposure TUNIT26 = ' ' / TFORM26 = 'E ' / TTYPE27 = 'LAT ' / observatory latitude TUNIT27 = 'degrees ' / TFORM27 = 'E ' / TTYPE28 = 'RA_CENTER' / image center TUNIT28 = 'degrees ' / TFORM28 = 'E ' / TTYPE29 = 'DEC_CENTER' / image center TUNIT29 = 'degrees ' / TFORM29 = 'E ' / TTYPE30 = 'RADIUS ' / image radius TUNIT30 = 'degrees ' / TFORM30 = 'E ' / TTYPE31 = 'REF_COLOR_BLUE' / median astrometry ref color TUNIT31 = ' ' / TFORM31 = 'E ' / TTYPE32 = 'REF_COLOR_RED' / median astrometry ref color TUNIT32 = ' ' / TFORM32 = 'E ' / TTYPE33 = 'NAME ' / name of original image TUNIT33 = ' ' / TFORM33 = '117A ' / TTYPE34 = 'DETECTION_LIMIT' / detection limit TUNIT34 = '10*mag ' / TFORM34 = 'B ' / TTYPE35 = 'SATURATION_LIMIT' / saturation limit TUNIT35 = '10*mag ' / TFORM35 = 'B ' / TTYPE36 = 'CERROR ' / astrometric error TUNIT36 = '50*arcsec' / TFORM36 = 'B ' / TTYPE37 = 'FWHM_X ' / PSF x width TUNIT37 = '25*arcsec' / TFORM37 = 'B ' / TTYPE38 = 'FWHM_Y ' / PSF y width TUNIT38 = '25*arcsec' / TFORM38 = 'B ' / TTYPE39 = 'TRATE ' / scan rate TUNIT39 = '100 usec/pixel' / TFORM39 = 'B ' / TTYPE40 = 'CCDNUM ' / CCD ID number TUNIT40 = ' ' / TFORM40 = 'B ' / TTYPE41 = 'FLAGS ' / image quality flags TUNIT41 = ' ' / TFORM41 = 'J ' / TTYPE42 = 'IMAGE_ID' / internal image ID TUNIT42 = ' ' / TFORM42 = 'J ' / TTYPE43 = 'PARENT_ID' / associated ref image TUNIT43 = ' ' / TFORM43 = 'J ' / TTYPE44 = 'EXTERN_ID' / external image ID TUNIT44 = ' ' / TFORM44 = 'J ' / TTYPE45 = 'SOURCE_ID' / analysis source ID TUNIT45 = ' ' / TFORM45 = 'I ' / TTYPE46 = 'NLINK_ASTROM' / mean number of matched measurements for astrometry TUNIT46 = ' ' / TFORM46 = 'I ' / TTYPE47 = 'NLINK_PHOTOM' / mean number of matched measurements for astrometry TUNIT47 = ' ' / TFORM47 = 'I ' / TTYPE48 = 'UBERCAL_DIST' / distance to nearest ubercal image TUNIT48 = ' ' / TFORM48 = 'I ' / TTYPE49 = 'XPIX_SYS_ERR' / systematic astrometry error in X TUNIT49 = ' ' / TFORM49 = 'E ' / TTYPE50 = 'YPIX_SYS_ERR' / systematic astrometry error in Y TUNIT50 = ' ' / TFORM50 = 'E ' / TTYPE51 = 'MAG_SYS_ERR' / systematic photometry error TUNIT51 = ' ' / TFORM51 = 'E ' / TTYPE52 = 'N_FIT_ASTROM' / number of stars used for astrometry cal TUNIT52 = ' ' / TFORM52 = 'I ' / TTYPE53 = 'N_FIT_PHOTOM' / number of stars used for photometry cal TUNIT53 = ' ' / TFORM53 = 'I ' / TTYPE54 = 'PHOTOM_MAP_ID' / reference to 2D zero point map TUNIT54 = ' ' / TFORM54 = 'J ' / TTYPE55 = 'ASTROM_MAP_ID' / reference to 2D astrometry map TUNIT55 = ' ' / TFORM55 = 'J ' / END @tÝ{ځ•¿ñ{ð'ÝÀý> +ª*8”`ì8”hóŸv™?xj¿xRŸvð.·²Û +žW#+7Åøž9ž +JF7Š@Ž9ÁhU:äYEºõÃU: ÷9ߘf:ʬºæ/!9äӕDEC--DISNk[œr¡uÀÀo5745g0516o.356887.cm.982631.smf[PHU]ÆeˆÄÆ +^q?†3,? +çó?üt<*Ô$=/Š’?ƵåÇ +µö,øµÉQº¶ +¶3­"‚µªÕœ-íÔÌ-è¶Ó. +Œ.a +-ÁDa-õàß­>2†-ÓÖ;DEC--WRPNk[q?Œ1îÀÿ'ÙB4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY01.hdr]œœF !ÆgÁÅ€~ò?ƒ‡¡?‚z?ÿò:š;¡ö“?ÿ0µ•mŽãMX1Ñ+Žä®‹²ÙnWŠ-Ç2-– +» +-2=O,Ðj—-=Î)¬†U9.‰kDEC--WRPNk[ ?Œ1îÀÿ'ÚB4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY02.hdr]¿œF!!Æi +Ã>/?‚Cj?‚ + +?ÿ6»¡‡Œ0O|?ü2µ,X5/ŽÞ5 +1µVn•ŽŒœN³  +É-†ãt­c‘p­¡©ê-ä(®-wl­{Z-6åg. +ØDEC--WRPNk[ +Ö?Œ1îÀÿ'ÛB4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY03.hdr]ÀœF"!Æj8þEœ©ì?€‘5?‹ ?ÿï:»ö8Žê?€3³Ûy±f'öŽ")”Ž”(Ü3-‚Ó4èŸ,žv*ÓÞz,¥‰æ- uo- +žð¬OÂ=-Š›Y.Ò4DEC--WRPNk[ /?Œ1îÀÿ'ÜB4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY04.hdr]ÀœF#!Æf æF +؈?†Œ?‚‹?ÿÙ;W*»l8>?ÿŒ¶ +wi3NÏ>4³Ž›f4“Ÿä5Ëë.'²Œ«š‹S­‹]Û,G Q­áA­›¹s,ÀwN-mívDEC--WRPNk[ À?Œ1îÀÿ'ÝB4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY05.hdr]¿œF$!ÆdøÞFjŒ?‡o?†~^?õ°Œ‘âŒßàë?çU¶åÂ5†tµ5ŒP3¶ûϵo¶î61S.Av&­„4š­h‘C-ßÖ±. +Q®Ro'­Kî +. +ÃDEC--WRPNk[f?Œ1îÀÿ'ÞB4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY06.hdr]ŸœF%!ÆûŸÆjûF?ƒÐà?‹1w?íÕ<¶ªŽ=¥Ò?⍵¯ +hµh¥­µ/Áø¶%4Þµš¬¶]“Ð-äþ©-SÈô-àƒ#.ŵ,Ÿ +•.jmÀ-ÓÄ .cédDEC--WRPNk[?Œ1îÀÿ'âB4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY10.hdr]Œ +œF&!Æ +ÕÆ#Žž?‚~ +?$ÿ?ÿìºËD?»Y +Ù?ÿ¢µ%HÒ²ïèBµî44`à·4¯ ‚Ž6²2-’[, +±-ælT¬Ž™N¬Lºt+Ø +¬ì‘«,žaóDEC--WRPNk[ó?Œ1îÀÿ'ãB4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY11.hdr]¿ + œF'!Æ,hÅ£• +?‚Vi?‚h8?ÿõ:˜™K<î¶?ýȵ4?3µa˜²Z +µC{îŽCqµtŒ-ƒ2-%§o«Æ_†-Æ[|,×f×-rQ,8-›ü DEC--WRPNk[Ù?Œ1îÀÿ'äB4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY12.hdr]¿ + +œF(!Æ + +óÃ'™?€W?KY?€9nú»"À?ÿÍ5 +é3zº+4Qás2óD”3Ì×öŽÉIu­V0ò«à/­²Ë¬ ˆ +¬ŸÔ`- +šŠ¬­ßÜ.X^DEC--WRPNk[ +A?Œ1îÀÿ'åB4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY13.hdr]¿ + +œF)!Æ +f4Eœ0?€ô?£?ÿ×»ԋ»=•˜?ÿ»5;(y4RÛî4qÚSŽ€SŽ.Ìö5늭€_Ϭœ8¬Ø•, +Š-Hç +­Ñmh+ðM-ÂzªDEC--WRPNk[7?Œ1îÀÿ'æB4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY14.hdr]À +œF*!Æ­üF +ÐJ?ƒ-1?‚2[?þ‘»ÚmœŒJw9?úìµ +‡j5dOÕ5Uûµ=|[Ží­€ +-%O­OË +­-.c­õ-È +'DEC--WRPNk[h?Œ1îÀÿ'èB4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY16.hdr]ŸœF,!Æ +ZFÂX?„y?œî?èیᣬœÐ§2?~‘Œµ» +:6‰ + 5}”ü¶Ýò]µ±E +¶@’Ô.7þP®Š÷[-ÂÈ.ìž«.ç3¯à +­€^ +®±»öDEC--WRPNk[û?Œ1îÀÿ'éB4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY17.hdr]œœF-!śòíÆp`Ž?ƒºÆ? +‰ý?ÿ»«Ð<¢ñ?þ µ[æ-³üPµÆ])µk!å4õí +µãƒ7-` ΩÇå3. (-1×ã-üèz-Ä¿ +­7C.L”DEC--WRPNk[?Œ1îÀÿ'ìB4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY20.hdr]ŸœF.!ŞþÆ"*"?€ç?‚Ge?ü±Œ"|M;І)?þµ2ˆ Ž0šá5T Iµ<å4ä15Žý -<ƒS©{/€­ÛÕ-aK€­˜¡-²$h­Œº- +ýDEC--WRPNk[ž?Œ1îÀÿ'íB4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY21.hdr]¿œF/!Ŝ=bÅ€§É?­r?—ò?€9AÖÏ:íyï?ÿäŽÚý³Ö|„Ž;ÁM³FþƒŽƒøŽÃb[-ŽE,û¹,ºœR«€OL,s_¹,Ö~^-G +-\ÞÿDEC--WRPNk[ °?Œ1îÀÿ'îB4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY22.hdr]¿œF0!Ŝ)Ðà +U2?‚f?4?ÿû:B¢àº“÷›?ÿõŽH,³ŠÜôŽ/Œ4Ðú4{¢ë3Õ:ã«àEÔ,íŠ}-©³)­ÙÖY­òoã­Œ‘k­HÜ +-:ODEC--WRPNk[ ƒ?Œ1îÀÿ'ïB4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY23.hdr]¿œF1!śø?E5Ð?·â?põ?ÿÛ» +-»ÖÆ?ÿÝŽ°là4•E +4‹^£ŽŒ³Æ,Z³Âú,,ÓÇ»­5<­˜óå-M¥5¬ö*° +@­!Å­­Œ×,DEC--WRPNk[ t?Œ1îÀÿ'ðB4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY24.hdr]¿œF2!Ŝe)F +ß¹?§“?€Ÿ +?ÿž;@Œ„:î?ÿýŽÑ.ó³.!Ö2ÏÆ4‹ìM46™/Žœ'A-fÀS¬W† +¬äšã¬^h-fN,«n­žDEC--WRPNk[ e?Œ1îÀÿ'ñB4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY25.hdr]ŸœF3!Ŝ +£FnåØ?©¶?Tœ?ÿýº&»|³ +?ÿƒ³˜É4®jS4Öd4³õ +гï`3_n"¬šŒÇ¬Ëú¿­Kõ±-™·±,å~‹,¢ œªŸlˆ,$>DEC--WRPNk[ =?Œ1îÀÿ'òB4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY26.hdr]ŸœF4!ŚºFŸry?‚Ž³?‚«å?þÙ»Âô¹ñÎ,?ÿþµVÉ4–"Â5có4³™©ïŽWŽ5Pµ-Ìä­‚«²W-Òô-Ôú{¬=åô+ +ys-–ÛDEC--WRPNk[3?Œ1îÀÿ'óB4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY27.hdr]ŒœF5!ÀêÆróž?šk?Ãm?ðÏŒ°&¥;{ÐÒ?ÿ„³jZ 3œ¥CŽ³¹£ŽÈjŸ5ž§Ž€ê-¡_-&ØZ¬`J­kˆ-$é|-npó­Ìœ,žn*DEC--WRPNk[ t?Œ1îÀÿ'öB4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY30.hdr]¿ +œF6!øžÆ#]?Ý¶?€G?ú’ŒUί:æ©Ã?ÿ厝–Ž=]€µ(S3³OÏ65}Y5/ù-kRÆ+ÓÝT-™ +ý-[O®-žh«=k­Þ”­’ÓÚDEC--WRPNk[^?Œ1îÀÿ'÷B4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY31.hdr]¿œF7!Á¬°þÅ¥>U?}N?ÇÉ?ÿo»‡æ¶žßˆŠ?€³/”³œª|ŽËÛŽœáÒ5Xñ1ŽÔ=c-”Ä«µŒÎ­wmî,¢z-Á +n-‰qŽ®+9-- +ªÃDEC--WRPNk[ v?Œ1îÀÿ'øB4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY32.hdr]Ÿ œF8!Áªg=ÃŽè?/H?8­?€¹Úæ¹íK?ÿþ4[‚£3Àš:4UÆm²xƒ 3c 3Kȱ­†¯ý­‡Ü¿®Šê, +)ž­1¬ëO¹¬¶‚ +,Øÿ&DEC--WRPNk[ +?Œ1îÀÿ'ùB4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY33.hdr]œ!œF9!Á™U +Eœ^Æ??\?,™?ÿý:©>¹ +KÉ?€3¥gáŽfJ³6 +Ð4‰C’3¯Ú1ŽX!X«œp- +”",óuÌ­ +­ìå-Ðk-2÷Ú­ñÅDEC--WRPNk[Å?Œ1îÀÿ'úB4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY34.hdr]œ" +œF:!Áƒ±¶F +“Ñ?[?€+?ÿçºâA—ºú?ë?ÿâ3k*³¥jM40¡µ +–Ž&4ƒÜ,; +,Á^o-U!Û­.ו,Ïq®+"ý¬ë + +,J[DEC--WRPNk[ o?Œ1îÀÿ'ûB4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY35.hdr]Ÿ# +œF;!@®wwFo©)?v?€XÄ?ý±Œ +°òºUM?ÿöŽ;‰24ÇB4Ë9Ž +&µ1|aŽù›¬ÀŸý¬×­`ï.-5þ*®ð¬Ê˜þ­–Ÿï­G·DEC--WRPNk[ €?Œ1îÀÿ'üB4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY36.hdr]Ÿ$ +œF<!ÁÔÅFŸYÃ?•j?‚Ìn?€8ŠÕ¢ºÚ?ÿ鮕šà4ð4^QìŽ%û3ÓS5bÖ-“Xg+›MŽ­ÛK-%üò,4÷[¬±‹+Œ0Ù-¢!ðDEC--WRPNk[s?Œ1îÀÿ'ýB4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY37.hdr]œ%œF=!@³#òFžÍa?‚ +R?ƒÈ¿ø<|LZ;mzÔ¿ÿ•³•c;ŽTLðµFC5í5\ +µ–Î$¬ÍÎU­LÀŽ­ÐÞ­šûÙ­±õ%-|‰ +-‡F7­»’ùDEC--WRPNk[ ‰?Œ1îÀÿ(B4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY40.hdr]œ( œF>!°nFnâ ?\Â?wÏ¿ÿ™;ec¹² +,¿ÿÿ4<ï1ÉQ]³ŠÌ•Ž=€o4ÀxÕŽ‚šªb{¬ç„ã-UÈÉ«À¥N+ÀŒ¬ô%_-;#)­ ]]DEC--WRPNk[ï?Œ1îÀÿ(B4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY41.hdr]¿)!œF?!®§ÑF +l?$¶?és¿ÿÿ9Û2&·ÃŸ[¿€Žžc³ +pÎ4(‹4ªžz³áÚµÀŸ,né¬$ž[­Äœ¬ƒÆb-yҙ-—C€,‡2~­¡UNDEC--WRPNk[ u?Œ1îÀÿ(B4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY42.hdr]Ÿ*"œF@!·#EœMž?G€?3û¿ÿù:nœº”W¿ÿõ2d+°3ŒÌè² +<µŽ +ͳž{4 +/k¬5Þ*¬ÅKc,‚ +)Ô++ÉUE­Ä4¬Í­E,“72DEC--WRPNk[Æ?Œ1îÀÿ(B4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY43.hdr]Ÿ+#œFA!³¢åÃb$?Há?>é¿ÿý:'0K¹r©×¿€2(%2Ï0,ŽWc²/Xã3•E<2ð‰@,Ž¬t5­5Vl¬Žæ”. +ŠÌ,.† +­i +­3òDEC--WRPNk[¿?Œ1îÀÿ(B4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY44.hdr]œ,$œFB!ºIûÅ€Ía?pP?‚«¿ÿúº`%"»š¿ÿÞ4fUþ3T\˜4y 4}Sö4csò58Ù­œw@,;m®­@¿¬Öÿº­./Ô­7Ži¬è—=­ŒèžDEC--WRPNk[æ?Œ1îÀÿ(B4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY45.hdr]Ÿ-%œFC!Âå¯ØÆ#TR?o?€è¿ÿ;t::ÎŒ¿ÿëŽ#s³­¹®4„‰5ŽÜ8™Ž!Ï,³ýý­.H–,ɍ-œ,—'Œ­.qï-j„,“Ý?,HÑDEC--WRPNk[ I?Œ1îÀÿ(B4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY46.hdr]Ÿ.&œFD!òBÆo˜i?‰Î?‡"ã¿ÿˆ;mLò¹÷ÕV¿ÿþ4€|4$q¯4‹¯ ³Ý?·³*"6+À­e¹¬lÿ‰¬1šÄ­‘¬ð¯ö,{>*Ӌ®WWŒDEC--WRPNk[k?Œ1îÀÿ(B4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY47.hdr]Œ/'œFE!ŝ×ÏF$Ý?5Ð?‡÷þ¿ÿ¥;M?*<‡Å¬¿÷ß3ΏLµ5D}³Ã'ù5Ézš4]:ƶ'ž­<Ó+N‡­CÕD®0°i­!ûS-ð^÷,œà!®=5°DEC--WRPNk[ ²?Œ1îÀÿ( +B4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY50.hdr]œ2(œFF!Ş>ýFnkr?‚M=?‚S¿ÿ!»šÉŸ:·í-¿ÿï5u‰1²¿*yŽc4è4¶صúûµ%:®o\«`=-}†&¬«_¬, +Eß-Xµ­‚݆­vnÆDEC--WRPNk[ +\?Œ1îÀÿ( +B4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY51.hdr]Ÿ3)œFG!ŞdàF +ZH?ÐÌ?øc¿ÿ‹»tJÏ:ýP"¿ÿá5îaŽ†gfŽ~ŠÏ4'\ÚµsSµ9XH­ƒ^Ê,ý +H-/­^w@+É"ƪ«Æ­ÇŽ{­õdfDEC--WRPNk[Õ?Œ1îÀÿ( +B4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY52.hdr]Ÿ4*œFH!ŞUE›ê»? +ü?m^¿ÿòº©ÊH;d\¿ÿ›Žn¶ÖŽÝÙø4R4$§“Žæ‘ÎŽ2U-‚ +,-žì«}$)fÒB.BU-t¢­Ïš×­ÄXDEC--WRPNk[«?Œ1îÀÿ( +B4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY53.hdr]¿5+œFI!ŝàþÃÛ¿?òu?L¶¿ÿûºHQv»XŠ +¿ÿ£5" +[4ÔÙ4â¢2Ze +³\®52dve­±í“Nà­cÌ"+³®À­«[l-ðÚ.2u5­ýÎDEC--WRPNk[c?Œ1îÀÿ(B4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY54.hdr]¿6,œFJ!ŞÈÅ€æ6?–ö?‚Jà¿ÿŠ;U«»7ê¿ÿ¿4Éá‚4€W]³Pš4ü[ŽÑQY5ZF­¥ +8¬š«-_‡ +­™ +¬þ×譋 +5-¡<­åð:DEC--WRPNk[ +?Œ1îÀÿ(B4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY55.hdr]¿7-œFK!Ŝ +Æ"ÐÆ?Ã?¯ +¿ùŒm +Û»q{"¿ÿŽ5q1Õ©Ž¯>5/1]5®|^3ËHo­Š³‚,€R­ +D£­*Ÿ°-ŀ­³Æê®'ý×,”ŠDEC--WRPNk[º?Œ1îÀÿ(B4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY56.hdr]œ8.œFL!şŸèÆoÁ?–Œ?†l¶¿ÿ;f§ì»¯>¿ÿ!51Î32a­³­0Š4ÿr§2vŒ#6%;­Mc+þD¢­Œ±­Pc,ë"z­x,Ù¬Ú +®0äÐDEC--WRPNk[#?Œ1îÀÿ(B4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY57.hdr]œ9/œFM!ÆKcFžÀ¹? +íC?ƒÛ›¿þž;Ð*jµšX­×Q«- +V".!•`­å3œ®E£ +-š­.wÞ­ÑwªDEC--WRPNk[¬?Œ1îÀÿ(B4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY60.hdr]µ<0œFN!Æ +’'Fp†À?€ŒÄ?€1T¿þ/»ô‘gŒ"á9¿üŸ³^>.4\?S5þ±µ6áVŽð4VÖœ*5{à+q&;«/èW-xåu-¶‚z­ \߬ÂT-+ŸFÂDEC--WRPNk[n?Œ1îÀÿ(B4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY61.hdr]Œ=1œFO!ÆØÊF +(=?‚U–?‚ +{¿ýñ<Â;ûÕ¿þ5}°µ5D]µ0ŽÄ4Úý +41CŽíC­.E²- +‚`-‡ ­br-­†yë,µ|]¬ +‰­q +HDEC--WRPNk[„?Œ1îÀÿ(B4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY62.hdr]Ÿ>2œFP!Æ +ÁÒE›¿F?Ÿ?^¿ÿ;s'ö;4% ¿ÿÀ4|‡jŽ +ŒÂµ&ã%Ž7·\4Q³ +3KÇ^¬Õ%É,zÔÄ-›-,Œãy®B>­J®Š®Œ(­rÉ·DEC--WRPNk[Š?Œ1îÀÿ(B4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY63.hdr]¿?3œFQ!Æ +(Ã3 +-?€{?R8¿ÿ²»Gz;žl¿ÿ>Ž÷éñŽô&4£Ý‰2Ðñ’3 +4cÔû-aZ-2­ +_S«Å'ð­ +¶ž,c—Z-)@­ÛjDEC--WRPNk[V?Œ1îÀÿ(B4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY64.hdr]¿@4œFR!Æ +oDÅ¥Ür?‚[Þ?'¿ÿý:pÓ;6ä¿ÿŸ500wŽ* +4g +ŽÌ +Œ³› ‚4 +"$­~]*ï9 +­Höò-e6y,T™{,-¥, +z±¬ÃjÝDEC--WRPNk[ +9?Œ1îÀÿ(B4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY65.hdr]¿A5œFS!Æ +“lÆ"$‘??îª¿ý˜Œ +vŒT¿ýÁ³sІ4ÿÛ3ÜÂ4}&5W~U4ä< +«šëP­.™Þ,Ã3­K¥ð­ZÏK« ŠQ­ŠT×­[§dDEC--WRPNk[ v?Œ1îÀÿ(B4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY66.hdr]¿B6œFT!Æ؜Ænª +?„} ?‡¿üÕŒ +ŽŒžÓ¿ïü5nò56+Ÿ5óëx5ü4è4& +¿6­M­–'Û­GZ­»/­é¡X®/®:2,gRé® +† DEC--WRPNk[?Œ1îÀÿ(B4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY67.hdr]»C7œFU!ÆdFkj‡?ˆdÞ?„Ý¿ی= +9Ö<±ë)¿ï¶67åµ?g¶'²b5æ’÷5–ÁLµšÔ®ÔD-(jÀ.]sk­Ì‹ª®#1.4-xmŽ•F*5[¯4¶è3Žjw†­¿y¬Ÿ¯»,Ö® ­Ý8W­ˆh¬×Ê-@ú­é[wDEC--WRPNk[þ?Œ1îÀÿ(!B4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY73.hdr]ŸI:œFX!ÆiËZÃ¥]?‚nC?<¿þm;ä‹uº'`È¿ÿý50ÆB3»EÛµ>3L4r +}³ÁL5w8­Ú³¬<ýß-  +–¬¯>W,ÛoS­€z<­"é­® +ÅDEC--WRPNk[÷?Œ1îÀÿ("B4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY74.hdr]¿J;œFY!ÆgJ€Å£óA? +|Ð?‚Ù>¿ýCŒÏ +»‚õš¿ÿu5çÀp4Ù5f³A5`Û24ÒÔ59§‰® sg«òý¯­¡«€­’åO­Ÿ£ï­«öù¬/ú­~–ôDEC--WRPNk[?Œ1îÀÿ(#B4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY75.hdr]ŸK<œFZ!ÆcâlÆ"V±?ˆw?‚HÁ¿Ô朓µ»˜ÜŽ¿ÿ86 +_Ò4ð×J6[ +ÉŽ+U5—³M5ƒ +¹® +À­m®ˆ;­-^L®i8­U§"­&kÖ­¿ DEC--WRPNk[x?Œ1îÀÿ($B4A¬©IA¥š$ÀÀo5745g0516o.356887.cm.982631.smf[XY76.hdr]œL=œF[!@tÝ{ÖeK¿ñ{óAi4ÀŽ×²?#G8”^P8”fžŸv€€?xjæ¿x#Ÿvڄ·¥Þžqê7 ¿â·žê?ž+üW6žï­9Å#©:å º÷59÷i_9ヵ:ÌÏ«ºåú9ä]'DEC--DISNmŽœt¡„ÀÀo5745g0526o.356899.cm.982643.smf[PHU]>ÆeÇÃÆ +c?†„d?†L¡?ÿy;ƒâ4=&ÕE?Ép¶—]µå¿µ€ôͶ +³l4ásµÆšö. +-Õ©.V.\žð-S +U.[ð­ +‡ˆ-ù×õDEC--WRPNmŽe?‹"ÑîÀÿ'ÙB4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY01.hdr]Œ?>œHð!Æg5Å€öž?„Nê?‚©ª?ÿ®;OÅ:òX?ÿ㵺xŽ{ÙšŽGÒŽýƳ6†¹µ +>¡-îmÌ,±˜Ù,¶s³-A--Ã-«-¬’Z².P†DEC--WRPNmŽÃ?‹"ÑîÀÿ'ÚB4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY02.hdr]œ@>œHñ!ÆhêrÃAØ ?‚ +)?‚ ;?ÿ»Ž öŒ6Hþ?ûêµIíI53vY52ù}µF<Ž¬ä€³AB-™Ê­­hXy­°ÿ“-ØÛ2-Noˆ­2]-KÉŠ.’–DEC--WRPNmŽ +u?‹"ÑîÀÿ'ÛB4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY03.hdr]ŸA>œHò!Æi§ EŽ?~ +?k|?ÿÀ;5 +Ü;ªah?ÿ +ŽŽ\ËŽ¢5ŽˆÈŽ®ï4ný4±Ë-2ìŠ,Ò +V,ÉYW,Ž¬æL¬ÌX-Þ. +$œHó!ÆfÏ9F +ƒV? +Ç}?‚:L?ÿö:•îö»Ð4?þ›¶ +k46×r4á¶çŽ×“84sõz5(õÖ. +þ¬Jzg­ +ה,Ø8r«›ä×­˜Åö-èM-‹h§DEC--WRPNmŽ f?‹"ÑîÀÿ'ÝB4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY05.hdr]ŸC>œHô!ÆdžFj\±?‡5?†Üð?ðŒµ&$ŒßäN?çe¶Š-5}þ‘5Ÿ€0¶ +úŸµ‚îÎ6ˆÁ.3q­j'ô­Á\-úÑi.ŋ®Y㉭_!.ÆDEC--WRPNmŽ)?‹"ÑîÀÿ'ÞB4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY06.hdr]œD>œHõ!ÆÙñÆk›‡?ƒÆ®?ŠŠ¯?éö<Éåž<á®8?鈵žê·µaš&µJ˜?¶ +e µµ]7¶Y¡L-Ê.Ú-hE-×ç»-ëgF- uÔ.G#™-Þ³@.g.&DEC--WRPNmŽå?‹"ÑîÀÿ'âB4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY10.hdr]Œ +E>œHö!ÆÚQÆ#¬j?‚æÁ?€÷ +?€8j±»Xh?ÿ¢µZSx³¿Âëµ%4§¯94yÚ`³ßœ‚-gÝ&,Þ + +-æai¬ïß_«†q*× +Ž¬°h,, + ®DEC--WRPNmŽ“?‹"ÑîÀÿ'ãB4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY11.hdr]Ÿ +F>œH÷!Æù+Å£tä?‚È6?‚–»?ÿÛ; +CÒ< +ñ?ý£µv}IµóÀ³uC§µNSþŽ‡µ3ï-Çà€-%ƒ\«”’Ÿ-¿Éb-«Õ-KsÜ,†&B-Ž wDEC--WRPNmŽƒ?‹"ÑîÀÿ'äB4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY12.hdr]Ÿ +G>œHø!Æ +ÛÃ-°ó?€_>?E +?€9GÆ»x8Ü?ÿ‰4÷ÒB4(ž4X¿Ü3~K3ÙdºŽÀžL­BKR¬”`ä­U¬ds¬Ûy,÷²“¬«Ô°.DDEC--WRPNmŽ +??‹"ÑîÀÿ'åB4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY13.hdr]Ÿ +H>œHù!Æ +r +E›Ïð?€n?±p?ÿÊ»$€¡»³?ÿ +57wK4Òõ4’ygŽ k[Ži&4ÿEè­ J}­&Ր¬þ¿ÿ,¥KÆ-€’­Ÿó,Ø1-Ÿh_DEC--WRPNmŽ$?‹"ÑîÀÿ'æB4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY14.hdr]¿I>œHú!ÆlUF +å?ƒË5?‚&?þm»åû!Œ66ð?û×µŽÕ75G¶s5bù#µ9ĒŽP/ã4ºÁ. [÷­{C­Ë|W-×YÕ-•Y +­¡ä««Ç--FÀDEC--WRPNmŽb?‹"ÑîÀÿ'çB4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY15.hdr]ŸJ>œHû!ÆtLFn-?[¥?ƒaú?þh»á +Ü»Zš?ÿߎœí§3Òc€ŽlirŽåfŒµ86N5”Ì-2‘µ«JgĬã[ +-(†]­VOÄ­;”á­Ò5Ü-Þ]üDEC--WRPNmŽ6?‹"ÑîÀÿ'èB4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY16.hdr]ŸK>œHü!ÆY÷F› +w?ƒK²?‚(?ÃØœ0þ•ŸEð³?{µ«Hx7 EI3Ã\é·Ls}¶,4‡¶â–. +„Û®ódÛ­4eµ/w.¬ +®ž¯^.ˆ®e+¯%ÊDEC--WRPNmŽê?‹"ÑîÀÿ'éB4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY17.hdr]œL>œHý!Ŝ +èÆoÿú?ƒñj?†,?ýžŒ˜ù<) +ß?üžµ„ ¢Ž†ƒµÉ'Œµ;”5(4Õµ÷Bœ-~zw,XZ. +éŒ-mÙ-êjâ-ä•D­g@§. +·ÛDEC--WRPNmŽ²?‹"ÑîÀÿ'ìB4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY20.hdr]œM>œHþ!ŝ„Æ" Š?€üt?‚a‹?þ•»Õ&õ;ß0?þƒŽójŽT±4øïTµKœë4ŠŽ)µ ë-Mʪَ୮®-5€¬ßÞô-Žª5¬µú-&QÅDEC--WRPNmŽ³?‹"ÑîÀÿ'íB4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY21.hdr]ŸN>œHÿ!ŜRyÅ€¿?«?¿?ÿý:‡^:²ûŸ?ÿðŽÈeè²Êú–ŽG»ª³Á†Ž”ŒÒŽžª–-y­5,kó,ÍÌ#)Ÿœé,o·°-ØÜ-1˜ª- +fnDEC--WRPNmŽ b?‹"ÑîÀÿ'îB4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY22.hdr]ŸO>œI!Ŝ:mà +QS?‘1?"?ÿõ:”€cº7/?ÿüŽLÙ³ôòŽ]ÅÍ4›+4@›Ã3z­ ++.'-8V-¯ŸÃ­®D­ÃHݎ>Ÿ­$ð +-úçDEC--WRPNmŽ F?‹"ÑîÀÿ'ïB4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY23.hdr]ŸP>œI!Ŝ% +E!8?—ø?f»?ÿåºë +»+:g?ÿÇŽN+o4ŸÓæ4_Y‘ŽK +m³†ØkŽ +óµ*5 &­vØ)­|P«- +8w­(™f,`1 ­ +ڋ­Ê†QDEC--WRPNmŽ '?‹"ÑîÀÿ'ðB4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY24.hdr]ŸQ>œI!Ŝ‡F +ÔH?¶¥?€Ô5?ÿ +;|vù9ñ§Ô?ÿþŽå ³x +ù3 +XÂ4gŽó4“‚_ŽnN-€`xª9Ž¬F˜ÃªƒN%«Â^Ž-?}Ô-/Bj¬š`XDEC--WRPNmŽþ?‹"ÑîÀÿ'ñB4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY25.hdr]œR>œI!Ŝ-ŠFn˜J?Éä?»Þ?€9I›»ÈKç?þÆŽ5Ôl4óx4æX,ŽŸøŒ2҄Ý4õô«/ø[­ƒ¢­T1s-Ç©æ-J‰+EoG+¿T*,Õ5DEC--WRPNmŽÞ?‹"ÑîÀÿ'òB4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY26.hdr]œS>œI!ŚeFŸr–?‚vc?‚¿F?þL»ëÄ€:¥—F?ÿóŽÿ¢4s€g5T+ 3©lŽ[`5`››-ÒOå¬çW8+<ތ-®¹-Ïr»+°> +«{Ã[-€nvDEC--WRPNmŽØ?‹"ÑîÀÿ'óB4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY27.hdr]ŒT>œI!ÃCâÆr?­õ?‚H¡?òmŒ¥ñä;ŠÂÆ?ÿk²ÔpJ4,;ŸŽÕ#ŠŽì5•;TŽóƃ-]ù,Éoè¬Iԋ­F ^-BÍI-Š +­Àdj-Ý DEC--WRPNmŽó?‹"ÑîÀÿ'öB4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY30.hdr]œ +U>œI!šìeÆ#xu?Ì??€0m?û¥Œ?E:áú–?ÿ掏à-Ž&€Uµ™1³‚*5_P5 š-®>ª€¶æ-Né»-Q%J-\ª'¥P­Ão­‡còDEC--WRPNmŽ?‹"ÑîÀÿ'÷B4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY31.hdr]ŸV>œI!ÁÄ(ËÅ¥S5?€¶?œC?ÿ}»qn¹Nn?€³‘·³4pŽÆpŽ¬D@5S²+Ž·tN-‹ü¬U,M­ =±,*® +-±mà-§öÞ®&,ឯDEC--WRPNmŽ +?‹"ÑîÀÿ'øB4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY32.hdr]œ W>œI!ÁÅ-)ÃV›?0t?:Ä?€¹;|Nº°y?ÿý4H̓3âKq4w\Ù²"R3ÈÍð3W¿0­;UŒ­Ñ¹®õd*6û”­%û/¬–v­š,žiDEC--WRPNmŽW?‹"ÑîÀÿ'ùB4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY33.hdr]Œ!X>œI !Á±Ñ0EœQ?Dþ?3Ñ?ÿý:àk¹æ}?ÿÿ3eš©Ž +zF²Í4uî3Qg†Ž%ì²)Ïá-/3/,‚˜Ð­[ϓ¬¡“y-¶?Ž,ôu£­X‹4DEC--WRPNmŽ”?‹"ÑîÀÿ'úB4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY34.hdr]Œ"Y>œI +!Ášw€F +†Ñ?cz?ŒZ?ÿèºÛ%§»ƒ?ÿÝ31 +³ÎÞ48p&µ&<Ž- +4° J,ä`õ,¿› +-t(E­Zªú-/]K®66z¬Á“,Ä@vDEC--WRPNmŽ ?‹"ÑîÀÿ'ûB4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY35.hdr]œ#Z>œI +!A +“ëFo«ö?‹?€Hñ?üàŒ!ŠPºòkê?ÿãŽz>œ4&¬o4fì“ŽvǵO<Áµ­î¬Èhš¬ +Ó­œ›-,m,RŠž­'Ð;­®9€­\qpDEC--WRPNmŽ {?‹"ÑîÀÿ'üB4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY36.hdr]œ$[>œI +!ÁrtÌFŸ‹Ó?›Æ?‚>š?ÿÔ»žº©ÇÑ?ÿòŽÇƒ3–¹ý4Nàg³ôá³Üut5&5-†Ñ ,¬Ø3­w1k-ˆo+âÝI¬‘vQ«Íp*-€ÈDEC--WRPNmŽb?‹"ÑîÀÿ'ýB4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY37.hdr]œ%\>œI +!Á?@Fž¿W?êÎ?„©¿ûœI!»œOFnÔù?K˜?šß¿ÿð:·'ž¶ +ñ9¿€3¥ÕÃ2S“æ°D5<Ž +&k4]ÑFŽ¢ö·-$wi¬¯±S-5íŒ*6,€Æ䬧ž,ójL­ +,†DEC--WRPNmŽÏ?‹"ÑîÀÿ(B4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY41.hdr]œ)^>œI!ªú¥F +rd? w?í>¿€¹=*ž¶ +N¿€Ž1­²•vº4Ž4¡aR³ÍՔµ%`,®£x¬¯7­¢yº«ÌœÝ-|í -—† ++£¯Î­€v»DEC--WRPNmŽÖ?‹"ÑîÀÿ(B4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY42.hdr]œ*_>œI!°À?EœYä?J@?9¿ÿú:`sºŠ`õ¿ÿ÷2Ó&’3¢^²˜lH³ÜÍL³±03ÉV8¬–Qž­ƒÎ+(€ç«+‹¬Œ­Š‘]­ѺªÎ]DEC--WRPNmŽ•?‹"ÑîÀÿ(B4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY43.hdr]Œ+`>œI!­/<ÃÎ?L?A¡¿ÿþ:³œž³ùm¿€2Ñh1±v'ŽEù²¥yM3¶§Ý3bF,n +§«Li­SZ-¬Wê².ÿû,VO­|Á­F–XDEC--WRPNmŽ{?‹"ÑîÀÿ(B4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY44.hdr]Œ,a>œI!¯")Å€¶Œ?pÀ?‚&û¿ÿâº÷óÿ»Ž~¿ÿÙ4w:3ZáZ4_ÒÞ4œ¹>4³~+5D‘­dú,D.­`§¬Õ_ +¬ðœ‘­qý”­Q§á­Éþ¿DEC--WRPNmŽ˜?‹"ÑîÀÿ(B4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY45.hdr]œ-b>œI!Â㧀Æ#_0?aÞ?€É€¿ÿy;ƒë†;µ¿ÿߎKŸ³œŸŸ4XIŒŽø +!Ž9ŠËŽN‹¬­^£+§N+-Ì'Å,ß²ü­$w-‚ââ,  +Q,Ž0 +DEC--WRPNmŽã?‹"ÑîÀÿ(B4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY46.hdr]œ.c>œI!Âø +cÆoÄ÷?‡›?†àß¿ÿ¶;:ð:EÊ\¿ÿü4j4‘Ý4„ +æŽióc³‘D6$é­CŒH¬-pÇ«_­\¬åè ¬ê‚Ý,ã‹pª«^®PòýDEC--WRPNmŽz?‹"ÑîÀÿ(B4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY47.hdr]Œ/d>œI!Ş€FœÝå?Yƒ?ˆŒû¿ÿö:„OZ<•À^¿ö34"ºCµPŸÉŽ&à5ÖuÜ3¹+ð¶8 +1­›Ÿø,\ñ­2.¬®7à ­;ܚ.Ä +,ŽW®ONdDEC--WRPNmŽ ˆ?‹"ÑîÀÿ( +B4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY50.hdr]œ2e>œI!ŝèHFnbû?‚RÞ?‚r ¿ÿ—»gR_:Ñeª¿ÿë5n'†² +˜EŽŽŸV4Ìæ`Žëе3üù®€Ú¬f +F-­'A¬{R»-Tœ­fŽ­‡§ôDEC--WRPNmŽ +ö?‹"ÑîÀÿ( +B4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY51.hdr]œ3f>œI!Ş7KF +]Ž?ï?‚O¿ÿ‹»tÀ2:ž?·¿ÿï5 ׶Ž:­“Ž„>ñ4%GJµdíµUJÖ­¥.¯,š-žNù­KÕ,(‘Û*ãdu­Åö7® +ä›DEC--WRPNmŽŽ?‹"ÑîÀÿ( +B4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY52.hdr]œ4g>œI!ŞBŸEœþ?¯?o®¿ÿðºŽþÌ;?m +¿ÿ¹Ž‘n·Ž®„ç4f#Û41|ÄŽÙR„Ž*fe-‘¥Ì-xo¬#g +«lQ¡.>Ž-mê +­Ž¡ +¬ûW{DEC--WRPNmŽn?‹"ÑîÀÿ( +B4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY53.hdr]Ÿ5h>œI!ŝѮùº?äÍ?U¿ÿúº[§X»\% +¿ÿ¡5•Ø4֝54”M­3FÓ³;/03@­¥÷í{5­CcJªØ0×­°‘Ö-k‡.2]o­ËšDEC--WRPNmŽ??‹"ÑîÀÿ(B4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY54.hdr]Ÿ6i>œI!ŝî©Å€Øo?¬¬?‚PS¿ÿ°;I’B»6"¿ÿÀ5dË4‹]“³ÎŠ}4òh˜ŽŒµ˜5a×+­Ùõ©¬³ +à-‹®e­ MQ­Vé­}Eª-ÞÄ­õ©¥DEC--WRPNmŽÞ?‹"ÑîÀÿ(B4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY55.hdr]Ÿ7j>œI!ś¹7Æ"Ó?±/?‚i¿÷Œƒ% +»±·¿ÿ +5 š¬4F$œŽ*RË5C +¶5·”W4mH­ˆë†¬w/ý­O–­G¹*-š‹z­Ähá®,&G¬9DEC--WRPNmŽš?‹"ÑîÀÿ(B4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY56.hdr]Œ8k>œI +!şF­Æoú?ëŒ?† +¿ÿ©;L@U»•9¿ÿ]5Fƒq3kÛ03\ +—4ÌÖ 1ÂÅ6 +Rú­¡,Œ+3•*­•\ä¬ñ3f,€|ñ­N£Ò«Æz2®*› +DEC--WRPNmŽ%?‹"ÑîÀÿ(B4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY57.hdr]Œ9l>œI +!ÆŸœFŸj? +M +?‚¿L¿þë»ÀŸ:íê[¿ÿã5Ò⊎ѿžµÕŒ^4'‘†µ8Mµd°*®Gš-„--ô +­ïS®hž,ø@­°ê­¢à4DEC--WRPNmŽb?‹"ÑîÀÿ(B4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY60.hdr]º<m>œI +!Æ +~ëFpâœI!Æž +F +ù;?‚˜š?‚;”¿ýk<Ӛ<#r¿üž5$@ïµPè{µDÕ5$Âr4ytŽòrü­XÍR-Ž­Ÿ-”8U­°©:­—±,Ø +Àª¬t­r8ðDEC--WRPNmŽ‰?‹"ÑîÀÿ(B4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY62.hdr]œ>o>œI !Æ +5ºE›Íì?€•á?Œ¿ÿÙ; + ;:ŠÄ¿ÿŒŽ»¥ŸŽvߎÇë뮀°ñ4bS³3—:i-!$,ŒKo-f!-ÇÌ® +p2­sI|®&.­|H5DEC--WRPNmŽ?‹"ÑîÀÿ(B4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY63.hdr]Ÿ?p>œI!!Æ +(¶Ã*ž?€`Î?k ¿ÿ“»jª;ay +¿ÿžµúfŽQ +Ç4ǔ€3ìï3 +oÅ4›E-€ˆ.,ž@­8v¬‰Ž­ŠI«üë*,à +è­çXrDEC--WRPNmŽ?‹"ÑîÀÿ(B4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY64.hdr]Ÿ@q>œI"!Æ +$Å¥³@?‚Çœ?&Ç¿ÿúºc®u;å¿ÿÖ5kƒ³ÝH4Ä\°ŽªW–²ÏÚ¡4Qß­°Œ)"ÿ¶­ÿã-K—+Ë<+€ +ß+Wk„¬’éÜDEC--WRPNmŽ +6?‹"ÑîÀÿ(B4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY65.hdr]ŸAr>œI#!Æ +ŠþÆ"Óê?~ß?,&¿ÿ»² +ï»;ÝÁ¿ÿ»4“<ô4Îg³qw2çLK5'ª^3jËÄ­I>Ï«Ç¿V-!w®­)¶ù­5Ô,œ +*­RNˆ¬šà¹DEC--WRPNmŽ ?‹"ÑîÀÿ(B4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY66.hdr]ŸBs>œI$!ƐÓÆoy? +oå?†©$¿↌ópõŒœö¿ô05 +®5 !6¶×5ê|ð5•±p5ÿ:ƒ­»ÈR¬àxf­ÀÒ8­µ>G®jð÷®:z"­M£¢®¹±DEC--WRPNmŽÅ?‹"ÑîÀÿ(B4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY67.hdr]»Ct>œI%!ÆcñSFjY?ˆ±3? +_X¿Ý·=Å<ó=^¿áŠ6 ¿ðµ”#5¶-8n5þê +5ƒíâµ­6S® +'í-â›.t‰_®b®óF."ôg-8ÝŽ­°†DEC--WRPNmŽŠ?‹"ÑîÀÿ(B4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY71.hdr]ŒGu>œI&!Æg}ÔF +á? +, ?‚ + ¿úvœI'!ÆhèÄEœí?ƒžŸ?‚>h¿ÿD;œ±ŒºÜ>¿ÿö5“ù”3¢G·Ž¶€5fá4Ëꎔ€'­À[è¬Oäþ-ȗ­å^R­ž]x*£ÃÃ-?[¢­÷ÈÑDEC--WRPNmŽù?‹"ÑîÀÿ(!B4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY73.hdr]œIw>œI(!ÆikÀà +Œ|?‚å?Ž›¿þL;Å@·¿ÿí5cï@4ŽµNS®46ܐ²C)5«Ä­Ÿ4†¬~¶p-²¬€à,³ó5­\€×­U⛮,™þDEC--WRPNmŽð?‹"ÑîÀÿ("B4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY74.hdr]ŸJx>œI)!ÆgKæÅ£Ÿ:? +yÅ?‚É¿ý­Œ +×滺­Z¿þä5è¬4Ø5^nã5^6.4š€š5*¹~® +e¬xã¿­§Š­•»&­›2ý­€à«Ý²Î­\†«DEC--WRPNmŽ +?‹"ÑîÀÿ(#B4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY75.hdr]ŸKy>œI*!ÆcZÆ#hÚ?‰{Ò?Šú¿۔œÊ;pŽÓ¿ÿ6)éw3æG6P +MŽÞ=@5Š'þ57ö£®(ƒ¬Æ×®zOl-P­®mñ?¬ü?Z¬Þ熭}·DEC--WRPNmŽ{?‹"ÑîÀÿ($B4A® A¥š$ÀÀo5745g0526o.356899.cm.982643.smf[XY76.hdr]ŒLz>œI+!@t×'è$ǿڌ蒞ËäÀ„ÛN?îé8”_Ú8”m +>˜pº?tN¿tqÚ>˜¢÷72 +ž5ö7‰¹¡¶—搞O¢é·£Ín¹ïÐ':à‡Kºí1Qº +æcºÓ°:ÏýDºàœ§º +),DEC--DISN[9œd¡ÀÀo5748g0436o.358811.cm.982690.smf[PHU]{Æg’9Æ3–?†Ñ(? +]»?õ~Œ”Iþ<Îcº?ꟶç뵁‹µ +í,µñ §5ϞXµ³}†.9m-;œ)-à`è.7KÚ«†"?-±àv®ÿ€-ô#IDEC--WRPN[9ß?Œ(öîÀÿ'ÙB4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY01.hdr]œ|{œSô!ÆgŒbţЌ?„i€?ƒ +K?þ8»óýÎ;ô*„?þ%µÍçŽÝV4ç-µb’V5Cöµýð.^æ,­Äð,®Ñ-À9ž­ŒÞ-Z( +®íê. +µbDEC--WRPN[9@?Œ(öîÀÿ'ÚB4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY02.hdr]Ÿ}{œSõ!Æj!à +Ε??€Ž!?ýŒŒK(ºë5?ÿåŽ{*3 O5ˆÉý4öU:3 +Ã×4`Ð!,çÖ6«žPíÿè%­$ÅK,–­uRç¬óº-)JŸDEC--WRPN[9 ¹?Œ(öîÀÿ'ÛB4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY03.hdr]¿~{œSö!Æiå‹Eœqu?=Á?Òœ?þ=;ï.ë:?ÿýŽOöɲš‘µ\F=ŽÈ™VŽ"Ÿ§5Œ’,îG-ª—-ÌÌú,®Ÿa-`e®,?á+Ú\Š¬Ú0DEC--WRPN[9î?Œ(öîÀÿ'ÜB4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY04.hdr]¿{œS÷!Æf ¬F +o?†“Z?–“?ýŒ!X'Œ:nè?ûm¶ +ki4õ'\5š‚Žºù׎¡ˆÏ3`$-."¢¬¬ò‚® +óÿ-—Üb-~ˆ-{W¯«Œd¬-u‹EDEC--WRPN[9 Š?Œ(öîÀÿ'ÝB4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY05.hdr]Ÿ€{œSø!ÆevFiŸm?†”{? +Ž|?õŒ–Âœ+õ?ÅU¶v¯5ûž¡5Ÿiô¶ +ωµp?Z5²Êû.3]® Xÿ­ž˜Y.=E.5­®Nâa­4t +-šsÃDEC--WRPN[9?Œ(öîÀÿ'ÞB4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY06.hdr]»{œSù!ƺÚÆpyÁ?ƒ¹?„—3?óŒ<ž€göŽ +»µ1]ü-@“æ-ÏÒ-<ôˆ¬Qà2-ŸnP+²„Ù+ +hy-ŽˆDEC--WRPN[9Õ?Œ(öîÀÿ'ãB4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY11.hdr]Ÿ +ƒ{œSû!ÆDïÅ€»?‚”¶?•G?ÿÛ; +¡2;‰œÿ?ÿjµo!Žˆ±3qí厏ߎž +ѳˆt-×#ê,˜ +X¬æñ-0±-'Ö¬ugV-㒪fmDEC--WRPN[9ž?Œ(öîÀÿ'äB4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY12.hdr]¿ +„{œSü!ÆöÊÃ&#?€é°?”X?ÿ€;X,iº_Eù?ÿú4Ýø²25:ŽlÁiŽ”ûà³'ÏŽòŸ¬qÕo«¥¹{-(¿,ʗè+Á›t-ÍtÅ,e÷O-©¹DEC--WRPN[9q?Œ(öîÀÿ'åB4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY13.hdr]¿ + +{œSý!ÆøŸEœïÞ?rá?Æ?ÿ³;FØQ;‹Ž?ÿhŽ€ŸÔŽ€ÎÀŽÆÙ4?Ôž2õM=Ž(ÀW-¡µ-$ø-V5Ȭª³­!æ†,á9ª­_Ó­Š{ DEC--WRPN[9Ô?Œ(öîÀÿ'æB4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY14.hdr]¿†{œSþ!ÆQ +F + +ù?‚,?žL?ÿs»‡$‰º÷?ÿâŽä*u4-Þñ4Ò‹ÙŽÑš(Ž$43Óñ,ðo§¬Ý!­^³«Î¢T,‘±T­Û>C¬Qœ“­_DEC--WRPN[9ì?Œ(öîÀÿ'çB4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY15.hdr]Ÿ‡{œSÿ!ÆaFoT‹?‚Ù?WÈ?ÿœ»;¢Ç:˜³æ?ÿôµƒD©4Œ³üc93ÛÍŽ¡V4nAN-Èg>­ +ʜ­IyÚ­A’­€ýÆ­BÿR­mc.,d°ÚDEC--WRPN[9 "?Œ(öîÀÿ'èB4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY16.hdr]Ÿˆ{œT!Æ™ÖF˜Ë3?„°^?† +.?ízŒÀÉôŸiG?ygÁµ£046Þ°€5ÒvÛ·—.µƒdþ¶«oB-‘Œì®iíq®(õ¿/™0-ܮ߯Â5Å­„AO¯?GVDEC--WRPN[9ñ?Œ(öîÀÿ'éB4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY17.hdr]œ‰{œT!řbBÆoå?‚* ?†ä?÷q<|‡;*Êû?ÿËŽ®ÒâŽ^]wµ#¿Žw¹mµTŸH¶"R©,c£,æÀ-düÈ- +þ,-fh_,âaÕ-oÖ7.ImMDEC--WRPN[9?Œ(öîÀÿ'ìB4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY20.hdr]œŠ{œT!řIéÆ!àV?ƒA?ƒ€?õ»<þ”;»n5?þîµ;_«Žw2¬µ®Y~µyšµƒ +µ„¡t-¹æ+Æ¢M-¹‚f-€k×. +è-H·Â-•±z-ÈÈ@DEC--WRPN[9?Œ(öîÀÿ'íB4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY21.hdr]Ÿ‹{œT!ŜÞ8Å€øe?U +?R!?ÿþ9ô"Ù;•è?ÿÜŽMŽù³h™ö3ÃÅø.¿.2ŽŽÔ³ËBu-[&-ªµí­žY@,×ìò-€Š¬ - +ùj,cÒzDEC--WRPN[9§?Œ(öîÀÿ'îB4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY22.hdr]ŸŒ{œT!Ŝ\LÃdy?Òê?n‘?ÿâ:÷ {;'/ +?ÿÉŽéú™ŽÍã,Ž§ +ÛŽî4FBÔŽPú-,ýI-x‡á-Ô#Ò,•£­¬Ö,Ÿž¬BX,VÎtDEC--WRPN[9Æ?Œ(öîÀÿ'ïB4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY23.hdr]¿{œT!Ŝj1E3>?Í?`a?ÿìºÌš':šu?ÿýŽVŒ2I4ãšU³–èa³'d+4ZÅž«–3à«z)x­ï !¬€¶—¬úÓ­8 +£­Sz +-‹MDEC--WRPN[9 ^?Œ(öîÀÿ'ðB4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY24.hdr]¿Ž{œT!ŜÊ[F +Z?È ?ŽJ?ÿ@;œç/º¿]Ž?ÿîµ%ƒ,1ƒoƒ³ 2ŽÐI4¥J +4ܝb-úG¡,LG>-EL-5ÕK-î«Âm”-eh-75DEC--WRPN[9 +$?Œ(öîÀÿ'ñB4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY25.hdr]Ÿ{œT!ŜÔFn?‚Ž?‚µX?ÿD;š +û»  +?ÿ:Ž¥24€k³5 +=Ž¢@”4ŸÒf5T¿,Yaõ­ns݊õY«‰¡J- +­)­oH0-(–j-ÆeDEC--WRPN[9Þ?Œ(öîÀÿ'òB4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY26.hdr]Ÿ{œT!śE©Fž)>?‚iý? +~º?ÿY»ŽÊûŒ?O_?ûŒµ;ûm5¬’5 +ÕµŽšaŽM¹05â -{á5¬ +–$­¶-.+Y,󬜭·bÚ¬‡.ùDEC--WRPN[9?Œ(öîÀÿ'óB4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY27.hdr]Œ‘{œT !¹X_Æpû?‚,j?„õ†?ýü»ûµ@<×?þ +Ž}=Ž’Û +µE`!µ0Åô4ÑfµÝ^­-ù±ñ*•N.+år +- +^Š-žr=-™ó<­Åð. +þâDEC--WRPN[9œ?Œ(öîÀÿ'öB4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY30.hdr]œ +’{œT +!ŸÎÆ#ä?€6©?+;?þ»ýyy:ÙP(?ÿé4i³V%}5—8³JV«4ø°Ê³—3«¡¥Ùª‰«­ºŽ,ÑÀ®2ÓÏ,˜ +­4æ+®ÿDEC--WRPN[9 ?Œ(öîÀÿ'÷B4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY31.hdr]Ÿ“{œT +!ÀišáÅ¥žp?{Ò?_³?ý$<$P8|v¡?€Ž³–¶³Õa=³ºŠŽª˜¡µ}~/²º¬à-×Cu«í +-2Má-.òŸ­Žø-”Ú .)óˆ¬ÃìúDEC--WRPN[9î?Œ(öîÀÿ'øB4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY32.hdr]œ ”{œT +!‡èà +" +?PŠ?;a?ÿô:ÿƹ’bO?ÿÿ³éVº3‹ä³\r^3í3šŽjÃ4ie‚-_ 9¬ÜÝ$­uÇ­;_Í-ç+®¬ò-Ϛ£­ÕªØDEC--WRPN[9??Œ(öîÀÿ'ùB4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY33.hdr]Œ!•{œT +!Áþs³Eœ&Æ?D?'Ü?ÿã:ó_¹ÅQ?ÿÿ4©Õ³ÃÀ⎊4Šz¯4L?ÿŽæg­šªŠ-cŸ-.ÿ•­BU€­ìá¡.'xÿ-5:Ä+s{’DEC--WRPN[9ú?Œ(öîÀÿ'úB4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY34.hdr]Œ"–{œT!ÁÅîF +GŽ?WÄ?‚P’?ÿž»>þ­ºiûÙ?ÿù³6C +34¢Ï3Œ—ŽŸ/ŽÎV45€] , +le«gï,',Ô+8Z0­ž€l­ˆÒu. +×DEC--WRPN[9 ;?Œ(öîÀÿ'ûB4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY35.hdr]Ÿ#—{œT!ÂZÚeFnÔG?œü?Í°?ÿ:;Ÿ +·; +í@?ÿϳø÷c4§åŠ5Ž$5‚4ä +4ñî-Š­±¬ü +^,(+¬-³íŽ-¹C‹-Ÿâ™-k ‚-c(ÂDEC--WRPN[9[?Œ(öîÀÿ'üB4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY36.hdr]Ÿ$˜{œT!ÃèZFŸ€p?Š¢?ªÔ?ô†<™Iv;.ӎ?ÿij“¢³ýò±4Âno4}ƒ5}µ 4ª¡-èa«ûŒà,„ +¬ÆÌ-;×k,“ãî-‘oà- +àDEC--WRPN[9M?Œ(öîÀÿ'ýB4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY37.hdr]Ÿ%™{œT!BKÜÉFž21?Fæ? +Àó¿òã<ža†:ÙÄ¿ÿê4›š=ŽP^ ²o4ö4ŸM5‰03¶^­²Úp*ðo,ã·ý­M{À«‹,-ŸIÓ®!>ŸDEC--WRPN[9Ç?Œ(öîÀÿ(B4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY40.hdr]Œ(š{œT!ÂÊÂÇFoI€??€Ûä¿ÿh»Œ-ðºlÃW¿ÿù4Ÿ”³¶€Ž5¡ÖŽƒš¿ŽŠ{ý4 +€¬þ­£¬ÈÓ-—4T¬ìè­+Ì,^­0I$¬äß¹,WHDEC--WRPN[9e?Œ(öîÀÿ(B4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY41.hdr]œ)›{œT!’DúF +¢?u{?—`¿ÿ÷º„ë¹fÏv¿€4:ÍùŽ ŽÄŽhÌô4ŒxNŽš +Ž—¢œ­Ëö×-˜Us­ +Å*œ»S­DŽ§-e¿ +­ƒŒ$­ +DEC--WRPN[9d?Œ(öîÀÿ(B4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY42.hdr]œ*œ{œT!€£Eœ%?M™?;_¿ÿÿ9˟_º_ ÷¿ÿú3rQ(³2[Í3M-z²Äj»³I‚ +43{ô­#G+«@þT¬¶]p­:êh-lŽu­±›Ö¬3 Í-ZyDEC--WRPN[9 á?Œ(öîÀÿ(B4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY43.hdr]œ+{œT!™â)à žø?K?3‡¿ÿþºAÛžX°›¿€2.Æ +ŽGµ2¿ß›ŽPÉ2Ö¬3Q䬙¿c-Òí +-Rœ•,{û3­,ˆF-|h-KխϬDEC--WRPN[9ú?Œ(öîÀÿ(B4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY44.hdr]Œ,ž{œT!²&NÅ€f‘?yÆ#¹¬?€šR?îv¿úŠ³Æp¡8?ûD? +W¿þ•»Ò2»AŽ¿ÿº1^‚Š4ñû58 +¡4›b4Á.g5ô`‰+ۏ¬ßÀ +«±pù­qx§­Šß–¬@ã¬Îüs®Ž}DEC--WRPN[9?Œ(öîÀÿ(B4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY47.hdr]Œ/¡{œT!Ŝê F¶Ÿ?‚n?†ú¿ÿöº‹ÙQ<(–‘¿üÀ5ÕÀŽ€ÿ¬µZÎ35™º˜Žl—í¶ž­&o,Ç$Ë-à­l +/­š-úxˬÌs]®1âÊDEC--WRPN[9 ?Œ(öîÀÿ( +B4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY50.hdr]Œ2¢{œT!śÜFnŒ?Äl-E€’DEC--WRPN[9 +Ã?Œ(öîÀÿ(B4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY54.hdr]¿6Š{œT +!ŝÚ=Å¥5?Fc?íj¿ÿÝ;»|8áv`¿€3oô4ÜI±‹Õ³V +ÜŽŽÄ$5S +,Ž©¬ð_F­A%,–Uâ-À©¬ª&ð-`’‡®ÏrDEC--WRPN[9 ”?Œ(öîÀÿ(B4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY55.hdr]Ÿ7§{œT!ŜpæÆ!»ê?íã?ƒi¿¿þ_»ä‚à»·Pî¿þÿ4ˆÀM3çˋ5|15=£^4ÁüE5¥p­ž-+ŒU¬›Åú­Ôä +­¡ +­kÈq¬Ï©5­ÏèzDEC--WRPN[9Ž¶Q5ŠŽ_,è +I­ŸŸ­5"µ*‰;­5L7­Œ'¿,ÐZõ­«"DEC--WRPN[9›?Œ(öîÀÿ(B4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY57.hdr]Œ9©{œT!!Æ/8Fh? +?ˆ +Ï¿çn<Ûû“<ÃD;¿î5ºUµ3 ¶1ç5÷€o5€Ñš¶&¬ô­‰ì-*øá.^è2®$'­Ô +.4-„Å®<ÈDEC--WRPN[9P?Œ(öîÀÿ(B4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY60.hdr]º<ª{œT"!Æ» +Fo„Ö?‚l?Q¿ÿ±;Jï%» +™#¿ÿØ5$Ûx³š7ŽÐmŽ³Áb²z­ ŽhA›­šÞÁ-1¬#(b-‡É÷­ªê-o +M¬äT¬xQ|DEC--WRPN[9 +?Œ(öîÀÿ(B4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY61.hdr]œ=«{œT#!ÆõÔF +i%?ƒ”¯?÷À¿ý‡<g;·‹¿þô5}¶Ž÷¡»µŽ›t5%a3ÅicŽŠ™x­·ís-÷*-ôˆy­ŠŠ|­Äþ-5šÚ¬ûKý¬šcÿDEC--WRPN[9 8?Œ(öîÀÿ(B4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY62.hdr]¿>¬{œT$!Æ +àLEœL·?ì?LP¿ÿœ»8oTžÔË⿀2與3Œñ +4•2D³—æ +3ž¥.ŽU§û¬Uß«!„™­G-ÜÛ­ºÐÌ-}ß­õ„Î,Õj×DEC--WRPN[9 +5?Œ(öîÀÿ(B4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY63.hdr]¿?­{œT%!Æ +£6Ã+R?:ó?„¿ÿž»_;‹©Û¿ÿh3’qŽ‘U4/4Á}ó4Ý&’4G†v¬)™¡-ÂB¬Šå­r +­‰?ثόú­°$œ­Ü‰çDEC--WRPN[9 o?Œ(öîÀÿ(B4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY64.hdr]¿@®{œT&!ÆÚ„Å€þ?ƒq?‚z¯¿ÿèºß#‚»Sڟ¿ÿš5ˆôS3˜/ï5ët4×[bŽЈ5|Xü­åìH* +t¬z™D¬Ý +<®/­‘L--»Ð4®’iDEC--WRPN[9 + +?Œ(öîÀÿ(B4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY65.hdr]¿A¯{œT'!Æ +tÝÆ#8ñ?º›?}¿üN<.T–;üŒ¿þ5¹Öµ[ÓŽL +Ž«ë͵24Ã<Ø­‹U-Ri"«Òq1,Ñ).,×f,[õX-:Üe­7ÄqDEC--WRPN[9(?Œ(öîÀÿ(B4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY66.hdr]œB°{œT(!Æ!K°Æoò?|ÐZ?†ùò¿Þ<öŒEŒd˜©¿úhµ1J$4ýbµ²To5žü±µ‘OS6£^¬m쫯.' +å®Ïé,â*Э£ø|-“¯É®2á.DEC--WRPN[9Š?Œ(öîÀÿ(B4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY67.hdr]¹C±{œT)!ÆdÜFk£8?‡‹? +¥»¿ã‹<ôð<vü¿÷“6>iµ€d¶ +O·5¹Œ5†?–µõ)•­ü!Ô,¡7.Rs®ìT®­§-±?‹-]\®4m^DEC--WRPN[9Ô?Œ(öîÀÿ(B4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY71.hdr]ºG²{œT*!ÆhŽ«F +šË?ƒ ð?‚pe¿ÿ¬;Pžó;“ +Í¿ÿT5 +8ƳÅ#µ( x5%™T3œº]µ of­“ª2+Ñš-ËJ;­?Æ,è®-ž²-5‰Ì­t²çDEC--WRPN[9 +R?Œ(öîÀÿ( B4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY72.hdr]¿H³{œT+!ÆjÞEžÅÊ?€Â×?Øe¿ÿúºa0 +Œ!/;¿üâ³/7Œ5––3°Û5 +Ì厃ŽŽeH«ü +­#**åà­‚mÃ-8}š,å' +-!¿ÿ­žlRDEC--WRPN[9 æ?Œ(öîÀÿ(!B4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY73.hdr]¿IŽ{œT,!Æj=ÌÚ©?ƒà?p0¿ÿøºë9‡w÷¿ÿÿ4“ö]²œªq3¶Ò†4_ñg4\³ø3‹™­A~+“c–¬­M­a©­™‘î,‚†Ÿ,ù +†­X|íDEC--WRPN[9 ?Œ(öîÀÿ("B4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY74.hdr]ÀJµ{œT-!Æi<'Å¥o,?ƒ4ý?ˆ2¿þb;é0O:8x\¿ÿü5„4B²$ +³p³:4<7\µi]¿4‹%¥­£5 +y­{°|¬Ò. +K8¬Ž¬h-a2G­2=DEC--WRPN[9 %?Œ(öîÀÿ(#B4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY75.hdr]ŸK¶{œT.!Æb"MÆ!è;?‹Of?ƒ ¿á2œhÝ»‘?¹¿ÿF6[ÿÛ4Z26@cq4ÚGx5]å5†×®W¬Ïµ‹®•AÝ+{­ðWœ­ËÁ«­K'–­š»ñDEC--WRPN[9q?Œ(öîÀÿ($B4A«˜uA¥š$ÀÀo5748g0436o.358811.cm.982690.smf[XY76.hdr]ŒL·{œT/!@t×#õGŽ>¿Ú‘þúØìÀS<²>µö8”]8”d +Ÿ`i€?y» +¿yÍvŸ`Æ!6W·ÏžtÚ»·¡î§·Ç|Už9ÏϵŽfb9¡ $:â·ÛºóÝ9ãêŒ9Ú>Ý:Ð5àºæ—Š9ÉDDEC--DISN]¯œq¡rÀÀo5748g0447o.358840.cm.982715.smf[PHU]žÆeÙ +Æ +j? +ÝÂ? +~O?ú;5|É$Žþڈ41mø5‡œn-¿¥Ï-4&í® +žA-H9:«®lÉ­8 +Î,©±<.%ºtDEC--WRPN]¯ z?‹ +îÀÿ'ÝB4A¬ÿÖA¥š$ÀÀo5748g0447o.358840.cm.982715.smf[XY05.hdr]ŸœžœYÔ!ÆfËFkL?ƒÝ +?†€Ž?¹CCŒ¡ +.?óѵ„gX5 +5•û¶v9µ£¬þ5ð*+-¡vó¬S—[­ª3r.)7-ë_®>Êý­ÁÓ(.)ÜDEC--WRPN]¯z?‹ +îÀÿ'ÞB4A¬ÿÖA¥š$ÀÀo5748g0447o.358840.cm.982715.smf[XY06.hdr]¹ŸžœYÕ!Æ +Æi¶†?„Êc?Jò?@ÿ=’Ø<úÍ7?äÞµ\[µ”Ùϵùí߶ —i¶h•Ž¶‹9-nM-R!Š-Ï).;Ã<.3Ö.ù.nØÁ.™‚êDEC--WRPN]¯˜?‹ +îÀÿ'âB4A¬ÿÖA¥š$ÀÀo5748g0447o.358840.cm.982715.smf[XY10.hdr]œ +¿žœYÖ!Æ«¢Æ!å3?‚9?ƒ„÷?þ©;ÏñS;?M?ÿºµZÌŽ’³ŽôšŽ`æŽøTµµ¬ä +-°BA-#áx,úx¬”ßa+•.{-ýÐ-ULu-ó|ìDEC--WRPN]¯Í?‹ +îÀÿ'ãB4A¬ÿÖA¥š$ÀÀo5748g0447o.358840.cm.982715.smf[XY11.hdr]Ÿ +ÀžœY×!Æ0•Å€R?‚XÁ?åF?ÿÇ;+9˜;† +?ÿ=µ4/ÚŽà*“Žn=2&·ŽŒ®µ-mœ- m-%Ž·¬T–,ç~s-_ O­<·,²ñù.ÔýDEC--WRPN]¯è?‹ +îÀÿ'äB4A¬ÿÖA¥š$ÀÀo5748g0447o.358840.cm.982715.smf[XY12.hdr]Ÿ +ÁžœYØ!Æ + +ÚŽ?€f\?Ž1?ÿó:€[=º]ãª?ÿú4öcõ³¶Ùk3Ô(Ž°÷M³†3%µ@ñ­LYa,Z'’«%?¿,®;®)H +ö-Ž-J,’ŒŠ.ËÌDEC--WRPN]¯ 8?‹ +îÀÿ'åB4A¬ÿÖA¥š$ÀÀo5748g0447o.358840.cm.982715.smf[XY13.hdr]¿ +žœYÙ!ÆøçEœ¿d? ,?W‘?ÿä:îY«: +†?ÿý²Â q³ +ÅKŽ³Ìųo3%4‰4eãü, ìñ+iÌ6-h#Y,zÁ­+²+ÅT-,gè-×ÍæDEC--WRPN]¯a?‹ +îÀÿ'æB4A¬ÿÖA¥š$ÀÀo5748g0447o.358840.cm.982715.smf[XY14.hdr]¿ÞœYÚ!Æ5öF +!ã?‚N£?‚Õ +?ÿµ»B†»w2¡?ÿŠµ\ê4CÛ5£¥µ +K3¹•ä5Ü,òw଀g­·u§,~%k¬ +{® `†,ç²-ãZDEC--WRPN]¯Ð?‹ +îÀÿ'çB4A¬ÿÖA¥š$ÀÀo5748g0447o.358840.cm.982715.smf[XY15.hdr]ŸÄžœYÛ!ÆŸâFnÖì?[x?‚8b?ÿ}»€ÃšºšWÓ?ÿõŽ~ñ/4 +žÓ³‚f +Ž(É +ŽÏšô5œ-RÚ$ªšÿÊ-Mð·-kß-Vɹ,Š¡Ë­ +)-ÐŸDEC--WRPN]¯ "?‹ +îÀÿ'èB4A¬ÿÖA¥š$ÀÀo5748g0447o.358840.cm.982715.smf[XY16.hdr]œÅžœYÜ!ÆÊ=FŠfœ?|. +?aàð?ö)Œž}³<‘¯@?ó4Ïéµ{„{¶/–é5Îœ¶ÀÒ· +Éô«¥—.Ä +-£ ù/pD®˜C/ +F„®Ž +߯ +DEC--WRPN]¯ +?‹ +îÀÿ'éB4A¬ÿÖA¥š$ÀÀo5748g0447o.358840.cm.982715.smf[XY17.hdr]»ÆžœYÝ!řìHÆoJ +?€©¥?‡^e?ö]< +‡[;ŸV¥?ÿM±›ìz3}!4øuÿµ:㎵–ó*¶$>;,¢ +ä,* +묞Ch¬Ñ¿­1hx-Ãȕ-ϒ‘.>üDEC--WRPN]¯#?‹ +îÀÿ'ìB4A¬ÿÖA¥š$ÀÀo5748g0447o.358840.cm.982715.smf[XY20.hdr]œÇžœYÞ!Ŝ¡ãÆ!°V?Åù?ƒ*9?ÿ»0h9;Ç +G?þѵI𙎻už4oòÂŽÿ’Œ1Ä*]µƒvo. ê,þúM¬.T- +±š­ÈY-9Ÿ,0ûé-ŸÇUHµBê«5Ÿ³},ÂpÖ,úÌ­Œ4-t¢b-u +¶­™+µ­[0†-Ë +mDEC--WRPN]¯X?‹ +îÀÿ'óB4A¬ÿÖA¥š$ÀÀo5748g0447o.358840.cm.982715.smf[XY27.hdr]ŒΞœYå!ÂÆÿÆrØ?âf?ƒEÉ?üŒ16»—ÉÃ?ÿP4îËQ3{Q(µ5e‚5#zW5šmµ†#M­‡4Ž­d­‘H +,‡[ +-¯¯@­Ýí= +2-·o{DEC--WRPN]¯ ˜?‹ +îÀÿ'öB4A¬ÿÖA¥š$ÀÀo5748g0447o.358840.cm.982715.smf[XY30.hdr]Ÿ +ÏžœYæ!™õPÆ"Ÿ“?‚n?™g?ûÆŒ:Í;ƒà?ÿÒµT3µ~àµBj*²¹Ž/5RͧŽŽT-Ÿ,šÅ.².è-›ó‰¬Úѳ­®Yâ,ü‚SDEC--WRPN]¯·?‹ +îÀÿ'÷B4A¬ÿÖA¥š$ÀÀo5748g0447o.358840.cm.982715.smf[XY31.hdr]ŸОœYç!>ӕùÅ¥Ê?lA?‚, +?þœ;Ê?•º™’?ÿõ³ +4(ÿŽ˜Å€ŽOÎ1µ0Nµ^)Ž,åçT¬ñ۬ы”­šÖ +-©Å-‡ÖV-ˆ¢.{iDEC--WRPN]¯œ?‹ +îÀÿ'øB4A¬ÿÖA¥š$ÀÀo5748g0447o.358840.cm.982715.smf[XY32.hdr]Ÿ ÑžœYè!ÁÁó +Úp?;1?NÉ?ÿú:\Š&º!?ÿþ3ŽôÅ3¡Ž%4$m?³ÀsŽ™Ž×±úP,‰ƒ­p»õ­ú}"-ˆ«ÕÐw¬)Óv-鍎+ÞõDEC--WRPN]¯?‹ +îÀÿ'ùB4A¬ÿÖA¥š$ÀÀo5748g0447o.358840.cm.982715.smf[XY33.hdr]Œ!ÒžœYé!Á·‹EœY?/þ?<†?ÿî:œIIº [?ÿý4j¯ŸŽ Žµ +õ3Àôê4 Ս2›rÁ­?0-›¹9.0eªÌœb­ÿí«-Ž‡ô¬L–-•0ÖDEC--WRPN]¯ À?‹ +îÀÿ'úB4A¬ÿÖA¥š$ÀÀo5748g0447o.358840.cm.982715.smf[XY34.hdr]œ"ÓžœYê!Á(©%F +gµ?„v?‚@Á?ÿ»ªw®žþ˜v?€³€“G²<çh4Æ@K³*A¹µ¹ +5ip +,ŸL.,€ßk,¬‚ ¬ 'Ð-¶C«åM5­É»2.ïDEC--WRPN]¯ Æ?‹ +îÀÿ'ûB4A¬ÿÖA¥š$ÀÀo5748g0447o.358840.cm.982715.smf[XY35.hdr]œ#ÔžœYë!ÂÎrFp.q?€d:¢Ïs¿ÿó³÷ê³à:ÖµÀ¥ŽÙ[54Ÿí#µƒî +,þR/¬¡(-Ã×-s°Ù-¹_-}à­(Ž‰-ÖDEC--WRPN]¯ ¶?‹ +îÀÿ(B4A¬ÿÖA¥š$ÀÀo5748g0447o.358840.cm.982715.smf[XY46.hdr]œ.ÝžœYô!ÂǍtÆpÇ?Ž€?†8p¿€¹;8c:•ªÁ¿ÿö4~ß4§0F4i⎈]ú3©ùP6 +Õ­rÝÿ¬åïí¬áH­ +p笙-ž”«Ëú®1¬øDEC--WRPN]¯ +?‹ +îÀÿ(B4A¬ÿÖA¥š$ÀÀo5748g0447o.358840.cm.982715.smf[XY47.hdr]Œ/ÞžœYõ!ŝ||Fž”€?‚AQ? +ò¿€7ãa€»#‚Ô¿ÿÎ4æ247uµ +c3©i³ÑuZµéx*­À~P­wôu¬—+­€šn+ƒ¬:›®-™DEC--WRPN]¯=?‹ +îÀÿ( +B4A¬ÿÖA¥š$ÀÀo5748g0447o.358840.cm.982715.smf[XY50.hdr]œ2ßžœYö!ŝ«YFn5±?„>?‚šK¿ÿÒ;Ú;†®S¿ÿu³ŽÌžŽz1Žë-4ý·2)ŵG6-^--XŸ-‚[¿¬ÙŸå­^-„üª«Ðap­•EDEC--WRPN]¯ Û?‹ +îÀÿ( +B4A¬ÿÖA¥š$ÀÀo5748g0447o.358840.cm.982715.smf[XY51.hdr]œ3àžœY÷!ŝ»òF +pM?jL?ÞŸ¿þÁ;É4:›Æ"¿ÿô2†é +¯Íª\ŽÀ>Í4ùÞ:5#GŽón<,µˆ«!–,|G¬Z·ø­¶ÕÌ.G-›¡Î­1ž[DEC--WRPN]¯À?‹ +îÀÿ( +B4A¬ÿÖA¥š$ÀÀo5748g0447o.358840.cm.982715.smf[XY52.hdr]œ4ážœYø!ŝíkEœ-V?Œp?rq¿ÿà:þ;^: +–„¿ÿý3cà²Ê·;µXD3ÏNݲ®À(Ž›mÖ-Ù+t œ.!»?+±.í-r +-y_ôªÆ]­2DEC--WRPN]¯»?‹ +îÀÿ( +B4A¬ÿÖA¥š$ÀÀo5748g0447o.358840.cm.982715.smf[XY53.hdr]Ÿ5âžœYù!ŝÇmà +!f?ù¯?cµ¿€¹^lbºØοÿé5]4Y‘–4vÖ636¯3ŽÑä4wº¬­˜üÇ­ +‹ +­š +(,zS¡-á7­1‘*È„­vÀÎDEC--WRPN]¯ ÷?‹ +îÀÿ(B4A¬ÿÖA¥š$ÀÀo5748g0447o.358840.cm.982715.smf[XY54.hdr]¿6ãžœYú!ŞÛÅ€ïo? +G?‚RZ¿ÿº;;‰» +¿ÿÙ4êàÇ4žŽ¢eê4ŒœŽNø²5'"­W•­> +v­j°ž­*>ð.$¬þF«'º¡® +FRDEC--WRPN]¯ Ò?‹ +îÀÿ(B4A¬ÿÖA¥š$ÀÀo5748g0447o.358840.cm.982715.smf[XY55.hdr]¿7äžœYû!ŝ +YÆ!QÚ?‚ +Î?„ÖÇ¿ÿ‡»sÀë»!ቿÿÏ4*Ü4s®ü5kó4wpÚ3ÊnÞ6 Õ\­°@w­ +Tý-šMU¬Õ'®>XݬÑx,šÅ¿®_Y”DEC--WRPN]¯ %?‹ +îÀÿ(B4A¬ÿÖA¥š$ÀÀo5748g0447o.358840.cm.982715.smf[XY56.hdr]Œ8åžœYü!ŝ +ÆoìŸ?‚;÷? +‰ö¿þ[»âpŒuË +¿øü5,Žo5 +¥µ4¿ç>5•‹4ör5ߔT­n®P­#¶q­€ýé­«u–¬F$ÿ­ß-“­Öã®>ýDEC--WRPN]¯ o?‹ +îÀÿ(B4A¬ÿÖA¥š$ÀÀo5748g0447o.358840.cm.982715.smf[XY57.hdr]œ9æžœYý!ÆŒµFÇ(?ƒž©?†{â¿ës<ÈÆç<‹žÃ¿öâ5taèµ[ݵ”‡€5¯t­5”³ +¶;&­ +.Š-/Âf-îÒq­ý÷­›(à-áYI-¯žx®òDEC--WRPN]¯È?‹ +îÀÿ(B4A¬ÿÖA¥š$ÀÀo5748g0447o.358840.cm.982715.smf[XY60.hdr]Œ<çžœYþ!ÆúéFnBr?‚€0?ƒ›¿ÿö:N;Uð’¿ÿ§5B<\Ž¥m^µlš4®;ªŽ†éžœZ!Æ +æŒEœ]I?~}ô?–¿ÿû:NðZº·¥¡¿ÿðµ ÞÝ48Ž(Žª wŽŠ‰"²Œže²ìI. +PˬxÂs-w?Ó-qÙ- +«,”9±,©s,,ŸæDEC--WRPN]¯ Ž?‹ +îÀÿ(B4A¬ÿÖA¥š$ÀÀo5748g0447o.358840.cm.982715.smf[XY63.hdr]Ÿ?êžœZ!Æ +úœÃ#*ì?€ï©?’¿ÿÿ¹›> +:À¡è¿ÿî³Ç‹³'0ÉŽTlG4ĹÆ4«¬(3°K÷+ÈÄn+‘X-F6­‹Çœ­{”-=œy¬3™•­Æ"ºDEC--WRPN]¯ d?‹ +îÀÿ(B4A¬ÿÖA¥š$ÀÀo5748g0447o.358840.cm.982715.smf[XY64.hdr]¿@ëžœZ!Æ +΀ť‰N?å?_m¿ÿß;w;”õ¿ÿà4vØü³lÏ+³W„üŽöcýŽ@f4ò¿Ï¬ÓÄF«Èê+žTõ-e»l,|'-$ +ï-9^á­¶‘ DEC--WRPN]¯ +M?‹ +îÀÿ(B4A¬ÿÖA¥š$ÀÀo5748g0447o.358840.cm.982715.smf[XY65.hdr]¿AìžœZ!Æ +nÆ"?ƒõ?‚¹I¿ÿ¥;XzŸ»F +t¿ÿ³5£óW4{5z2¬Ï£³«™m3ºk5‹e­÷ÂĬåº+­utC¬`Û-¿+ˆ’|¬Â +Ë­öbâDEC--WRPN]¯ ?‹ +îÀÿ(B4A¬ÿÖA¥š$ÀÀo5748g0447o.358840.cm.982715.smf[XY66.hdr]¿BížœZ!ÆwRÆmNA?ƒ?Œ?ˆë¿ÜÿœOYŒ¯J[¿ò64ŠŒk5,55àèm5èŽ85º^l65t8,/f8¬¹ìâ­íŸl®ï|® +]®~­°A³®E„óDEC--WRPN]¯ä?‹ +îÀÿ(B4A¬ÿÖA¥š$ÀÀo5748g0447o.358840.cm.982715.smf[XY67.hdr]ºCîžœZ!Æg}Fma.? +Ñ?ƒŽL¿òù<¥8,<;x¿ýÀ5ÄŽ·7 +µ›€J5E +ù5‡ +ûµ”œ8­æÉè,³0×-žÄ_­‡Eº­ÞÊ-‹'-”­ß"ôDEC--WRPN]¯u?‹ +îÀÿ(B4A¬ÿÖA¥š$ÀÀo5748g0447o.358840.cm.982715.smf[XY71.hdr]»GïžœZ!ÆfŸF +#?†mÜ?‚¯ +¿û